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

Fix indexing error in ImplicitSurfChem::eval

Browse source 
This method was not working in the case where there was more than one surface
phase, since 'loc' was always set to zero.
This commit is contained in:
Ray Speth 2015-11-24 13:31:45 -05:00
parent e805d4d9ad
commit 3043a5d872
1 changed files with 1 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
src/kinetics/ImplicitSurfChem.cpp
View file

@ -155,13 +155,12 @@ void ImplicitSurfChem::eval(doublereal time, doublereal* y,
{
updateState(y); // synchronize the surface state(s) with y
doublereal rs0, sum;
size_t loc, kstart;
size_t loc = 0, kstart;
for (size_t n = 0; n < m_nsurf; n++) {
rs0 = 1.0/m_surf[n]->siteDensity();
m_vecKinPtrs[n]->getNetProductionRates(DATA_PTR(m_work));
kstart = m_vecKinPtrs[n]->kineticsSpeciesIndex(0,m_surfindex[n]);
sum = 0.0;
loc = 0;
for (size_t k = 1; k < m_nsp[n]; k++) {
ydot[k + loc] = m_work[kstart + k] * rs0 * m_surf[n]->size(k);
sum -= ydot[k];