Skip to content
Snippets Groups Projects
Commit e90c22b2 authored by Sten Delos's avatar Sten Delos
Browse files

minor fix; no need for growth stats to account for radiation because they are...

minor fix; no need for growth stats to account for radiation because they are still inaccurate in that context
parent fd9ca4ab
No related branches found
No related tags found
No related merge requests found
......@@ -119,12 +119,8 @@ class ngenic : public pm_mpi_fft
static double ngenic_growth_int(double a, void *param)
{
#if defined(RADIATION) && !defined(SMOOTHMATTER)
return pow(a / (All.Omega0 + (1 - All.Omega0 - All.OmegaLambda) * a + All.OmegaLambda * a * a * a + All.OmegaR / a), 1.5);
#elif !defined(RADIATION) && defined(SMOOTHMATTER)
#ifdef SMOOTHMATTER
return pow(a / (All.Omega0 + All.OmegaSmooth + (1 - All.Omega0 - All.OmegaLambda - All.OmegaSmooth) * a + All.OmegaLambda * a * a * a), 1.5);
#elif defined(RADIATION) && defined(SMOOTHMATTER)
return pow(a / (All.Omega0 + All.OmegaSmooth + (1 - All.Omega0 - All.OmegaLambda - All.OmegaSmooth - All.OmegaR) * a + All.OmegaLambda * a * a * a + All.OmegaR / a), 1.5);
#else
return pow(a / (All.Omega0 + (1 - All.Omega0 - All.OmegaLambda) * a + All.OmegaLambda * a * a * a), 1.5);
#endif
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment