CCL
Data Structures | Macros | Functions
ccl_correlation.c File Reference
#include "gsl/gsl_integration.h"
#include "ccl_cls.h"
#include "gsl/gsl_errno.h"
#include "gsl/gsl_roots.h"
#include "gsl/gsl_spline.h"
#include "gsl/gsl_sf_bessel.h"
#include "gsl/gsl_sf_legendre.h"
#include "ccl_error.h"
#include "ccl_utils.h"
#include "ccl_correlation.h"
#include <stdlib.h>
#include <math.h>
#include <stdio.h>
#include <string.h>
#include "ccl_power.h"
#include "ccl.h"
#include "fftlog.h"

Data Structures

struct  corr_int_par
 

Macros

#define ELL_MIN_FFTLOG   0.01
 
#define ELL_MAX_FFTLOG   60000
 
#define N_ELL_FFTLOG   5000
 

Functions

void ccl_correlation (ccl_cosmology *cosmo, int n_ell, double *ell, double *cls, int n_theta, double *theta, double *wtheta, int corr_type, int do_taper_cl, double *taper_cl_limits, int flag_method, int *status)
 

Macro Definition Documentation

#define ELL_MAX_FFTLOG   60000
#define ELL_MIN_FFTLOG   0.01
#define N_ELL_FFTLOG   5000

Function Documentation

void ccl_correlation ( ccl_cosmology cosmo,
int  n_ell,
double *  ell,
double *  cls,
int  n_theta,
double *  theta,
double *  wtheta,
int  corr_type,
int  do_taper_cl,
double *  taper_cl_limits,
int  flag_method,
int *  status 
)

Computes the correlation function (wrapper)

Parameters
cosmo:Cosmological parameters
n_ell: number of multipoles in the input power spectrum
ell: multipoles at which the power spectrum is evaluated
cls: input power spectrum
n_theta: number of output values of the separation angle (theta)
theta: values of the separation angle in degrees.
wtheta: the values of the correlation function at the angles above will be returned in this array, which should be pre-allocated
do_taper_cl:
taper_cl_limits
flag_method: method to compute the correlation function. Choose between:
  • CCL_CORR_FFTLOG : fast integration with FFTLog
  • CCL_CORR_BESSEL : direct integration over the Bessel function
  • CCL_CORR_LGNDRE : brute-force sum over legendre polynomials
corr_type: type of correlation function. Choose between:
  • CCL_CORR_GG : spin0-spin0
  • CCL_CORR_GL : spin0-spin2
  • CCL_CORR_LP : spin2-spin2 (xi+)
  • CCL_CORR_LM : spin2-spin2 (xi-) Currently supported spin-0 fields are number counts and CMB lensing. The only spin-2 is currently shear.