leftover commit

mlmc/conf/mlmc.conf

@@ -15,9 +15,9 @@ flux_alpha = 1  #Upwind: 1, Central: 0
 #Problem = StochasticLaplace2D
 #Problem = StochasticPollution2D
 #Problem = DeterministicPollution2D
-#Problem = StochasticPollutionCosHat2D
+Problem = StochasticPollutionCosHat2D
 #Problem = DeterministicPollutionCosHat2D
-Problem = DeterministicPollutionCosHatSquare500
+#Problem = DeterministicPollutionCosHatSquare500
 
 StochasticField = LogNormal
 #StochasticField = Gaussian
@@ -26,14 +26,14 @@ StochasticField = LogNormal
 Experiment = MLMCExperiment
 #Experiment = MLMCOverEpsilon
 
-initLevels = 2, 3
-initSampleAmount = 1, 0
+initLevels = 4, 5, 6
+initSampleAmount = 0, 0, 1
 
 epsilon = 0.01
 mcOnly = true
 epsilon_lst = 0.05, 0.01, 0.005
 
-maxLevel = 4
+maxLevel = 6
 plevel = 1
 degree = 2
 
@@ -87,7 +87,7 @@ NewtonVerbose = -1
 DebugLevel = -1
 
 #Transport
-scaling = 16
+scaling = 8
 TimeSeries = uniform
 startTime = 0.0
 endTime = 1.0

tools/plot_mlmc.py

@@ -15,7 +15,7 @@ L2 = "--------------------------------------------------------------------------
 plt.rcParams.update({'font.size': 14})
 
 
-def read_log(_file_name_):
+def read_log(_file_name_=None):
     log_file = '../build/log/log' if _file_name_ is None else _file_name_
     with open(log_file) as file:
         content = file.readlines()
@@ -113,15 +113,5 @@ def plot_statistics(dfs, labels, plot_name=None):
 if __name__ == "__main__":
     argv = sys.argv
     argv1 = argv[1] if len(argv) == 2 else None
-
-    # results14 = read_log('../results/mlmc_alpha/14')
-    # plot_statistics(results14, None, argv1)
-    # results18 = read_log('../results/mlmc_alpha/18')
-    # plot_statistics(results18, None, argv1)
-
-    resultslin = read_log('../results/mlmc_disc/LinearFEM')
-    resultsquad = read_log('../results/mlmc_disc/QuadraticFEM')
-    resultshybr = read_log('../results/mlmc_disc/HybridFEM')
-    results = [resultslin, resultsquad, resultshybr]
-    labels = ['Linear FEM', 'Quadratic FEM', 'Hybrid FEM']
-    plot_statistics(results, labels, argv1)
+    results = read_log()
+    plot_content(results)