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;
}
