Fast Fourier solvers for the tensor product high-order FEM for a Poisson type equation

DOI10.1134/S096554252002013XMaRDI QIDQ2206483zbMATH OpenOpenAlexFDO

Authors I. A. Zlotnik, Alexander Zlotnik

Publication date 22 October 2020

Published in Computational Mathematics and Mathematical Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1701.03967

zbMATH Keywords

fast direct algorithmFFThigh-order finite element methodPoisson equation

Mathematics Subject Classification ID

Boundary eigenvalue problems for ordinary differential equations (34B09) Eigenfunctions, eigenfunction expansions, completeness of eigenfunctions of ordinary differential operators (34L10) Boundary value problems for second-order elliptic equations (35J25) Numerical quadrature and cubature formulas (65D32) Numerical solutions to equations with nonlinear operators (65J15) Numerical methods for eigenvalue problems for boundary value problems involving PDEs (65N25) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Numerical methods for discrete and fast Fourier transforms (65T50)

Abstract: We present direct logarithmically optimal in theory and fast in practice algorithms to implement the tensor product high order finite element method on multi-dimensional rectangular parallelepipeds for solving PDEs of the Poisson kind. They are based on the well-known Fourier approaches. The key new points are the fast direct and inverse FFT-based algorithms for expansion in eigenvectors of the 1D eigenvalue problems for the high order FEM. The algorithms can further be used for numerous applications, in particular, to implement the tensor product high order finite element methods for various time-dependent PDEs. Results of numerical experiments in 2D and 3D cases are presented.

Recommendations

• A fast direct algorithm for implementing a high-order finite element method on rectangles as applied to boundary value problems for the Poisson equation
• Tensor product based fast spectral element solvers
• A fast algorithm for solving the tensor product collocation equations
• Fast tensor-product solvers: partially deformed three-dimensional domains
• scientific article; zbMATH DE number 749520

Cites work

• A fast direct algorithm for implementing a high-order finite element method on rectangles as applied to boundary value problems for the Poisson equation
• Accelerating MATLAB performance: 1001 tips to speed up MATLAB programs
• An efficient direct parallel spectral-element solver for separable elliptic problems
• Finite element method with discrete transparent boundary conditions for the time-dependent 1D Schrödinger equation
• scientific article; zbMATH DE number 193948 (Why is no real title available?)
• scientific article; zbMATH DE number 919117 (Why is no real title available?)
• Matrix decomposition algorithms for arbitrary order \(C^0\) tensor product finite element systems
• Matrix decomposition algorithms for elliptic boundary value problems: A survey
• Matrix decomposition algorithms for the \(C^{0}\)-quadratic finite element Galerkin method
• The Methods of Cyclic Reduction, Fourier Analysis and the FACR Algorithm for the Discrete Solution of Poisson’s Equation on a Rectangle
• The splitting in potential Crank–Nicolson scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip

Cited in

(5)

• A fast direct algorithm for implementing a high-order finite element method on rectangles as applied to boundary value problems for the Poisson equation
• On compact 4th order finite-difference schemes for the wave equation
• Tensor product based fast spectral element solvers
• A High-Order Spectral Element Fast Fourier Transform for the Poisson Equation
• Fast tensor-product solvers: partially deformed three-dimensional domains

Describes a project that uses

Uses Software

• FFTW
• LTFAT

This page was built for publication: Fast Fourier solvers for the tensor product high-order FEM for a Poisson type equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2206483)