From d7463cc22665acf4b5a869f492a5d2b2c6d73eed Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Fri, 20 Jul 2012 18:26:50 +0000 Subject: [PATCH] Removed unimplemented TST and LandauTeller kinetic rates --- include/cantera/kinetics/RxnRates.h | 73 ----------------------------- 1 file changed, 73 deletions(-) diff --git a/include/cantera/kinetics/RxnRates.h b/include/cantera/kinetics/RxnRates.h index 4a4772490..ce951c6a7 100644 --- a/include/cantera/kinetics/RxnRates.h +++ b/include/cantera/kinetics/RxnRates.h @@ -246,52 +246,6 @@ protected: }; -#ifdef INCL_TST - -class TST -{ - -public: - static int type() { - return TSTRATE; - } - TST() {} - TST(const vector_fp& c) { - m_b.resize(10); - copy(c.begin(), c.begin() + 10, m_b.begin()); - m_k = int(c[10]); - } - - void update_C(const vector_fp& c) { - doublereal ck = c[m_k]; - delta_s0 = m_b[0] + m_b[1]*ck + m_b[2]*ck*ck; - delta_e0 = m_b[5] + m_b[6]*ck + m_b[7]*ck*ck; - } - - doublereal update(doublereal logT, doublereal recipT) const { - doublereal delta_s = delta_s0*(1.0 + m_b[3]*logT + m_b[4]*recipT); - doublereal delta_E = delta_e0*(1.0 + m_b[8]*logT + m_b[9]*recipT); - return logBoltz_Planck + logT + delta_s - delta_E*recipT; - } - - doublereal updateRC(doublereal logT, doublereal recipT) const { - doublereal lres = update(logT, recipT); - return exp(lres); - } - - void writeUpdateRHS(std::ostream& s) const {} - -protected: - doublereal delta_s0, delta_e0; - int m_k; - vector_fp m_b; -}; - -#endif - - - - //! Arrhenius reaction rate type depends only on temperature /** * A reaction rate coefficient of the following form. @@ -675,33 +629,6 @@ protected: vector_fp dotProd_; //!< dot product of chebCoeffs with the reduced pressure polynomial }; -// class LandauTeller { - -// public: -// static int type(){ return LANDAUTELLER; } -// LandauTeller(){} -// LandauTeller( const vector_fp& c ) : m_c(c) { m_c[0] = log(c[0]); } - -// doublereal update(doublereal logT, doublereal recipT) const { -// return m_c[0] + m_c[1]*tt[1] - m_c[2]*tt[2] -// + m_c[3]*tt[3] + m_c[4]*tt[4]; -// } - -// //void writeUpdateRHS(ostream& s) const { -// // s << exp(m_logA); -// // s << " * exp("; -// // if (m_b != 0.0) s << m_b << " * tlog"; -// // if (m_E != 0.0) s << " - " << m_E << " * rt"; -// // if (m_E != 0.0) s << " - " << m_E << " * rt"; -// // s << ");" << endl; -// // } -// //} - -// protected: -// doublereal m_logA, m_b, m_E; -// }; - -//} } #endif