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

https://hal.inria.fr/hal-01914257
Contributor : Sébastien Imperiale <>
Submitted on : Wednesday, September 11, 2019 - 9:24:42 AM
Last modification on : Friday, April 30, 2021 - 9:54:53 AM
Long-term archiving on: : Friday, February 7, 2020 - 11:34:39 PM

File

HOFFT_vHal.pdf
Files produced by the author(s)

Identifiers

Citation

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⟩

Share

Metrics

Record views

128

Files downloads

994