worked on problem

mlmc/src/problem/StochasticEllipticProblem.cpp

@@ -3,12 +3,20 @@
 
 template<>
 Tensor StochasticLaplace1DT<Vector>::Permeability(const cell &c) const {
-    mout << *this->fineSample << endl;
-    mout << this->fineSample->size() << endl;
-    return this->fineSample->operator()(c(), 0) * One;
+    if (fineSample && !coarseSample)
+        return this->fineSample->operator()(c(), 0) * One;
+    else if (coarseSample && !fineSample)
+        return this->coarseSample->operator()(c(), 0) * One;
+    else
+        Exit("Check generators")
 }
 
 template<>
 Tensor StochasticLaplace2DT<Vector>::Permeability(const cell &c) const {
-    return this->fineSample->operator()(c(), 0) * One;
+    if (fineSample && !coarseSample)
+        return this->fineSample->operator()(c(), 0) * One;
+    else if (coarseSample && !fineSample)
+        return this->coarseSample->operator()(c(), 0) * One;
+    else
+    Exit("Check generators")
 }
\ No newline at end of file