Quantcast
Channel: CosmoCoffee
Viewing all articles
Browse latest Browse all 171

Computers and software • Cobaya error about CamSpec

$
0
0
Hi,
I'm trying to run Cobaya with CAMB using the CamSpec data.


At first, I installed Cobaya and theory code as:

Code:

python -m pip install cobaya --upgradecobaya-install cosmo -p /my/path/to/packages
Next, I ran the analysis without Camspec and using the following data:

Code:

likelihood:  H0.riess2020: null  bao.desi_2024_bao_all: null
and I successfully obtained the minimized points and the MCMC chains.

However, when I include Camspec,

Code:

likelihood:  H0.riess2020: null  bao.desi_2024_bao_all: null  planck_2018_lowl.TT: null  planck_2018_lowl.EE: null  planck_NPIPE_highl_CamSpec.TTTEEE: null  planckpr4lensing:    package_install:      github_repository: carronj/planck_PR4_lensing      min_version: 1.0.2
we obtained the following error:

Code:

[root@hpc02 work]# cobaya-run -M CosmoTest.yaml[output] Output to be read-from/written-into folder 'chains', with prefix 'CosmoTest1'[camb] `camb` module loaded successfully from /root/work/packages/code/CAMB/camb[bao.desi_2024_bao_all] Initialized.[planck_npipe_highl_camspec.ttteee] L-range for 143x143: 30 2000[planck_npipe_highl_camspec.ttteee] L-range for 217x217: 500 2500[planck_npipe_highl_camspec.ttteee] L-range for 143x217: 500 2500[planck_npipe_highl_camspec.ttteee] L-range for TE: 30 2000[planck_npipe_highl_camspec.ttteee] L-range for EE: 30 2000[planck_npipe_highl_camspec.ttteee] Number of data points: 9915[minimize] Initializing[minimize] Run 1/2 will start from random initial point:[minimize] {'logA': 3.049692689364472, 'ns': 0.9620239278781086, 'theta_MC_100': 1.0407515608913094, 'ombh2': 0.022275311383953073, 'omch2': 0.12081426211442803, 'tau': 0.04487196507983622, 'A_planck': 0.9994165227203023, 'amp_143': 8.010976957793789, 'amp_217': 19.316919371486325, 'amp_143x217': 10.208012448402823, 'n_143': 0.9686371541402705, 'n_217': 0.7650130660513249, 'n_143x217': 0.8268599187620077, 'calTE': 0.9901009524048164, 'calEE': 1.0142806131456061}[minimize] Run 2/2 will start from random initial point:[minimize] {'logA': 3.0508032390094764, 'ns': 0.9594274212556125, 'theta_MC_100': 1.0405115430120404, 'ombh2': 0.02248236826803105, 'omch2': 0.11946044683479742, 'tau': 0.05227357491133455, 'A_planck': 0.9960570237424191, 'amp_143': 11.01004234862195, 'amp_217': 19.57638904791919, 'amp_143x217': 12.164165697936852, 'n_143': 1.0039038929391748, 'n_217': 1.0729755026155867, 'n_143x217': 0.984968304353735, 'calTE': 1.017725075267489, 'calEE': 0.9943596793081703}[prior] *WARNING* There are unbounded parameters (['calTE', 'A_planck', 'calEE']). Prior bounds are given at 0.9999995 confidence level. Beware of likelihood modes at the edge of the prior[minimize] Starting run 1/2[minimize] Run 1/2 converged.[minimize] Starting run 2/2[minimize] Run 2/2 converged.[minimize] Finished successfully![minimize] -log(posterior) minimized to 222.226[minimize] *ERROR* Cannot reproduce log minimum to within 0.01. Maybe your likelihood is stochastic or large numerical error? Recomputed min: nan (was -222.226) at array([ 3.06194569,  1.00648616,  1.04562493,  0.02762708,  0.12227564,        0.05471409,  0.99981794, 16.4099343 , 17.06925873, 12.02435367,        0.34590346,  2.0576749 ,  0.26125204,  0.9996899 ,  0.99974252])[exception handler] ---------------------------------------Traceback (most recent call last):  File "/usr/local/bin/cobaya-run", line 8, in <module>    sys.exit(run_script())  File "/usr/local/lib/python3.9/site-packages/cobaya/run.py", line 198, in run_script    run(info, **arguments.__dict__)  File "/usr/local/lib/python3.9/site-packages/cobaya/run.py", line 146, in run    sampler.run()  File "/usr/local/lib/python3.9/site-packages/cobaya/samplers/minimize/minimize.py", line 342, in run    self.process_results(*mpi.zip_gather(  File "/usr/local/lib/python3.9/site-packages/cobaya/mpi.py", line 274, in wrapper    result = method(self, *args, **kwargs)  File "/usr/local/lib/python3.9/site-packages/cobaya/samplers/minimize/minimize.py", line 389, in process_results    raise LoggedError(cobaya.log.LoggedError: Cannot reproduce log minimum to within 0.01. Maybe your likelihood is stochastic or large numerical error? Recomputed min: nan (was -222.226) at array([ 3.06194569,  1.00648616,  1.04562493,  0.02762708,  0.12227564,        0.05471409,  0.99981794, 16.4099343 , 17.06925873, 12.02435367,        0.34590346,  2.0576749 ,  0.26125204,  0.9996899 ,  0.99974252])-------------------------------------------------------------
Could you tell me why I got an error and how to fix it.
I'm using
  • almalinux:8
    Python 3.9.13
Thank you for your attention and future advice.

Statistics: Posted by Yo Toda — February 13 2025



Viewing all articles
Browse latest Browse all 171

Trending Articles