[Kinetics] Restore old handling of repeated species in path diagrams

Fixes an error introduced in 37f71bd9 which caused these reactions to be
ignored. However, flux calculations for reactions such as H + HO2 -> 2 OH are
still incorrect.
This commit is contained in:
Ray Speth 2016-11-24 11:42:17 -05:00
parent 0d982c8f58
commit 3093e6e6d4

View file

@ -566,10 +566,10 @@ int ReactionPathBuilder::init(ostream& logfile, Kinetics& kin)
vector<vector<size_t> > allReactants(m_nr);
for (size_t i = 0; i < m_nr; i++) {
for (size_t k = 0; k < m_ns; k++) {
if (kin.reactantStoichCoeff(k, i)) {
for (int n = 0; n < kin.reactantStoichCoeff(k, i); n++) {
allReactants[i].push_back(k);
}
if (kin.productStoichCoeff(k, i)) {
for (int n = 0; n < kin.productStoichCoeff(k, i); n++) {
allProducts[i].push_back(k);
}
}