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

Fixed nan problem if m_normResid0 = 0.0 at line 965

Browse source
This commit is contained in:
Christopher Lueth 2010-10-13 22:06:03 +00:00
parent ae989985eb
commit 9630216f3e
1 changed files with 6 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
Cantera/src/numerics/NonlinearSolver.cpp
View file

@ -961,7 +961,12 @@ namespace Cantera {
* We aren't going to solve the system if we don't need to. Therefore, return an estimate
* of the next solution update based on the ratio of the residual reduction.
*/
s1 = s0 * m_normResidTrial / m_normResid0;
if (m_normResid0 > 0.0) {
s1 = s0 * m_normResidTrial / m_normResid0;
}
else {
s1 = 0;
}
if (m_normResidTrial < 1.0) {
retnTrial = 3;
} else {