From c72470024668c4defdc8701231e8f388f79887ef Mon Sep 17 00:00:00 2001
From: Guilhem Lavaux <lavaux@iap.fr>
Date: Thu, 27 Feb 2014 16:03:42 +0100
Subject: [PATCH] Fixed calls to get_K_p

---
 src/fourier/details/euclidian_spectrum_1d.hpp | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/src/fourier/details/euclidian_spectrum_1d.hpp b/src/fourier/details/euclidian_spectrum_1d.hpp
index 36b524a..4d6fb9b 100644
--- a/src/fourier/details/euclidian_spectrum_1d.hpp
+++ b/src/fourier/details/euclidian_spectrum_1d.hpp
@@ -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