-
Bispectrum_ReadParams(B, outroot, Ini)
-
TBispectrumParams :: B
-
character(LEN=*) intent(in) :: outroot
-
TIniFile :: Ini
Bispectrum_SetDefParams(B)
TBispectrumParams :: B
GetBispectrum(State, CTrans)
Note: may need high maxetak to make sure oscillatory k integrals cancel correctly
for accurate alpha(r), beta(r), e.g. 8000; not so important for bispectrum
increase accuracy_boost
CAMBdata target :: State
ClTransferData :: CTrans
InitBesselDerivs(CTrans)
j_l' array for interpolation if needed; not needed for local fnl
ClTransferData :: CTrans
NonGauss_deriv_l_r(CP, CTrans, dJl, dddJl, indP, resP, r)
As above, but integral against derivative of bessel function to get derivative of function
CAMBParams :: CP
ClTransferData :: CTrans
real(dl) intent(in) :: dJl(BessRanges%npoints,CTrans%ls%nl)
real(dl) intent(in) :: dddJl(BessRanges%npoints,CTrans%ls%nl)
integer intent(in) :: indP(:)
real(dl) :: resP(CTrans%ls%nl,size(indP)
real(dl) intent(in) :: r
NonGauss_l_r(CP, CTrans, nfields, ind, indP, res, r)
functions of the form int dk k^2 k^i j_l(kr) Delta_l(k) [P]
ind and indP are arrays of required k^i powers
res and resP are the results without and with the power spectrum P in the integrand
Output of P scaled by 1d10 (so bispectrum by 1d20)
CAMBParams :: CP
ClTransferData :: CTrans
integer :: nfields
integer intent(in) :: ind(:)
integer intent(in) :: indP(:)
real(dl) :: res(CTrans%ls%nl,size(ind)
real(dl) intent(in) :: r
NonGauss_l_r_localOpt(CP, CTrans, ind, indP, nfields, res, r)
functions of the form int dk k^2 k^i j_l(kr) Delta_l(k) [P]
ind and indP are arrays of required k^i powers
res and resP are the results without and with the power spectrum P in the integrand
CAMBPArams :: CP
ClTransferData :: CTrans
integer intent(in) :: ind(:)
integer intent(in) :: indP(:)
integer :: nfields
real(dl) :: res(CTrans%ls%nl,size(ind)
real(dl) intent(in) :: r