We study the performance of adaptive Fourier-Galerkin methods in a periodic box in R^d with dimension d>1. These methods offer unlimited approximation power only restricted by solution and data regularity. They are of intrinsic interest but are also a first step towards understanding adaptivity for the hp-FEM. We examine two nonlinear approximation classes, one classical corresponding to algebraic decay of Fourier coefficients and another associated with exponential decay typical of spectral approximation. We investigate the natural sparsity class for the operator range and find that the exponential class is not preserved, thus in contrast with the algebraic class. This entails a striking different behavior of the feasible residuals that lead to practical algorithms, influencing the overall optimality. The sparsity degradation for the exponential class is partially compensated with coarsening. We present several feasible adaptive Fourier algorithms, prove their contraction properties, and examine the cardinality of the activated sets. The Galerkin approximations at the end of each iteration are quasi-optimal for both classes, but inner loops or intermediate approximations are sub-optimal for the exponential class.

Adaptive Fourier-Galerkin methods / Canuto, Claudio; Nochetto, R. H.; Verani, M.. - In: MATHEMATICS OF COMPUTATION. - ISSN 0025-5718. - STAMPA. - 83:(2014), pp. 1645-1687. [10.1090/S0025-5718-2013-02781-0]

Adaptive Fourier-Galerkin methods

CANUTO, CLAUDIO;
2014

Abstract

We study the performance of adaptive Fourier-Galerkin methods in a periodic box in R^d with dimension d>1. These methods offer unlimited approximation power only restricted by solution and data regularity. They are of intrinsic interest but are also a first step towards understanding adaptivity for the hp-FEM. We examine two nonlinear approximation classes, one classical corresponding to algebraic decay of Fourier coefficients and another associated with exponential decay typical of spectral approximation. We investigate the natural sparsity class for the operator range and find that the exponential class is not preserved, thus in contrast with the algebraic class. This entails a striking different behavior of the feasible residuals that lead to practical algorithms, influencing the overall optimality. The sparsity degradation for the exponential class is partially compensated with coarsening. We present several feasible adaptive Fourier algorithms, prove their contraction properties, and examine the cardinality of the activated sets. The Galerkin approximations at the end of each iteration are quasi-optimal for both classes, but inner loops or intermediate approximations are sub-optimal for the exponential class.
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/2506031
 Attenzione

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