A Non-Iterative Matrix Method for Constraint Molecular Dynamics Simulations | Synapse