From 7c611bab70fb683eb5dc5c28ddac1d7f86be89eb Mon Sep 17 00:00:00 2001
From: David van der Spoel <spoel@xray.bmc.uu.se>
Date: Thu, 21 Jun 2012 18:14:33 +0200
Subject: [PATCH] Fixes #961 - Shake reports wrong error

The constraint index reported by shake was off by 1.

Change-Id: I6f4d0c013b30f5963fb312b41fb87e90e6314619
---
 src/mdlib/shakef.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/mdlib/shakef.c b/src/mdlib/shakef.c
index f6fd555c4e..7f003fe670 100644
--- a/src/mdlib/shakef.c
+++ b/src/mdlib/shakef.c
@@ -150,7 +150,7 @@ void cshake(atom_id iatom[],int ncon,int *nnit,int maxnit,
           rrpr    = rijx*tx+rijy*ty+rijz*tz;
 	
           if (rrpr < toler*mytol) 
-              error=ll;
+              error=ll+1;
           else {
               acor      = omega*diff*m2[ll]/rrpr;
               lagr[ll] += acor;
-- 
2.22.0