*potential = quadrFac - 3 * (linFac - constFac);
- RealType lambdaFacRevInv = gmx::maskzInv(1.0 - lambdaFac, dvdlMask);
+ RealType lambdaFacRevInv = gmx::maskzInv(1 - lambdaFac, dvdlMask);
*dvdl = dLambdaFac * 0.5_real * (lambdaFac * lambdaFacRevInv) * (quadrFac - 2 * linFac + constFac);
}
RealType forceQuad = -quadrFac + linearFac;
RealType potentialQuad = c_half * quadrFac - linearFac + constFac;
RealType dvdlQuad = dLambdaFac * 28 * (lambdaFac * lambdaFacRevInv)
- * ((6.5_real * rInv14C - rInv8C) - (13 * rInv13C - 2. * rInv7C)
+ * ((6.5_real * rInv14C - rInv8C) - (13 * rInv13C - 2 * rInv7C)
+ (6.5_real * rInv12C - rInv6C));
potentialQuad = potentialQuad
fscal_vdw[i] = -quadrFac + linearFac;
fscal_vdw[i] *= rpinv;
- vvdw[i] = 0.5 * quadrFac - linearFac + constFac;
+ vvdw[i] = 0.5_real * quadrFac - linearFac + constFac;
dvdl_vdw[i] += DLF[i] * 28 * (LFV[i] / (one - LFV[i]))
* ((6.5_real * rInv14 - rInv8) - (13 * rInv13 - 2 * rInv7)