Multigrid solver implemented, convergence is there, parameters not optimized