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

Added the correction from Jennings (2012) for the velocity field

Browse Source
This commit is contained in:
Guilhem Lavaux 2012-12-12 17:41:14 -05:00
parent a3ef55c318
commit 975156d1f3
2 changed files with 16 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/powerSpectrum.cpp
View File

@ -1,3 +1,4 @@
#include <cassert>
#include <vector> #include <vector>
#include <algorithm> #include <algorithm>
#include <iostream> #include <iostream>
@ -651,4 +652,17 @@ double computeCorrel2(double powNorm, double topHatRad1, double topHatRad2)
#endif #endif
} }
double vvCorrection(double P_deltadelta, double k)
{
static const double alpha0 = -12480.5, alpha1 = 1.824, alpha2 = 2165.87, alpha3=1.796;
if (k > 0.3)
return 0;
double r =(alpha0*sqrt(P_deltadelta) + alpha1*P_deltadelta*P_deltadelta)/(alpha2 + alpha3*P_deltadelta);
assert(P_deltadelta > 0);
if (r < 0)
return 0;
return r;
}
}; };

2
src/powerSpectrum.hpp
View File

@ -30,6 +30,8 @@ namespace Cosmology {
double computeVarianceZero(double powNorm); double computeVarianceZero(double powNorm);
double computeCorrel(double powNorm, double topHatRad1); double computeCorrel(double powNorm, double topHatRad1);
double computeCorrel2(double powNorm, double topHatRad1, double topHatRad2); double computeCorrel2(double powNorm, double topHatRad1, double topHatRad2);
double vvCorrection(double P_deltadelta, double k);
}; };
#endif #endif