Fix LINCS triangle constraining with OpenMP

[alexxy/gromacs.git] / src / gromacs / mdlib / lincs.cpp

diff --git a/src/gromacs/mdlib/lincs.cpp b/src/gromacs/mdlib/lincs.cpp
index 0e15e5c8ede23a5cc9e5be756951908e2bcf0609..543c1848dd2a2239f1848676de5e3aa04f22f474 100644 (file)
--- a/src/gromacs/mdlib/lincs.cpp
+++ b/src/gromacs/mdlib/lincs.cpp
@@ -1377,7 +1377,7 @@ static void set_lincs_matrix(Lincs *li, real *invmass, real lambda)
         {
             set_lincs_matrix_task(li, &li->task[th], invmass,
                                   &ncc_triangle, &nCrossTaskTriangles);
-            ntriangle = li->task[th].ntriangle;
+            ntriangle += li->task[th].ntriangle;
         }
         GMX_CATCH_ALL_AND_EXIT_WITH_FATAL_ERROR;
     }