Montepython issues
GCspectro likelihood
Run technicalities in [[MontePython Oneloop Runs]]
- Installation issues:
- Class Oneloop C-code not compatible with older
gcc 9.4, it is compatible withgcc~>11.0 - Classy does not work with
numpy >=2.0. - Montepython requires
numpy>1.23.5because of ascipy.optimize import fsolve
- Class Oneloop C-code not compatible with older
- w0, wa priors are not important in (Analytical Marginalized) AM likelihoods.
- w0, wa priors are important in full likelihoods, getting constantly errors like:
desinomarg.err:48:Error in Class: perturbations_init(L:796) :condition (w_fld_ini >= 0.) is true; The fluid is meant to be negligible at early time, and unimportant for defining the initial conditions of other species. You are using parameters for which this assumption may break down, since at early times you have w_fld(a--->0) = 1.406016e-03 >= 0
- AM likelihoods do not generate their own fiducials, one has to generate them using the full likelihoods and be very careful with fiducial parameters.
- MP likelihoods when reaching max number of Nsteps, finish (on the cluster at least) with
FAILED, ExitCode 137giving the wrong impression that the run has failed or there was a problem. They should exit withExitCode 0 - When generating the fiducial, the MP likelihood returns :
256 warnings.warn(
257 "Writing fiducial model in %s, for %s likelihood\n" % (
258 self.fid_filename, self.name))
259 return 1j
This is a strange complex number return, because in the case of using the option --display-each-chi2 at the same time as -f0 to generate the fiducial, this return value gets added to the other chi2 of other likelihoods and one gets the wrong impression that something is going on. The logic here should be clarified.
- The public (Casas et al 2023) ISTF-like
euclid_spectroscopiclikelihood does not work form_nu=0because it automatically requestsP_cb,f_cband othercbquantities, which are not available in this case from class. - The oneloop spectroscopic likelihoods do not contain
P_cbexplicitly, working withP_mmonly. - The runs of the MP validation paper were done in the
Omega_m, sigma8and camb basis. The oneloop runs have been done using the class basis andlogAsmaking a comparison more difficult.