Fix sign in 14-interactions with sc-gapsys.
authorSebastian Kehl <sebastian.kehl@mpcdf.mpg.de>
Fri, 24 Sep 2021 09:49:24 +0000 (11:49 +0200)
committerMagnus Lundborg <magnus.lundborg@scilifelab.se>
Tue, 5 Oct 2021 11:38:46 +0000 (11:38 +0000)
src/gromacs/listed_forces/pairs.cpp

index c97718459e7bab513f94d8dcb7811e1437c5d406..59fb748c553092a971292b97d41e6bc1b0f58dc3 100644 (file)
@@ -339,7 +339,7 @@ static real free_energy_evaluate_single(real
                 real quadrFac = linFac * r * rInvQ;
 
                 /* Computing Coulomb force and potential energy */
-                fscal_elec[i] = 2 * quadrFac - 3 * linFac;
+                fscal_elec[i] = -2 * quadrFac + 3 * linFac;
                 fscal_elec[i] *= rpinv;
 
                 velec[i] = quadrFac - 3 * (linFac - constFac);
@@ -416,7 +416,7 @@ static real free_energy_evaluate_single(real
                 constFac  = 91 * rInv12 - 28 * rInv6;
 
                 /* Computing LJ force and potential energy*/
-                fscal_vdw[i] = quadrFac - linearFac;
+                fscal_vdw[i] = -quadrFac + linearFac;
                 fscal_vdw[i] *= rpinv;
 
                 vvdw[i] = 0.5 * quadrFac - linearFac + constFac;