mirror of
https://github.com/cosmo-sims/MUSIC.git
synced 2024-09-19 17:03:46 +02:00
Fixed a bug in Eisenstein&Hu fitting functions, affecting power spectra at a few per cents. The TFset_parameters function takes Omega_b/Omega_m, *not* Omega_b/Omega_c as wrongly stated in the Eisenstein&Hu source code.
This commit is contained in:
parent
1b3386d7c7
commit
31d2dbf2f2
1 changed files with 2 additions and 1 deletions
|
@ -162,7 +162,7 @@ struct eisenstein_transfer
|
|||
{
|
||||
m_h0 = cosmo.H0*0.01;
|
||||
TFset_parameters( (cosmo.Omega_m)*cosmo.H0*cosmo.H0*(0.01*0.01),
|
||||
cosmo.Omega_b/(cosmo.Omega_m-cosmo.Omega_b), Tcmb);
|
||||
cosmo.Omega_b/cosmo.Omega_m, Tcmb);
|
||||
|
||||
fb_ = cosmo.Omega_b/(cosmo.Omega_m);
|
||||
fc_ = (cosmo.Omega_m-cosmo.Omega_b)/(cosmo.Omega_m) ;
|
||||
|
@ -299,6 +299,7 @@ public:
|
|||
*pow(1.5/wdmgx_,0.29);
|
||||
break;
|
||||
}
|
||||
std::cerr << "WDM alpha = " << m_WDMalpha << std::endl;
|
||||
}
|
||||
|
||||
inline double compute( double k, tf_type type )
|
||||
|
|
Loading…
Reference in a new issue