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

moved declaration of nCutoff to avoid jumping over it, causing a compile failure with gcc 4.0.1 on linuxChemEquil.cpp

Browse source
This commit is contained in:
Dave Goodwin 2006-10-23 00:56:31 +00:00
parent ab15545caf
commit d4672515bd
2 changed files with 7 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
Cantera/src/ChemEquil.cpp
View file

@ -1237,6 +1237,8 @@ namespace Cantera {
}
#endif
double nCutoff;
bool normalStep = true;
/*
* Decide if we are to do a normal step or a modified step
@ -1304,8 +1306,8 @@ namespace Cantera {
* The logic below looks for obvious cases where the current element
* coefficient matrix is rank deficient.
*
* The way around rank-deficiency is to lump-sum the corresponding row of the
* matrix. Note, lump-summing seems to work very well in terms of
* The way around rank-deficiency is to lump-sum the corresponding row
* of the matrix. Note, lump-summing seems to work very well in terms of
* its stability properties, i.e., it heads in the right direction,
* albeit with lousy convergence rates.
*
@ -1317,7 +1319,8 @@ namespace Cantera {
for (m = 0; m < m_mm; m++) {
lumpSum[m] = 1;
}
double nCutoff = 1.0E-4 * n_t_calc;
nCutoff = 1.0E-4 * n_t_calc;
#ifdef DEBUG_HKM_EPEQUIL
if (debug_prnt_lvl > 0) {
printf(" Lump Sum Elements Calculation: \n");

1
Cantera/src/MultiPhaseEquil.cpp
View file

@ -730,6 +730,7 @@ namespace Cantera {
for (k = 0; k < m_nsp; k++) {
kc = m_species[k];
if (m_mix->speciesPhaseIndex(kc) == ip) {
// bug fixed 7/12/06 DGG
stoich = nu[k]; // nu[kc];
psum += stoich * stoich;
}