ccl_corr_sample.c File Reference

#include "ccl.h"
#include "ccl_correlation.h"
#include "ccl_utils.h"
#include "ctest.h"
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <time.h>
#include <string.h>

Include dependency graph for ccl_corr_sample.c:

Go to the source code of this file.

Macros
#define	OC   0.25
#define	OB   0.05
#define	OK   0.00
#define	ON   0.00
#define	HH   0.70
#define	W0   -1.0
#define	WA   0.00
#define	NS   0.96
#define	NORMPS   0.80
#define	ZD   0.5
#define	NZ   128
#define	Z0_GC   0.50
#define	SZ_GC   0.05
#define	NL   512
#define	PS   0.1
#define	NEFF   3.046
#define	ELL_MAX_CL   20000

Functions
int	main (int argc, char **argv)

Macro Definition Documentation

#define ELL_MAX_CL   20000

Definition at line 28 of file ccl_corr_sample.c.

Referenced by main().

#define HH   0.70

Definition at line 16 of file ccl_corr_sample.c.

Referenced by main().

#define NEFF   3.046

Definition at line 27 of file ccl_corr_sample.c.

Referenced by main().

#define NL   512

Definition at line 25 of file ccl_corr_sample.c.

#define NORMPS   0.80

Definition at line 20 of file ccl_corr_sample.c.

Referenced by main().

#define NS   0.96

Definition at line 19 of file ccl_corr_sample.c.

Referenced by main().

#define NZ   128

Definition at line 22 of file ccl_corr_sample.c.

Referenced by main().

#define OB   0.05

Definition at line 13 of file ccl_corr_sample.c.

Referenced by main().

#define OC   0.25

Definition at line 12 of file ccl_corr_sample.c.

Referenced by main().

#define OK   0.00

Definition at line 14 of file ccl_corr_sample.c.

Referenced by main().

#define ON   0.00

Definition at line 15 of file ccl_corr_sample.c.

#define PS   0.1

Definition at line 26 of file ccl_corr_sample.c.

#define SZ_GC   0.05

Definition at line 24 of file ccl_corr_sample.c.

Referenced by main().

#define W0   -1.0

Definition at line 17 of file ccl_corr_sample.c.

Referenced by main().

#define WA   0.00

Definition at line 18 of file ccl_corr_sample.c.

Referenced by main().

#define Z0_GC   0.50

Definition at line 23 of file ccl_corr_sample.c.

Referenced by main().

#define ZD   0.5

Definition at line 21 of file ccl_corr_sample.c.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 31 of file ccl_corr_sample.c.

References ccl_angular_cls(), ccl_cl_tracer_free(), ccl_cl_tracer_number_counts_simple(), ccl_cl_workspace_free(), ccl_cl_workspace_new_limber(), CCL_CORR_FFTLOG, CCL_CORR_GG, ccl_correlation(), ccl_cosmology_create(), ccl_cosmology_free(), ccl_log_spacing(), ccl_mnu_sum, ccl_parameters_create(), cl_cmbl_bm::cosmo, default_config, ELL_MAX_CL, HH, ccl_test_distances::mnu, NEFF, NORMPS, NS, NZ, OB, OC, OK, halomod_bm::params, pow(), SZ_GC, w, W0, WA, and Z0_GC.

{
  // Use the default configuration, plus the cosmological parameters that were 
  // defined above
  int status=0;
  
  double mnuval=0.;
  double* mnu;
  mnu = &mnuval;
  ccl_mnu_convention mnu_type = ccl_mnu_sum;
  
  ccl_configuration config = default_config;
  ccl_parameters params = ccl_parameters_create(OC, OB, OK, NEFF, mnu, mnu_type,
                                                W0, WA, HH, NORMPS, NS, -1, -1, -1, -1,
                        NULL, NULL, &status);
  ccl_cosmology *cosmo = ccl_cosmology_create(params,config);

  // Create example number density and bias for tracer
  double z_arr_gc[NZ], nz_arr_gc[NZ], bz_arr[NZ];
  for(int i=0; i < NZ; i++) {
    z_arr_gc[i] = Z0_GC - 5.*SZ_GC + 10.*SZ_GC*(i + 0.5)/NZ;
    nz_arr_gc[i] = exp(-0.5*pow((z_arr_gc[i] - Z0_GC) / SZ_GC, 2));
    bz_arr[i] = 1. + z_arr_gc[i];
  }

  // Define a galaxy clustering tracer and calculate C_ell's
  double l_logstep = 1.05;
  double l_linstep = 5.;
  CCL_ClWorkspace *w=ccl_cl_workspace_new_limber(ELL_MAX_CL+1,l_logstep,l_linstep,&status);
  //CCL_ClTracer *tr_nc_1=ccl_cl_tracer_number_counts_simple(cosmo,nz,zarr_1,pzarr_1,nz,zarr_1,bzarr,&status);
  CCL_ClTracer *ct_gc = ccl_cl_tracer_number_counts_simple(cosmo, 
                               NZ, z_arr_gc, nz_arr_gc,
                               NZ, z_arr_gc, bz_arr,
                               &status);
  int il;
  double *clarr = malloc(ELL_MAX_CL*sizeof(double));
  int *ilarr = malloc(ELL_MAX_CL*sizeof(int));
  double *larr = malloc(ELL_MAX_CL*sizeof(double));
  for(il=0; il < ELL_MAX_CL; il++)
    // Calculate auto-correlation of this tracer
    larr[il] = il;
  ccl_angular_cls(cosmo, w, ct_gc, ct_gc, ELL_MAX_CL, ilarr, clarr, &status);
  
  // Define cosine tapering, to reduce ringing. The first two numbers are 
  // [lmin, lmax] for the low-ell taper, and the last two are [lmin, lmax] for 
  // the high-ell taper.
  double taper_cl_limits[4] = {1, 2, 10000, 15000};
  
  double *clustering_corr, *theta;
  int ntheta = 15;
  theta = ccl_log_spacing(0.01, 5., ntheta); // New array with log spacing
  clustering_corr = malloc(ntheta*sizeof(double));
  
  // Calculate correlation function from angular power spectrum
  ccl_correlation(cosmo, ELL_MAX_CL, larr, clarr, 
                  ntheta, theta, clustering_corr, 
                  CCL_CORR_GG, 0, taper_cl_limits, 
                  CCL_CORR_FFTLOG, &status);
  
  // Print results
  for(int it=0; it < ntheta; it++)
    printf("%le %le\n", theta[it], clustering_corr[it]);
  
  // Free tracers and other allocated memory
  ccl_cl_workspace_free(w);
  ccl_cl_tracer_free(ct_gc);
  ccl_cosmology_free(cosmo);
  free(clustering_corr);
  free(ilarr);
  free(larr);
  free(clarr);

  return 0;
}