diff options
-rw-r--r-- | log-fourier_integrator.cpp | 4 |
1 files changed, 1 insertions, 3 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index 5ff27a4..f02aebc 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -156,7 +156,7 @@ int main(int argc, char* argv[]) { if (μ₋ > 0 && μ₊ > 0) { μₜ = (μ₊ + μ₋) / 2; } else { - μₜ *= C₀; + μₜ *= pow(tanh(C₀-1)+1, x); } } @@ -212,8 +212,6 @@ int main(int argc, char* argv[]) { } } - μₜ *= pow(tanh(C₀-1)+1, x); - ΔCₜ = 0; for (unsigned i = 0; i < N; i++) { ΔCₜ += std::norm(Cₜ[i] - Cₜ₊₁[i]); |