From 5125a858781fdeea091d9a8189b6cbf5e4758d22 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Jean-S=C3=A9bastien?= Date: Thu, 9 Dec 2021 20:21:04 +0100 Subject: [PATCH] Change momentum factor for catalogue --- src/EXECS/LiebLin_Catalogue_Fixed_c_k_Nscaling.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/EXECS/LiebLin_Catalogue_Fixed_c_k_Nscaling.cc b/src/EXECS/LiebLin_Catalogue_Fixed_c_k_Nscaling.cc index d0701a1..8884d6f 100644 --- a/src/EXECS/LiebLin_Catalogue_Fixed_c_k_Nscaling.cc +++ b/src/EXECS/LiebLin_Catalogue_Fixed_c_k_Nscaling.cc @@ -30,7 +30,7 @@ int main(int argc, char* argv[]) cout << "char whichDSF \t\t Which structure factor should be calculated ? Options are: " "d for rho rho, g for psi psi{dagger}, o for psi{dagger} psi" << endl; cout << "DP c_int \t\t Value of the interaction parameter: use positive real values only" << endl; - cout << "int kfact \t\t momentum factor: momemntum will be set to kfact * kF/16" << endl; + cout << "int kfact \t\t momentum factor: momemntum will be set to kfact * kF/4" << endl; cout << "DP kBT \t\t Temperature (positive only of course)" << endl; cout << "DP target_sumrule \t sumrule saturation you're satisfied with" << endl; cout << "int Hrs \t\t Allowed computational time (hours)" << endl; @@ -67,7 +67,7 @@ int main(int argc, char* argv[]) int N = Nv[iN]; DP L = N; - int iKmin = (kfact * N)/32; + int iKmin = (kfact * N)/8; int iKmax = iKmin; DP srsat = 0.0; bool refine = false;