Assertion error in ForwardProblem::getSimulationStateTimepoint

[dweindl]

I got the following assertion error in ForwardProblem::getSimulationStateTimepoint when simulating Weber_BMC2015 with AMICI 0.26.1:

AMICI/include/amici/forwardproblem.h:200: const amici::SimulationState& amici::ForwardProblem::getSimulationStateTimepoint(int) const: Assertion `map_iter != timepoint_states_.end()' failed.

Needs further investigation.

Seems to be non-deterministic, but happens repeatedly.

Might be some memory corruption issue, potentially caused outside amici.

[dweindl]

This happens if the solver fails right at an output timepoint (here: root-after-reinitialization error). In this case, the model state from that timepoint does not get stored in timepoint_states_ (because handleDataPoint doesn't get called anymore), and the reported problem occurs. In non-debug builds, this usually segfaults when the returned invalid SimulationState is used).

AMICI/src/forwardproblem.cpp

Lines 169 to 196 in 91054b3

	int status = solver->run(next_t_stop);
	/* sx will be copied from solver on demand if sensitivities
	are computed */
	solver->writeSolution(&t_, x_, dx_, sx_, dx_);
	if (status == AMICI_ILL_INPUT) {
	/* clustering of roots => turn off root-finding */
	solver->turnOffRootFinding();
	} else if (status == AMICI_ROOT_RETURN || t_ == next_t_event) {
	// solver-tracked or time-triggered event
	solver->getRootInfo(roots_found_.data());
	// check if we are at a trigger timepoint.
	// if so, set the root-found flag
	if (t_ == next_t_event) {
	for (auto ie : model->state_independent_events_[t_]) {
	// determine direction of root crossing from
	// root function value at the previous event
	roots_found_[ie] = std::copysign(1, -rootvals_[ie]);
	}
	++it_trigger_timepoints;
	}
	handleEvent(&tlastroot_, false, false);
	}
	}
	}
	handleDataPoint(next_t_out);