research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[1D/Python] Fix spurious 'divide by zero' warning

Browse source
This commit is contained in:
Ray Speth 2014-09-02 22:34:21 +00:00
parent 46231379c1
commit 1e2ba1f0d2
1 changed files with 2 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
interfaces/cython/cantera/onedim.py
View file

@ -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))