diff --git a/include/cantera/kinetics/InterfaceKinetics.h b/include/cantera/kinetics/InterfaceKinetics.h index 331514ca5..18aa1064c 100644 --- a/include/cantera/kinetics/InterfaceKinetics.h +++ b/include/cantera/kinetics/InterfaceKinetics.h @@ -133,31 +133,21 @@ public: * @param fwdROP vector of rates of progress. * length = number of reactions, Units are kmol m-2 s-1. */ - virtual void getFwdRatesOfProgress(doublereal* fwdROP) { - updateROP(); - std::copy(m_ropf.begin(), m_ropf.end(), fwdROP); - } + virtual void getFwdRatesOfProgress(doublereal* fwdROP); //! Return the reverse rates of progress for each reaction /*! * @param revROP vector of rates of progress. * length = number of reactions, Units are kmol m-2 s-1. */ - virtual void getRevRatesOfProgress(doublereal* revROP) { - updateROP(); - std::copy(m_ropr.begin(), m_ropr.end(), revROP); - } + virtual void getRevRatesOfProgress(doublereal* revROP); //! Return the net rates of progress for each reaction /*! * @param netROP vector of rates of progress. * length = number of reactions, Units are kmol m-2 s-1. */ - virtual void getNetRatesOfProgress(doublereal* netROP) { - updateROP(); - std::copy(m_ropnet.begin(), m_ropnet.end(), netROP); - } - + virtual void getNetRatesOfProgress(doublereal* netROP); //! Get the equilibrium constants of all reactions, whether //! the reaction is reversible or not. diff --git a/src/kinetics/InterfaceKinetics.cpp b/src/kinetics/InterfaceKinetics.cpp index cc49853c0..82867bbe3 100644 --- a/src/kinetics/InterfaceKinetics.cpp +++ b/src/kinetics/InterfaceKinetics.cpp @@ -215,6 +215,7 @@ void InterfaceKinetics::setElectricPotential(int n, doublereal V) thermo(n).setElectricPotential(V); m_redo_rates = true; } + //==================================================================================================================== // Update properties that depend on temperature /* @@ -393,6 +394,23 @@ void InterfaceKinetics::checkPartialEquil() } } +void InterfaceKinetics::getFwdRatesOfProgress(doublereal* fwdROP) +{ + updateROP(); + std::copy(m_ropf.begin(), m_ropf.end(), fwdROP); +} + +void InterfaceKinetics::getRevRatesOfProgress(doublereal* revROP) +{ + updateROP(); + std::copy(m_ropr.begin(), m_ropr.end(), revROP); +} + +void InterfaceKinetics::getNetRatesOfProgress(doublereal* netROP) +{ + updateROP(); + std::copy(m_ropnet.begin(), m_ropnet.end(), netROP); +} /** * Get the equilibrium constants of all reactions, whether