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

Buf fixes for handling boundaries.

Browse source
This commit is contained in:
Harry Moffat 2011-08-15 03:57:40 +00:00
parent d617f43013
commit 36a56093d1
1 changed files with 2 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
Cantera/src/numerics/RootFind.cpp
View file

@ -772,14 +772,14 @@ namespace Cantera {
if (xnew < xmin) {
bottomBump++;
if (bottomBump < 3) {
xnew = x2 + (x2 - xmin) / 2.0;
xnew = x2 - (x2 - xmin) / 2.0;
} else {
if (x2 == xmin || x1 == xmin) {
// we are here when we are bumping against the bottom limit.
// No further action is possible
if (xnew < xmin) {
slope = (f2 - f1) / delXtmp;
xmin = x2 - f2 / slope;
xnew = x2 - f2 / slope;
if (xnew < xmin) {
retn = ROOTFIND_SOLNLOWERTHANXMIN;
*xbest = xnew;