From ae6faa03fd2f3e6746024a089db90a19b50f4dbd Mon Sep 17 00:00:00 2001 From: damonge Date: Fri, 10 Apr 2026 10:33:16 +0100 Subject: [PATCH 1/2] implemented, not tested! --- Makefile | 4 ++-- src/kappa.c | 2 +- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/Makefile b/Makefile index fa7e3ff..4a1225c 100644 --- a/Makefile +++ b/Makefile @@ -33,8 +33,8 @@ USE_MPI = yes #GSL #GSL_INC = -I/add/path #GSL_LIB = -L/add/path -GSL_INC = -I/home/alonso/include -GSL_LIB = -L/home/alonso/lib +GSL_INC = -I/home/damonge/include +GSL_LIB = -L/home/damonge/lib #FFTW FFTW_INC = FFTW_LIB = diff --git a/src/kappa.c b/src/kappa.c index bb961fd..9527710 100644 --- a/src/kappa.c +++ b/src/kappa.c @@ -98,7 +98,7 @@ void kappa_get_beam_properties(ParamCoLoRe *par) double *inv_r_max=my_malloc(kmap->nr*sizeof(double)); for(i_r=0;i_rnr;i_r++) { int i_r_here=(int)(kmap->rf[i_r]*idr+0.5); - inv_r_max[i_r]=1./(i_r_here*dr); + inv_r_max[i_r]=1./kmap->rf[i_r]; i_r_max_arr[i_r]=MIN(i_r_here,nr-1); } i_r_min_arr[0]=0; From a2e4e07dd4d786364fbb645f4a4879cd47ee4880 Mon Sep 17 00:00:00 2001 From: damonge Date: Fri, 10 Apr 2026 10:44:10 +0100 Subject: [PATCH 2/2] tested --- param_example.cfg | 4 ++-- src/cosmo.c | 6 ++++-- 2 files changed, 6 insertions(+), 4 deletions(-) diff --git a/param_example.cfg b/param_example.cfg index 9935fc5..31dfe05 100644 --- a/param_example.cfg +++ b/param_example.cfg @@ -122,7 +122,7 @@ imap1: kappa: { #Redshifts at which maps should be generated - z_out= [0.4] + z_out= [1100.] nside= 64 } @@ -130,7 +130,7 @@ kappa: #of the ISW effect (time derivative of the gravitational potential). isw: { - z_out= [0.4] + z_out= [1100.] nside= 64 } diff --git a/src/cosmo.c b/src/cosmo.c index 2211849..63b9d8a 100644 --- a/src/cosmo.c +++ b/src/cosmo.c @@ -803,7 +803,8 @@ void cosmo_set(ParamCoLoRe *par) fprintf(par->f_dbg,"\n"); #endif //_DEBUG #else //_ADD_EXTRA_KAPPA - report_error(1,"Source plane %d outside redshift range\n",ii+1); + print_info("Lensing source plane %d outside redshift range. Will only fill up to z=%.3lf\n", + ii+1,par->z_max); #endif //_ADD_EXTRA_KAPPA } } @@ -843,7 +844,8 @@ void cosmo_set(ParamCoLoRe *par) fprintf(par->f_dbg,"\n"); #endif //_DEBUG #else //_ADD_EXTRA_KAPPA - report_error(1,"Source plane %d outside redshift range\n",ii+1); + print_info("ISW source plane %d outside redshift range. Will only fill up to z=%.3lf\n", + ii+1,par->z_max); #endif //_ADD_EXTRA_KAPPA } }