static gdouble _nc_window_tophat_eval_fourier (const NcWindow * wp, const gdouble k, const gdouble R) { gdouble kR = k * R; NCM_UNUSED (wp); if (kR == 0.0) return 1.0; return 3.0 * gsl_sf_bessel_j1 (kR) / kR; }
/// Integrand for var_TH: 1/k \Delta^2(k) (3.*j1(kR)/(kR))^2 k in hMpc^{-1} double dvar_TH(double x, void * params) { cvar_params c1 = *(cvar_params *) params; // std::cout<<1<<std::endl; cosmology *c2; double *R; double *z; bool *psi; c2=c1.cptr; R=c1.R; z=c1.z; psi=c1.psinit; double arg=(*R)*x; double d2; if(*psi){ d2=(*c2).Delta2_L_num(x,*z); }else{ d2=(*c2).Delta2_L(x,*z); } double bes= pow(3.*gsl_sf_bessel_j1(arg)/(arg) ,2.); //std::cout<<1./x*d2*bes<<"Is this quick"<<std::endl; return 1./x*d2*bes; }