ccl_3dcorr_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_3dcorr_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	K_MIN   0.001
#define	K_MAX   100
#define	N_MAX_PK   20000

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

Macro Definition Documentation

#define HH   0.70

Definition at line 16 of file ccl_3dcorr_sample.c.

Referenced by main().

#define K_MAX   100

Definition at line 29 of file ccl_3dcorr_sample.c.

Referenced by main().

#define K_MIN   0.001

Definition at line 28 of file ccl_3dcorr_sample.c.

#define N_MAX_PK   20000

Definition at line 30 of file ccl_3dcorr_sample.c.

#define NEFF   3.046

Definition at line 27 of file ccl_3dcorr_sample.c.

Referenced by main().

#define NL   512

Definition at line 25 of file ccl_3dcorr_sample.c.

#define NORMPS   0.80

Definition at line 20 of file ccl_3dcorr_sample.c.

Referenced by main().

#define NS   0.96

Definition at line 19 of file ccl_3dcorr_sample.c.

Referenced by main().

#define NZ   128

Definition at line 22 of file ccl_3dcorr_sample.c.

#define OB   0.05

Definition at line 13 of file ccl_3dcorr_sample.c.

Referenced by main().

#define OC   0.25

Definition at line 12 of file ccl_3dcorr_sample.c.

Referenced by main().

#define OK   0.00

Definition at line 14 of file ccl_3dcorr_sample.c.

Referenced by main().

#define ON   0.00

Definition at line 15 of file ccl_3dcorr_sample.c.

#define PS   0.1

Definition at line 26 of file ccl_3dcorr_sample.c.

#define SZ_GC   0.05

Definition at line 24 of file ccl_3dcorr_sample.c.

#define W0   -1.0

Definition at line 17 of file ccl_3dcorr_sample.c.

Referenced by main().

#define WA   0.00

Definition at line 18 of file ccl_3dcorr_sample.c.

Referenced by main().

#define Z0_GC   0.50

Definition at line 23 of file ccl_3dcorr_sample.c.

#define ZD   0.5

Definition at line 21 of file ccl_3dcorr_sample.c.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 33 of file ccl_3dcorr_sample.c.

References ccl_correlation_3d(), ccl_cosmology_create(), ccl_cosmology_free(), ccl_halofit, ccl_log_spacing(), ccl_mnu_sum, ccl_parameters_create(), cl_cmbl_bm::cosmo, default_config, HH, ccl_configuration::matter_power_spectrum_method, ccl_test_distances::mnu, NEFF, NORMPS, NS, OB, OC, OK, halomod_bm::params, W0, and WA.

{
  // Use the default configuration, plus the cosmological parameters that were 
  // defined above
  int status=0;
  ccl_configuration config = default_config;
  config.matter_power_spectrum_method=ccl_halofit;
  
  double mnuval = 0.;
  double* mnu;
  mnu = &mnuval;
  ccl_mnu_convention mnu_type = ccl_mnu_sum;

  ccl_parameters params = ccl_parameters_create(OC, OB, OK, NEFF, mnu, mnu_type,
                                                W0, WA, HH, NORMPS, NS,
                        14.079181246047625, 0.5, 55, 0, NULL, NULL, &status);
  ccl_cosmology *cosmo = ccl_cosmology_create(params,config);
  
  // Define cosine tapering, to reduce ringing. The first two numbers are 
  // [kmin, kmax] for the low-k taper, and the last two are [kmin, kmax] for 
  // the high-k taper.
  double taper_pk_limits[4] = {0.001, 0.002, 50, 75};
  
  double *xi, *r;
  int n_r = 20;
  r = ccl_log_spacing(0.1, 50., n_r); // New array with log spacing
  xi = malloc(n_r*sizeof(double));

  // Calculate 3dcorrelation function from power spectrum
  ccl_correlation_3d(cosmo, 1.0, 
                  n_r, r, xi,
                  0, taper_pk_limits, 
                  &status);
  
  // Print results
  for(int it=0; it < n_r; it++)
    printf("%le %le\n", r[it], xi[it]);
  
  // Free allocated memory
  ccl_cosmology_free(cosmo);
  free(xi);

  return 0;
}