@@ -167,21 +167,21 @@ def test_compute_flux(self):
167167
168168 # Check that we're computing the species fluxes correctly
169169 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 (
171171 1e-6 * reaction_rates [i , 0 ]
172172 )
173- assert abs (reaction_rates [i , 0 ] - - species_rates [i , 1 ]) < abs (
173+ assert abs (reaction_rates [i , 0 ] - - species_rates [i , 1 ]) <= abs (
174174 1e-6 * reaction_rates [i , 0 ]
175175 )
176- assert abs (reaction_rates [i , 0 ] - - species_rates [i , 2 ]) < abs (
176+ assert abs (reaction_rates [i , 0 ] - - species_rates [i , 2 ]) <= abs (
177177 1e-6 * reaction_rates [i , 0 ]
178178 )
179- assert abs (reaction_rates [i , 0 ] - species_rates [i , 3 ]) < abs (
179+ assert abs (reaction_rates [i , 0 ] - species_rates [i , 3 ]) <= abs (
180180 1e-6 * reaction_rates [i , 0 ]
181181 )
182182
183183 # 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
185185
186186 def test_jacobian (self ):
187187 """
0 commit comments