Add Solver::get_current_time

Prerequisite for #2246.

include/amici/solver.h

@@ -964,6 +964,12 @@ class Solver {
      */
     int getMaxConvFails() const;
 
+    /**
+     * @brief Get the current time of the solver.
+     * @return internal time of the solver
+     */
+    virtual realtype get_current_time() const = 0;
+
     /**
      * @brief Serialize Solver (see boost::serialization::serialize)
      * @param ar Archive to serialize to

include/amici/solver_cvodes.h

@@ -91,6 +91,8 @@ class CVodeSolver : public Solver {
 
     void setNonLinearSolverB(int which) const override;
 
+    realtype get_current_time() const override;
+
   protected:
     void calcIC(realtype tout1) const override;
 

include/amici/solver_idas.h

@@ -106,6 +106,8 @@ class IDASolver : public Solver {
 
     void setNonLinearSolverB(int which) const override;
 
+    realtype get_current_time() const override;
+
   protected:
     /**
      * @brief Postprocessing of the solver memory after a discontinuity

src/solver_cvodes.cpp

@@ -386,6 +386,17 @@ void CVodeSolver::setNonLinearSolverB(int const which) const {
         throw CvodeException(status, "CVodeSetNonlinearSolverB");
 }
 
+realtype CVodeSolver::get_current_time() const {
+    if (!solver_memory_)
+        return std::numeric_limits<realtype>::quiet_NaN();
+
+    realtype time;
+    int status = CVodeGetCurrentTime(solver_memory_.get(), &time);
+    if (status != CV_SUCCESS)
+        throw CvodeException(status, "CVodeGetCurrentTime");
+    return time;
+}
+
 void CVodeSolver::setErrHandlerFn() const {
     int status = CVodeSetErrHandlerFn(
         solver_memory_.get(), wrapErrHandlerFn,

src/solver_idas.cpp

@@ -876,6 +876,17 @@ void IDASolver::setNonLinearSolverB(int which) const {
         throw CvodeException(status, "CVodeSetNonlinearSolverB");
 }
 
+realtype IDASolver::get_current_time() const {
+    if (!solver_memory_)
+        return std::numeric_limits<realtype>::quiet_NaN();
+
+    realtype time;
+    int status = IDAGetCurrentTime(solver_memory_.get(), &time);
+    if (status != IDA_SUCCESS)
+        throw CvodeException(status, "IDAGetCurrentTime");
+    return time;
+}
+
 /**
  * @brief Jacobian of xdot with respect to states x
  * @param N number of state variables