serializer->doInt(reinterpret_cast<int*>(&fepvals->softcoreFunction));
serializer->doReal(&fepvals->scGapsysScaleLinpointLJ);
serializer->doReal(&fepvals->scGapsysScaleLinpointQ);
- serializer->doReal(&fepvals->scSigmaLJGapsys);
+ serializer->doReal(&fepvals->scGapsysSigmaLJ);
}
else
{
fepvals->softcoreFunction = SoftcoreType::Beutler;
fepvals->scGapsysScaleLinpointLJ = 0.85;
fepvals->scGapsysScaleLinpointQ = 0.3;
- fepvals->scSigmaLJGapsys = 0.3;
+ fepvals->scGapsysSigmaLJ = 0.3;
}
/* handle lambda_neighbors */
const real gmx_unused gapsysScaleLinpointCoul = scParams.gapsysScaleLinpointCoul;
const real gmx_unused gapsysScaleLinpointVdW = scParams.gapsysScaleLinpointVdW;
- const real gmx_unused sigma6VdWGapsys = scParams.sigma6VdWGapsys;
+ const real gmx_unused gapsysSigma6VdW = scParams.gapsysSigma6VdW;
const bool gmx_unused doShiftForces = ((flags & GMX_NONBONDED_DO_SHIFTFORCE) != 0);
const bool doPotential = ((flags & GMX_NONBONDED_DO_POTENTIAL) != 0);
alignas(GMX_SIMD_ALIGNMENT)
real gmx_unused preloadGapsysScaleLinpointCoul[DataTypes::simdRealWidth];
alignas(GMX_SIMD_ALIGNMENT)
- real gmx_unused preloadSigma6VdWGapsys[NSTATES][DataTypes::simdRealWidth];
+ real gmx_unused preloadGapsysSigma6VdW[NSTATES][DataTypes::simdRealWidth];
alignas(GMX_SIMD_ALIGNMENT) real preloadLjPmeC6Grid[NSTATES][DataTypes::simdRealWidth];
#else
real preloadPairIsValid[DataTypes::simdRealWidth];
real gmx_unused preloadAlphaCoulEff[DataTypes::simdRealWidth];
real gmx_unused preloadGapsysScaleLinpointVdW[DataTypes::simdRealWidth];
real gmx_unused preloadGapsysScaleLinpointCoul[DataTypes::simdRealWidth];
- real gmx_unused preloadSigma6VdWGapsys[NSTATES][DataTypes::simdRealWidth];
+ real gmx_unused preloadGapsysSigma6VdW[NSTATES][DataTypes::simdRealWidth];
real preloadLjPmeC6Grid[NSTATES][DataTypes::simdRealWidth];
#endif
for (int s = 0; s < DataTypes::simdRealWidth; s++)
if (c6 > 0 && c12 > 0)
{
/* c12 is stored scaled with 12.0 and c6 is scaled with 6.0 - correct for this */
- preloadSigma6VdWGapsys[i][s] = 0.5_real * c12 / c6;
+ preloadGapsysSigma6VdW[i][s] = 0.5_real * c12 / c6;
}
else
{
- preloadSigma6VdWGapsys[i][s] = sigma6VdWGapsys;
+ preloadGapsysSigma6VdW[i][s] = gapsysSigma6VdW;
}
}
}
preloadLjPmeC6Grid[i][s] = 0;
preloadQq[i][s] = 0;
preloadSigma6[i][s] = 0;
- preloadSigma6VdWGapsys[i][s] = 0;
+ preloadGapsysSigma6VdW[i][s] = 0;
}
}
}
RealType gmx_unused alphaCoulEff;
RealType gmx_unused gapsysScaleLinpointVdWEff;
RealType gmx_unused gapsysScaleLinpointCoulEff;
- RealType gmx_unused sigmaVdWGapsysEff[NSTATES];
+ RealType gmx_unused gapsysSigma6VdWEff[NSTATES];
for (int i = 0; i < NSTATES; i++)
{
gmx::gatherLoadTranspose<2>(nbfp.data(), typeIndices[i], &c6[i], &c12[i]);
}
if constexpr (softcoreType == KernelSoftcoreType::Gapsys)
{
- sigmaVdWGapsysEff[i] = gmx::load<RealType>(preloadSigma6VdWGapsys[i]);
+ gapsysSigma6VdWEff[i] = gmx::load<RealType>(preloadGapsysSigma6VdW[i]);
}
}
if constexpr (softcoreType == KernelSoftcoreType::Beutler)
rSq,
LFV[i],
DLF[i],
- sigmaVdWGapsysEff[i],
+ gapsysSigma6VdWEff[i],
gapsysScaleLinpointVdWEff,
repulsionShift,
dispersionShift,
fepVals_.bScCoul = true;
fepVals_.scGapsysScaleLinpointLJ = 0.3;
fepVals_.scGapsysScaleLinpointQ = 0.3;
- fepVals_.scSigmaLJGapsys = 0.3;
+ fepVals_.scGapsysSigmaLJ = 0.3;
fepVals_.softcoreFunction = SoftcoreType::Beutler;
}
fep->bScCoul = (getEnum<Boolean>(&inp, "sc-coul", wi) != Boolean::No);
fep->scGapsysScaleLinpointLJ = get_ereal(&inp, "sc-gapsys-scale-linpoint-LJ", 0.85, wi);
fep->scGapsysScaleLinpointQ = get_ereal(&inp, "sc-gapsys-scale-linpoint-Q", 0.3, wi);
- fep->scSigmaLJGapsys = get_ereal(&inp, "sc-sigma-LJ-gapsys", 0.3, wi);
+ fep->scGapsysSigmaLJ = get_ereal(&inp, "sc-gapsys-sigma-LJ", 0.3, wi);
fep->dh_hist_size = get_eint(&inp, "dh_hist_size", 0, wi);
fep->dh_hist_spacing = get_ereal(&inp, "dh_hist_spacing", 0.1, wi);
fep->separate_dhdl_file = getEnum<SeparateDhdlFile>(&inp, "separate-dhdl-file", wi);
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
sc-coul = no
sc-gapsys-scale-linpoint-LJ = 0.85
sc-gapsys-scale-linpoint-Q = 0.3
-sc-sigma-LJ-gapsys = 0.3
+sc-gapsys-sigma-LJ = 0.3
separate-dhdl-file = yes
dhdl-derivatives = yes
dh_hist_size = 0
real fscal_vdw[2], fscal_elec[2];
real velec[2], vvdw[2];
real dvdl_elec[2], dvdl_vdw[2];
- real gapsysScaleLinpointCoul, gapsysScaleLinpointVdW, sigma6VdWGapsys[2];
+ real gapsysScaleLinpointCoul, gapsysScaleLinpointVdW, gapsysSigma6VdW[2];
real rQ, rLJ;
real scaleDvdlRCoul;
int i, ntab;
/* The c6 & c12 coefficients now contain the constants 6.0 and 12.0, respectively.
* Correct for this by multiplying with (1/12.0)/(1/6.0)=6.0/12.0=0.5.
*/
- sigma6VdWGapsys[i] = half * c12[i] / c6[i];
+ gapsysSigma6VdW[i] = half * c12[i] / c6[i];
}
else
{
- sigma6VdWGapsys[i] = scParams.sigma6VdWGapsys;
+ gapsysSigma6VdW[i] = scParams.gapsysSigma6VdW;
}
}
}
if (LFV[i] < 1)
{
- rLJ = gmx::sixthroot(c_twentySixSeventh * sigma6VdWGapsys[i] * (one - LFV[i]));
+ rLJ = gmx::sixthroot(c_twentySixSeventh * gapsysSigma6VdW[i] * (one - LFV[i]));
rLJ *= gapsysScaleLinpointVdW;
}
else
fepVals_.bScCoul = true;
fepVals_.scGapsysScaleLinpointLJ = 0.3;
fepVals_.scGapsysScaleLinpointQ = 0.3;
- fepVals_.scSigmaLJGapsys = 0.3;
+ fepVals_.scGapsysSigmaLJ = 0.3;
fepVals_.softcoreFunction = SoftcoreType::Beutler;
fr_.ic = std::make_unique<interaction_const_t>();
PS("sc-function", enumValueToString(fep->softcoreFunction));
PR("sc-gapsys-scale-linpoint-LJ", fep->scGapsysScaleLinpointLJ);
PR("sc-gapsys-scale-linpoint-Q", fep->scGapsysScaleLinpointQ);
- PR("sc-sigma-LJ-gapsys", fep->scSigmaLJGapsys);
+ PR("sc-gapsys-sigma-LJ", fep->scGapsysSigmaLJ);
};
static void pr_pull(FILE* fp, int indent, const pull_params_t& pull)
fep2->scGapsysScaleLinpointQ,
ftol,
abstol);
- cmp_real(fp, "inputrec->fepvals->scSigmaLJGapsys", -1, fep1->scSigmaLJGapsys, fep2->scSigmaLJGapsys, ftol, abstol);
+ cmp_real(fp, "inputrec->fepvals->scGapsysSigmaLJ", -1, fep1->scGapsysSigmaLJ, fep2->scGapsysSigmaLJ, ftol, abstol);
}
void cmp_inputrec(FILE* fp, const t_inputrec* ir1, const t_inputrec* ir2, real ftol, real abstol)
//! scale for the linearization point for the coulomb interaction with gapsys soft-core
real scGapsysScaleLinpointQ = 0.3;
//! lower bound for c12/c6 in gapsys soft-core
- real scSigmaLJGapsys = 0.3;
+ real scGapsysSigmaLJ = 0.3;
//! Whether to print the dvdl term associated with this term; if it is not specified as separate, it is lumped with the FEP term
gmx::EnumerationArray<FreeEnergyPerturbationCouplingType, bool> separate_dvdl;
//! Whether to write a separate dhdl.xvg file note: NOT a gmx_bool, but an enum
softcoreType(fepvals.softcoreFunction),
gapsysScaleLinpointVdW(fepvals.scGapsysScaleLinpointLJ),
gapsysScaleLinpointCoul(fepvals.scGapsysScaleLinpointQ),
- sigma6VdWGapsys(gmx::power6(fepvals.scSigmaLJGapsys))
+ gapsysSigma6VdW(gmx::power6(fepvals.scGapsysSigmaLJ))
{
// This is checked during tpr reading, so we can assert here
GMX_RELEASE_ASSERT(fepvals.sc_r_power == 6.0, "We only support soft-core r-power 6");
// (gapsys sc) linearization point scaling for Coulomb interactions
real gapsysScaleLinpointCoul;
// (gapsys sc) lower bound/replacement for c12/c6 in vdw interactions
- real sigma6VdWGapsys;
+ real gapsysSigma6VdW;
};
/* VdW */