We derive a new formula for the asymptotic eigenvalue distribution of stiffness matrices obtained by applying Fi finite elements with standard mesh refinement to the semielliptic PDE of second order in divergence form -▽(Κ▽Tu) = f on Ω, u = g on ∂Ω. Here Ω ⊂ ℝ2, and K is supposed to be piecewise continuous and point wise symmetric semipositive definite. The symbol describing this asymptotic eigenvalue distribution depends on the PDE, but also both on the numerical scheme for approaching the underlying bilinear form and on the geometry of triangulation of the domain. Our work is motivated by recent results on the superlinear convergence behavior of the conjugate gradient method, which requires the knowledge of such asymptotic eigenvalue distributions for sequences of matrices depending on a discretization parameter h when h → 0. We compare our findings with similar results for the finite difference method which were published in recent years. In particular we observe that our sequence of stiffness matrices is part of the class of generalized locally Toeplitz sequences for which many theoretical tools are available. This enables us to derive some results on the conditioning and preconditioning of such stiffness matrices.

On the asymptotic spectrum of finite element matrix sequences

SERRA CAPIZZANO, STEFANO
2007-01-01

Abstract

We derive a new formula for the asymptotic eigenvalue distribution of stiffness matrices obtained by applying Fi finite elements with standard mesh refinement to the semielliptic PDE of second order in divergence form -▽(Κ▽Tu) = f on Ω, u = g on ∂Ω. Here Ω ⊂ ℝ2, and K is supposed to be piecewise continuous and point wise symmetric semipositive definite. The symbol describing this asymptotic eigenvalue distribution depends on the PDE, but also both on the numerical scheme for approaching the underlying bilinear form and on the geometry of triangulation of the domain. Our work is motivated by recent results on the superlinear convergence behavior of the conjugate gradient method, which requires the knowledge of such asymptotic eigenvalue distributions for sequences of matrices depending on a discretization parameter h when h → 0. We compare our findings with similar results for the finite difference method which were published in recent years. In particular we observe that our sequence of stiffness matrices is part of the class of generalized locally Toeplitz sequences for which many theoretical tools are available. This enables us to derive some results on the conditioning and preconditioning of such stiffness matrices.
2007
finite element methods, matrix sequence, asymptotic eigenvalue distribution
Beckermann, B.; SERRA CAPIZZANO, Stefano
File in questo prodotto:
File Dimensione Formato  
beckermann_finite_elem.pdf

accesso aperto

Descrizione: PDF editoriale
Tipologia: Altro materiale allegato
Licenza: Creative commons
Dimensione 271.05 kB
Formato Adobe PDF
271.05 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11383/1671187
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 40
  • ???jsp.display-item.citation.isi??? 37
social impact