High-order discrete fourier transform for the solution of the Poisson equation - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue SIAM Journal on Scientific Computing Année : 2019

High-order discrete fourier transform for the solution of the Poisson equation

Résumé

The aim of this work is to propose a novel, fast, matrix-free solver for the Poisson problem discretised with High-Order Spectral Element Methods (HO-SEM). This method is based on the use of the Discrete Fourier Transform to reduce the problem to the inversion of the symbol of the operator in frequency space. The solver proposed is endowed with several properties. First, it preserves the efficiency of standard FFT algorithm; then, the matrix storage is minimised; a pseudo-explicit Singular Value Decomposition (SVD) is used for the inversion of the symbols; finally, it can be easily extended to multiple dimensions and non-periodic boundary conditions. In particular, due to the underlying HO-SEM discretisation, the multi-dimensional symbol of the operator can be efficiently computed from the one-dimensional symbol by tensorisation.
Fichier principal
Vignette du fichier
HOFFT_vHal.pdf (700.54 Ko) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

hal-01914257 , version 1 (06-11-2018)
hal-01914257 , version 2 (11-09-2019)

Identifiants

Citer

Federica Caforio, Sébastien Imperiale. High-order discrete fourier transform for the solution of the Poisson equation. SIAM Journal on Scientific Computing, 2019, 41 (5), pp.A2747-A2771. ⟨10.1137/18M1225410⟩. ⟨hal-01914257v2⟩
189 Consultations
738 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More