Fix: Don't access AMICI rdata->x with non-full reporting mode

ICB-DCM · Jul 27, 2020 · 27fe478 · 27fe478
1 parent b2a0035
commit 27fe478
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/parpeamici/multiConditionProblem.cpp b/src/parpeamici/multiConditionProblem.cpp
@@ -415,7 +415,7 @@ AmiciSimulationRunner::AmiciResultPackageSimple runAndLogSimulation(
              * better exception handling, we check those fields to deduce where
              * the error occurred
              */
-            bool forwardFailed = std::isnan(rdata->x[rdata->x.size() - 1]);
+            bool forwardFailed = std::isnan(rdata->llh);
             bool backwardFailed = solver->getSensitivityOrder() >= amici::SensitivityOrder::first
                     && solver->getSensitivityMethod() == amici::SensitivityMethod::adjoint
                     && !rdata->sllh.empty() && std::isnan(rdata->sllh[0]);