From 3baa0eae9439f213420bcfe9690cf5fd9b660ce0 Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Sun, 30 Dec 2012 00:16:56 +0000 Subject: [PATCH] [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. --- src/oneD/Domain1D.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/oneD/Domain1D.cpp b/src/oneD/Domain1D.cpp index ee6026f76..0629905d4 100644 --- a/src/oneD/Domain1D.cpp +++ b/src/oneD/Domain1D.cpp @@ -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(jpt-1, 0); + jmin = std::max(jpt, 1) - 1; jmax = std::min(jpt+1,m_points-1); }