Eddy

class maelstrom.eddy.Eddy(time: numpy.ndarray, tau: numpy.ndarray, period_guess=None, asini_guess=None)

A PyMC3 model for modelling light travel time variations in a binary system using the subdividing method.

Parameters

• time (array) – Time observations

• tau (array) – Flux observations

• period_guess (float, optional) – Initial guess of the orbital period,

• units of days. Defaults to None. (in) –

• asini_guess (float, optional) – Initial guess of the projected

• axis, in units of seconds. Defaults to None. (semi-major) –

optimize(vars=None)

Optimises the model.

Parameters

• vars (list, optional) – List of parameters in the model to optimize.

• none, will optimize all parameters at once. Defaults to None. (If) –

Returns

Optimisation results

Return type

dict

sample(tune=3000, draws=3000, start=None, target_accept=0.9, **kwargs)

Samples the midel using the exoplanet PyMC3 sampler. By default, this will sample from 2 chains over 2 cores simultaneously.

Parameters

• tune (int, optional) – Number of tuning steps. Defaults to 3000.

• draws (int, optional) – Number of samples. Defaults to 3000.

• start (dict, optional) – Starting location of the sampler. If none is

• the sampler will first optimize the model. Defaults to None. (supplied,) –

• target_accept (float, optional) – Target acceptance ratio of the NUTS

• Defaults to 0.9. (sampler.) –

Returns

A PyMC3 trace of the posterior.

Return type

trace