Add fix for boundaries outside of theta values

nlmixr2 · Nov 12, 2024 · f026e1d · f026e1d
1 parent 30969e7
commit f026e1d
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/src/inner.cpp b/src/inner.cpp
@@ -5278,11 +5278,15 @@ NumericMatrix foceiCalcCov(Environment e){
               // Upper only
               checkUpperBound = true;
             }
-            if (checkLowerBound && (std::fabs((cur-op_focei.lower[k])/cur) < op_focei.boundTol)) {
+            if (checkLowerBound &&
+                j < thetaNames.size() &&
+                (std::fabs((cur-op_focei.lower[k])/cur) < op_focei.boundTol)) {
               boundary = true;
               boundStr += "\"" + thetaNames[j] + "\" ";
             }
-            if (checkUpperBound && (std::fabs((op_focei.upper[k]-cur)/cur) < op_focei.boundTol)) {
+            if (checkUpperBound &&
+                j < thetaNames.size() &&
+                (std::fabs((op_focei.upper[k]-cur)/cur) < op_focei.boundTol)) {
               boundary = true;
               boundStr += "\"" + thetaNames[j] + "\" ";
             }