As a new user to Cosmomc I have a few very basic questions:
1) I intend to use WMAP5yr + Acbar + SN data. For this, I have increased the lmax to 3000 but running cosmomc still gives the warning:
Should I still increase the lmax?Warning: Window function non-zero outside l_max: data/windows/acbar2007_v3_corr_21 assuming contribution is negligible.
2) My first model should use seven parameters, which are already defined in params.ini and elsewhere, so I only need to reduce the nbr of already defined parameters accordingly. Should I do this by setting all the values (max, min etc) of the excess parameters to "0" in the params.ini or comment them away, i.e.:
Code: Select all
param[w] = 0 0 0 0 0
Code: Select all
#param[w] = -1 -1 -1 0 0
Code: Select all
cov_matrix_dimension = 7
3) If I change the pivot scale in params.ini to e.g. 0.01, do I need to change it elsewhere also, e.g. in cmbtypes.f90, where the pivot scale also seems to be given?
4) After re-compiling with my changes, the program starts to run but stops pretty soon with a segmentation fault, entire output:
_____________________/cosmomc$ ./cosmomc params.ini
Random seeds: 16082, 6581 rand_inst: 0
Reading BCP data set: data/acbar2007_v3_corr.newdat
Reading: acbar2007_v3_corr_
Using selected band ranges
Warning: Window function non-zero outside l_max: data/windows/acbar2007_v3_corr_21
assuming contribution is negligible
Warning: Window function non-zero outside l_max: data/windows/acbar2007_v3_corr_22
assuming contribution is negligible
Warning: Window function non-zero outside l_max: data/windows/acbar2007_v3_corr_23
assuming contribution is negligible
Warning: Window function non-zero outside l_max: data/windows/acbar2007_v3_corr_24
assuming contribution is negligible
Warning: Window function non-zero outside l_max: data/windows/acbar2007_v3_corr_25
assuming contribution is negligible
Computing tensors: F
Doing CMB lensing: T
lmax = 3000
Number of C_ls = 3
Varying 7 parameters ( 0 fast)
reading: params_CMB.covmat
starting Monte-Carlo
Reading: supernovae data (SDSS).
Reading: data/MLCS.FITRES
Using 288 SNe from subset SDSS + ESSENCE + SNLS + LOWZ + HST.
reading WMAP5 data
Initializing WMAP likelihood, version v3
Segmentation fault
Where can this segmentation fault come from? I have only changed lines in the files params.ini, distparams.ini and cmbtypes.f90. In addition to the abovementioned changes, I have set:
Without my changes described above, cosmomc seems to compile and run just fine.use_HST = F
...
use_Age_Tophat_Prior = F
use_SN = T
...
CMB_lensing = T
I will be very grateful if somebody can take the time to guide me a bit forward.
Br,
Marianne