[1D] Fixed an indexing error that affected Jacobian evaluations

Not all components of the Jacobian were being computed for the first point in
the first flow Domain. Cherry pick of r2018.
This commit is contained in:
Ray Speth 2013-01-22 20:28:18 +00:00
parent 9656f3aa07
commit 7e065242d8

View file

@ -102,7 +102,7 @@ eval(size_t jg, doublereal* xg, doublereal* rg,
jmin = 0;
jmax = m_points - 1;
} else { // evaluate points for Jacobian
jmin = std::max<size_t>(jpt-1, 0);
jmin = std::max<size_t>(jpt, 1) - 1;
jmax = std::min(jpt+1,m_points-1);
}