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

[Equil] Check index bounds for known problems in VCS_SOLVE

Browse source 
See #187
This commit is contained in:
Ray Speth 2016-05-02 18:22:01 -04:00
parent aa5f26bda6
commit 2561ab82e3
1 changed files with 8 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
src/equil/vcs_phaseStability.cpp
View file

@ -223,6 +223,10 @@ size_t VCS_SOLVE::vcs_popPhaseID(std::vector<size_t> & phasePopPhaseIDs)
// Single Phase Stability Resolution
size_t kspec = Vphase->spGlobalIndexVCS(0);
size_t irxn = kspec - m_numComponents;
if (irxn > m_deltaGRxn_old.size()) {
throw CanteraError("VCS_SOLVE::vcs_popPhaseID",
"Index out of bounds due to logic error.");
}
doublereal deltaGRxn = m_deltaGRxn_old[irxn];
Fephase = exp(-deltaGRxn) - 1.0;
if (Fephase > 0.0) {
@ -370,6 +374,10 @@ int VCS_SOLVE::vcs_popPhaseRxnStepSizes(const size_t iphasePop)
double delmol = deltaMolNumPhase * X_est[k];
if (kspec >= m_numComponents) {
irxn = kspec - m_numComponents;
if (irxn > m_stoichCoeffRxnMatrix.nColumns()) {
throw CanteraError("VCS_SOLVE::vcs_popPhaseRxnStepSizes",
"Index out of bounds due to logic error.");
}
for (size_t j = 0; j < m_numComponents; ++j) {
double stoicC = m_stoichCoeffRxnMatrix(j,irxn);
if (stoicC != 0.0 && m_elType[j] == VCS_ELEM_TYPE_ABSPOS) {