Change to convention normalizing Γ₀ by τ₀

1 files changed, 4 insertions, 2 deletions

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index fb8664b..51f5326 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -61,11 +61,13 @@ int main(int argc, char* argv[]) {
 
   LogarithmicFourierTransform fft(N, k, Δτ, 4);
 
-  Real Γ₀ = 1.0;
-  Real μ = Γ₀;
+  Real Γ₀ = 1.0 + τ₀;
+  Real μ = 1.0;
+/*  Real μ = Γ₀;
   if (τ₀ > 0) {
     μ = (sqrt(1+4*Γ₀*τ₀)-1)/(2*τ₀);
   }
+*/
 
   std::vector<Real> Cₜ₋₁(N);
   std::vector<Real> Rₜ₋₁(N);