From e774c063a46f8ffb4d05c41d65afbbcacad3f207 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jean-S=C3=A9bastien?= <git@jscaux.org>
Date: Wed, 15 Dec 2021 08:33:24 +0100
Subject: [PATCH] Tweak script parameters

---
 scripts/catalogue/c_fixed_N_prog_dsfs.sh  | 2 +-
 scripts/catalogue/c_fixed_N_prog_plots.sh | 4 ----
 scripts/catalogue/c_scan_N_fixed_data.sh  | 2 +-
 scripts/catalogue/c_scan_N_fixed_dsfs.sh  | 4 ++--
 scripts/catalogue/c_scan_N_fixed_plot.sh  | 6 ------
 5 files changed, 4 insertions(+), 14 deletions(-)

diff --git a/scripts/catalogue/c_fixed_N_prog_dsfs.sh b/scripts/catalogue/c_fixed_N_prog_dsfs.sh
index 3aa8f99..9e164bc 100755
--- a/scripts/catalogue/c_fixed_N_prog_dsfs.sh
+++ b/scripts/catalogue/c_fixed_N_prog_dsfs.sh
@@ -53,7 +53,7 @@ clist=(0.015625 0.03125 0.0625 0.125 0.25 0.5 1 2 4 8 16 32 64 128 256 512 1024)
 basedir="$(pwd)"
 
 ommin=0
-ommax=$(($nkmax*($nkmax - 8))) # use lower branch of 1ph continuum
+ommax=$((0.5 * $nkmax*($nkmax + 8))) # use (half of) TG Type 1 dispersion
 Nom=1000
 
 for c in $clist
diff --git a/scripts/catalogue/c_fixed_N_prog_plots.sh b/scripts/catalogue/c_fixed_N_prog_plots.sh
index 7d560f9..3fd89b4 100755
--- a/scripts/catalogue/c_fixed_N_prog_plots.sh
+++ b/scripts/catalogue/c_fixed_N_prog_plots.sh
@@ -52,10 +52,6 @@ clist=(0.015625 0.03125 0.0625 0.125 0.25 0.5 1 2 4 8 16 32 64 128 256 512 1024)
 
 basedir="$(pwd)"
 
-ommin=0
-ommax=$(($nkmax*($nkmax - 8))) # use lower branch of 1ph continuum
-Nom=1000
-
 for c in $clist
 do
     echo 'Starting plotting of DSFs for c =' $c
diff --git a/scripts/catalogue/c_scan_N_fixed_data.sh b/scripts/catalogue/c_scan_N_fixed_data.sh
index cb3233f..9e26418 100755
--- a/scripts/catalogue/c_scan_N_fixed_data.sh
+++ b/scripts/catalogue/c_scan_N_fixed_data.sh
@@ -64,7 +64,7 @@ refine=0
 clist_raw=()
 for nc in {-160..160}
 do
-    clist_raw=($clist_raw $(( 2**($nc/16.) )))
+    clist_raw=($clist_raw $(( 4 * 2**($nc/16.) )))
 done
 # Now cast the integer values to true integers
 zmodload zsh/mathfunc
diff --git a/scripts/catalogue/c_scan_N_fixed_dsfs.sh b/scripts/catalogue/c_scan_N_fixed_dsfs.sh
index c5c3201..417064a 100755
--- a/scripts/catalogue/c_scan_N_fixed_dsfs.sh
+++ b/scripts/catalogue/c_scan_N_fixed_dsfs.sh
@@ -57,7 +57,7 @@ fi
 iKmax=$(($nkmax * $N/8))
 
 ommin=0
-ommax=$(($nkmax*$nkmax)) # use mid-point between upper and lower branch of 1ph continuum
+ommax=$((0.5 * $nkmax*($nkmax + 8))) # use (half of) TG Type 1 dispersion
 Nom=1000
 
 # store location for these runs
@@ -67,7 +67,7 @@ dirstore='Lieb-Liniger/'$correlator'/T_'$kBT'/c_scan_N_fixed/k_0_to_'${(l:2::0:)
 clist_raw=()
 for nc in {-160..160}
 do
-    clist_raw=($clist_raw $(( 2**($nc/16.) )))
+    clist_raw=($clist_raw $((4 * 2**($nc/16.) )))
 done
 # Now cast the integer values to true integers
 zmodload zsh/mathfunc
diff --git a/scripts/catalogue/c_scan_N_fixed_plot.sh b/scripts/catalogue/c_scan_N_fixed_plot.sh
index b42f962..c4074fd 100755
--- a/scripts/catalogue/c_scan_N_fixed_plot.sh
+++ b/scripts/catalogue/c_scan_N_fixed_plot.sh
@@ -54,12 +54,6 @@ elif [[ $whichDSF == 'g' ]]; then
     correlator='psi-psidag'
 fi
 
-iKmax=$(($nkmax * $N/8))
-
-ommin=0
-ommax=$(($nkmax*$nkmax)) # use mid-point between upper and lower branch of 1ph continuum
-Nom=1000
-
 # store location for these runs
 dirstore='Lieb-Liniger/'$correlator'/T_'$kBT'/c_scan_N_fixed/k_0_to_'${(l:2::0:)nkmax}'kFo4/sr_'$target_sumrule'/store'