In this paper, we address the problem of preconditioning sequences of large sparse non-symmetric systems of linear equations and present two new strategies to construct approximate updates of factorized preconditioners. Both updates are based on the availability of an incomplete LU (ILU) factorization for one matrix of the sequence and differ in the approximation of the so-called ideal updates. The first strategy is an approximate diagonal update of the ILU factorization; the second strategy relies on banded approximations of the factors in the ideal update. The efficiency and reliability of the proposed preconditioners are shown in the solution of nonlinear systems of equations by preconditioned inexact Newton-Krylov methods. Matrix-free implementations of the updating strategy are provided and numerical experiments are carried out on application problems.
New updates of incomplete LU factorizations and applications to large nonlinear systems / S. Bellavia; B. Morini; M. Porcelli. - In: OPTIMIZATION METHODS & SOFTWARE. - ISSN 1029-4937. - ELETTRONICO. - 29:(2014), pp. 321-340. [10.1080/10556788.2012.762517]
New updates of incomplete LU factorizations and applications to large nonlinear systems
BELLAVIA, STEFANIA;MORINI, BENEDETTA;PORCELLI, MARGHERITA
2014
Abstract
In this paper, we address the problem of preconditioning sequences of large sparse non-symmetric systems of linear equations and present two new strategies to construct approximate updates of factorized preconditioners. Both updates are based on the availability of an incomplete LU (ILU) factorization for one matrix of the sequence and differ in the approximation of the so-called ideal updates. The first strategy is an approximate diagonal update of the ILU factorization; the second strategy relies on banded approximations of the factors in the ideal update. The efficiency and reliability of the proposed preconditioners are shown in the solution of nonlinear systems of equations by preconditioned inexact Newton-Krylov methods. Matrix-free implementations of the updating strategy are provided and numerical experiments are carried out on application problems.I documenti in FLORE sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.