From 7264db9324c64936008c8978cab864a77bd6ac19 Mon Sep 17 00:00:00 2001
From: David van der Spoel <spoel@xray.bmc.uu.se>
Date: Sat, 30 Jun 2012 15:18:31 +0200
Subject: [PATCH] Fixed bug in free energy calculation for Morse potentials

Introduced recently when FEP was added to the Morse potential.

Change-Id: I90edfd13835637c872af31f2c7abf900ce8a8e5f
---
 src/gmxlib/bondfree.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/gmxlib/bondfree.c b/src/gmxlib/bondfree.c
index 789702a9ca..0b1c6799d6 100644
--- a/src/gmxlib/bondfree.c
+++ b/src/gmxlib/bondfree.c
@@ -159,7 +159,7 @@ real morse_bonds(int nbonds,
     if (temp == one)
     {
         /* bonds are constrainted. This may _not_ include bond constraints if they are lambda dependent */
-        *dvdlambda = cbB-cbA;
+        *dvdlambda += cbB-cbA;
         continue;
     }
 
-- 
2.22.0