fix bug in Reference CustomCompoundBondForce

c52a3d25 · Mike Garrahan · 8369f91d · c52a3d25
Commit c52a3d25 authored Mar 04, 2014 by Mike Garrahan
Show whitespace changes
Inline Side-by-side

Showing with 1 addition and 1 deletion

platforms/reference/src/SimTKReference/ReferenceCustomCompoundBondIxn.cpp ...nce/src/SimTKReference/ReferenceCustomCompoundBondIxn.cpp +1 -1

No files found.
--- a/platforms/reference/src/SimTKReference/ReferenceCustomCompoundBondIxn.cpp
+++ b/platforms/reference/src/SimTKReference/ReferenceCustomCompoundBondIxn.cpp
@@ -124,7 +124,7 @@ void ReferenceCustomCompoundBondIxn::calculateOneIxn(int bond, vector<RealVec>&

    // Compute all of the variables the energy can depend on.

-    const vector<int>& atoms = bondAtoms[0];
+    const vector<int>& atoms = bondAtoms[bond];
    for (int i = 0; i < (int) particleTerms.size(); i++) {
        const ParticleTermInfo& term = particleTerms[i];
        variables[term.name] = atomCoordinates[term.atom][term.component];