Clean up base cases (and also writeHeatMapBMP)

Base case for multi grid is a computeS.
base case for fullMG calls multiGrid instead of approxOp