A Simple Proof of Gustafsson’s Conjecture in Case of Poisson Equation on Rectangular Domains
- 1 Institute of Mathematical Sciences, Ewha Womans University, Seoul, South Korea
- 2 Department of Mathematics, Ewha Womans University, Seoul, South Korea
Abstract
We consider the standard five-point finite difference method for solving the Poisson equation with the Dirichlet boundary condition. Its associated matrix is a typical ill-conditioned matrix whose size of the condition number is as big as . Among ILU, SGS, modified ILU (MILU) and other ILU-type preconditioners, Gustafson shows that only MILU achieves an enhancement of the condition number in different order as . His seminal work, however, is not for the MILU but for a perturbed version of MILU and he observes that without the perurbation, it seems to reach the same result in practice. In this work, we give a simple proof of Gustafsson's conjecture on the unnecessity of perturbation in case of Poisson equation on rectangular domains. Using the Cuthill-Mckee ordering, we simplify the recursive equation in two dimensional grid nodes into a recursive one in the level that is one-dimensional. Due to the simplification, our proof is easy to follow and very short.
- Dupont, T., Kendall, R. and Rachford, H. (1968) An Approximate Factorization Procedure for Solving Self-Adjoint Elliptic Difference Equations. SIAM Journal on Numerical Analysis, 5, 559-573. http://dx.doi.org/10.1137/0705045
- Gustafsson, I. (1978) A Class of First Order Factorization Methods. BIT Numerical Mathematics, 18, 142-156. http://dx.doi.org/10.1007/BF01931691
- Greenbaum, A. (1997) Iterative Methods for Solving Linear Systems. SIAM, Philadelphia. http://dx.doi.org/10.1137/1.9781611970937
- Greenbaum, A. and Rodrigue, G.H. (1989) Optimal Preconditioners of a Given Sparsity Pattern. BIT Numerical Mathematics, 29, 610-634. http://dx.doi.org/10.1007/BF01932737
- Axelsson, O. (1972) A Generalized SSOR Method. BIT Numerical Mathematics, 13, 443-467. http://dx.doi.org/10.1007/BF01932955
- Axelsson, O. (1994) Iterative Solution Methods. Cambridge University Press, Cambridge.
- Chan, T.F. and van der Vorst, H.A. (1997) Approximate and Incomplete Factorizations. In: Keyes, D.E., Sameh, A. and Venkatakrishnan, V., Eds., Parallel Numerical Algorithms, ICASE/LaRC Interdisciplinary Series in Science and Engineering, 4, Kluwer Academic, Dordrecht, 167-202.
- Hackbusch, W. (1994) Iterative Solution of Large Sparse Linear Systems of Equations. Applied Math. Sci. Series, No 95, Springer-Verlag, New York. http://dx.doi.org/10.1007/978-1-4612-4288-8
- Notay, Y. (1992) Upper Eigenvalue Bounds and Related Modified Incomplete Factorization Strategies. In: Beauwens, R. and de Groen, P., Eds., Iterative Methods in Linear Algebra, Amsterdam, North-Holland, 551-562.
- Bruaset, A.M. and Tveito, A. (1992) RILU Preconditioning; A Computational Study. Journal of Computational and Applied Mathematics, 39, 259-275. http://dx.doi.org/10.1016/0377-0427(92)90203-A
- Gustafsson, I. (1979) Stability and Rate of Convergence of Modified Incomplete Cholesky Factorization Method. Research Report 79.02R, Department of Computer Sciences, Chalmers University of Technology and University of Göteborg, Göteborg, Sweden.
- Beauwens, R. (1984) Upper Eigenvalue Bounds for Pencils of Matrices. Linear Algebra and Its Applications, 62, 87-104. http://dx.doi.org/10.1016/0024-3795(84)90088-0
- Beauwens, R. (1985) On Axelsson’s Perturbations. Linear Algebra and Its Applications, 68, 221-242. http://dx.doi.org/10.1016/0024-3795(85)90214-9