A fourth-order iterative solver for the singular poisson equation

Abstract

A compact fourth-order finite difference scheme solver devoted to the singular-Poisson equation is proposed and verified. The solver is based on a mixed formulation: the Poisson equation is splitted into a system of partial differential equations of the first order. This system is then discretized using a fourth-order compact scheme. This leads to a sparse linear system but introduces new variables related to the gradient of an unknow function. The Schur factorization allows us to work on a linear sub-problem for which a conjugated-gradient preconditioned by an algebraic multigrid method is proposed.Numerical results show that the new proposed Poisson solver is efficient while retaining the fourth-order compact accuracy. © 2013 Springer-Verlag.