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

Fixed calls to get_K_p

Browse Source
This commit is contained in:
Guilhem Lavaux 2014-02-27 16:03:42 +01:00
parent aed2920be7
commit c724700246
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/fourier/details/euclidian_spectrum_1d.hpp
View File

@ -115,7 +115,7 @@ namespace CosmoTool
for (long p = 1; p < rand_map.size(); p++)
{
double A_k = std::sqrt(0.5*V*f(rand_map.get_K(p)));
double A_k = std::sqrt(0.5*V*f(rand_map.get_K_p(p)));
d[p] = std::complex<T>(gsl_ran_gaussian(rng, A_k),
gsl_ran_gaussian(rng, A_k));
}
@ -164,7 +164,7 @@ namespace CosmoTool
std::complex<T> *d = m.data();
for (long p = 0; p < m_c.size(); p++)
d[p] *= f(m_c.get_K(p));
d[p] *= f(m_c.get_K_p(p));
}
template<typename T>
@ -174,7 +174,7 @@ namespace CosmoTool
std::complex<T> *d = m.data();
for (long p = 0; p < m_c.size(); p++)
d[p] *= std::sqrt(f(m_c.get_K(p)));
d[p] *= std::sqrt(f(m_c.get_K_p(p)));
}
template<typename T>
@ -185,7 +185,7 @@ namespace CosmoTool
for (long p = 0; p < m_c.size(); p++)
{
T A = f(m_c.get_K(p));
T A = f(m_c.get_K_p(p));
if (A==0)
d[p] = 0;
else
@ -201,7 +201,7 @@ namespace CosmoTool
for (long p = 0; p < m_c.size(); p++)
{
T A = std::sqrt(f(m_c.get_K(p)));
T A = std::sqrt(f(m_c.get_K_p(p)));
if (A == 0)
d[p] = 0;
else