diff options
| -rw-r--r-- | log-fourier_integrator.cpp | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index 55f358a..017def7 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -140,7 +140,8 @@ int main(int argc, char* argv[]) { std::vector<Real> Rₜ₊₁ = fft.inverse(Ȓₜ₊₁); std::vector<Real> Cₜ₊₁ = fft.inverse(Ĉₜ₊₁); - μₜ *= pow(tanh(Cₜ₊₁[0]-1)+1, x); + Real C₀ = (Cₜ₊₁[0] + Cₜ₊₁[1]) / 2; + μₜ *= pow(tanh(C₀-1)+1, x); ΔCₜ = 0; for (unsigned i = 0; i < N; i++) { |