Analyses on the finite difference method by Gibou et al. for Poisson equation

Gibou et al. in [4] introduced a finite difference method for solving the Poisson equation in irregular domains with the Dirichlet boundary condition. Contrary to its great importance, its properties have not been mathematically analyzed, but have just been numerically observed. In this article, we present two analyses for the method. One proves that its solution is second order accurate, and the other estimates the condition number of its linear system. According to our estimation, the condition number of the unpreconditioned linear system is of size O(1/(h⋅hmin)), and each of Jacobi, SGS, and ILU preconditioned systems is of size O(h−2). Furthermore, our analysis shows that the condition number of MILU is of size O(h−1), the most successful one.