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

added statements to getEquilibriumConstants to force a T update on the next call. Array m_rkcn was being left containing incorrect data. Problem is now fixed.

Browse source
This commit is contained in:
Dave Goodwin 2005-12-15 17:07:06 +00:00
parent 33fbbf0b93
commit 9a67ae84c4
1 changed files with 5 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
Cantera/src/GasKinetics.cpp
View file

@ -139,6 +139,10 @@ namespace Cantera {
for (i = 0; i < m_ii; i++) {
kc[i] = exp(-rkc[i]*rrt + m_dn[i]*logStandConc);
}
// force an update of T-dependent properties, so that m_rkcn will
// be updated before it is used next.
m_kdata->m_temp = 0.0;
}
/**
@ -440,7 +444,7 @@ namespace Cantera {
*/
const vector_fp& m_rkc = m_kdata->m_rkcn;
for (int i = 0; i < m_ii; i++) {
krev[i] *= m_rkc[i];
krev[i] *= m_rkc[i];
}
}
}