Python/numpy implementation of Gauss-Sidel, based on Ed Bueler's Obstacle Problem solver at: http://www.dms.uaf.edu/~bueler/obstacle.m