Newton–Krylov methods, a combination of Newton-like methods and Krylov sub- space methods for solving the Newton equations, often need adequate preconditioning in order to be successful. Approximations of the Jacobian matrices are required to form preconditioners, and this step is very often the dominant cost of Newton–Krylov methods. Therefore, working with precon- ditioners may destroy the “Jacobian-free” (or matrix-free) setting where the single Jacobian-vector product can be provided without forming and storing the element of the true Jacobian. In this paper, we propose and analyze a preconditioning technique for sequences of nonsymmetric Jacobian matri- ces based on the update of an earlier preconditioner. The proposed strategy can be implemented in a matrix-free manner. Numerical experiments on popular test problems confirm the effectiveness of the approach in comparison with the standard ILU-preconditioned Newton–Krylov approaches.
Bellavia, S., Bertaccini, D., Morini, B. (2011). Nonsymmetric Preconditioner Updates in Newton–Krylov Methods for Nonlinear Systems. SIAM JOURNAL ON SCIENTIFIC COMPUTING, 33(5), 2595-2619 [10.1137/100789786].
Nonsymmetric Preconditioner Updates in Newton–Krylov Methods for Nonlinear Systems
BERTACCINI, DANIELE;
2011-01-01
Abstract
Newton–Krylov methods, a combination of Newton-like methods and Krylov sub- space methods for solving the Newton equations, often need adequate preconditioning in order to be successful. Approximations of the Jacobian matrices are required to form preconditioners, and this step is very often the dominant cost of Newton–Krylov methods. Therefore, working with precon- ditioners may destroy the “Jacobian-free” (or matrix-free) setting where the single Jacobian-vector product can be provided without forming and storing the element of the true Jacobian. In this paper, we propose and analyze a preconditioning technique for sequences of nonsymmetric Jacobian matri- ces based on the update of an earlier preconditioner. The proposed strategy can be implemented in a matrix-free manner. Numerical experiments on popular test problems confirm the effectiveness of the approach in comparison with the standard ILU-preconditioned Newton–Krylov approaches.File | Dimensione | Formato | |
---|---|---|---|
78978.pdf
solo utenti autorizzati
Dimensione
290.76 kB
Formato
Adobe PDF
|
290.76 kB | Adobe PDF | Visualizza/Apri Richiedi una copia |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.