prototype implementation

AMICI-dev · May 17, 2024 · 601aa75 · 601aa75
1 parent cc620e1
commit 601aa75
Show file tree

Hide file tree

Showing 3 changed files with 17 additions and 0 deletions.
diff --git a/include/amici/backwardproblem.h b/include/amici/backwardproblem.h
@@ -70,6 +70,12 @@ class BackwardProblem {
      */
     AmiVector const& getAdjointState() const { return xB_; }
 
+    /**
+     * @brief Accessor for dxB
+     * @return xB
+     */
+    AmiVector const& getAdjointStateDerivative() const { return dxB_; }
+
     /**
      * @brief Accessor for xQB
      * @return xQB

diff --git a/include/amici/edata.h b/include/amici/edata.h
@@ -445,6 +445,12 @@ class ExpData : public SimulationParameters {
         std::vector<realtype> const& input, char const* fieldname
     ) const;
 
+    /**
+     * @brief check whether this instance implements only equilibration
+     * without simulation
+     */
+    bool isEquilibration() const;
+
     /** @brief number of observables */
     int nytrue_{0};
 

diff --git a/src/edata.cpp b/src/edata.cpp
@@ -386,6 +386,11 @@ void ExpData::checkEventsDimension(
         );
 }
 
+bool ExpData::isEquilibration() const {
+    return ts_.size() == 1 && std::isinf(ts_[0]);
+}
+
+
 void checkSigmaPositivity(
     std::vector<realtype> const& sigmaVector, char const* vectorName
 ) {