adapted MatrixGraphs constructor

mlmc/src/generators/algorithms/CirculantEmbedding.hpp

@@ -48,7 +48,7 @@ public:
     explicit CirculantEmbedding(Meshes &meshes) :
         meshes(meshes),
         normalDist(meshes),
-        cellMGraphs(meshes, dof(new LagrangeDoF(0))) {
+        cellMGraphs(meshes, new LagrangeDoF(0)) {
 
         config.get("evtol", evtol);
         config.get("StochasticField", fieldType);
@@ -220,4 +220,4 @@ public:
     string Name() const override { return "Circulant Embedding"; }
 };
 
-#endif //M_CIRCULANTEMBEDDING_H
\ No newline at end of file
+#endif //M_CIRCULANTEMBEDDING_H

mlmc/src/generators/algorithms/HybridFluxGenerator.hpp

@@ -67,8 +67,8 @@ public:
 
 
 
-//    MatrixGraphs cellMGraphs;         cellMGraphs(MatrixGraphs(meshes, dof(new CellDoF(3)))),
-//    MatrixGraphs faceMGraphs;         faceMGraphs(MatrixGraphs(meshes, dof(new FaceDoF(1)))),
+//    MatrixGraphs cellMGraphs;         cellMGraphs(MatrixGraphs(meshes, new CellDoF(3))),
+//    MatrixGraphs faceMGraphs;         faceMGraphs(MatrixGraphs(meshes, new FaceDoF(1))),
 //
 //    Vector *faceFlux = nullptr;
 //    Vector *faceValues = nullptr;

mlmc/src/pdesolver/PDESolver.cpp

@@ -12,7 +12,7 @@ MatrixGraphs *PDESolver::CreateSolutionMatrixGraphs(Meshes &meshes) {
         auto assemble = this->GetAssemble();
         auto assembleType = typeid(*assemble).name();
         if (assembleType == typeid(HybridEllipticAssemble).name())
-            return new MatrixGraphs(meshes, dof(new FaceDoF(1)));
+            return new MatrixGraphs(meshes, new FaceDoF(1));
         if (assembleType == typeid(MixedEllipticAssemble).name())
             return new MatrixGraphs(meshes, *disc);
     }
@@ -65,11 +65,11 @@ void EllipticPDESolver::plotSolution(SampleSolution &solution) {
 
 void EllipticPDESolver::createOtherMatrixGraphs(Meshes &meshes) {
     if (plotting) {
-        fluxMGraphs = new MatrixGraphs(meshes, dof(new LagrangeDoF(0, 3)));
+        fluxMGraphs = new MatrixGraphs(meshes, new LagrangeDoF(0, 3));
         if (typeid(*assemble).name() == typeid(HybridEllipticAssemble).name())
-            pressureMGraphs = new MatrixGraphs(meshes, dof(new LagrangeDoF(0, 1)));
+            pressureMGraphs = new MatrixGraphs(meshes, new LagrangeDoF(0, 1));
         else
-            pressureMGraphs = new MatrixGraphs(meshes, dof(new LagrangeDoF(1, 1)));
+            pressureMGraphs = new MatrixGraphs(meshes, new LagrangeDoF(1, 1));
     } else {
         fluxMGraphs = nullptr;
         pressureMGraphs = nullptr;

mlmc/tests/basics/TestPlotMap.cpp

@@ -42,7 +42,7 @@ protected:
             new MatrixGraphs(*meshes, *(new LagrangeDiscretization(*meshes, 3))),
 
             new MatrixGraphs(*meshes, *(new RTLagrangeDiscretization(*meshes, 0, 0))),
-            new MatrixGraphs(*meshes, dof(new FaceDoF(1))),
+            new MatrixGraphs(*meshes, new FaceDoF(1)),
 
             new CellMatrixGraphs(*meshes, *(new DGDiscretization(*meshes, 0))),
             new CellMatrixGraphs(*meshes, *(new DGDiscretization(*meshes, 1))),