Allow value for mui that is different from PiT?
For the upcoming project, where we want to simulate \mu_i
values that are different from \pi\cdot T
, I adjusted the local version of BaHaMAS to allow such values. The value of mui would then be given in the path like for the other parameters. The value of \mu_i
is given as a fraction of of \pi\cdot T
and in the path it would look like .../Nf5.0/muiPiT2700/mass...
, where the value for \mu_i
is then \mu_i=0.27\cdot \pi\cdot T
.
The implementation was quite straight-forward in the file ProduceInputFile.bash, where \pi\cdot T
fractions are allowed that are smaller 1. git diff
on my local adjusted BaHaMAS gives:
diff --git a/LQCD_Software_Dependent_Code/CL2QCD/ProduceInputFile.bash b/LQCD_Software_Dependent_Code/CL2QCD/ProduceInputFile.bash
index bc4d78c..5002a99 100644
--- a/LQCD_Software_Dependent_Code/CL2QCD/ProduceInputFile.bash
+++ b/LQCD_Software_Dependent_Code/CL2QCD/ProduceInputFile.bash
@@ -79,6 +79,10 @@ function ProduceInputFile_CL2QCD()
__static__AddToInputFile "useChemicalPotentialIm=1"
if [[ ${BHMAS_chempot} = "PiT" ]]; then
__static__AddToInputFile "chemicalPotentialIm=$(awk -v ntime="${BHMAS_ntime}" 'BEGIN{printf "%.15f\n", atan2(0, -1)/ntime}')"
+ elif [[ ${BHMAS_chempot} =~ ^PiT[0-9]{4}$ ]]; then
+ local PiT_fraction
+ PiT_fraction="${BHMAS_chempot#PiT}"
+ __static__AddToInputFile "chemicalPotentialIm=$(awk -v ntime="${BHMAS_ntime}" -v f=0."${PiT_fraction}" 'BEGIN{printf "%.15f\n", atan2(0, -1)/ntime*f}')"
else
Fatal ${BHMAS_fatalValueError} "Unknown value " emph "${BHMAS_chempot}" " of imaginary chemical potential for input file!"
fi
So my question is, Alessandro, would you adopt this change for BaHaMAS? I think that Alfredo also needs to simulate at \mu_i
different from 0 and \pi\cdot T
.