Commit aba4c054 authored by Jason Swails's avatar Jason Swails
Browse files

Fix old behavior.

parent eb215a4d
...@@ -783,9 +783,8 @@ def readAmberSystem(topology, prmtop_filename=None, prmtop_loader=None, shake=No ...@@ -783,9 +783,8 @@ def readAmberSystem(topology, prmtop_filename=None, prmtop_loader=None, shake=No
atomI = topatoms[iAtom] atomI = topatoms[iAtom]
atomJ = topatoms[jAtom] atomJ = topatoms[jAtom]
atomK = topatoms[kAtom] atomK = topatoms[kAtom]
numH = ((atomI.element.atomic_number == 1) + (atomJ.element.atomic_number == 1) + numH = ((atomI.element.atomic_number == 1) + (atomK.element.atomic_number == 1))
(atomK.element.atomic_number == 1)) constrained = (numH == 2 or (numH == 1 and atomJ.element is elem.oxygen))
constrained = (numH == 2 or (numH == 1 and atomK.element is elem.oxygen))
else: else:
constrained = False constrained = False
if constrained: if constrained:
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment