implement kappa plotting in pde solver

mlmc/src/pdesolver/PDESolver.cpp

@@ -60,10 +60,10 @@ void EllipticPDESolver::plotSolution(SampleSolution &solution) {
   assemble->SetFlux(solution.U, flux.U);
   mpp::plot(flux.IdString()) << flux.U << mpp::endp;
 
-//  LagrangeDiscretization kappaDisc(meshes, 0, 0);
-//  SampleSolution kappa(&kappaDisc, solution.Id(), "Kappa");
-//  assemble->GetProblem()->Permeability(kappa.U);
-//  mpp::plot(kappa.IdString()) << kappa.U << mpp::endp;
+  LagrangeDiscretization kappaDisc(meshes, 0, 1);
+  SampleSolution kappa(&kappaDisc, solution.Id(), "Kappa");
+  assemble->GetProblem()->Permeability(kappa.U);
+  mpp::plot(kappa.IdString()) << kappa.U << mpp::endp;
 
   delete pressureDisc;
 }