CosmoMC error adding the scalar field cosmological model's calculation

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Koichi Hirano
Posts: 2
Joined: January 24 2007
Affiliation: Tsuru University

CosmoMC error adding the scalar field cosmological model'

Post by Koichi Hirano » January 26 2007

Hello,

I added the calculation of the scalar field cosmological model in subroutine INIT_BACKGROUND of CAMB.

I started CosmoMC with this CAMB,
then, This running without probrem untill the ten times calling CAMB,
that time, I wrote the value of Params%P in the function GetLogLike on DISPLAY
--------------------------------------------------------------------------------------------------------------------------------
Params%P= 2.7000001E-02 8.9349999E-02 1.040000 7.9999998E-02
0.0000000E+00 0.0000000E+00 -1.000000 1.000000 0.0000000E+00
-0.2800000 3.200000 0.0000000E+00 0.0000000E+00
--------------------------------------------------------------------------------------------------------------------------------

but, about ten times calling CAMB, CosmoMC stop.
That time, I wrote the value of Params%P in the function GetLogLike on DISPLAY
--------------------------------------------------------------------------------------------------------------------------------
Params%P= NaN 8.8799998E-02 1.040000 7.9999998E-02
0.0000000E+00 0.0000000E+00 -1.000000 1.000000 0.0000000E+00
-0.2029928 3.034900 0.0000000E+00 0.0000000E+00
--------------------------------------------------------------------------------------------------------------------------------
Probrem is that the Value of Baryon density is NaN.

I not modify the CosmoMC and CAMB program except adding model calculation in subroutine INIT_BACKGROUND of CAMB.

When I delete this model calculation line, CosmoMC run without probrem.

In this cosmological model calculation, I use large number of array.

Is there someone who faced the same probrem?

(I use intel ifort compiler. dell's Desktop)

Pascal Vaudrevange
Posts: 50
Joined: March 26 2006
Affiliation: DESY
Contact:

CosmoMC error adding the scalar field cosmological model'

Post by Pascal Vaudrevange » January 27 2007

Maybe you wrote outside the bounds of an array in the computation that you added?

Try compiling with compiler option '-CB'.

Pascal

Post Reply