LINCS: a linear constraint solver
This code shows how to implement LINCS in C++. LINCS is a LINear Constraint Solver for molecular simulations with bond constraints.
LINCS is described by: Hess B, Bekker H, Berendsen HJC and Fraaije JGEM (1997) The original paper can be found here
This implementation was created by Marnix Kok and Wieger Hofstra and used in a project that simulates the coils of a rope that falls to the ground.
Download the algorithm implementation here: lincs.cpp