Log In Sign Up

A Hermite-like basis for faster matrix-free evaluation of interior penalty discontinuous Galerkin operators

by   Martin Kronbichler, et al.

This work proposes a basis for improved throughput of matrix-free evaluation of discontinuous Galerkin symmetric interior penalty discretizations on hexahedral elements. The basis relies on ideas of Hermite polynomials. It is used in a fully discontinuous setting not for higher order continuity but to minimize the effective stencil width, namely to limit the neighbor access of an element to one data point for the function value and one for the derivative. The basis is extended to higher orders with nodal contributions derived from roots of Jacobi polynomials and extended to multiple dimensions with tensor products, which enable the use of sum factorization. The beneficial effect of the reduced data access on modern processors is shown. Furthermore, the viability of the basis in the context of multigrid solvers is analyzed. While a plain point-Jacobi approach is less efficient than with the best nodal polynomials, a basis change via sum-factorization techniques enables the combination of the fast matrix-vector products with effective multigrid constituents. The basis change is essentially for free on modern hardware because these computations can be hidden behind the cost of the data access.


page 1

page 2

page 3

page 4


Efficient low-order refined preconditioners for high-order matrix-free continuous and discontinuous Galerkin methods

In this paper, we design preconditioners for the matrix-free solution of...

Efficient Explicit Time Stepping of High Order Discontinuous Galerkin Schemes for Waves

This work presents algorithms for the efficient implementation of discon...

Linearizing the hybridizable discontinuous Galerkin method: A linearly scaling operator

This paper proposes a matrix-free residual evaluation technique for the ...

Fast matrix-free evaluation of discontinuous Galerkin finite element operators

We present an algorithmic framework for matrix-free evaluation of discon...

Constraint preserving discontinuous Galerkin method for ideal compressible MHD on 2-D Cartesian grids

We propose a constraint preserving discontinuous Galerkin method for ide...

Alias-free, matrix-free, and quadrature-free discontinuous Galerkin algorithms for (plasma) kinetic equations

Understanding fundamental kinetic processes is important for many proble...

Change of basis for m-primary ideals in one and two variables

Following recent work by van der Hoeven and Lecerf (ISSAC 2017), we disc...