AI summary · arXiv: 2408.10494v1 · math.NA

Tensor-Product Split-Simplex Summation-By-Parts Operators

Zelalem Arega Worku, Jason E. Hicken, David W. Zingg
Original paper Submitted to arXiv on 20 Aug. 2024 28 pages, 15 figures License: NONEXCLUSIVE-DISTRIB 1.0

How to read this summary? Illustrated guide

The essentials

  • Authors introduce innovative approach to constructing efficient sparse summation-by-parts (SBP) operators on triangles and tetrahedra with tensor-product structure
  • Method involves splitting simplices into quadrilateral or hexahedral subdomains, mapping tensor-product SBP operators onto subdomains, and assembling them using continuous-Galerkin-type procedure
  • Tensor-product split-simplex operators do not have repeated degrees of freedom at interior interfaces, ensuring SBP property by construction for stable discretizations
  • Accuracy and sparsity of operators enhance efficiency of SBP discretizations on simplicial meshes, crucial for entropy-stable discretizations based on two-point flux functions
  • Numerical experiments show proposed operators outperform existing dense multidimensional SBP operators by more than an order of magnitude in many cases in terms of accuracy per degree of freedom and computational time
  • Work represents significant advancement in numerical methods for solving partial differential equations on complex geometries by leveraging tensor-product structure and innovative splitting techniques

Key points written by AI: they may contain errors, check them against the paper.

AI-generated keywords: Tensor-Product Split-Simplex Summation-By-Parts Operators

Read further

  • Comprehensive Summary
  • Layman's Summary
  • Blog article
  • Original abstract
In their paper titled "Tensor-Product Split-Simplex Summation-By-Parts Operators," authors Zelalem Arega Worku, Jason E. Hicken, and David W. Zingg introduce an innovative approach to constructing efficient sparse summation-by-parts (SBP) operators on triangles and tetrahedra with a tensor-product structure. The key idea behind their method is to split the simplices into quadrilateral or hexahedral subdomains, map tensor-product SBP operators onto these subdomains, and then assemble them back together using a continuous-Galerkin-type procedure. One of the distinguishing features of the tensor-product split-simplex operators developed by the authors is that they do not have repeated degrees of freedom at the interior interfaces between the split subdomains. This property ensures that the operators satisfy the SBP property by construction, leading to stable discretizations. Moreover, the accuracy and sparsity of these operators significantly enhance the efficiency of SBP discretizations on simplicial meshes. The sparsity aspect of these operators is particularly crucial for entropy-stable discretizations based on two-point flux functions as it reduces the number of two-point flux computations required. Through a series of numerical experiments, the authors demonstrate that their proposed operators outperform existing dense multidimensional SBP operators by more than an order of magnitude in many cases. This superiority is evident in both accuracy per degree of freedom and computational time required to achieve a specified error threshold. Overall, this work represents a significant advancement in the field of numerical methods for solving partial differential equations on complex geometries. By leveraging a tensor-product structure and innovative splitting techniques, the authors have developed efficient SBP operators that offer improved performance and stability compared to traditional approaches.

Was this synthesis useful?

Score: 0

Why do we need votes?

Votes are used to determine whether we need to re-run our summarizing tools. If the count reaches -10, our tools can be restarted.

Figures and plots

Ask about a figure by number or use its question button to request visual analysis. Up to two figures and four image panels per question. Values read from plots are approximate; check the original data.

AI assistant instructions?

Tools for this paper

Explain it differently 3 free rewrites left today
Save to a collection Group papers in collections to compare them, chat with them and export them. Create a free account Log in

Advanced tools Premium

For researchers who follow a field closely: citation alerts, journal-club slides, version comparison and Zotero export. Discover Premium

Zotero
Citation alerts Citation alerts guide

Papers around this one

See them on a map

Ranked by how close their title and abstract are to this paper’s. Select some to read them together.

Closest papers

  1. Close
    Summation-by-parts operators for general function spaces: Flux differencing, entropy stability, and accuracy
    H. Haase, Andrew R. Winters, J. Glaubitz · 2026 · Not summarized yet
    Abstract

    Summation-by-parts (SBP) operators are a cornerstone of provably stable high-order methods for hyperbolic conservation laws. For nonlinear problems, the flux-differencing framework of Fisher and Carpenter combines such SBP operators with entropy-conservative two-point fluxes to obtain semi-discretizations that are conservative in the sense of Lax--Wendroff and discretely entropy-conservative or stable. So far, this framework has been developed exclusively for operators that are exact for polynomials. Function-space SBP (FSBP) operators lift this restriction by being exact on a general finite-dimensional function space that can be tailored to the problem at hand. In this work, we extend flux …

  2. Close
    Entropy Stable Element-Local Continuous-Discontinuous Residual Correction Summation-By-Parts Discretization
    Zelalem Arega Worku, D. C. Del Rey Fernández · 2026 · cited 1 time · Not summarized yet
    Abstract

    Residual correction provides an efficient alternative to two-point volume flux differencing for enforcing entropy stability in high-order discretizations of hyperbolic partial differential equations. We develop an element-local residual-correction summation-by-parts (SBP) method that combines continuous-SBP (C-SBP) and discontinuous-SBP (D-SBP) coupling on conforming meshes. The construction retains continuous assembly without simultaneous approximation terms (SATs) wherever suitable, restricts SATs to selected discontinuous facets, and permits the coupling to change in time as a foundation for adaptive discretizations. The scalar entropy defect and a graph-based distribution vector acting o…

  3. Close
    TensorSpectra: Bridging the Spectral and Tensor Worlds for High-Dimensional PDEs
    Bing-Cheng Hu, Tian-Chi Yu, Zhao-Xiang Li et al. · 2026 · Not summarized yet
    Abstract

    High-dimensional partial differential equations (PDEs) remain challenging for spectral methods due to the exponential growth of multidimensional coefficient spaces. We propose TensorSpectra, a low-rank spectral framework that directly parameterizes spectral coefficient tensors in canonical polyadic (CP) form, reducing representation complexity from $\mathcal{O}(N^d)$ to $\mathcal{O}(dNr)$. TensorSpectra combines structure-preserving spectral calculus, de-aliased nonlinear interactions, operator-aware residual weighting, and progressive low-rank corrections to enable efficient PDEs solving in compressed coefficient space. We further establish a theoretical framework connecting spectral trunca…

  4. Related
    Adaptively-refinable polar-spline discrete differential forms: hierarchical construction, exactness, and applications
    Diogo C. Cabanas, Deepesh Toshniwal, Rafael Vázquez · 2026 · Not summarized yet
    Abstract

    A common approach to representing disk-like and sphere-like geometries in isogeometric analysis is to use collapsed-edge singularities in the underlying parameterization. The resulting tensor-product B-spline spaces on such singular parameterizations lack the required smoothness to be used in isogeometric analysis. Polar splines (Toshniwal et al., CMAME 2017) rectify this issue by extracting a smoother subspace, and have been successfully used to discretize high-order partial differential equations, as well as to perform structure-preserving discretizations of the de Rham complex.The latter is particularly relevant for mixed formulations that appear in applications such as electromagnetism a…

  5. Related
    Simultaneous Matrix Multiplication by Tensor Contraction Aggregation
    Edgar Solomonik, Markus Bläser · 2026 · Not summarized yet
    Abstract

    We present an interpretation of the fast matrix multiplication algorithms of Coppersmith and Winograd as aggregations of tensor contractions. This point of view provides a direct connection to older constructions by Pan and to more recent work on bilinear algorithms for symmetric tensor contraction. Further, we provide a simple explicit approach for exact aggregation of tensor contractions, each one yielding an independent matrix-matrix product. These constructions and border rank adaptations thereof do not come close to improving the best known bounds on $\omega$ (the matrix multiplication exponent), but provide explicit moderate-size fast matrix multiplication algorithms that achieve a bet…

  6. Related
    Matrix Product Operators In The Age of Block Encoding
    Eugene Dumitrescu · 2026 · Quantum Physics · cited 2 times · AI summary available
    Abstract

    We develop a block-encoding compiler that speeds up linear combination of unitaries Hamiltonian simulation programs by treating matrix product operators as compressed, virtual-path LCU programs. In showing how these new conditional PREP and SELECT stages are compiled in terms of a parent matrix product operator, we go beyond typical operator splitting product formulas and illustrate how tensor networks are a natural language and valid intermediate representation for quantum circuits. Our results are numerically verified for two important cases, namely, Heisenberg and perturbed Heisenberg-adjacent chain real-time evolution, and highlight polynomial speedups. Specifically, we highlight a polyn…

  7. Related
    A Polynomial-Scaling PDE Solver with Entanglement-Basis Tensor Networks
    Abhijatmedhi Chotrattanapituk, Michael J. Landry, Chu-Liang Fu et al. · 2026 · Not summarized yet
    Abstract

    We develop a finite element method (FEM) for partial differential equation (PDE) solver based on the entanglement-basis representation introduced in our companion work. By lifting non-linear finite-element equations into an augmented coefficient space, the governing PDE together with boundary, initial, and inter-element constraints can be expressed through a unified quadratic residual minimization. Although this augmented space grows exponentially with the number of elements, its tensor-product structure allows it to be represented efficiently using tensor networks. Using the matrix product state (MPS) as a concrete example, we show that density matrix renormalization group (DMRG) sweeps ena…

  8. Related
    Efficient Construction of T-Mesh Bases via a Constraint Transformation Matrix
    C. Provatidis · 2026 · Axioms · Not summarized yet
    Abstract

    This paper addresses the treatment of large T-elements, based on Lagrange interpolation, through the application of constraints along the horizontal and vertical mesh lines of a quadrilateral patch. Instead of resorting to manual or symbolic factorization to determine explicit closed-form expressions of the shape functions, we propose constructing a transformation matrix and applying it directly to the background tensor-product grid, thereby establishing a linear relationship between the active (primary or master) and the auxiliary (secondary or slave) degrees of freedom. Software implementation is greatly simplified when all data (boundary conditions, slave-node definitions, etc.) refer dir…

Papers that cite it

  1. Related
    A collocation scheme that is equivalent to discontinuous Galerkin discretizations
    Jason E. Hicken · 2026 · arXiv.org · preprint, not peer reviewed · cited 1 time · Not summarized yet
    Abstract

    A spectral collocation operator with the summation-by-parts property was introduced by Chan to develop entropy-stable discontinuous Galerkin (DG) semi-discretizations (https://doi.org/10.1016/j.jcp.2018.02.033). The present work shows that semi-discretizations based on this collocation operator produce solutions that are equivalent to solutions of a DG semi-discretization using the same underlying quadrature. The equivalence holds regardless of the number of degrees of freedom in the collocation scheme and when the quadrature is not strictly positive. Extraneous degrees of freedom in the collocation scheme are associated with the nullspace of the operator and remain zero throughout an unstea…

  2. Related
    Finite-difference compatible entropy-conserving schemes for the compressible Euler equations
    C. De Michele, A. Edoh, G. Coppola · 2024 · Journal of Computational Physics · cited 10 times · Not summarized yet
    Abstract

    This paper introduces a family of entropy-conserving finite-difference discretizations for the compressible flow equations. In addition to conserving the primary quantities of mass, momentum, and total energy, the methods also preserve kinetic energy and pressure equilibrium. The schemes are based on finite-difference (FD) representations of the logarithmic mean, establishing and leveraging a broader link between linear and nonlinear two-point averages and FD forms. The schemes are locally conservative due to the summation-by-parts property and therefore admit a local flux form, making them applicable also in finite-volume and finite-element settings. The effectiveness of these schemes is va…

How are these papers found?

We compare the title and abstract of this paper with those of every paper summarized on this site, and with the papers Semantic Scholar recommends, those it cites and those citing it (for a much-cited paper, its most recent citations and, from OpenAlex, the most cited ones). Among the prior work and the follow-ups, the papers that use each other in depth, a method it uses and the much-cited follow-ups come a little higher. The levels (very close, close, related) come from the similarity of the texts, not from a judgment of quality; * marks a paper without an abstract, compared by its title. The sentence under a paper is written by an AI from the two abstracts: check it in the papers. A paper marked not relevant by three readers leaves the list.

Updated on 1 November 2026.

The previous summary was created more than a year ago and can be re-run (if necessary) by clicking on the Run button below.

Disclaimer: The AI-based summarization tool and virtual assistant provided on this website may not always provide accurate and complete summaries or responses. We encourage you to carefully review and evaluate the generated content to ensure its quality and relevance to your needs.