From 975156d1f3fa15c2c4139ae7281773499e6de94e Mon Sep 17 00:00:00 2001
From: Guilhem Lavaux <guilhem.lavaux@laposte.net>
Date: Wed, 12 Dec 2012 17:41:14 -0500
Subject: [PATCH] Added the correction from Jennings (2012) for the velocity
 field

---
 src/powerSpectrum.cpp | 14 ++++++++++++++
 src/powerSpectrum.hpp |  2 ++
 2 files changed, 16 insertions(+)

diff --git a/src/powerSpectrum.cpp b/src/powerSpectrum.cpp
index 6ad7c89..02549d6 100644
--- a/src/powerSpectrum.cpp
+++ b/src/powerSpectrum.cpp
@@ -1,3 +1,4 @@
+#include <cassert>
 #include <vector>
 #include <algorithm>
 #include <iostream>
@@ -651,4 +652,17 @@ double computeCorrel2(double powNorm,  double topHatRad1, double topHatRad2)
 #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;
+  }
+
 };
diff --git a/src/powerSpectrum.hpp b/src/powerSpectrum.hpp
index ad0b0a9..de7ce53 100644
--- a/src/powerSpectrum.hpp
+++ b/src/powerSpectrum.hpp
@@ -30,6 +30,8 @@ namespace Cosmology {
   double computeVarianceZero(double powNorm);
   double computeCorrel(double powNorm,  double topHatRad1);
   double computeCorrel2(double powNorm,  double topHatRad1, double topHatRad2);
+
+  double vvCorrection(double P_deltadelta, double k);
 };
 
 #endif