A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations