void VCS_SOLVE::vcs_SSPhase()
{
vector_int numPhSpecies(m_numPhases, 0);
for (size_t kspec = 0; kspec < m_nsp; ++kspec) {
numPhSpecies[m_phaseID[kspec]]++;
}
// Handle the special case of a single species in a phase that has been
// earmarked as a multispecies phase. Treat that species as a single-species
// phase
for (size_t iph = 0; iph < m_numPhases; iph++) {
vcs_VolPhase* Vphase = m_VolPhaseList[iph].get();
Vphase->m_singleSpecies = false;
if (TPhInertMoles[iph] > 0.0) {
Vphase->setExistence(2);
}
if (numPhSpecies[iph] <= 1 && TPhInertMoles[iph] == 0.0) {
Vphase->m_singleSpecies = true;
}
}
// Fill in some useful arrays here that have to do with the static
// information concerning the phase ID of species. SSPhase = Boolean
// indicating whether a species is in a single species phase or not.
for (size_t kspec = 0; kspec < m_nsp; kspec++) {
size_t iph = m_phaseID[kspec];
vcs_VolPhase* Vphase = m_VolPhaseList[iph].get();
if (Vphase->m_singleSpecies) {
m_SSPhase[kspec] = true;
} else {
m_SSPhase[kspec] = false;
}
}
}
void VCS_SOLVE::vcs_SSPhase()
{
std::vector<int> numPhSpecies(m_numPhases, 0);
for (size_t kspec = 0; kspec < m_numSpeciesTot; ++kspec) {
numPhSpecies[m_phaseID[kspec]]++;
}
/*
* Handle the special case of a single species in a phase that
* has been earmarked as a multispecies phase.
* Treat that species as a single-species phase
*/
for (size_t iph = 0; iph < m_numPhases; iph++) {
vcs_VolPhase* Vphase = m_VolPhaseList[iph];
Vphase->m_singleSpecies = false;
if (TPhInertMoles[iph] > 0.0) {
Vphase->setExistence(2);
}
if (numPhSpecies[iph] <= 1) {
if (TPhInertMoles[iph] == 0.0) {
Vphase->m_singleSpecies = true;
}
}
}
/*
* Fill in some useful arrays here that have to do with the
* static information concerning the phase ID of species.
* SSPhase = Boolean indicating whether a species is in a
* single species phase or not.
*/
for (size_t kspec = 0; kspec < m_numSpeciesTot; kspec++) {
size_t iph = m_phaseID[kspec];
vcs_VolPhase* Vphase = m_VolPhaseList[iph];
if (Vphase->m_singleSpecies) {
m_SSPhase[kspec] = true;
} else {
m_SSPhase[kspec] = false;
}
}
}
