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

Fixed illegal access to empty vector in MultiPhaseEquil constructor

Browse source
This commit is contained in:
Ray Speth 2012-01-09 17:35:02 +00:00
parent 0bd3a247bf
commit 53eaae006a
1 changed files with 4 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
Cantera/src/equil/MultiPhaseEquil.cpp
View file

@ -181,8 +181,10 @@ namespace Cantera {
// Take a very small step in composition space, so that no
// species has precisely zero moles.
vector_fp dxi(m_nsp - m_nel, 1.0e-20);
multiply(m_N, DATA_PTR(dxi), DATA_PTR(m_work));
unsort(m_work);
if (!dxi.empty()) {
multiply(m_N, DATA_PTR(dxi), DATA_PTR(m_work));
unsort(m_work);
}
for (k = 0; k < m_nsp; k++) {
m_moles[k] += m_work[k];