using HybridRTLagrangeDiscretization

mlmc/src/pdesolver/PDESolverCreator.cpp

@@ -21,7 +21,7 @@ PDESolver *PDESolverCreator::Create(Meshes &meshes) {
     if (_model == "HybridElliptic")
         return new EllipticPDESolver(
             new HybridEllipticAssemble(
-                new RTLagrangeDiscretization(meshes, 0, 0),
+                new HybridRTLagrangeDiscretization(meshes, 0, 0),
                 CreateStochasticEllipticProblem(_problem, meshes)
             ), meshes, _quantity, _costMeasure
         );