removed empiric measures

mlmc/src/montecarlo/MonteCarlo.cpp

@@ -6,15 +6,6 @@ void MonteCarlo::Method() {
   vout(1) << "Start with: " << aggregate;
   method();
   aggregate.UpdateParallel();
-
-  vout(2) << sums;
-  avgs.Update(sums, aggregate.ctr.M);
-  vars.Update(avgs);
-  kurtosis.Update(avgs, vars);
-
-  aggregate.SVarQ = vars.Q;
-  aggregate.SVarY = vars.Y;
-
   vout(1) << "End with: " << aggregate;
   mout.EndBlock(verbose == 0, "Took");
 }
@@ -22,14 +13,12 @@ void MonteCarlo::Method() {
 void MonteCarlo::method() {
   SampleSolution cSolution(pdeSolver->GetDisc(), cId);
   SampleSolution fSolution(pdeSolver->GetDisc(), fId);
+
   while (aggregate.ctr.dMcomm != 0) {
     computeSampleSolution(aggregate.index(), fId, fSolution);
     if (onlyFine) cSolution.Init();
     else computeSampleSolution(aggregate.index(), cId, cSolution);
     aggregate.Update(fSolution, cSolution);
-
-    sums.Update(fSolution, cSolution); // todo rmv
-    vout(3) << sums; // todo rmv
   }
 }
 

mlmc/src/montecarlo/MonteCarlo.hpp

@@ -14,31 +14,23 @@ protected:
 
   int plotting = 0;
 
+  int level;
+
   bool onlyFine;
 
   bool parallel;
 
-  int level;
-
   Meshes *meshes;
 
   PDESolver *pdeSolver;
 
 public:
-  WelfordAggregate aggregate;
-
-  Sums sums;
-
-  Averages avgs;
-
-  Variances vars;
-
-  Kurtosis kurtosis;
-
   PDESolverCreator pdeSolverCreator;
 
   MeshesCreator meshesCreator;
 
+  WelfordAggregate aggregate;
+
   SampleID cId;
 
   SampleID fId;