From d422676a24bc3967ab3f319d60f7798081dab4e5 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 30 Apr 2025 09:33:16 -0300
Subject: Made the exponent for the μ adjustment adjustable
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 log-fourier_integrator.cpp | 8 ++++++--
 1 file 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++) {
-- 
cgit v1.2.3-70-g09d2