start adding code for C^{-1} = [S+N]^{-1} on flat-sky maps.