thejoker
thejoker copied to clipboard
adrn
A custom Monte Carlo sampler for the (gravitational) two-body problem
This is a rebase of #107 using the new infrastructure from v1.3. Right now, this only works if the velocity measurements can be identified with the correct binary component. Ideally,...
See note and sample data from @AlexLaroche7 @astroludwig
the function plot_phase_fold raise error. In the plot.py file the funciton plot_phase_fold are defined as: def plot_phase_fold(sample, data=None, ax=None, with_time_unit=False, n_phase_samples=4096, add_labels=True, show_s_errorbar=True, residual=False, remove_trend=True, plot_kwargs=None, data_plot_kwargs=None) this function will...
In paper, the orbital paramters are v0, K, e, φ0, P and omega, but in code, there is M0 but no φ0.
Hi, I want to fix the period, which should how to set up? Thanks!
Disabled because pytables would not install on windows / github actions properly.
See email by Song Wang: ```python with schwimmbad.MultiPool() as pool: joker_mcmc = tj.TheJoker(prior_mcmc, pool=pool, random_state=rnd) mcmc_init = joker_mcmc.setup_mcmc(data, samples) with schwimmbad.MultiPool() as pool: joker = tj.TheJoker(prior, pool=pool, random_state=rnd) prior_samples =...
As originally asked by @astrosong: If the sigma_K0 and sigma_v (and P_min, P_max) are different for each object, which means a different prior for each object, then can we still...
I am running the sampler using MultiPool(), and keeping an eye on my activity monitor. I see eight python processes run for a few minutes, and then they stop --...
