Commit 6aea38c9 authored by pgrinaway's avatar pgrinaway
Browse files

make clearer

parent 405607e0
...@@ -856,13 +856,13 @@ class ForceField(object): ...@@ -856,13 +856,13 @@ class ForceField(object):
uniquePropers = set() uniquePropers = set()
for angle in data.angles: for angle in data.angles:
for atom in bondedToAtom[angle[0]]: for atom in bondedToAtom[angle[0]]:
if atom != angle[1] and atom != angle[2]: if atom not in angle:
if atom < angle[2]: if atom < angle[2]:
uniquePropers.add((atom, angle[0], angle[1], angle[2])) uniquePropers.add((atom, angle[0], angle[1], angle[2]))
else: else:
uniquePropers.add((angle[2], angle[1], angle[0], atom)) uniquePropers.add((angle[2], angle[1], angle[0], atom))
for atom in bondedToAtom[angle[2]]: for atom in bondedToAtom[angle[2]]:
if atom != angle[1] and atom != angle[0]: if atom not in angle:
if atom > angle[0]: if atom > angle[0]:
uniquePropers.add((angle[0], angle[1], angle[2], atom)) uniquePropers.add((angle[0], angle[1], angle[2], atom))
else: else:
......
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