Several old and new finite-element preconditioners for nodal-based spectral discretizations of -Delta u = f in the domain Omega = (-1, 1)(d) (d = 2 or 3), with Dirichlet or Neumann boundary conditions, are considered and compared in terms of both condition number and computational efficiency. The computational domain covers the case of classical single-domain spectral approximations (see [C. Canuto et al., Spectral Methods. Fundamentals in Single Domains, Springer, Heidelberg, 2006]), as well as that of more general spectral-element methods in which the preconditioners are expressed in terms of local (upon every element) algebraic solvers. The primal spectral approximation is based on the Galerkin approach with numerical integration (G-NI) at the Legendre-Gauss-Lobatto (LGL) nodes in the domain. The preconditioning matrices rely on either P-1, Q(1), or Q(1), (NI) (i.e., with numerical integration) finite elements on meshes whose vertices coincide with the LGL nodes used for the spectral approximation. The analysis highlights certain preconditioners, which yield the solution at an overall cost proportional to Nd+1, where N denotes the polynomial degree in each direction.

Finite-Element Preconditioning of G-NI Spectral Methods / Canuto, Claudio; P., Gervasio; A., Quarteroni. - In: SIAM JOURNAL ON SCIENTIFIC COMPUTING. - ISSN 1064-8275. - 31:(2010), pp. 4422-4451. [10.1137/090746367]

Finite-Element Preconditioning of G-NI Spectral Methods

CANUTO, CLAUDIO;
2010

Abstract

Several old and new finite-element preconditioners for nodal-based spectral discretizations of -Delta u = f in the domain Omega = (-1, 1)(d) (d = 2 or 3), with Dirichlet or Neumann boundary conditions, are considered and compared in terms of both condition number and computational efficiency. The computational domain covers the case of classical single-domain spectral approximations (see [C. Canuto et al., Spectral Methods. Fundamentals in Single Domains, Springer, Heidelberg, 2006]), as well as that of more general spectral-element methods in which the preconditioners are expressed in terms of local (upon every element) algebraic solvers. The primal spectral approximation is based on the Galerkin approach with numerical integration (G-NI) at the Legendre-Gauss-Lobatto (LGL) nodes in the domain. The preconditioning matrices rely on either P-1, Q(1), or Q(1), (NI) (i.e., with numerical integration) finite elements on meshes whose vertices coincide with the LGL nodes used for the spectral approximation. The analysis highlights certain preconditioners, which yield the solution at an overall cost proportional to Nd+1, where N denotes the polynomial degree in each direction.
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

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/11583/2302863
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo