Multigrid and preconditioning strategies for implicit PDE solvers for degenerate parabolic equations