Hi,
For testing purposes I am trying to run cosmomc (std version) with Omega_k = 0.1.
Somewhere for a similar problem Anthony suggested :
" Check the hard coded H0 priors? (or give θ a larger starting range) "
I've been focusing on enlargening \theta range.
(I did not try yet to directly modify hard coded priors on H0)
And it still fails to run.
I am setting : use_HST = F ;
param[theta] = 1.04 0.5 10 0.002 0.002
param[omegak] = 0.1 0.1 0.1 0 0
I keep getting :
MCMC.f90:Couldn't start after 1000 tries - check starting ranges
However , for the same inputs , when I set 0 < Omega_k <= 0.6,
this basic cosmomc works fine.
Is there a subtle conflict that I am not aware of ?
Thank you for any hint on how to overcome this,
Patrice
Struggling to run comomc for Omega_k = 0.1 & slightly above
-
- Posts: 19
- Joined: November 05 2009
- Affiliation: University of Cape Town
-
- Posts: 19
- Joined: November 05 2009
- Affiliation: University of Cape Town
Struggling to run comomc for Omega_k = 0.1 & slightly ab
Apparently, the best option is to change hard coded priors on H0 accordingly.
By allowing a larger range of values for H0 + right priors on \theta in params.ini,
everything seems to be running now.
Patrice
By allowing a larger range of values for H0 + right priors on \theta in params.ini,
everything seems to be running now.
Patrice