From 1e2ba1f0d258332b2728f4e779fba05e56b76553 Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Tue, 2 Sep 2014 22:34:21 +0000 Subject: [PATCH] [1D/Python] Fix spurious 'divide by zero' warning --- interfaces/cython/cantera/onedim.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/interfaces/cython/cantera/onedim.py b/interfaces/cython/cantera/onedim.py index be3421539..98ea9e759 100644 --- a/interfaces/cython/cantera/onedim.py +++ b/interfaces/cython/cantera/onedim.py @@ -543,8 +543,8 @@ class CounterflowDiffusionFlame(FlameBase): stoich += self.gas.n_atoms(fuel, 'C') d_u_d_z = np.gradient(self.u) / np.gradient(self.grid) - phi = (self.X[self.gas.species_index(fuel)] - / self.X[self.gas.species_index(oxidizer)] * stoich) + phi = (self.X[self.gas.species_index(fuel)] * stoich / + np.maximum(self.X[self.gas.species_index(oxidizer)], 1e-20)) z_stoich = np.interp(-1., -phi, self.grid) return np.abs(np.interp(z_stoich, self.grid, d_u_d_z))