Cloudy & Associates

Commit ce3fb9de authored by Robin Williams's avatar Robin Williams
Browse files

Tighten threshold on negative abundances in chemistry solver

Previous implementation had been intended to ignore values close to
underflow, but the number chosen was large enough to be a significant
source for low abundance species.  Tightening the criterion to cut the
step length if any negative abundance appears works fine after wide
testing, so seems like a better approach.
parent 1e690db4
......@@ -140,10 +140,7 @@ double mole_solve()
newMolsTmp = abs(newmols[mol]);
iworst = mol;
}
if( newmols[mol] < -SMALLABUND)
{
++nBad;
}
++nBad;
newmols[mol] = 0.;
}
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment