added PDESolver to MLMC

mlmc/src/montecarlo/MultilevelMonteCarlo.hpp

@@ -19,6 +19,7 @@ public:
 
     Meshes *meshes;
     IStochasticAssemble *assemble;
+    PDESolver *pdeSolver;
 
     vector<int> levels = {};
     vector<int> samples = {};
@@ -41,11 +42,13 @@ public:
     MultilevelMonteCarlo(vector<int> &initLevels,
                          vector<int> &initSampleAmount,
                          Meshes *meshes,
-                         IStochasticAssemble *assemble) :
+                         IStochasticAssemble *assemble,
+                         PDESolver *pdeSolver) :
         initLevels(initLevels),
         initSampleAmount(initSampleAmount),
         meshes(meshes),
-        assemble(assemble) {
+        assemble(assemble),
+        pdeSolver(pdeSolver) {
 
         config.get("mcOnly", mcOnly);
         config.get("MLMCVerbose", verbose);