Fixed id coarse-to-fine conversion, added eye-ball norm, and added comparison of project_vector to interpolation operator
in the case of a uniform mesh. For the latter, there must be at least 3 coarsenings performed.
Added in manual Petsc solve. It does not yet check the solution against the actual solution.
Cleaned up the code to give a better performance show. Note that rediscretization is not working as of yet, and so MG uses Galerkin.
Fixed MPI bug where nz,nnz entries were both zero in some places.
In progress of fixing issue of penalty-based boundary conditions in matrix setup.