pith. sign in

arxiv: 2604.04613 · v1 · submitted 2026-04-06 · 🧮 math.NA · cs.NA

A Convergent Hybridizable Discontinuous Galerkin Method for Einstein--Scalar Equations

Mukul Dwivedi , Andreas Rupp This is my paper

Pith reviewed 2026-05-10 19:29 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords hybridizable discontinuous GalerkinEinstein-scalar equationsBondi gaugespherical symmetrynumerical relativityerror estimatesstabilitywell-posedness
0
0 comments X

The pith

A hybridized discontinuous Galerkin scheme for the radial Einstein-scalar system in Bondi gauge is locally well-posed, globally L2-stable, and converges optimally for the evolution variable when the polynomial degree is at least one.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a hybridizable discontinuous Galerkin method for the spherically symmetric Einstein-scalar equations written in Bondi gauge. The system is first recast as a local first-order PDE-ODE problem through scaled auxiliary variables. Element unknowns are then solved locally while interface traces on the radial mesh provide the coupling; these traces can be eliminated recursively, leaving only the main evolution variable to be advanced in time. The remaining metric quantities are recovered afterward from discrete versions of the constraint equations. Analysis establishes local semidiscrete well-posedness, a global L2 stability bound, and optimal-order L2 error estimates for the evolution variable together with reconstruction bounds for the metric components and mass functional.

Core claim

The central claim is that the proposed HDG discretization of the radial Einstein-scalar system permits recursive elimination of all interface traces, thereby reducing the time-stepping problem to a single unknown while still delivering local well-posedness, a global L2 stability estimate, optimal L2 convergence of order k+1 for the evolution variable when polynomials of degree k greater than or equal to 1 are used, and controlled reconstruction errors for the metric variables and the associated mass functional.

What carries the argument

Hybridizable discontinuous Galerkin discretization with recursive trace elimination on the one-dimensional radial mesh skeleton, which decouples local element solves from the global evolution and recovers metric variables from discrete constraints.

If this is right

  • Only the principal evolution variable requires time advancement; all metric quantities are reconstructed afterward from algebraic constraints.
  • Optimal L2 error bounds hold for the evolution variable and extend to the reconstructed metric components and mass functional.
  • The scheme remains stable for large-data collapse profiles as well as for smooth-pulse evolution.
  • Numerical experiments confirm the predicted spatial convergence rates on radial meshes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same recursive-elimination idea could be tested on other one-dimensional reductions of hyperbolic systems that arise in general relativity.
  • Guaranteed L2 stability may allow longer integration times for critical phenomena without added artificial viscosity.
  • If similar trace elimination proves feasible in angular directions, the method might extend beyond pure radial symmetry.

Load-bearing premise

The one-dimensional radial geometry permits complete recursive elimination of interface traces without loss of stability or accuracy in the hybridized scheme.

What would settle it

A sequence of refined radial meshes applied to a known smooth exact solution of the Einstein-scalar system in which the computed L2 error for the main evolution variable fails to decrease at the optimal rate once the polynomial degree is at least one.

Figures

Figures reproduced from arXiv: 2604.04613 by Andreas Rupp, Mukul Dwivedi.

Figure 1
Figure 1. Figure 1: Convergence history for Example 5.1. Panel (a) shows the error of uh with guide slopes k + 1, while panel (b) shows the error of gh with guide slopes k + 2 [PITH_FULL_IMAGE:figures/full_fig_p015_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Solutions of the Einstein–scalar system for Example 5.1. In each row, from left to right, the three panels show u˜h, g˜h, and gh. 5.2. Example 2. (Long-time evolution for a lower-energy large data set [8, Ex. 2]). In this example, we consider a large data set with a smaller effective concentration, so that the collapse develops more slowly and the intermediate evolution can be observed more clearly. We tak… view at source ↗
Figure 3
Figure 3. Figure 3: Solutions of the Einstein–scalar system for Example 5.2. In each row, from left to right, the three panels show u˜h, g˜h, and gh. rapidly. A natural reference point is the Bondi-coordinate study of Pürrer et al. [43, Sec. 4], who considered Gaussian-like initial data in their analysis of critical and subcritical Einstein–scalar evolutions. Motivated by this literature, we consider the localized pulse on th… view at source ↗
Figure 4
Figure 4. Figure 4: Solutions of the Einstein–scalar system for the Gaussian-pulse test of Example 5.3. In each row, from left to right, the three panels show u˜h, g˜h, and gh. [10] D. Christodoulou. The problem of a self-gravitating scalar field. Comm. Math. Phys., 105(3):337–361, 1986. [11] D. Christodoulou. A mathematical theory of gravitational collapse. Comm. Math. Phys., 109(4):613–647, 1987. [12] D. Christodoulou. The … view at source ↗
Figure 5
Figure 5. Figure 5: Time history of the numerical Bondi-mass proxy for Example 5.3. [19] B. Cockburn and K. Shi. Superconvergent HDG methods for linear elasticity with weakly symmetric stresses. IMA J. Numer. Anal., 33(3):747–770, 2013. [20] B. Cockburn and C.-W. Shu. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. II. General framework. Math. Comp., 52(186):411–435, 1989. … view at source ↗
read the original abstract

We propose and analyze a hybridized discontinuous Galerkin (HDG) method for the spherically symmetric Einstein--scalar system in Bondi gauge. After rewriting the model as a local first-order PDE--ODE system by introducing suitable scaled variables, we construct a semidiscrete scheme in which the element unknowns are computed locally and the coupling is carried by traces on the mesh skeleton. In the present radial setting, these traces can be eliminated recursively, so that only the main evolution variable is advanced in time, while the metric variables are recovered from discrete constraint relations. We prove local semidiscrete well-posedness, derive a global \(L^2\)--stability estimate, establish an optimal order \(L^2\) error bound for the main evolution variable for polynomial degree \(k\ge 1\), and obtain reconstruction error estimates for the metric variables and the associated mass functional. Numerical experiments verify the predicted spatial convergence rate and illustrate qualitative features of the Einstein--scalar dynamics, including large-data collapse profiles and smooth-pulse evolution.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 3 minor

Summary. The manuscript proposes a hybridized discontinuous Galerkin (HDG) method for the spherically symmetric Einstein-scalar system in Bondi gauge. The Einstein-scalar equations are first rewritten as a local first-order PDE-ODE system via scaled variables. A semidiscrete HDG scheme is constructed with element-local unknowns coupled through skeleton traces. In the radial (one-dimensional) setting, these traces are eliminated recursively, allowing the scheme to advance only the main evolution variable while recovering metric variables from discrete constraints. The paper proves local semidiscrete well-posedness, a global L²-stability estimate, an optimal-order L² error bound for the evolution variable when the polynomial degree k ≥ 1, and reconstruction error estimates for the metric variables and mass functional. Numerical experiments confirm the predicted convergence rates and illustrate qualitative features such as large-data collapse and smooth-pulse evolution.

Significance. If the central claims hold, the work supplies a convergent, stable, and computationally efficient discretization for a nonlinear hyperbolic system in general relativity. The recursive trace elimination in the radial setting is a notable technical device that reduces the global system to a single evolution variable while preserving the HDG structure. The combination of local well-posedness, global L² stability, optimal error estimates, and reconstruction bounds for the metric and mass functional would constitute a solid contribution to the numerical analysis of GR models, provided the elimination step does not compromise the underlying energy estimates or constraint consistency.

major comments (3)
  1. [§3.2] §3.2 (recursive elimination): The claim that recursive elimination of all interface traces preserves the HDG stability and accuracy properties is load-bearing for every subsequent result. The manuscript must supply an explicit lemma showing that the eliminated system retains the same energy identity used in the unreduced HDG formulation; without it, the global L²-stability estimate in §4 and the optimal error bound in §5 rest on an unverified reduction.
  2. [Theorem 4.1] Theorem 4.1 (local well-posedness): The local solvability argument appears to invoke the eliminated trace relations directly. It is unclear whether the discrete constraint relations used for metric recovery remain consistent with the original first-order system after elimination; a separate verification that the eliminated scheme satisfies the same algebraic constraints as the unreduced HDG scheme is required.
  3. [Theorem 5.3] Theorem 5.3 (optimal L² error for k ≥ 1): The error analysis relies on projection estimates that must account for both the PDE and ODE components after trace elimination. The manuscript should clarify whether the approximation properties at r = 0 (where spherical symmetry introduces singular coefficients) are handled by the same projection operators used for the interior elements.
minor comments (3)
  1. [Introduction] The abstract and introduction cite standard HDG references but omit recent works on HDG for hyperbolic systems with constraints; adding one or two such references would clarify the novelty of the recursive elimination technique.
  2. [Figure 1] Figure 1 (mesh and variable layout) would benefit from an explicit indication of which variables are eliminated versus retained after the recursive step.
  3. [§2.1] Notation for the scaled variables introduced in §2.1 is used inconsistently in the error analysis; a short table summarizing the variable definitions and their discrete counterparts would improve readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the detailed and insightful report. The comments highlight important points regarding the recursive elimination procedure and its implications for the stability and error analysis. We address each major comment below, providing clarifications and indicating the revisions we plan to make to strengthen the manuscript.

read point-by-point responses

  1. Referee: [§3.2] §3.2 (recursive elimination): The claim that recursive elimination of all interface traces preserves the HDG stability and accuracy properties is load-bearing for every subsequent result. The manuscript must supply an explicit lemma showing that the eliminated system retains the same energy identity used in the unreduced HDG formulation; without it, the global L²-stability estimate in §4 and the optimal error bound in §5 rest on an unverified reduction.

    Authors: We agree that an explicit statement is necessary to make the argument fully rigorous. In the revised version, we will insert a new lemma (Lemma 3.2) immediately after the description of the recursive elimination in §3.2. The lemma states that the energy identity for the reduced system is identical to that of the unreduced HDG scheme because the elimination is performed by solving the local algebraic problems for the traces and substituting back, which preserves the structure of the bilinear form and the resulting energy estimate. The proof is by induction over the mesh elements, starting from the outermost element and proceeding inward, confirming that no additional terms are introduced that would affect the L² stability. revision: yes

  2. Referee: [Theorem 4.1] Theorem 4.1 (local well-posedness): The local solvability argument appears to invoke the eliminated trace relations directly. It is unclear whether the discrete constraint relations used for metric recovery remain consistent with the original first-order system after elimination; a separate verification that the eliminated scheme satisfies the same algebraic constraints as the unreduced HDG scheme is required.

    Authors: The local well-posedness proof in Theorem 4.1 is carried out on the reduced system, but the consistency with the original constraints follows directly from the way the traces are eliminated using the exact HDG numerical flux definitions. To address the concern, we will add a corollary to Theorem 4.1 that explicitly verifies the satisfaction of the discrete constraints post-elimination. This verification shows that the metric variables recovered via the discrete constraints match those that would be obtained from the unreduced scheme, ensuring algebraic consistency with the first-order system. revision: yes

  3. Referee: [Theorem 5.3] Theorem 5.3 (optimal L² error for k ≥ 1): The error analysis relies on projection estimates that must account for both the PDE and ODE components after trace elimination. The manuscript should clarify whether the approximation properties at r = 0 (where spherical symmetry introduces singular coefficients) are handled by the same projection operators used for the interior elements.

    Authors: The projection operators used in the error analysis (defined in §5.1) are constructed to incorporate the radial weighting and the singular coefficients at r=0 through appropriate scaling in the definition of the scaled variables. These operators satisfy the same approximation properties uniformly across all elements, including the one adjacent to the origin, because the spherical symmetry is built into the formulation from the outset. We will add a clarifying paragraph in the proof of Theorem 5.3 referencing the weighted approximation results from Lemma 5.1, which explicitly handle the behavior at r=0. revision: yes

Circularity Check

0 steps flagged

No circularity: standard HDG construction and independent analysis in radial 1D setting

full rationale

The paper rewrites the Einstein-scalar system as a first-order PDE-ODE, builds a standard HDG semidiscrete scheme with local unknowns and skeleton traces, then exploits the radial (1D) structure to eliminate traces recursively. It subsequently proves local well-posedness, global L2 stability, optimal L2 error bounds for k≥1, and reconstruction estimates. These proofs are presented as direct consequences of the HDG formulation and the elimination step; no result is obtained by fitting parameters to data, renaming known patterns, or invoking self-citations whose validity depends on the current paper. The derivation chain is self-contained against external benchmarks of HDG theory.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard numerical-analysis assumptions (smoothness of solutions, well-posedness of the continuous problem) and the domain-specific choice of Bondi gauge and scaled variables; no free parameters or invented entities are introduced.

axioms (2)
  • domain assumption The spherically symmetric Einstein-scalar system in Bondi gauge admits a local first-order PDE-ODE reformulation via suitable scaled variables.
    Explicitly invoked in the abstract as the starting point for constructing the semidiscrete scheme.
  • domain assumption The radial mesh allows exact recursive elimination of all skeleton traces while preserving the local well-posedness and stability properties.
    Required for reducing the global system to a single evolution variable.

pith-pipeline@v0.9.0 · 5477 in / 1353 out tokens · 44621 ms · 2026-05-10T19:29:57.388341+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages

  1. [1]

    Arnowitt, S

    R. Arnowitt, S. Deser, and C. W. Misner. The dynamics of general relativity. InGravitation: An introduction to current research, pages 227–265. Wiley, New York-London, 1962

    work page 1962

  2. [2]

    J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and J. Van Meter. Gravitational-wave extraction from an inspiraling configuration of merging black holes.Phys. Rev. Lett., 96(11):111102, 2006

    work page 2006

  3. [3]

    R. Bartnik. Einstein equations in the null quasispherical gauge.Classical Quantum Gravity, 14(8):2185–2194, 1997

    work page 1997

  4. [4]

    Bondi, M

    H. Bondi, M. G. J. van der Burg, and A. W. K. Metzner. Gravitational waves in general relativity. VII. Waves from axi-symmetric isolated systems.Proc. Roy.Soc. London Ser. A, 269:21–52, 1962

    work page 1962

  5. [5]

    Campanelli, C

    M. Campanelli, C. O. Lousto, P. Marronetti, and Y. Zlochower. Accurate evolutions of orbiting black-hole binaries without excision.Phys. Rev. Lett., 96(11):111101, 2006

    work page 2006

  6. [6]

    H. Chen, P. Lu, and X. Xu. A hybridizable discontinuous Galerkin method for the Helmholtz equation with high wave number.SIAM J. Numer. Anal., 51(4):2166–2188, 2013

    work page 2013

  7. [7]

    Y. Chen, B. Dong, and J. Jiang. Optimally convergent hybridizable discontinuous Galerkin method for fifth- order Korteweg–de Vries type equations.ESAIM Math. Model. Numer. Anal., 52(6):2283–2306, 2018

    work page 2018

  8. [8]

    Chen, C.-W

    Y. Chen, C.-W. Shu, and S.-T. Yau. On the convergence of the discontinuous Galerkin scheme for Einstein- scalar equations.Math. Comp., 95(358):557–587, 2026

    work page 2026

  9. [9]

    Christodoulou

    D. Christodoulou. Global existence of generalized solutions of the spherically symmetric Einstein-scalar equa- tions in the large.Comm. Math. Phys., 106(4):587–621, 1986. HDG FOR EINSTEIN SCALAR EQUATION 19 r0 5 10 15 20 ~uh 0 0.2 0.4 0.6 ~uh r0 5 10 15 20 ~gh 0 0.2 0.4 0.6 ~gh r0 5 10 15 20 gh 0 0.2 0.4 0.6 0.8 1 gh (a)t=0;fromlefttoright: ~uh,~gh,gh. r0...

    work page 1986

  10. [10]

    Christodoulou

    D. Christodoulou. The problem of a self-gravitating scalar field.Comm. Math. Phys., 105(3):337–361, 1986

    work page 1986

  11. [11]

    Christodoulou

    D. Christodoulou. A mathematical theory of gravitational collapse.Comm. Math. Phys., 109(4):613–647, 1987

    work page 1987

  12. [12]

    Christodoulou

    D. Christodoulou. The structure and uniqueness of generalized solutions of the spherically symmetric Einstein- scalar equations.Comm. Math. Phys., 109(4):591–611, 1987

    work page 1987

  13. [13]

    Christodoulou

    D. Christodoulou. The formation of black holes and singularities in spherically symmetric gravitational collapse. Comm. Pure Appl. Math., 44(3):339–373, 1991

    work page 1991

  14. [14]

    P. G. Ciarlet.The finite element method for elliptic problems, volume Vol. 4 ofStudies in Mathematics and its Applications. North-Holland Publishing Co., Amsterdam-New York-Oxford, 1978

    work page 1978

  15. [15]

    Cockburn, B

    B. Cockburn, B. Dong, J. Guzmán, M. Restelli, and R. Sacco. A hybridizable discontinuous Galerkin method for steady-state convection-diffusion-reaction problems.SIAM J. Sci. Comput., 31(5):3827–3846, 2009

    work page 2009

  16. [16]

    Cockburn and J

    B. Cockburn and J. Gopalakrishnan. The derivation of hybridizable discontinuous Galerkin methods for Stokes flow. SIAM J. Numer. Anal., 47(2):1092–1125, 2009

    work page 2009

  17. [17]

    Cockburn, J

    B. Cockburn, J. Gopalakrishnan, and R. Lazarov. Unified hybridization of discontinuous Galerkin, mixed, and continuous Galerkin methods for second order elliptic problems.SIAM J. Numer. Anal., 47(2):1319–1365, 2009

    work page 2009

  18. [18]

    Cockburn, J

    B. Cockburn, J. Gopalakrishnan, N. C. Nguyen, J. Peraire, and F.-J. Sayas. Analysis of HDG methods for Stokes flow.Math. Comp., 80(274):723–760, 2011. 20 MUKUL DWIVEDI AND ANDREAS RUPP t0 10 20 30 40 50 60 Mh(t) 3.919 3.92 3.921 3.922 3.923 3.924 3.925 Mh(t) Figure 5.Time history of the numerical Bondi-mass proxy for Example 5.3

    work page 2011

  19. [19]

    Cockburn and K

    B. Cockburn and K. Shi. Superconvergent HDG methods for linear elasticity with weakly symmetric stresses. IMA J. Numer. Anal., 33(3):747–770, 2013

    work page 2013

  20. [20]

    Cockburn and C.-W

    B. Cockburn and C.-W. Shu. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. II. General framework.Math. Comp., 52(186):411–435, 1989

    work page 1989

  21. [21]

    Cockburn and C.-W

    B. Cockburn and C.-W. Shu. Runge-Kutta discontinuous Galerkin methods for convection-dominated problems. J. Sci. Comput., 16(3):173–261, 2001

    work page 2001

  22. [22]

    Cui and W

    J. Cui and W. Zhang. An analysis of HDG methods for the Helmholtz equation. IMA J. Numer. Anal., 34(1):279–295, 2014

    work page 2014

  23. [23]

    R. A. d’Inverno, M. R. Dubal, and E. A. Sarkies. Cauchy-characteristic matching for a family of cylindrical solutions possessing both gravitational degrees of freedom.Classical QuantumGravity, 17(16):3157–3170, 2000

    work page 2000

  24. [24]

    B. Dong. Optimally convergent HDG method for third-order Korteweg–de Vries type equations. J. Sci. Comput., 73(2-3):712–735, 2017

    work page 2017

  25. [25]

    Dwivedi, R

    M. Dwivedi, R. Gutendorf, and A. Rupp. A hybridizable discontinuous Galerkin method for the non–local Camassa–Holm–Kadomtsev–Petviashvili equation.arXiv preprint arXiv:2601.13800, 2026

    work page arXiv 2026

  26. [26]

    Dwivedi and A

    M. Dwivedi and A. Rupp. A hybridizable discontinuous Galerkin method for the Ostrovsky equation.arXiv preprint arXiv:2602.13786, 2026

    work page arXiv 2026

  27. [27]

    Dwivedi and T

    M. Dwivedi and T. Sarkar. Stability of fully discrete local discontinuous Galerkin method for the generalized Benjamin–Ono equation.arXiv preprint arXiv:2405.08360, 2024

    work page arXiv 2024

  28. [28]

    Gómez, R

    R. Gómez, R. A. Isaacson, J. S. Welling, and J. Winicour. Numerical relativity on null cones. InDynamical spacetimes and numerical relativity (Philadelphia, PA, 1985), pages 236–255. Cambridge Univ. Press, Cam- bridge, 1986

    work page 1985

  29. [29]

    Gómez, L

    R. Gómez, L. Lehner, R. L. Marsa, and J. Winicour. Moving black holes in 3D.Phys. Rev. D (3), 57(8):4778– 4788, 1998

    work page 1998

  30. [30]

    Griesmaier and P

    R. Griesmaier and P. Monk. Error analysis for a hybridizable discontinuous Galerkin method for the Helmholtz equation. J. Sci. Comput., 49(3):291–310, 2011

    work page 2011

  31. [31]

    S. G. Hahn and R. W. Lindquist. The two-body problem in geometrodynamics.Ann. Physics, 29:304–331, 1964

    work page 1964

  32. [32]

    Kabaria, A

    H. Kabaria, A. J. Lew, and B. Cockburn. A hybridizable discontinuous Galerkin formulation for non-linear elasticity.Comput. Methods Appl. Mech. Engrg., 283:303–329, 2015

    work page 2015

  33. [33]

    P. Lu, R. Maier, and A. Rupp. A localized orthogonal decomposition strategy for hybrid discontinuous Galerkin methods. ESAIM Math. Model. Numer. Anal., 59(2):1213–1237, 2025

    work page 2025

  34. [34]

    P. Lu, A. Rupp, and G. Kanschat. Analysis of injection operators in geometric multigrid solvers for HDG methods. SIAM J. Numer. Anal., 60(4):2293–2317, 2022

    work page 2022

  35. [35]

    P. Lu, A. Rupp, and G. Kanschat. Homogeneous multigrid for HDG.IMA J. Numer. Anal., 42(4):3135–3153, 2022

    work page 2022

  36. [36]

    P. Lu, W. Wang, G. Kanschat, and A. Rupp. Homogeneous multigrid for HDG applied to the Stokes equation. IMA J. Numer. Anal., 44(5):3124–3152, 2024

    work page 2024

  37. [37]

    N. C. Nguyen and J. Peraire. Hybridizable discontinuous Galerkin methods for the two-dimensional Monge- Ampère equation.J. Sci. Comput., 100(2):Paper No. 44, 32, 2024. HDG FOR EINSTEIN SCALAR EQUATION 21

    work page 2024

  38. [38]

    N. C. Nguyen, J. Peraire, and B. Cockburn. An implicit high-order hybridizable discontinuous Galerkin method for nonlinear convection-diffusion equations.J. Comput. Phys., 228(23):8841–8855, 2009

    work page 2009

  39. [39]

    N. C. Nguyen, J. Peraire, and B. Cockburn. A hybridizable discontinuous Galerkin method for Stokes flow. Comput. Methods Appl. Mech. Engrg., 199(9-12):582–597, 2010

    work page 2010

  40. [40]

    N. C. Nguyen, J. Peraire, and B. Cockburn. An implicit high-order hybridizable discontinuous Galerkin method for the incompressible Navier-Stokes equations.J. Comput. Phys., 230(4):1147–1170, 2011

    work page 2011

  41. [41]

    Pretorius

    F. Pretorius. Evolution of binary black-hole spacetimes.Phys. Rev. Lett., 95(12):121101, 4, 2005

    work page 2005

  42. [42]

    Pretorius

    F. Pretorius. Numerical relativity using a generalized harmonic decomposition.Classical Quantum Gravity, 22(2):425–451, 2005

    work page 2005

  43. [43]

    Pürrer, S

    M. Pürrer, S. Husa, and P. C. Aichelburg. News from critical collapse: Bondi mass, tails, and quasinormal modes. Phys. Rev. D, 71(10):104005, 2005

    work page 2005

  44. [44]

    W. Qiu, J. Shen, and K. Shi. An HDG method for linear elasticity with strong symmetric stresses.Math. Comp., 87(309):69–93, 2018

    work page 2018

  45. [45]

    Rhebergen and B

    S. Rhebergen and B. Cockburn. A space-time hybridizable discontinuous Galerkin method for incompressible flows on deforming domains.J. Comput. Phys., 231(11):4185–4204, 2012

    work page 2012

  46. [46]

    G. R. Richter. An optimal-order error estimate for the discontinuous Galerkin method. Math. Comp., 50(181):75–88, 1988

    work page 1988

  47. [47]

    R. K. Sachs. Gravitational waves in general relativity. VIII. Waves in asymptotically flat space-time.Proc. Roy.Soc. London Ser. A, 270:103–126, 1962

    work page 1962

  48. [48]

    A.Samii, N.Panda, C.Michoski, andC.Dawson.AhybridizeddiscontinuousGalerkinmethodforthenonlinear Korteweg–de Vries equation.J. Sci. Comput., 68(1):191–212, 2016

    work page 2016

  49. [49]

    Sun and C.-W

    Z. Sun and C.-W. Shu. Strong stability of explicit Runge-Kutta time discretizations.SIAM J. Numer. Anal., 57(3):1158–1182, 2019. Department of Mathematics, Saarland University, Saarbrücken, Germany Email address:{mukul.dwivedi;andreas.rupp}@uni-saarland.de

    work page 2019