Method

NumCosmoWLEllipticitySeriesTraceDeteval

Declaration [src]

void
nc_wl_ellipticity_series_trace_det_eval (
  NcWLEllipticitySeriesTraceDet* ser,
  gdouble rho
)

Description [src]

Refills sers chi and jac series (see nc_wl_ellipticity_series_trace_det_get_chi()/ _get_jac()) at $\chi_L=\rho\,w$.

eps: $\chi_I(\chi_L,g)=(\chi_L-g)/(1-g\chi_L)$, closed form, no recursion: $c_0=\chi_L={1{:}\rho}$; $c_k=\chi_L^{k-1}(\chi_L^2-1)= \rho^{k+1}w^{k+1}-\rho^{k-1}w^{k-1}$ for $k\ge1$ — a 2-term Laurent series at every order, written directly into sers own stable slots (no scratch needed at all).

eps: $f_\mathrm{eps}’(\chi_L,-g)=(1-g^2)/(1-g\chi_L)^2$ Taylor-in-$g$ coefficients — closed form: $1/(1-g\chi_L)^2=\sum_n(n+1)\chi_L^n g^n$ (literal binomial series, single-term at harmonic $n$ each), times $(1-g^2)$. $\mathrm{Jac}=|\mathrm{that}|^2$ is this series convolved with its own conjugate. Also fills nc_wl_ellipticity_series_trace_det_get_abs_sq() ($|\chi_I(\chi_L,g)|^2=\chi_I\overline{\chi_I}$, population-independent, exactly the $\rho^2(g)$ input nc_galaxy_shape_pop_eval_p_rho2_g_series() needs), since it’s a pure function of chi already computed here.

Parameters

rho

Type: gdouble

$\rho=|\chi_L|$.