PETSc DM/multigrid memory management is broken

[Prof Garth Wells]

There is something wrong the memory management of PETSc DM objects, which are used in the geometric multigrid code.

A symptom is that the test test_mg_solver_laplace in test/python/la/test_mg_solver.py frequently segfaults.