mirror/MUSIC
1
0
Fork 0
mirror of https://github.com/cosmo-sims/MUSIC.git synced 2024-09-19 17:03:46 +02:00
Code Releases Activity

fixed a wrong sign for quintessence models (undocumented feature still, not affecting public version)

Browse source
This commit is contained in:
Oliver Hahn 2015-06-23 15:42:44 +02:00
parent 7684780175
commit b7e6ed4a3d
1 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
cosmology.hh
View file

@ -66,7 +66,7 @@ public:
Cosmology *cosm = (Cosmology*)Params;
double a2 = a*a;
double Ha = sqrt(cosm->Omega_m/(a2*a) + cosm->Omega_k/a2
+ cosm->Omega_DE * pow(a,3.+3.*cosm->w_0) * exp(-3.*(a-1.0)*cosm->w_a) );
+ cosm->Omega_DE * pow(a,3.+3.*cosm->w_0) * exp(3.*(a-1.0)*cosm->w_a) );
return Ha;
}
@ -75,9 +75,9 @@ public:
Cosmology *cosm = (Cosmology*)Params;
double a2 = a*a;
double H = H_of_a( a, Params );
double Hprime = -1.5 * cosm->Omega_m / (a2*a2*H) - cosm->Omega_k / (a2*a*H)
+ 0.5 * cosm->Omega_DE * pow( a, 3.+3.*cosm->w_0 ) * exp( -3.*(a-1.)*cosm->w_a )
* 1.5 / a * ( 1. + cosm->w_0 - a * cosm->w_a );
double Hprime = 1/(a*H) * ( -1.5 * cosm->Omega_m / (a2*a) - cosm->Omega_k / a2
+ 1.5 * cosm->Omega_DE * pow( a, 3.+3.*cosm->w_0 ) * exp( 3.*(a-1.)*cosm->w_a )
* ( 1. + cosm->w_0 + a * cosm->w_a ) );
return Hprime;
}