Make stateUpdaters cloneable

src/state/dieterichstateupdater.hh

@@ -16,6 +16,8 @@ class DieterichStateUpdater
   virtual void solve(VectorType const &velocity_field);
   virtual void extractState(SingletonVectorType &state);
 
+  virtual StateUpdater<SingletonVectorType, VectorType> *clone();
+
 private:
   SingletonVectorType alpha_old;
   SingletonVectorType alpha;
@@ -57,4 +59,10 @@ void DieterichStateUpdater<SingletonVectorType, VectorType>::extractState(
   state = alpha;
 }
 
+template <class SingletonVectorType, class VectorType>
+StateUpdater<SingletonVectorType, VectorType> *
+DieterichStateUpdater<SingletonVectorType, VectorType>::clone() {
+  return new DieterichStateUpdater<SingletonVectorType, VectorType>(*this);
+}
+
 #endif

src/state/ruinastateupdater.hh

@@ -15,6 +15,8 @@ class RuinaStateUpdater : public StateUpdater<SingletonVectorType, VectorType> {
   virtual void solve(VectorType const &velocity_field);
   virtual void extractState(SingletonVectorType &state);
 
+  virtual StateUpdater<SingletonVectorType, VectorType> *clone();
+
 private:
   SingletonVectorType alpha_old;
   SingletonVectorType alpha;
@@ -55,4 +57,10 @@ void RuinaStateUpdater<SingletonVectorType, VectorType>::extractState(
   state = alpha;
 }
 
+template <class SingletonVectorType, class VectorType>
+StateUpdater<SingletonVectorType, VectorType> *
+RuinaStateUpdater<SingletonVectorType, VectorType>::clone() {
+  return new RuinaStateUpdater<SingletonVectorType, VectorType>(*this);
+}
+
 #endif

src/state/stateupdater.hh

@@ -7,6 +7,8 @@ template <class SingletonVectorType, class VectorType> class StateUpdater {
   virtual void setup(double _tau) = 0;
   virtual void solve(VectorType const &velocity_field) = 0;
   virtual void extractState(SingletonVectorType &state) = 0;
+
+  virtual StateUpdater<SingletonVectorType, VectorType> *clone() = 0;
 };
 
 #endif