fixed PME forces with FE, fixes #929 introduced in recent merge

Change-Id: I4f3bc59ed8c0351c4c59feceb3c2116b2f063daf

diff --git a/src/mdlib/pme.c b/src/mdlib/pme.c
index 1337a0511cfe907d74ed9db8c5f4309b01d8f182..ed45fbfb4310a7feec78d608989c24071303c15d 100644 (file)
--- a/src/mdlib/pme.c
+++ b/src/mdlib/pme.c
@@ -2147,10 +2147,10 @@ for(ithx=0; (ithx<order); ithx++)              \
 }
 
 
-void gather_f_bsplines(gmx_pme_t pme,real *grid,
-                       gmx_bool bClearF,pme_atomcomm_t *atc,
-                       splinedata_t *spline,
-                       real scale)
+static void gather_f_bsplines(gmx_pme_t pme,real *grid,
+                              gmx_bool bClearF,pme_atomcomm_t *atc,
+                              splinedata_t *spline,
+                              real scale)
 {
     /* sum forces for local particles */
     int     nn,n,ithx,ithy,ithz,i0,j0,k0;
@@ -2193,7 +2193,7 @@ void gather_f_bsplines(gmx_pme_t pme,real *grid,
     for(nn=0; nn<spline->n; nn++)
     {
         n  = spline->ind[nn];
-        qn = atc->q[n];
+        qn = scale*atc->q[n];
 
         if (bClearF)
         {