Scalable Evaluation of Hadamard Products with Tensor Product Basis for Entropy-Stable High-Order Methods

arXiv2306.11665MaRDI QIDQ6440887FDOQ6440887

Authors: Alexander Cicchino, S. K. Nadarajah 

Publication date: 20 June 2023

Abstract: A sum-factorization form for the evaluation of Hadamard products with a tensor product basis is derived in this work. The proposed algorithm allows for Hadamard products to be computed in $mathcalOleft(n^{d+1}ight)$ flops rather than $mathcalOleft(n^{2d}ight)$, where $d$ is the dimension of the problem. With this improvement, entropy conserving and stable schemes, that require a dense Hadamard product in the general modal case, become computationally competitive with the modal discontinuous Galerkin (DG) scheme. We numerically demonstrate the application of the sum-factorized Hadamard product in our in-house partial differential equation solver PHiLiP based on the Nonlinearly Stable Flux Reconstruction scheme. We demonstrate that the entropy conserving flow solver scales at $mathcalOleft(n^{d+1}ight)$ for three-dimensional compressible flow in curvilinear coordinates, along with a computational cost comparison with the modal DG and over-integrated DG schemes.

This page was built for publication: Scalable Evaluation of Hadamard Products with Tensor Product Basis for Entropy-Stable High-Order Methods

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