diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-30 09:33:16 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-30 09:33:16 -0300 |
| commit | d422676a24bc3967ab3f319d60f7798081dab4e5 (patch) | |
| tree | d181cce9e9b0f19a42e0f66bed0109eadd8967ec | |
| parent | e108c33ab66d71c53fa7e9d81b95f6a00f61a114 (diff) | |
| download | code-d422676a24bc3967ab3f319d60f7798081dab4e5.tar.gz code-d422676a24bc3967ab3f319d60f7798081dab4e5.tar.bz2 code-d422676a24bc3967ab3f319d60f7798081dab4e5.zip | |
Made the exponent for the μ adjustment adjustable
| -rw-r--r-- | log-fourier_integrator.cpp | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp index 4cd18ee..9bbf554 100644 --- a/log-fourier_integrator.cpp +++ b/log-fourier_integrator.cpp @@ -17,6 +17,7 @@ int main(int argc, char* argv[]) { /* Iteration parameters */ Real ε = 1e-14; Real γ₀ = 1; + Real x = 0.5; Real β₀ = 0; Real βₘₐₓ = 0.7; Real Δβ = 0.01; @@ -25,7 +26,7 @@ int main(int argc, char* argv[]) { int opt; - while ((opt = getopt(argc, argv, "p:s:2:T:t:b:d:g:k:D:e:0:lS:")) != -1) { + while ((opt = getopt(argc, argv, "p:s:2:T:t:b:d:g:k:D:e:0:lS:x:")) != -1) { switch (opt) { case 'p': p = atoi(optarg); @@ -60,6 +61,9 @@ int main(int argc, char* argv[]) { case '0': β₀ = atof(optarg); break; + case 'x': + x = atof(optarg); + break; case 'l': loadData = true; break; @@ -132,7 +136,7 @@ int main(int argc, char* argv[]) { std::vector<Real> Rₜ₊₁ = fft.inverse(Ȓₜ₊₁); std::vector<Real> Cₜ₊₁ = fft.inverse(Ĉₜ₊₁); - μₜ *= pow(tanh(Cₜ₊₁[0]-1)+1, 0.5); + μₜ *= pow(tanh(Cₜ₊₁[0]-1)+1, x); ΔCₜ = 0; for (unsigned i = 0; i < N; i++) { |