Still in debugging in multigrid.  About to do some major rewriting to
calculating the coefficients and to dealing with rho in relax and resid.
Lots of print statements are included for debugging purposes.