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

Took the line search algorithm out , and put it into an ifdef block.

Browse source 
It's not needed.
This commit is contained in:
Harry Moffat 2008-06-18 22:23:05 +00:00
parent 7a76bf5444
commit ac0aa34b32
1 changed files with 4 additions and 17 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -915,6 +915,8 @@ namespace VCSnonideal {
}
}
}
#ifdef VCS_LINE_SEARCH
/*********************************************************************/
/*** LINE SEARCH ALGORITHM FOR MAJOR SPECIES IN NON-IDEAL PHASES *****/
/*********************************************************************/
@ -926,33 +928,18 @@ namespace VCSnonideal {
(doPhaseDeleteIph == -1) &&
(m_speciesUnknownType[kspec] != VCS_SPECIES_TYPE_INTERFACIALVOLTAGE)) {
double dx_old = dx;
#ifdef DEBUG_HKM_NOT
std::vector<double> feSpecies_tmp(m_numSpeciesTot);
feSpecies_tmp = m_feSpecies_old;
std::vector<double> molNumSpecies_tmp(m_numSpeciesTot);
molNumSpecies_tmp = m_molNumSpecies_old;
#endif
#ifdef DEBUG_MODE
dx = vcs_line_search(irxn, dx_old, ANOTE);
#else
dx = vcs_line_search(irxn, dx_old);
#endif
#ifdef DEBUG_HKM_NOT
for (int kk = 0; kk < m_numSpeciesTot; kk++) {
if (feSpecies_tmp[kk] != m_feSpecies_old[kk]) {
printf("we are here\n");
}
if (molNumSpecies_tmp[kk] != m_molNumSpecies_old[kk]) {
printf("we are here\n");
}
}
#endif
vcs_setFlagsVolPhases(false, VCS_STATECALC_NEW);
}
m_deltaMolNumSpecies[kspec] = dx;
#endif
} /* End of Loop on ic[irxn] -> the type of species */
/***********************************************************************/
/****** CALCULATE KMOLE NUMBER CHANGE FOR THE COMPONENT BASIS **********/
/***********************************************************************/