[Kinetics] Buffer rate index for MultiRate objects

Cantera · Sep 22, 2021 · dda597f · dda597f
1 parent bbc27b4
commit dda597f
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/kinetics/BulkKinetics.cpp b/src/kinetics/BulkKinetics.cpp
@@ -140,6 +140,9 @@ bool BulkKinetics::addReaction(shared_ptr<Reaction> r, bool resize)
             m_bulk_rates.back()->resize(m_kk, nReactions());
         }
 
+        // Set index of rate to number of reaction within kinetics
+        rate->setRateIndex(nReactions() - 1);
+
         // Add reaction rate to evaluator
         size_t index = m_bulk_types[rate->type()];
         m_bulk_rates[index]->add(nReactions() - 1, *rate);