Poisson equation discretized into block tridiagonal