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)
CosmoMC error adding the scalar field cosmological model's calculation
-
- Posts: 2
- Joined: January 24 2007
- Affiliation: Tsuru University
-
- Posts: 50
- Joined: March 26 2006
- Affiliation: DESY
- Contact:
CosmoMC error adding the scalar field cosmological model'
Maybe you wrote outside the bounds of an array in the computation that you added?
Try compiling with compiler option '-CB'.
Pascal
Try compiling with compiler option '-CB'.
Pascal