ccl_test_correlation_3d.c File Reference

#include "ccl.h"
#include "ctest.h"
#include <stdlib.h>
#include <stdio.h>
#include <math.h>

Include dependency graph for ccl_test_correlation_3d.c:

Go to the source code of this file.

Classes
struct	corrs_3d_data

Functions
void	__attribute__ ((weak))
static int	linecount (FILE *f)
static void	compare_correlation_3d (int i_model, struct corrs_3d_data *data)
void	__ctest_corrs_3d_model_1_run (struct corrs_3d_data *data)
void	__ctest_corrs_3d_model_2_run (struct corrs_3d_data *data)
void	__ctest_corrs_3d_model_3_run (struct corrs_3d_data *data)

Variables
static struct corrs_3d_data	__ctest_corrs_3d_data

Function Documentation

void __attribute__

(

(weak)

)

Definition at line 23 of file ccl_test_correlation_3d.c.

References ccl_mnu_sum.

                      {
  data->Omega_c = 0.25;
  data->Omega_b = 0.05;
  data->h = 0.7;
  data->A_s = 0.8;
  data->n_s = 0.96;
  data->sigma8=0.8;
  data->Neff=3.046;
  double mnuval = 0.;
  data->mnu= &mnuval;
  data->mnu_type = ccl_mnu_sum;

  double Omega_v[5]={0.7, 0.7, 0.7, 0.65, 0.75};
  double w_0[5] = {-1.0, -0.9, -0.9, -0.9, -0.9};
  double w_a[5] = {0.0, 0.0, 0.1, 0.1, 0.1};
  for(int i=0;i<5;i++) {
    data->Omega_v[i] = Omega_v[i];
    data->w_0[i] = w_0[i];
    data->w_a[i] = w_a[i];
    data->Omega_k[i] = 1.0 - data->Omega_c - data->Omega_b - data->Omega_v[i];
  }
}

void __ctest_corrs_3d_model_1_run

(

struct corrs_3d_data *

data

)

Definition at line 154 of file ccl_test_correlation_3d.c.

References compare_correlation_3d().

                         {
  int model=1;
  compare_correlation_3d(model,data);
}

void __ctest_corrs_3d_model_2_run

(

struct corrs_3d_data *

data

)

Definition at line 159 of file ccl_test_correlation_3d.c.

References compare_correlation_3d().

                         {
  int model=2;
  compare_correlation_3d(model,data);
}

void __ctest_corrs_3d_model_3_run

(

struct corrs_3d_data *

data

)

Definition at line 164 of file ccl_test_correlation_3d.c.

References compare_correlation_3d().

                         {
  int model=3;
  compare_correlation_3d(model,data);
}

static void compare_correlation_3d ( int  i_model, struct corrs_3d_data *  data  )

static

Definition at line 58 of file ccl_test_correlation_3d.c.

References corrs_3d_data::A_s, ASSERT_DBL_NEAR_TOL, ASSERT_NOT_NULL, ccl_boltzmann, ccl_correlation_3d(), ccl_cosmology_create(), ccl_cosmology_free(), ccl_halofit, ccl_parameters_create(), cl_cmbl_bm::cosmo, default_config, corrs_3d_data::h, linecount(), ccl_configuration::matter_power_spectrum_method, corrs_3d_data::mnu, corrs_3d_data::mnu_type, corrs_3d_data::n_s, corrs_3d_data::Neff, corrs_3d_data::Omega_b, corrs_3d_data::Omega_c, ccl_parameters::Omega_g, corrs_3d_data::Omega_k, ccl_parameters::Omega_l, corrs_3d_data::Omega_v, halomod_bm::params, corrs_3d_data::sigma8, ccl_parameters::sigma8, ccl_cosmology::status_message, ccl_configuration::transfer_function_method, corrs_3d_data::w_0, corrs_3d_data::w_a, and z.

Referenced by __ctest_corrs_3d_model_1_run(), __ctest_corrs_3d_model_2_run(), and __ctest_corrs_3d_model_3_run().

{
  int nk,nr,i,j;
  int status=0;
  char fname[256],str[1024];
  char* rtn;
  FILE *f;
  ccl_configuration config = default_config;
  config.matter_power_spectrum_method= ccl_halofit;
  config.transfer_function_method = ccl_boltzmann;
  ccl_parameters params = ccl_parameters_create(data->Omega_c,data->Omega_b,data->Omega_k[i_model-1],
        data->Neff, data->mnu, data->mnu_type, data->w_0[i_model-1],data->w_a[i_model-1],
        data->h,data->A_s,data->n_s,-1, -1, -1, -1,NULL,NULL, &status);
  params.Omega_g=0.0;
  params.Omega_l=data->Omega_v[i_model-1];
  params.sigma8=data->sigma8;
  ccl_cosmology * cosmo = ccl_cosmology_create(params, config);
  ASSERT_NOT_NULL(cosmo);      
  
  sprintf(fname,"./tests/benchmark/model%d_xi.txt",i_model);
  f=fopen(fname,"r");
  if(f==NULL) {
    fprintf(stderr,"Error opening file %s\n",fname);
    exit(1);
  }
  nr=linecount(f)-1; rewind(f);

  // tolerence on abs difference in r^2 xi(r) for the range r = 0.1 - 100 Mpc (40 points in r) for z=0,1,2,3,4,5
  double CORR_TOLERANCE1[6] = {2.84e-2, 3.57e-3, 1.37e-3, 1.26e-3, 2.0e-4, 1.8e-4};
  // tolerence on abs difference in r^2 xi(r) for the range r = 50 - 250 Mpc (100 points in r) for z=0,1,2,3,4,5
  double CORR_TOLERANCE2[6] = {3.15e-2, 9.87e-3, 5.04e-3, 3.05e-3, 9.98e-4, 7.14e-4};

  int N1=40;
  double *r_arr1=malloc(N1*sizeof(double));
  double *r_arr2=malloc((nr-N1)*sizeof(double));
  double (*ximm_bench_arr)[6]=malloc(6*nr*sizeof(double));
  
  rtn = fgets(str, 1024, f);
  for(i=0;i<nr;i++) {
    double r_h;
    int stat;
    stat=fscanf(f,"%lf",&r_h);
    if(stat!=1) {
      fprintf(stderr,"Error reading file %s, line %d\n",fname,i+2);
      exit(1);
    }

    if(i<N1)
      r_arr1[i]=r_h;
    else
      r_arr2[i-N1]=r_h;

    for(j=0;j<6;j++) {
      double ximm_bench;      
      stat=fscanf(f,"%lf",&ximm_bench);
      if(stat!=1) {
    fprintf(stderr,"Error reading file %s, line %d\n",fname,i+2);
    exit(1);
      }

      ximm_bench_arr[i][j]=ximm_bench;
    }
  }

  double *ximm_ccl_out1=malloc(N1*sizeof(double));
  double *ximm_ccl_out2=malloc((nr-N1)*sizeof(double));

   for(j=0;j<6;j++) {
      double z = j+0.;  
      ccl_correlation_3d(cosmo,1.0/(j+1),N1,r_arr1,ximm_ccl_out1,0,NULL,&status);
      ccl_correlation_3d(cosmo,1.0/(j+1),nr-N1,r_arr2,ximm_ccl_out2,0,NULL,&status);

      if (status) printf("%s\n",cosmo->status_message);
      for(i=0;i<nr;i++){     
      double err;
      if(i<N1){
      err=fabs(r_arr1[i]*r_arr1[i]*(ximm_ccl_out1[i]-ximm_bench_arr[i][j])); 
      ASSERT_DBL_NEAR_TOL(0.,err,CORR_TOLERANCE1[j]);
      }
      else{
      err=fabs(r_arr2[i-N1]*r_arr2[i-N1]*(ximm_ccl_out2[i-N1]-ximm_bench_arr[i][j])); 
      ASSERT_DBL_NEAR_TOL(0.,err,CORR_TOLERANCE2[j]);
      }
      }   
  }
  fclose(f);
  
  free(r_arr1);
  free(r_arr2);
  free(ximm_bench_arr);
  free(ximm_ccl_out1);
  free(ximm_ccl_out2); 

  ccl_cosmology_free(cosmo);
}

static int linecount ( FILE *  f )

static

Definition at line 46 of file ccl_test_correlation_3d.c.

Referenced by compare_correlation_3d().

{
  // Counts #lines from file
  int i0=0;
  char ch[1000];
  while((fgets(ch,sizeof(ch),f))!=NULL) {
    i0++;
  }
  return i0;
}

Variable Documentation

static struct corrs_3d_data __ctest_corrs_3d_data

static

Definition at line 154 of file ccl_test_correlation_3d.c.