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 metadatas

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, September 13, 2019 - 1:14:48 AM

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

12

Files downloads

133