ccl_test_bcm.c File Reference

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

Include dependency graph for ccl_test_bcm.c:

Go to the source code of this file.

Classes
struct	bcm_data

Macros
#define	BCM_TOLERANCE   1e-4

Functions
void	__attribute__ ((weak))
static int	linecount (FILE *f)
static void	compare_bcm (int i_model, struct bcm_data *data)
void	__ctest_bcm_model_1_run (struct bcm_data *data)

Variables
static struct bcm_data	__ctest_bcm_data

Macro Definition Documentation

#define BCM_TOLERANCE   1e-4

Definition at line 6 of file ccl_test_bcm.c.

Referenced by compare_bcm().

Function Documentation

void __attribute__

(

(weak)

)

Definition at line 24 of file ccl_test_bcm.c.

References ccl_mnu_sum.

                 {
  data->Omega_c = 0.25;
  data->Omega_b = 0.05;
  data->h = 0.7;
  data->A_s = 2.2e-9;
  data->n_s = 0.96;
  data->Neff = 3.046;
  double mnuval = 0.;
  data->m_nu= &mnuval;
  data-> mnu_type = ccl_mnu_sum;

  double Omega_v[1]={0.7};
  double w_0[1] = {-1.0};
  double w_a[1] = {0.0};

  for(int i=0;i<1;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_bcm_model_1_run

(

struct bcm_data *

data

)

Definition at line 150 of file ccl_test_bcm.c.

References compare_bcm().

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

static void compare_bcm ( int  i_model, struct bcm_data *  data  )

static

Definition at line 59 of file ccl_test_bcm.c.

References bcm_data::A_s, ASSERT_DBL_NEAR_TOL, ASSERT_NOT_NULL, ccl_configuration::baryons_power_spectrum_method, BCM_TOLERANCE, ccl_bcm, ccl_bcm_model_fka(), ccl_cosmology_create(), ccl_cosmology_free(), ccl_nonlin_matter_power(), ccl_parameters_create(), cl_cmbl_bm::cosmo, default_config, bcm_data::h, linecount(), bcm_data::m_nu, bcm_data::mnu_type, bcm_data::n_s, bcm_data::Neff, bcm_data::Omega_b, bcm_data::Omega_c, ccl_parameters::Omega_g, bcm_data::Omega_k, ccl_parameters::Omega_l, halomod_bm::params, ccl_cosmology::status_message, bcm_data::w_0, and bcm_data::w_a.

Referenced by __ctest_bcm_model_1_run().

{
  int nk,i,j;
  int status =0;
  char fname[256],str[1024];
  char fname2[256];
  char* rtn;
  FILE *f,*f2;
  ccl_configuration config = default_config;
  config.baryons_power_spectrum_method=ccl_bcm;
  ccl_parameters params = ccl_parameters_create(data->Omega_c,data->Omega_b,data->Omega_k[i_model-1],
                        data->Neff, data->m_nu, data-> mnu_type,
                        data->w_0[i_model-1],data->w_a[i_model-1],
                        data->h,data->A_s,data->n_s,14,-1,-1,-1,NULL,NULL, &status);
                                                                        
  params.Omega_l=params.Omega_l+params.Omega_g;
  params.Omega_g=0;
  ccl_cosmology * cosmo = ccl_cosmology_create(params, config);
  ASSERT_NOT_NULL(cosmo);
  ccl_configuration config_nobar = default_config;
  ccl_parameters params_nobar = ccl_parameters_create(data->Omega_c,data->Omega_b,data->Omega_k[i_model-1],
                        data->Neff, data->m_nu, 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.sigma_8=data->sigma_8;
  params_nobar.Omega_l=params_nobar.Omega_l+params_nobar.Omega_g;
  params_nobar.Omega_g=0;
  ccl_cosmology * cosmo_nobar = ccl_cosmology_create(params_nobar, config_nobar);
  ASSERT_NOT_NULL(cosmo_nobar);
  
  sprintf(fname,"./tests/benchmark/bcm/w_baryonspk_nl.dat");
  f=fopen(fname,"r");
  if(f==NULL) {
    fprintf(stderr,"Error opening file %s\n",fname);
    exit(1);
  }
  sprintf(fname2,"./tests/benchmark/bcm/wo_baryonspk_nl.dat");
  f2=fopen(fname2,"r");
  if(f2==NULL) {
    fprintf(stderr,"Error opening file %s\n",fname2);
    exit(1);
  }
  nk=linecount(f)-1; rewind(f);

  //Skip first 4 lines
  rtn = fgets(str, 1024, f);
  rtn = fgets(str, 1024, f);
  rtn = fgets(str, 1024, f);
  rtn = fgets(str, 1024, f);
  
  rtn = fgets(str, 1024, f2);
  rtn = fgets(str, 1024, f2);
  rtn = fgets(str, 1024, f2);
  rtn = fgets(str, 1024, f2);
  
  for(i=0;i<nk-4;i++) {
    double k_h,k;
    int stat;
    double psbar,psnobar,fbcm_bench,err;
    double psbar_bench,psnobar_bench;
    stat=fscanf(f,"%le %le",&k_h,&psbar);
    if(stat!=2) {
      fprintf(stderr,"Error reading file %s, line %d\n",fname,i+5);
      exit(1);
    }
    stat=fscanf(f2,"%*le %le",&psnobar);
    if(stat!=1) {
      fprintf(stderr,"Error reading file %s, line %d\n",fname2,i+5);
      exit(1);
    }
    k=k_h*data->h;
    //Check baryonic correction directly
    fbcm_bench=ccl_bcm_model_fka(cosmo,k,1.,&status);
    if (status) printf("%s\n",cosmo->status_message);
    err=fabs(psbar/psnobar/fbcm_bench-1);
    ASSERT_DBL_NEAR_TOL(err,0.,BCM_TOLERANCE);
    //And check the ratio between power spectra
    psbar_bench=ccl_nonlin_matter_power(cosmo,k,1.,&status);
    if (status) printf("%s\n",cosmo->status_message);
    psnobar_bench=ccl_nonlin_matter_power(cosmo_nobar,k,1.,&status);
    if (status) printf("%s\n",cosmo_nobar->status_message);
    err=fabs(psbar/psnobar/(psbar_bench/psnobar_bench)-1);
    ASSERT_DBL_NEAR_TOL(err,0.,BCM_TOLERANCE);
    
  }
  fclose(f);
  fclose(f2);
  
  ccl_cosmology_free(cosmo);
}

static int linecount ( FILE *  f )

static

Definition at line 47 of file ccl_test_bcm.c.

Referenced by compare_bcm().

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

Variable Documentation

struct bcm_data __ctest_bcm_data

static

Definition at line 150 of file ccl_test_bcm.c.