summaryrefslogtreecommitdiff
path: root/log-fourier_integrator.cpp
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2025-04-19 17:35:30 -0300
committerJaron Kent-Dobias <jaron@kent-dobias.com>2025-04-19 17:35:30 -0300
commit0c08fa2b2dc82cbbd2e6bb4ceb201f337961a75b (patch)
treedb2f218db6fc30cbc2053b9fa3d4960a68fd836d /log-fourier_integrator.cpp
parentbb76b20bd43b427c5fb100eb431ee5b8487a4613 (diff)
downloadcode-0c08fa2b2dc82cbbd2e6bb4ceb201f337961a75b.tar.gz
code-0c08fa2b2dc82cbbd2e6bb4ceb201f337961a75b.tar.bz2
code-0c08fa2b2dc82cbbd2e6bb4ceb201f337961a75b.zip
More tweaks
Diffstat (limited to 'log-fourier_integrator.cpp')
-rw-r--r--log-fourier_integrator.cpp4
1 files changed, 2 insertions, 2 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index 28dd6e7..3049bd5 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -148,7 +148,7 @@ int main(int argc, char* argv[]) {
ΔC₀ = ΔC;
}
- std::cerr << β << " " << μ << " " << ΔC << " " << γ;
+ std::cerr << β << " " << μ << " " << ΔC << " " << γ << " " << Cₜ[0];
std::cerr << "\r";
}
if (std::isnan(Cₜ[0])) {
@@ -172,7 +172,7 @@ int main(int argc, char* argv[]) {
/* We found an upper bound */
μ₂ = μ;
}
- μ *= Cₜ[0];
+ μ *= sqrt(Cₜ[0]);
} else {
/* Once the bounds are set, we can use bisection */
if (Cₜ[0] > 1) {