Skip to content

Commit 27bd1f3

Browse files
rwestssun30
authored andcommitted
Fix ArrheniusBM.get_activation_energy unit test.
The E0 had been changed, but the expected Ea not updated.
1 parent 4fad006 commit 27bd1f3

1 file changed

Lines changed: 7 additions & 1 deletion

File tree

rmgpy/kinetics/arrheniusTest.py

Lines changed: 7 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -528,7 +528,13 @@ def test_get_activation_energy(self):
528528
"""
529529
Hrxn = -44000 # J/mol
530530
Ea = self.arrhenius_bm.get_activation_energy(Hrxn)
531-
self.assertAlmostEqual(Ea, 95074, delta=1e1)
531+
532+
w = self.w0
533+
E0 = self.E0
534+
Vp = 2 * w * (w + E0)/(w - E0)
535+
Ea_exp = (w + Hrxn/2) * (Vp - 2*w + Hrxn)**2 / (Vp*Vp - 4*w*w + Hrxn*Hrxn)
536+
537+
self.assertAlmostEqual(Ea, Ea_exp, delta=1e1)
532538

533539

534540
################################################################################

0 commit comments

Comments
 (0)