Skip to Main content Skip to Navigation
Journal articles

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.
Document type :
Journal articles
Complete list of metadata

Cited literature [21 references]  Display  Hide  Download
Contributor : Sébastien Imperiale Connect in order to contact the contributor
Submitted on : Wednesday, September 11, 2019 - 9:24:42 AM
Last modification on : Friday, July 8, 2022 - 10:05:25 AM
Long-term archiving on: : Friday, February 7, 2020 - 11:34:39 PM


Files produced by the author(s)



Federica Caforio, Sébastien Imperiale. High-order discrete fourier transform for the solution of the Poisson equation. SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2019, 41 (5), pp.A2747-A2771. ⟨10.1137/18M1225410⟩. ⟨hal-01914257v2⟩



Record views


Files downloads