Merge pull request #1628 from peastman/variable

Fixed bug in reference implementation of VariableLangevinIntegrator

Merge pull request #1628 from peastman/variable
Fixed bug in reference implementation of VariableLangevinIntegrator
32d2707c · peastman · GitHub · 49cb1b43 · 38f8ffb5 · 32d2707c
Commit 32d2707c authored Oct 13, 2016 by peastman Committed by GitHub Oct 13, 2016
Hide whitespace changes
Inline Side-by-side

Showing with 3 additions and 3 deletions

platforms/reference/src/SimTKReference/ReferenceVariableStochasticDynamics.cpp ...rc/SimTKReference/ReferenceVariableStochasticDynamics.cpp +3 -3

No files found.
--- a/platforms/reference/src/SimTKReference/ReferenceVariableStochasticDynamics.cpp
+++ b/platforms/reference/src/SimTKReference/ReferenceVariableStochasticDynamics.cpp
@@ -266,11 +266,11 @@ void ReferenceVariableStochasticDynamics::update(const OpenMM::System& system, v

   // copy xPrime -> atomCoordinates

+   RealOpenMM invStepSize = 1.0/getDeltaT();
   for (int ii = 0; ii < numberOfAtoms; ii++) {
       if (masses[ii] != 0.0) {
-           atomCoordinates[ii][0] = xPrime[ii][0];
-           atomCoordinates[ii][1] = xPrime[ii][1];
-           atomCoordinates[ii][2] = xPrime[ii][2];
+           velocities[ii] = (xPrime[ii]-atomCoordinates[ii])*invStepSize;
+           atomCoordinates[ii] = xPrime[ii];
       }
   }