Technical Report 2015-005

Spectral Analysis and Spectral Symbol of Matrices in Isogeometric Galerkin Methods

Carlo Garoni, Carla Manni, Stefano Serra-Capizzano, Debora Sesana, and Hendrik Speleers

January 2015

A linear full elliptic second order Partial Differential Equation (PDE), defined on a d-dimensional domain Omega, is approximated by the isogeometric Galerkin method based on uniform tensor-product B-splines of degrees (p1,...,pd). The considered approximation process leads to a d-level stiffness matrix, banded in a multilevel sense. This matrix is close to a d-level Toeplitz structure when the PDE coefficients are constant and the physical domain Omega is just the hypercube (0,1)d without using any geometry map. In such a simplified case, a detailed spectral analysis of the stiffness matrices has been carried out in a previous work. In this paper, we complete the picture by considering non-constant PDE coefficients and an arbitrary domain Omega, parameterized with a non-trivial geometry map. We compute and study the spectral symbol of the related stiffness matrices. This symbol describes the asymptotic eigenvalue distribution when the fineness parameters tend to zero (so that the matrix-size tends to infinity). The mathematical technique used for computing the symbol is based on the theory of Generalized Locally Toeplitz (GLT) sequences.

Available as PDF (287 kB, no cover)

Download BibTeX entry.