Fast Barycentric-Based Evaluation Over Spectral/hp Elements

by   Edward Laughton, et al.

As the use of spectral/hp element methods, and high-order finite element methods in general, continues to spread, community efforts to create efficient, optimized algorithms associated with fundamental high-order operations have grown. Core tasks such as solution expansion evaluation at quadrature points, stiffness and mass matrix generation, and matrix assembly have received tremendousattention. With the expansion of the types of problems to which high-order methods are applied, and correspondingly the growth in types of numerical tasks accomplished through high-order methods, the number and types of these core operations broaden. This work focuses on solution expansion evaluation at arbitrary points within an element. This operation is core to many postprocessing applications such as evaluation of streamlines and pathlines, as well as to field projection techniques such as mortaring. We expand barycentric interpolation techniques developed on an interval to 2D (triangles and quadrilaterals) and 3D (tetrahedra, prisms, pyramids, and hexahedra) spectral/hp element methods. We provide efficient algorithms for their implementations, and demonstrate their effectiveness using the spectral/hp element library Nektar++.



There are no comments yet.


page 25


Low-order preconditioning for the high-order finite element de Rham complex

In this paper we present a unified framework for constructing spectrally...

Computational performance of a parallelized high-order spectral and mortar element toolbox

In this paper, a comprehensive performance review of a MPI-based high-or...

Local Fourier Analysis of P-Multigrid for High-Order Finite Element Operators

Multigrid methods are popular for solving linear systems derived from di...

A High-Order Lower-Triangular Pseudo-Mass Matrix for Explicit Time Advancement of hp Triangular Finite Element Methods

Explicit time advancement for continuous finite elements requires the in...

Exposing and exploiting structure: optimal code generation for high-order finite element methods

Code generation based software platforms, such as Firedrake, have become...

Enhancing data locality of the conjugate gradient method for high-order matrix-free finite-element implementations

This work investigates a variant of the conjugate gradient (CG) method a...

Portable high-order finite element kernels I: Streaming Operations

This paper is devoted to the development of highly efficient kernels per...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.