Robust and optimal multi-iterative techniques for IgA collocation linear systems

• DOI: 10.1016/j.cma.2014.11.036
• Publication Date: 2014-11-30T01:24:57Z
• Authors: Marco Donatelli; Carlo Garoni; Carla Manni; Stefano Serra-Cap...; Hendrik Speleers
• Keywords: Isogeometric analysis; B-splines; Collocation method; Toeplitz matrices; Symbol; P-GMRES, multigrid, and multi-iterative methods; 01 natural sciences; B-splines; Collocation method; Isogeometric analysis; Multi-iterative methods; Multigrid; P-GMRES; Symbol; Toeplitz matrices; Computer Science Applications1707 Computer Vision and Pattern Recognition; Computational Mechanics; Mechanics of Materials; Mechanical Engineering; Physics and Astronomy; P-GMRES; and multi-iterative methods; B-spline; Symbol; 0101 mathematics; Toeplitz matrice; Settore MAT/08 - ANALISI NUMERICA; multigrid; Collocation method; Isogeometric analysi

We consider fast solvers for the large linear systems coming from the Isogeometric Analysis (IgA) collocation approximation based on B-splines of full elliptic d-dimensional Partial Differential Equations (PDEs). We are interested in designing iterative algorithms which are optimal and robust. The former property implies that the computational cost is linear with respect to the number of degrees of freedom (i.e. the matrix size). The latter property means that the convergence rate is completely independent of (or only mildly dependent on) all the relevant parameters: in our setting, we can mention the coefficients of the PDE, the dimensionality d, the geometric map G describing the physical domain, the matrix size (related to the fineness parameters), and the spline degrees (associated with the IgA approximation order). Our approach is based on the spectral symbol, which describes the global eigenvalue behavior of the IgA collocation matrices. It is precisely the spectral information contained in the symbol that is exploited for the design of an optimal and robust multi-iterative solver of multigrid type. Several numerical experiments are presented and discussed in view of recent theoretical findings concerning the symbol and its properties.