parallel runnable

mlmc/conf/mlmc.conf

@@ -18,12 +18,12 @@ Problem = StochasticPollution2D
 StochasticField = LogNormal
 #StochasticField = Gaussian
 
-Experiment = ConvergenceTest
-#Experiment = MLMCExperiment
+#Experiment = ConvergenceTest
+Experiment = MLMCExperiment
 #Experiment = MLMCOverEpsilon
 
-initLevels = 3, 4
-initSampleAmount = 0, 1
+initLevels = 3, 4, 5
+initSampleAmount = 10, 5, 2
 
 #l_init = 3,  4,  5
 #M_init = 45, 15, 5

mlmc/src/stochastics/HybridFluxGenerator.C

@@ -25,7 +25,7 @@ void HybridFluxGenerator2D::GetFineSample(Vector &fineField) {
 
     assemble.problem->LoadNewSample(&perm);
 
-    Preconditioner *pc = GetPC(faceMatrixGraph, assemble, "LIB_PS");
+    Preconditioner *pc = GetPC(faceMatrixGraph, assemble, "SuperLU");
 
     Solver solver(pc, "GMRES");