Easy and Efficient preconditioning of the Isogeometric Mass Matrix

DOI10.1016/J.CAMWA.2020.12.009MaRDI QIDQ6342015arXivFDO

Authors Gabriele Loli, Giancarlo Sangalli, Mattia Tani

Publication date 3 June 2020

Numerical computation using splines (65D07) Preconditioners for iterative methods (65F08) Iterative numerical methods for linear systems (65F10) Numerical computation of matrix norms, conditioning, scaling (65F35) Spline approximation (41A15) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs (65M55) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Mesh generation, refinement, and adaptive methods for boundary value problems involving PDEs (65N50) Finite element methods applied to problems in solid mechanics (74S05) PDE constrained optimization (numerical aspects) (49M41)

Abstract: This paper deals with the fast solution of linear systems associated with the mass matrix, in the context of isogeometric analysis. We propose a preconditioner that is both efficient and easy to implement, based on a diagonal-scaled Kronecker product of univariate parametric mass matrices. Its application is faster than a matrix-vector product involving the mass matrix itself. We prove that the condition number of the preconditioned matrix converges to 1 as the mesh size is reduced, that is, the preconditioner is asymptotically equivalent to the exact inverse. Moreover, we give numerical evidence of its good behaviour with respect to the spline degree and the (possibly singular) geometry parametrization. We also extend the preconditioner to the multipatch case through an Additive Schwarz method.

This page was built for publication: Easy and Efficient preconditioning of the Isogeometric Mass Matrix

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6342015)