Lusin theorem, GLT sequences and matrix computations: An application to the spectral analysis of PDE discretization matrices

详细信息    查看全文

• 作者：Carlo Garoni^a ; ^b ; ^{garoni@mat.uniroma2.it" class="auth_mail" title="E-mail the corresponding author} ; ^{carlo.garoni@uninsubria.it" class="auth_mail" title="E-mail the corresponding author} ; Carla Manni^a ; ^{manni@mat.uniroma2.it" class="auth_mail" title="E-mail the corresponding author} ; Stefano Serra-Capizzano^b ; ^c ; ^{stefano.serrac@uninsubria.it" class="auth_mail" title="E-mail the corresponding author} ; ^{stefano.serra@it.uu.se" class="auth_mail" title="E-mail the corresponding author} ; Debora Sesana^b ; ^{debora.sesana@uninsubria.it" class="auth_mail" title="E-mail the corresponding author} ; Hendrik Speleers^a ; ^{speleers@mat.uniroma2.it" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Lusin theorem ; Generalized locally Toeplitz sequences ; Spectral distribution and symbol ; PDE discretization matrices ; Isogeometric analysis
• 刊名：Journal of Mathematical Analysis and Applications
• 出版年：2017
• 出版时间：1 February 2017
• 年：2017
• 卷：446
• 期：1
• 页码：365-382
• 全文大小：479 K

文摘

In this paper we consider a general d  -dimensional second-order elliptic Partial Differential Equation (PDE) with variable coefficients, and we extend previous results on the spectral distribution of discretization matrices arising from B-spline Isogeometric Analysis (IgA). First, we provide the spectral symbol of the Galerkin B-spline IgA stiffness matrices, under the assumption that the PDE coefficients only belong to ser=111111111&_pii=S0022247X16304498&_rdoc=1&_issn=0022247X&md5=9c2aa00a8a38161bdb4a8d23d641dc2a" title="Click to view the MathML source">L^∞$L^{\infty }$. This symbol describes the asymptotic spectral distribution when the fineness parameters tend to zero (so that the matrix-size tends to infinity). Second, we prove the positive semi-definiteness of the ser=111111111&_pii=S0022247X16304498&_rdoc=1&_issn=0022247X&md5=258d2f8f7cd594ddb01aa3e26ed03f4b" title="Click to view the MathML source">d×d$d\times d$ symmetric matrix in the Fourier variables ser=111111111&_pii=S0022247X16304498&_rdoc=1&_issn=0022247X&md5=0ba3f4b0a1991bcc5488f1c1e14a965b" title="Click to view the MathML source">(θ₁,…,θ_d), which appears in the expression of the symbol. This matrix is related to the discretization of the (negative) Hessian operator, and its positive semi-definiteness implies the non-negativity of the symbol. The mathematical arguments used in our derivation are based on the Lusin theorem, on the theory of Generalized Locally Toeplitz (GLT) sequences, and on careful Linear Algebra manipulations of matrix determinants. These arguments are very general and can also be applied to other (local) PDE discretization methods, different from B-spline IgA.