--- branches/development/src/nonbonded/GB.cpp	2012/02/16 15:59:20	1674
+++ branches/development/src/nonbonded/GB.cpp	2012/03/10 04:21:44	1686
@@ -374,10 +374,10 @@ namespace OpenMD {
     Vector3d uxu = cross(ul1, ul2);
 
     (*(idat.pot))[VANDERWAALS_FAMILY] += U *  *(idat.sw);
-    *(idat.f1) += dUdr * rhat + dUda * ul1 + dUdb * ul2;    
-    *(idat.t1) += dUda * rxu1 - dUdg * uxu;
-    *(idat.t2) += dUdb * rxu2 + dUdg * uxu;
-    *(idat.vpair) += U * *(idat.sw);
+    *(idat.f1) += (dUdr * rhat + dUda * ul1 + dUdb * ul2) * *(idat.sw);
+    *(idat.t1) += (dUda * rxu1 - dUdg * uxu) * *(idat.sw);
+    *(idat.t2) += (dUdb * rxu2 + dUdg * uxu) * *(idat.sw);
+    *(idat.vpair) += U;
 
     return;