guilhem_lavaux/cosmotool
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Moved s_tilde

Browse Source
This commit is contained in:
Guilhem Lavaux 2015-02-24 14:45:09 +01:00
parent c39f89d030
commit 535f25d710
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
src/cosmopower.cpp
View File

@ -172,9 +172,10 @@ double CosmoPower::powerHuWiggles(double k)
double beta_c = ehu.beta_c; double beta_c = ehu.beta_c;
double alpha_b = ehu.alpha_b; double alpha_b = ehu.alpha_b;
double beta_b = ehu.beta_b; double beta_b = ehu.beta_b;
double s_tilde = s * pow(1 + pow(ehu.beta_node / (k * s), 3), -1./3);
double xx = k * s; double xx = k * s;
double s_tilde = s * pow(1 + pow(ehu.beta_node / (xx), 3), -1./3);
double f = 1 / (1 + pow(xx / 5.4, 4)); double f = 1 / (1 + pow(xx / 5.4, 4));
double q = k / (13.41 * k_eq); double q = k / (13.41 * k_eq);
double T_c = f * T_tilde_0(q, 1, beta_c) + (1 - f) * T_tilde_0(q, alpha_c, beta_c); double T_c = f * T_tilde_0(q, 1, beta_c) + (1 - f) * T_tilde_0(q, alpha_c, beta_c);