CosmoCoffee

[Rajeev Kumar Jain]

Hi,
I've modified cosmomc to find out the best fit values of the parameters of the inflationary
scenario I am considerning. The modified camb produces the angular power spectrum
without any problem but the cosmomc run gives the following error message after some
time and then stops. I am running cosmomc on 4 processors of the same node.

forrtl: error (75): floating point exception
0: __FINI_00_remove_gp_range [0x3ff81a21488]
1: __FINI_00_remove_gp_range [0x3ff81a2a4d0]
2: __FINI_00_remove_gp_range [0x3ff800d9cc0]
3: wmap_likelihood₅yrwmap_likelihood_compute_ [WMAP_5yr_likelihood.F90: 317,
0x12004c470]
4: cmbdatamaplnlike_ [cmbdata.F90: 959, 0x1200b2c30]
5: cmbdatacmblnlike_ [cmbdata.F90: 918, 0x1200b1e80]
6: calclikegetloglikepost_ [calclike.f90: 121, 0x1200d6720]
7: calclikegetloglike_ [calclike.f90: 84, 0x1200d63b8]
8: montecarlomcmcsample_ [MCMC.f90: 535, 0x1200e9ab8]
9: solvecosmology_ [driver.F90: 320, 0x1200ed820]
10: main [for_main.c: 203, 0x1200edacc]
11: _{_start} [0x120037508]
MPI process 3583599 died from signal 6 (IOT/Abort trap)
MPI process 3583598 died from signal 9 (Killed)
MPI process 3583587 died from signal 9 (Killed)

I have no idea what this problem is due to. Please advise.

Regards,
Rajeev Jain

[Pascal Vaudrevange]

Did you check the C_l's that are coming out of CAMB for your model? Maybe something funky is going on, i.e. some C_l^TT<0 because of interpolation issues or something?