diff --git a/deps/AMICI/include/amici/forwardproblem.h b/deps/AMICI/include/amici/forwardproblem.h index 67658ebe..b2b33a21 100644 --- a/deps/AMICI/include/amici/forwardproblem.h +++ b/deps/AMICI/include/amici/forwardproblem.h @@ -9,6 +9,7 @@ #include #include +#include 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; };