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

Several fixes due to Solaris port

Browse source
This commit is contained in:
Harry Moffat 2009-02-05 17:10:03 +00:00
parent a1f3936c42
commit 69840fd235
2 changed files with 6 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
Cantera/src/equil/vcs_phaseStability.cpp
View file

@ -108,7 +108,9 @@ namespace VCSnonideal {
for (iph = 0; iph < m_numPhases; iph++) {
Vphase = m_VolPhaseList[iph];
int existence = Vphase->exists();
#ifdef DEBUG_MODE
strcpy(anote, "");
#endif
if (existence > 0) {
#ifdef DEBUG_MODE

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

@ -120,6 +120,8 @@ namespace VCSnonideal {
vcs_VolPhase *Vphase;
double *sc_irxn = NULL; /* Stoichiometric coefficients for cur rxn */
double *dnPhase_irxn;
double atomComp;
int iphasePop;
#ifdef DEBUG_MODE
char ANOTE[128];
/*
@ -450,7 +452,7 @@ namespace VCSnonideal {
* First step is a major branch in the algorithm.
* We first determine if a phase pops into existence.
*/
int iphasePop = vcs_popPhaseID();
iphasePop = vcs_popPhaseID();
/*
*
*/
@ -635,7 +637,7 @@ namespace VCSnonideal {
for (int j = 0; j < m_numElemConstraints; ++j) {
int elType = m_elType[j];
if (elType == VCS_ELEM_TYPE_ABSPOS) {
double atomComp = m_formulaMatrix[j][kspec];
atomComp = m_formulaMatrix[j][kspec];
if (atomComp > 0.0) {
double maxPermissible = m_elemAbundancesGoal[j] / atomComp;
if (maxPermissible < VCS_DELETE_MINORSPECIES_CUTOFF) {