LINCS constraint triangle accuracy issue with OpenMP
When multiple OpenMP threads are used, LINCS will only apply double the iterations for constraint triangles (which usually only occurs when replacing hydrogens with vsites in CH3/NH3) when the last thread has at least one such triangle.
Fix LINCS triangle constraining with OpenMP
LINCS would only apply double the iterations to constraints
in triangles when the last OpenMP task contained at least one such