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

Remove unusable KINETICS_WITH_INTERMEDIATE_ZEROED_PHASES macro

Browse source
This commit is contained in:
Ray Speth 2013-06-13 22:50:26 +00:00
parent 75e586da0f
commit 112d9b5093
2 changed files with 0 additions and 57 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

30
include/cantera/kinetics/InterfaceKinetics.h
View file

@ -613,36 +613,6 @@ protected:
*/
std::vector<std::vector<bool> > m_rxnPhaseIsProduct;
#ifdef KINETICS_WITH_INTERMEDIATE_ZEROED_PHASES
//! Vector of ints indicating whether zeroed phase is an intermediate for
//! the formation of another phase
/*!
* If a phase is zeroed out but it is an intermediate, then the phase can
* be formed whether it is stable or not, but the destruction rate of
* species in that phase can't exceed the formation rate for species in
* that phase.
*
* length = number of phases in the object. By default all phases are not
* intermediates
*/
std::vector<int> m_phaseIsIntermediate;
int m_numIntermediatePhases;
//! Reaction rate reduction factor for intermediates
/*!
* Individual reaction rates are reduced to accommodate the requirements of intermediate
* zero phases.
*
* length = number of reactions in the object. By default all phases are
* not intermediates.
*/
std::vector<doublereal> m_rxnRateFactorPhaseIntermediates;
//! Work vector having length number of species
std::vector<doublereal> m_speciesTmpP;
std::vector<doublereal> m_speciesTmpD;
#endif
int m_ioFlag;
};
}

27
src/kinetics/InterfaceKinetics.cpp
View file

@ -606,33 +606,6 @@ void InterfaceKinetics::updateROP()
m_ROP_ok = true;
}
#ifdef KINETICS_WITH_INTERMEDIATE_ZEROED_PHASES
//=================================================================================================
InterfaceKinetics::adjustRatesForIntermediatePhases()
{
doublereal sFac = 1.0;
getCreatingRates(DATA_PTR(m_speciestmpP));
getDestructionRates(DATA_PTR(m_speciestmpD));
for (iphase = 0; iphase < nphases; iphase++) {
if (m_intermediatePhases(iphase)) {
for (isp = 0; isp < nspecies; isp++) {
if (m_speciesTmpD[ispI] > m_speciesTmpP[I]) {
sFac = m_speciesTmpD[ispI]/ m_speciesTmpP[I];
}
// Loop over reactions that are reactants for the species in the phase
// reducing their rates.
}
}
}
}
#endif
void InterfaceKinetics::getDeltaGibbs(doublereal* deltaG)
{
/*