pyseobnr.eob.dynamics.postadiabatic_C_prec

Contains the functions needed for computing the precessing post-adiabatic dynamics

Functions

compute_adiabatic_solution(r, omega, H, ...)

Compute adiabatic solution for pphi by solving dH_dr =0 with prstar=0

compute_combined_dynamics_exp_v1(omega_ref, ...)

Compute the dynamics by using a postadiabatic scheme up to a certain separation from which an ODE integration is started.

compute_postadiabatic_dynamics(omega_ref, ...)

Compute postadiabatic dynamics from the spins specified at a certain reference frequency and with a certain starting frequency.

compute_postadiabatic_solution(r, pphi, ...)

Compute postadiabatic solution for prstar, pphi and omega up to a certain postadiabatic order.

compute_pphi(r, pr, pphi, omega, H, RR, ...)

Compute pphi at a given postadiabatic order by solving the equation dpr_dr * dH_dpr + dH_dr - flux[0] / csi, where flux[0] is the radial component of the radiation reaction force and csi is the tortoise coordinate transformation factor, for pphi.

compute_pr(r, pr, pphi, omega, H, RR, ...[, ...])

Compute prstar at a given postadiabatic order by solving the equation dpphi_dr * dH_dpr * csi - flux[1] = 0, where flux[1] is the azimuthal component of the radiation reaction force and csi is the tortoise coordinate transformation factor, for prstar.

compute_prec_cycles(r_final, t_pn, omega_pn, ...)

Estimate the number of precession cycles from LNhat, computed as the phase of the precession frequency at the final point of the radial grid divided by 2 pi.

j0_eqn(j0_sol, r, H, chi1_v, chi2_v, m_1, ...)

Equation dH_dr =0 evaluated at prstar=0 used to obtain the adiabatic solution for the orbital angular momentum (pphi)

pphi_eqn(pphi_sol, r, pr, dpr_dr, H, RR, ...)

Equation dpr_dr * dH_dpr + dH_dr - flux[0] / csi, where flux[0] is the radial component of the radiation reaction force and csi is the tortoise coordinate transformation factor, to solve for ppphi in the numerical postadiabatic routine.

pphi_eqn_analytic(pphi_sol, r, pr, dpr_dr, ...)

Equation dpr_dr * dH_dpr + dH_dr - flux[0] / csi, where flux[0] is the radial component of the radiation reaction force and csi is the tortoise coordinate transformation factor, to solve for pphi in the analytic postadiabatic routine.

pr_eqn(pr_sol, r, pphi, dpphi_dr, H, RR, ...)

Equation dpphi_dr * dH_dpr * csi - flux[1] = 0, where flux[1] is the azimuthal component of the radiation reaction force and csi is the tortoise coordinate transformation factor, to solve for prstar in the numerical postadiabatic routine.

pr_eqn_analytic(pr_sol, r, pphi, dpphi_dr, ...)

Equation dpphi_dr * dH_dpr * csi - flux[1] = 0, where flux[1] is the azimuthal component of the radiation reaction force and csi is the tortoise coordinate transformation factor, to solve for prstar in the analytic postadiabatic routine.

precessing_final_spin(chi1_LN, chi2_LN, ...)

Magnitude of the final spin computed from NR fits of precessing NR simulations.

spline_derivative(x, y)

Compute derivative of y(x) using cubic splines.