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

[Kinetics] Relax tolerance on reaction element balance

Browse source 
This allows use of some mechanisms with imprecise fractional product
stoichiometries.
This commit is contained in:
Ray Speth 2015-02-20 00:11:12 +00:00
parent 79cf2dcba7
commit e545d123c2
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/kinetics/Kinetics.cpp
View file

@ -329,7 +329,7 @@ void Kinetics::checkReactionBalance(const Reaction& R)
const string& elem = iter->first;
double elemsum = balr[elem] + balp[elem];
double elemdiff = fabs(balp[elem] - balr[elem]);
if (elemsum > 0.0 && elemdiff/elemsum > 1e-5) {
if (elemsum > 0.0 && elemdiff/elemsum > 1e-4) {
ok = false;
msg += " " + elem + " " + fp2str(balr[elem]) +
" " + fp2str(balp[elem]) + "\n";