Finite element block-factorized preconditioners

Robust preconditioners on block-triangular and block-factorized form for three types of linear systems of two-by-two block form are studied in this thesis. The first type of linear systems, which are dense, arise from a boundary element type of discretization of crack propagation problems. Numerical experiment show that simple algebraic preconditioning strategies results in iterative schemes that are highly competitive with a direct solution method. The second type of algebraic systems, which are sparse, indefinite and nonsymmetric, arise from a finite element (FE) discretization of the partial differential equations (PDE) that describe (visco)elastic glacial isostatic adjustment (GIA). The Schur complement approximation in the block preconditioners is constructed by assembly of local, exactly computed Schur matrices. The quality of the approximation is verified in numerical experiments. When the block preconditioners for the indefinite problem are combined with an inner iterative scheme preconditioned by a (nearly) optimal multilevel preconditioner, the resulting preconditioner is (nearly) optimal and robust with respect to problem size, material parameters, number of space dimensions, and coefficient jumps. Two approaches to mathematically formulate the PDEs for GIA are compared. In the first approach the equations are formulated in their full complexity, whereas in the second their formulation is confined to the features and restrictions of the employed FE package. Different solution methods for the algebraic problem are used in the two approaches. Analysis and numerical experiments reveal that the first strategy is more accurate and efficient than the latter. The block structure in the third type of algebraic systems is due to a fine-coarse splitting of the unknowns. The inverse of the pivot block is approximated by a sparse matrix which is assembled from local, exactly inverted matrices. Numerical experiments and analysis of the approximation show that it is robust with respect to problem size and coefficient jumps.

[1]  Yvan Notay Optimal V-cycle Algebraic Multilevel Preconditioning , 1998 .

[2]  Yvan Notay,et al.  Robust parameter‐free algebraic multilevel preconditioning , 2002, Numer. Linear Algebra Appl..

[3]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[4]  Owe Axelsson,et al.  Preconditioning methods for linear systems arising in constrained optimization problems , 2003, Numer. Linear Algebra Appl..

[5]  O. Axelsson,et al.  Algebraic multilevel preconditioning methods, II , 1990 .

[6]  Yousef Saad,et al.  Multilevel ILU With Reorderings for Diagonal Dominance , 2005, SIAM J. Sci. Comput..

[7]  Owe Axelsson Stabilization of algebraic multilevel iteration methods; additive methods , 2004, Numerical Algorithms.

[8]  P. Vassilevski,et al.  Algebraic multilevel preconditioning methods. I , 1989 .

[9]  Johannes K. Kraus,et al.  Algebraic multilevel preconditioning of finite element matrices using local Schur complements , 2006, Numer. Linear Algebra Appl..

[10]  Owe Axelsson,et al.  Preconditioning of Boundary Value Problems Using Elementwise Schur Complements , 2009, SIAM J. Matrix Anal. Appl..

[11]  Victor Eijkhout,et al.  The Nested Recursive Two-Level Factorization Method for Nine-Point Difference Matrices , 1991, SIAM J. Sci. Comput..

[12]  O. Axelsson,et al.  Finite element solution of boundary value problemes - theory and computation , 2001, Classics in applied mathematics.

[13]  I. Gustafsson,et al.  Preconditioning and two-level multigrid methods of arbitrary degree of approximation , 1983 .

[14]  曹志浩,et al.  ON ALGEBRAIC MULTILEVEL PRECONDITIONING METHODS , 1993 .

[15]  Edmond Chow,et al.  Approximate Inverse Techniques for Block-Partitioned Matrices , 1997, SIAM J. Sci. Comput..

[16]  Owe Axelsson,et al.  Variable-step multilevel preconditioning methods, I: Self-adjoint and positive definite elliptic problems , 1994, Numer. Linear Algebra Appl..

[17]  E. F. F. Botta,et al.  Matrix Renumbering ILU: An Effective Algebraic Multilevel ILU Preconditioner for Sparse Matrices , 1999, SIAM J. Matrix Anal. Appl..

[18]  Victor Eijkhout,et al.  The Role of the Strengthened Cauchy-Buniakowskii-Schwarz Inequality in Multilevel Methods , 1991, SIAM Rev..

[19]  Yvan Notay,et al.  Using approximate inverses in algebraic multilevel methods , 1998, Numerische Mathematik.

[20]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[21]  Yvan Notay,et al.  Optimal v-cycle algebraic multilevel preconditioning , 1998, Numer. Linear Algebra Appl..

[22]  Yousef Saad,et al.  ARMS: an algebraic recursive multilevel solver for general sparse linear systems , 2002, Numer. Linear Algebra Appl..

[23]  Panayot S. Vassilevski,et al.  On Two Ways of Stabilizing the Hierarchical Basis Multilevel Methods , 1997, SIAM Rev..

[24]  Owe Axelsson,et al.  Two Simple Derivations of Universal Bounds for the C.B.S. Inequality Constant , 2004 .