bug fix flux computation

mlmc/src/pdesolver/assembling/HybridEllipticAssemble.cpp

@@ -214,7 +214,7 @@ void HybridEllipticAssemble::SetFlux(const Vector &u, Vector &flux) {
             }
         }
         F *= (1 / area);
-        for (int d = 0; d < 3; ++d)
+        for (int d = 0; d < SpaceDimension; ++d)
             flux(c(), d) = F[d];
     }
 }

mlmc/src/pdesolver/assembling/LagrangeEllipticAssemble.cpp

@@ -329,7 +329,7 @@ void LagrangeEllipticAssemble::SetFlux(const Vector &u, Vector &flux) {
             area += w;
         }
         F *= (1 / area);
-        for (int d = 0; d < 3; ++d) {
+        for (int d = 0; d < SpaceDimension; ++d) {
             flux(c(), d) = F[d];
         }
     }

mlmc/src/pdesolver/assembling/MixedEllipticAssemble.cpp

@@ -355,7 +355,7 @@ void MixedEllipticAssemble::SetFlux(const Vector &u, Vector &flux) {
             area += w;
         }
         F *= (1 / area);
-        for (int d = 0; d < 3; ++d) {
+        for (int d = 0; d < SpaceDimension; ++d) {
             flux(c(), d) = F[d];
         }
     }