diff options
Diffstat (limited to 'log-fourier_integrator.cpp')
| -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++) { |