From 01f36b9be7d3ab0185b5451fb5c36fc04bb85b95 Mon Sep 17 00:00:00 2001
From: Guilhem Lavaux <guilhem.lavaux@gmail.com>
Date: Fri, 30 Nov 2012 17:11:21 -0500
Subject: [PATCH] Added a get_delta_k function to get the element phase volume

---
 src/fourier/details/euclidian_maps.hpp | 5 +++++
 src/powerSpectrum.cpp                  | 2 +-
 2 files changed, 6 insertions(+), 1 deletion(-)

diff --git a/src/fourier/details/euclidian_maps.hpp b/src/fourier/details/euclidian_maps.hpp
index cfadc72..6a650fe 100644
--- a/src/fourier/details/euclidian_maps.hpp
+++ b/src/fourier/details/euclidian_maps.hpp
@@ -113,6 +113,11 @@ namespace CosmoTool
 		   this->delta_k);
     }
 
+    const std::vector<double>& get_delta_k() const
+    {
+      return this->delta_k;
+    }
+
     template<typename Array, typename Array2>
     void get_Kvec(const Array& ik, Array2& kvec)
     {
diff --git a/src/powerSpectrum.cpp b/src/powerSpectrum.cpp
index 9f3d7c8..6ad7c89 100644
--- a/src/powerSpectrum.cpp
+++ b/src/powerSpectrum.cpp
@@ -21,7 +21,7 @@ using namespace std;
 #define POWER_BDM 7
 #define POWER_TEST 8
 
-#define POWER_SPECTRUM HU_WIGGLES
+#define POWER_SPECTRUM POWER_EFSTATHIOU
 
 namespace Cosmology {