NumCosmoMathCSQ1Deval_delta_theta_at

gdouble
ncm_csq1d_eval_delta_theta_at (
  NcmCSQ1D* csq1d,
  const gdouble t
)

Returns the accumulated residual phase shift $\delta\theta(t)$ computed during the last call to ncm_csq1d_prepare().

The full phase entering the mode functions is $\theta(t) = \int_{t_i}^{t} \nu(t’)\,\mathrm{d}t’ + \delta\theta(t)$.

For the full phase and residual-phase equations, including the numerically stable form used internally, see CSQ1D Formalism.

Note: You must call ncm_csq1d_prepare_phase_splines() before using this function, or it will throw an error.

t

Type: const gdouble

Time $t$.

Type: gdouble

$\delta\theta(t)$.