Skip to content

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.

Edited by Alessandro Sciarra
To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information