pyseobnr.eob.dynamics.integrate_ode_prec.compute_dynamics_prec_opt

pyseobnr.eob.dynamics.integrate_ode_prec.compute_dynamics_prec_opt(omega_ref: float, omega_start: float, omegaPN_f: float, H: Hamiltonian, RR: Callable, m_1: float, m_2: float, splines: dict, t_pn: array, dynamics_pn: array, params: EOBParams, rtol: float = 1e-12, atol: float = 1e-12, step_back: float = 250.0, y_init=None, initial_conditions: Literal['adiabatic', 'postadiabatic'] = 'adiabatic', initial_conditions_postadiabatic_type: Literal['analytic', 'numeric'] = 'analytic')[source]

Function to perform a non-precessing EOB evolution with the spins modified at every timestep according to the values from the precessing-spin PN evolution.

Parameters:

  • omega_ref (float) – Reference orbital frequency at which the spins are defined

  • omega_start (float) – Starting orbital frequency

  • omegaPN_f (float) – Final orbital frequency from the precessing-spin PN evolution

  • H (Hamiltonian) – Hamiltonian class

  • RR (Callable) – RR force

  • m_1 (float) – Mass component of the primary

  • m_2 (float) – Mass component of the secondary

  • splines (dict) – Dictionary containing the splines in orbital frequency of the vector components of the spins, LN and L as well as the spin projections onto LN and L

  • t_pn (np.array) – Time array of the PN evolution of the spins and Newtonian angular momentum.

  • dynamics_pn (np.array) – Array of the spin-precessing PN evolution. It contains the Newtonian angular momentum, the dimensionful spin vectors and the PN orbital frequency.

  • params (EOBParams) – Container of additional inputs

  • rtol (float, optional) – Relative tolerance for EOB integration. Defaults to 1e-12

  • atol (float, optional) – Absolute tolerance for EOB integration. Defaults to 1e-12

  • step_back (float, optional) – Amount of time to step back for fine interpolation. Defaults to 250.

  • y_init (np.ndarray, optional) – Initial condition vector (r,phi,pr,pphi).

  • initial_conditions (str, optional) – Type of initial conditions for the ODE evolution (‘adiabatic’ or ‘postadiabatic’).

  • initial_conditions_postadiabatic_type (str, optional) – Type of postadiabatic initial conditions for the ODE evolution (‘analytic’ or ‘numeric’).

Returns:

Low and high sampling rate dynamics, unit Newtonian orbital angular momentum,

assembled dynamics and the index splitting the low and high sampling rate dynamics

Return type:

(tuple)