get initial rootvals

dweindl · Dec 5, 2023 · 36f98c7 · 36f98c7
1 parent d8e140c
commit 36f98c7
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 1 deletion.
diff --git a/include/amici/forwardproblem.h b/include/amici/forwardproblem.h
@@ -7,7 +7,6 @@
 #include "amici/vector.h"
 #include <amici/amici.h>
 
-#include <memory>
 #include <sundials/sundials_direct.h>
 #include <vector>
 

diff --git a/src/forwardproblem.cpp b/src/forwardproblem.cpp
@@ -130,6 +130,8 @@ void ForwardProblem::workForwardProblem() {
 
     /* store initial state and sensitivity*/
     initial_state_ = getSimulationState();
+    // store root information at t0
+    model->froot(t_, x_, dx_, rootvals_);
 
     // get list of trigger timepoints for fixed-time triggered events
     auto trigger_timepoints = model->get_trigger_timepoints();