Change the way convergence is calculated.

 * Create a mask using the coarsest grid.
 * Calculate error at the points on the coarsest grid.
 * Use four grids, calculate convergence rate using the three finest.
Use PyAMG to solve the pressure system.