add initialization identifier

include/amici/edata.h

@@ -398,10 +398,15 @@ class ExpData : public SimulationParameters {
     void clear_observations();
 
     /**
-     * @brief Arbitrary (not necessarily unique) identifier.
+     * @brief Unique identifier.
      */
     std::string id;
 
+    /**
+     * @brief Unique Identifier of the ExpData to use for initialization.
+     */
+    std::string initialization_id;
+
   protected:
     /**
      * @brief resizes observedData, observedDataStdDev, observedEvents and
@@ -480,6 +485,7 @@ inline bool operator==(ExpData const& lhs, ExpData const& rhs) {
     return *dynamic_cast<SimulationParameters const*>(&lhs)
                == *dynamic_cast<SimulationParameters const*>(&rhs)
            && lhs.id == rhs.id && lhs.nytrue_ == rhs.nytrue_
+           && lhs.initialization_id == rhs.initialization_id
            && lhs.nztrue_ == rhs.nztrue_ && lhs.nmaxevent_ == rhs.nmaxevent_
            && is_equal(lhs.observed_data_, rhs.observed_data_)
            && is_equal(lhs.observed_data_std_dev_, rhs.observed_data_std_dev_)

include/amici/rdata.h

@@ -95,10 +95,15 @@ class ReturnData : public ModelDimensions {
         Model& model, Solver const& solver, ExpData const* edata
     );
     /**
-     * @brief Arbitrary (not necessarily unique) identifier.
+     * @brief Unique identifier.
      */
     std::string id;
 
+    /**
+     * @brief Unique Identifier of the ExpData to use for initialization.
+     */
+    std::string initialization_id;
+
     /**
      * timepoints (shape `nt`)
      */

src/edata.cpp

@@ -126,6 +126,7 @@ ExpData::ExpData(
     }
 
     id = rdata.id;
+    initialization_id = rdata.initialization_id;
 }
 
 void ExpData::setTimepoints(std::vector<realtype> const& ts) {