Skip to content

Update prec model

Serguei Ossokine requested to merge update_prec_model into main

This is an omnibus update to the precessing model. Main changes are:

  • Lowered the ODE tolerances from 1e-9 to 1e-8 for faster model evaluation.
  • Added in integrate_ode_prec.py an extra stopping condition of the ODE integration if omega becomes nan.
  • Remove the storage of LNhat (it is an unnecessary calculation).
  • Check if (omega_start)**(-2./3.) < 10.5 and if so use a omega_start corresponding to 10.5M.
  • Finalized version of the PA routine (postadiabatic_C_prec.pyx).
  • PA initial conditions which can be used by the non-PA model (although the default model is the PA model). Improved documentation.

Merge request reports