We consider fast solvers for large linear systems arising from the Galerkin approximation based on B-splines of classical d-dimensional elliptic problems, d≥. 1, in the context of isogeometric analysis. Our ultimate goal is to design iterative algorithms with the following two properties. First, their computational cost is optimal, that is linear with respect to the number of degrees of freedom, i.e.the resulting matrix size. Second, they are totally robust, i.e., their convergence speed is substantially independent of all the relevant parameters: in our case, these are the matrix size (related to the fineness parameter), the spline degree (associated to the approximation order), and the dimensionality d of the problem. We review several methods like PCG, multigrid, multi-iterative algorithms, and we show how their numerical behavior (in terms of convergence speed) can be understood through the notion of spectral distribution, i.e.through a compact symbol which describes the global eigenvalue behavior of the considered stiffness matrices. As a final step, we show how we can design an optimal and totally robust multi-iterative method, by taking into account the analytic features of the symbol. A wide variety of numerical experiments, few open problems and perspectives are presented and critically discussed.

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

DONATELLI, MARCO;SERRA CAPIZZANO, STEFANO;
2015

Abstract

We consider fast solvers for large linear systems arising from the Galerkin approximation based on B-splines of classical d-dimensional elliptic problems, d≥. 1, in the context of isogeometric analysis. Our ultimate goal is to design iterative algorithms with the following two properties. First, their computational cost is optimal, that is linear with respect to the number of degrees of freedom, i.e.the resulting matrix size. Second, they are totally robust, i.e., their convergence speed is substantially independent of all the relevant parameters: in our case, these are the matrix size (related to the fineness parameter), the spline degree (associated to the approximation order), and the dimensionality d of the problem. We review several methods like PCG, multigrid, multi-iterative algorithms, and we show how their numerical behavior (in terms of convergence speed) can be understood through the notion of spectral distribution, i.e.through a compact symbol which describes the global eigenvalue behavior of the considered stiffness matrices. As a final step, we show how we can design an optimal and totally robust multi-iterative method, by taking into account the analytic features of the symbol. A wide variety of numerical experiments, few open problems and perspectives are presented and critically discussed.
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING
http://www.sciencedirect.com/science/article/pii/S0045782514001844
B-splines; Galerkin method; Isogeometric analysis; PCG, multigrid, and multi-iterative methods; Symbol; Toeplitz matrices
File in questo prodotto:
Non ci sono file associati a questo prodotto.

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: http://hdl.handle.net/11383/1926921
 Attenzione

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

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 56
  • ???jsp.display-item.citation.isi??? 56
social impact