Skip to content

Commit

Permalink
fixes #5252 (#5271)
Browse files Browse the repository at this point in the history
  • Loading branch information
ptosco committed May 8, 2022
1 parent ab760b8 commit ba0a7c6
Showing 1 changed file with 12 additions and 12 deletions.
24 changes: 12 additions & 12 deletions Code/ForceField/Wrap/testConstraints.py
Expand Up @@ -148,23 +148,23 @@ def testUFFAngleConstraints(self):
self.assertTrue(r == 0)
conf = m.GetConformer()
angle = rdMolTransforms.GetAngleDeg(conf, 1, 3, 6)
self.assertTrue(int(angle) == 90)
self.assertAlmostEqual(angle, 90, delta=0.1)
ff = ChemicalForceFields.UFFGetMoleculeForceField(m)
self.assertTrue(ff)
ff.UFFAddAngleConstraint(1, 3, 6, True, -10.0, 10.0, 100.0)
r = ff.Minimize()
self.assertTrue(r == 0)
conf = m.GetConformer()
angle = rdMolTransforms.GetAngleDeg(conf, 1, 3, 6)
self.assertTrue(int(angle) == 100)
self.assertAlmostEqual(angle, 100, delta=0.1)
ff = ChemicalForceFields.UFFGetMoleculeForceField(m)
self.assertTrue(ff)
ff.UFFAddAngleConstraint(1, 3, 6, False, -10.0, 10.0, 100.0)
r = ff.Minimize()
self.assertTrue(r == 0)
conf = m.GetConformer()
angle = rdMolTransforms.GetAngleDeg(conf, 1, 3, 6)
self.assertEqual(int(angle), 10)
self.assertAlmostEqual(angle, 10, delta=0.1)

def testUFFTorsionConstraints(self):
m = Chem.MolFromMolBlock(self.molB, True, False)
Expand All @@ -176,7 +176,7 @@ def testUFFTorsionConstraints(self):
r = ff.Minimize()
self.assertTrue(r == 0)
dihedral = rdMolTransforms.GetDihedralDeg(conf, 1, 3, 6, 8)
self.assertTrue(int(dihedral) == 20)
self.assertAlmostEqual(dihedral, 20, delta=0.1)
rdMolTransforms.SetDihedralDeg(conf, 1, 3, 6, 8, -30.0)
ff = ChemicalForceFields.UFFGetMoleculeForceField(m)
self.assertTrue(ff)
Expand All @@ -185,7 +185,7 @@ def testUFFTorsionConstraints(self):
self.assertTrue(r == 0)
conf = m.GetConformer()
dihedral = rdMolTransforms.GetDihedralDeg(conf, 1, 3, 6, 8)
self.assertEquals(int(dihedral), -40)
self.assertAlmostEqual(dihedral, -40, delta=0.1)
rdMolTransforms.SetDihedralDeg(conf, 1, 3, 6, 8, -10.0)
ff = ChemicalForceFields.UFFGetMoleculeForceField(m)
self.assertTrue(ff)
Expand All @@ -194,7 +194,7 @@ def testUFFTorsionConstraints(self):
self.assertTrue(r == 0)
conf = m.GetConformer()
dihedral = rdMolTransforms.GetDihedralDeg(conf, 1, 3, 6, 8)
self.assertTrue(int(dihedral) == -10)
self.assertAlmostEqual(dihedral, -10, delta=0.1)

def testUFFPositionConstraints(self):
m = Chem.MolFromMolBlock(self.molB, True, False)
Expand Down Expand Up @@ -262,23 +262,23 @@ def testMMFFAngleConstraints(self):
self.assertTrue(r == 0)
conf = m.GetConformer()
angle = rdMolTransforms.GetAngleDeg(conf, 1, 3, 6)
self.assertTrue(int(angle) == 90)
self.assertAlmostEqual(angle, 90, delta=0.1)
ff = ChemicalForceFields.MMFFGetMoleculeForceField(m, mp)
self.assertTrue(ff)
ff.MMFFAddAngleConstraint(1, 3, 6, True, -10.0, 10.0, 100.0)
r = ff.Minimize()
self.assertTrue(r == 0)
conf = m.GetConformer()
angle = rdMolTransforms.GetAngleDeg(conf, 1, 3, 6)
self.assertTrue(int(angle) == 100)
self.assertAlmostEqual(angle, 100, delta=0.1)
ff = ChemicalForceFields.MMFFGetMoleculeForceField(m, mp)
self.assertTrue(ff)
ff.MMFFAddAngleConstraint(1, 3, 6, False, -10.0, 10.0, 100.0)
r = ff.Minimize()
self.assertTrue(r == 0)
conf = m.GetConformer()
angle = rdMolTransforms.GetAngleDeg(conf, 1, 3, 6)
self.assertEquals(int(angle), 10) #(int(angle) == 10)
self.assertAlmostEqual(angle, 10, delta=0.1)

def testMMFFTorsionConstraints(self):
m = Chem.MolFromMolBlock(self.molB, True, False)
Expand All @@ -291,7 +291,7 @@ def testMMFFTorsionConstraints(self):
r = ff.Minimize()
self.assertTrue(r == 0)
dihedral = rdMolTransforms.GetDihedralDeg(conf, 1, 3, 6, 8)
self.assertEquals(int(dihedral), 20)
self.assertAlmostEqual(dihedral, 20, delta=0.1)
rdMolTransforms.SetDihedralDeg(conf, 1, 3, 6, 8, -30.0)
ff = ChemicalForceFields.MMFFGetMoleculeForceField(m, mp)
self.assertTrue(ff)
Expand All @@ -300,7 +300,7 @@ def testMMFFTorsionConstraints(self):
self.assertTrue(r == 0)
conf = m.GetConformer()
dihedral = rdMolTransforms.GetDihedralDeg(conf, 1, 3, 6, 8)
self.assertTrue(int(dihedral) == -40)
self.assertAlmostEqual(dihedral, -40, delta=0.1)
rdMolTransforms.SetDihedralDeg(conf, 1, 3, 6, 8, -10.0)
ff = ChemicalForceFields.MMFFGetMoleculeForceField(m, mp)
self.assertTrue(ff)
Expand All @@ -309,7 +309,7 @@ def testMMFFTorsionConstraints(self):
self.assertTrue(r == 0)
conf = m.GetConformer()
dihedral = rdMolTransforms.GetDihedralDeg(conf, 1, 3, 6, 8)
self.assertTrue(int(dihedral) == -10)
self.assertAlmostEqual(dihedral, -10, delta=0.1)

def testMMFFPositionConstraints(self):
m = Chem.MolFromMolBlock(self.molB, True, False)
Expand Down

0 comments on commit ba0a7c6

Please sign in to comment.