On the direct solution of Poisson's equation on a non-uniform grid (Q1208847)

The direct solution of systems of linear equations arising from the discretization of Poisson's equation with Neumann boundary conditions is considered under the assumption that the nonuniform grid allows a one- dimensional transformation from the physical domain into a uniform grid in a computational rectangular domain. Hence the system matrix is of block tridiagonal form. Under the additional assumptions that the above mentioned transformation has the property of odd symmetry and that the order of the block matrices as well as their number are odd it is first shown how the total matrix can be symmetrized. In a second step an extension of the method of matrix decomposition is applied thus reducing the given problem to a set of tridiagonal systems of linear equations by means of the eigenproblem of a pair of symmetric matrices.