Test amici #2491

ICB-DCM · Sep 24, 2024 · b053270 · b053270
1 parent 77df243
commit b053270
Showing 1 changed file with 10 additions and 0 deletions.
diff --git a/deps/AMICI/include/amici/forwardproblem.h b/deps/AMICI/include/amici/forwardproblem.h
@@ -9,6 +9,7 @@
 
 #include <sundials/sundials_direct.h>
 #include <vector>
+#include <iostream>
 
 namespace amici {
 
@@ -197,6 +198,15 @@ class ForwardProblem {
         if (model->getTimepoint(it) == initial_state_.t)
             return getInitialSimulationState();
         auto map_iter = timepoint_states_.find(model->getTimepoint(it));
+        if (map_iter == timepoint_states_.end()) {
+            std::cerr << "Simulation state not found for timepoint #"
+                      << it << " " << model->getTimepoint(it) << std::endl;
+            std::cerr << "Available timepoints: ";
+            for (auto const& [t, state] : timepoint_states_)
+                std::cerr << t << " ";
+            std::cerr << std::endl;
+        }
+
         assert(map_iter != timepoint_states_.end());
         return map_iter->second;
     };