Dependencies

Subroutines
  • CorrFuncFullSkyImpl(lmax)
    Accurate curved sky correlation function method Uses non-perturbative isotropic term with 2nd order expansion in C_{gl,2} Neglects C_{gl}(theta) terms (very good approx)
    • integer intent(in) :: lmax
  • GetBessels(MaxArg)
    • real(dl) intent(in) :: MaxArg
Functions
  • double precision
    BESSI(N, X)
    from http://perso.orange.fr/jean-pierre.moreau/Fortran/tbessi_f90.txt This subroutine calculates the first kind modified Bessel function of integer order N, for any REAL X. We use here the classical recursion formula, when X > N. For X < N, the Miller's algorithm is used to avoid overflows. REFERENCE: C.W.CLENSHAW, CHEBYSHEV SERIES FOR MATHEMATICAL FUNCTIONS, MATHEMATICAL TABLES, VOL.5, 1962.
    • integer intent(in) :: N
    • double precision  :: X
  • double precision
    BESSI0(X)
    • double precision  :: X
  • double precision
    BESSI1(X)
    • double precision  :: X
  • real(dl)
    bessj(x, n)
    • real(dl)  :: x
    • integer  :: n
  • real(dl)
    bessj0(x)
    • real(dl)  :: x