<abstract xmlns="http://www.w3.org/1999/xhtml"><p>The problem of solving large-scale systems of linear algebraic equations arises in a wide range of applications. In many cases the preconditioned iterative method is a method of choice. This paper deals with the approximate inverse preconditioning AINV/SAINV based on the incomplete generalized Gram–Schmidt process. This type of the approximate inverse preconditioning has been repeatedly used for matrix diagonalization in computation of electronic structures but approximating inverses is of an interest in parallel computations in general. Our approach uses adaptive dropping of the matrix entries with the control based on the computed intermediate quantities. Strategy has been introduced as a way to solve di cult application problems and it is motivated by recent theoretical results on the loss of orthogonality in the generalized Gram– Schmidt process. Nevertheless, there are more aspects of the approach that need to be better understood. The diagonal pivoting based on a rough estimation of condition numbers of leading principal submatrices can sometimes provide inefficient preconditioners. This short study proposes another type of pivoting, namely the pivoting that exploits incremental condition estimation based on monitoring both direct and inverse factors of the approximate factorization. Such pivoting remains rather cheap and it can provide in many cases more reliable preconditioner. Numerical examples from real-world problems, small enough to enable a full analysis, are used to illustrate the potential gains of the new approach.</p></abstract>

The problem of solving large-scale systems of linear algebraic equations arises in a wide range of applications. In many cases the preconditioned iterative method is a method of choice. This paper deals with the approximate inverse preconditioning AINV/SAINV based on the incomplete generalized Gram–Schmidt process. This type of the approximate inverse preconditioning has been repeatedly used for matrix diagonalization in computation of electronic structures but approximating inverses is of an interest in parallel computations in general. Our approach uses adaptive dropping of the matrix entries with the control based on the computed intermediate quantities. Strategy has been introduced as a way to solve di cult application problems and it is motivated by recent theoretical results on the loss of orthogonality in the generalized Gram– Schmidt process. Nevertheless, there are more aspects of the approach that need to be better understood. The diagonal pivoting based on a rough estimation of condition numbers of leading principal submatrices can sometimes provide inefficient preconditioners. This short study proposes another type of pivoting, namely the pivoting that exploits incremental condition estimation based on monitoring both direct and inverse factors of the approximate factorization. Such pivoting remains rather cheap and it can provide in many cases more reliable preconditioner. Numerical examples from real-world problems, small enough to enable a full analysis, are used to illustrate the potential gains of the new approach.

A note on adaptivity in factorized approximate inverse preconditioning

Analele ştiinţifice ale Universităţii "Ovidius" Constanţa. Seria Matematică

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

{"article-title":"A note on adaptivity in factorized approximate inverse preconditioning"}

The problem of solving large-scale systems of linear algebraic equations arises in a wide range of applications. In many cases the preconditioned iterative...

A note on adaptivity in factorized approximate inverse preconditioning

Publicado en línea: 22 sept 2020

Páginas: 149 - 159

Recibido: 10 jul 2019

Aceptado: 16 dic 2019

DOI: https://doi.org/10.2478/auom-2020-0024

Palabras clave
sparse approximate inverse preconditioners, approximate factorization, generalized Gram-Schmidt process

© 2020 Jiří Kopal et al., published by Sciendo

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

A note on adaptivity in factorized approximate inverse preconditioning

Publicado en línea: 22 sept 2020

Páginas: 149 - 159

Recibido: 10 jul 2019

Aceptado: 16 dic 2019

DOI: https://doi.org/10.2478/auom-2020-0024

Palabras clavesparse approximate inverse preconditioners, approximate factorization, generalized Gram-Schmidt process

© 2020 Jiří Kopal et al., published by Sciendo

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

Palabras clave
sparse approximate inverse preconditioners, approximate factorization, generalized Gram-Schmidt process