Premium
Incomplete block‐matrix factorization of M ‐matrices using two‐step iterative method for matrix inversion and preconditioning
Author(s) -
Buranay S.C.,
Iyikal O.C.
Publication year - 2021
Publication title -
mathematical methods in the applied sciences
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.719
H-Index - 65
eISSN - 1099-1476
pISSN - 0170-4214
DOI - 10.1002/mma.6502
Subject(s) - mathematics , block matrix , iterative method , incomplete cholesky factorization , incomplete lu factorization , matrix splitting , matrix (chemical analysis) , factorization , discretization , matrix decomposition , square matrix , mathematical analysis , algorithm , symmetric matrix , eigenvalues and eigenvectors , physics , materials science , quantum mechanics , composite material
Using the general method of Owe Axelsson given in 1986 for incomplete factorization of M ‐matrices in block‐matrix form, we give a recursive approach to construct incomplete block‐matrix factorization of M ‐matrices by proposing a two‐step iterative method for the approximation of the inverse of diagonal pivoting block matrices at each stage of the recursion. For various predescribed tolerances in the accuracy of the approximation of the inverses, the obtained incomplete block‐matrix factorizations are used to precondition the iterative methods as one‐step stationary iterative (OSSI) method and biconjugate gradient stabilized method (BI‐CGSTAB). Certain applications are conducted on M ‐matrices occurring from the discretization of two Dirichlet boundary value problems of Laplace's equation on a rectangle using finite difference method. Numerical results justify that the given incomplete block‐matrix factorization of M ‐matrices using the two‐step iterative method to approximate the inverse of diagonal pivoting block matrices at each stage of the recursion give robust preconditioners. The obtained results are presented through tables and figures.