@@ -167,21 +167,21 @@ def test_compute_flux(self):
167
167
168
168
# Check that we're computing the species fluxes correctly
169
169
for i in range (t .shape [0 ]):
170
- assert abs (reaction_rates [i , 0 ] - species_rates [i , 0 ]) < abs (
170
+ assert abs (reaction_rates [i , 0 ] - species_rates [i , 0 ]) <= abs (
171
171
1e-6 * reaction_rates [i , 0 ]
172
172
)
173
- assert abs (reaction_rates [i , 0 ] - - species_rates [i , 1 ]) < abs (
173
+ assert abs (reaction_rates [i , 0 ] - - species_rates [i , 1 ]) <= abs (
174
174
1e-6 * reaction_rates [i , 0 ]
175
175
)
176
- assert abs (reaction_rates [i , 0 ] - - species_rates [i , 2 ]) < abs (
176
+ assert abs (reaction_rates [i , 0 ] - - species_rates [i , 2 ]) <= abs (
177
177
1e-6 * reaction_rates [i , 0 ]
178
178
)
179
- assert abs (reaction_rates [i , 0 ] - species_rates [i , 3 ]) < abs (
179
+ assert abs (reaction_rates [i , 0 ] - species_rates [i , 3 ]) <= abs (
180
180
1e-6 * reaction_rates [i , 0 ]
181
181
)
182
182
183
183
# Check that we've reached equilibrium
184
- assert abs (reaction_rates [- 1 , 0 ] - 0.0 ) < 1e-2
184
+ assert abs (reaction_rates [- 1 , 0 ] - 0.0 ) <= 1e-2
185
185
186
186
def test_jacobian (self ):
187
187
"""
0 commit comments