summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--log-fourier_integrator.cpp4
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]);