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

Abstract : 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.
Type de document :
Pré-publication, Document de travail
Liste complète des métadonnées

Contributeur : Sébastien Imperiale <>
Soumis le : mardi 6 novembre 2018 - 17:44:17
Dernière modification le : mardi 13 novembre 2018 - 10:12:37


Fichiers produits par l'(les) auteur(s)


  • HAL Id : hal-01914257, version 1


Federica Caforio, Sébastien Imperiale. High-order discrete fourier transform for the solution of the Poisson equation. 2018. 〈hal-01914257〉