workaround for dune-functions problems

src/pfc.cc

@@ -103,7 +103,8 @@ int main(int argc, char** argv)
     Mu.getVector() = 0.0;
     Nu.getVector() = 0.0;
     std::srand( time(0) );
-    Psi.interpol([psi_mean](auto const& x) { return (std::rand()/double(RAND_MAX) - 0.5) + psi_mean; });
+//     Psi.interpol([psi_mean](auto const& x) { return (std::rand()/double(RAND_MAX) - 0.5) + psi_mean; });
+    Psi = psi_mean;
     
 //     using PreconType = PfcPrecon<typename PfcProblem::SystemMatrixType::MultiMatrix, typename PfcProblem::SystemVectorType::MultiVector>;
 //     PreconType precon(prob.getSystemMatrix().getMatrix(), &tau, M0);