Pith. sign in

REVIEW 4 major objections 5 minor 20 references

DiscoTEX 1.0: Discontinuous collocation and implicit-turned-explicit (IMTEX) integration symplectic, symmetric numerical algorithms with high order jumps for differential equations II: extension to higher-orders of numerical convergence

T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that DiscoTEX's jump-corrected Hermite time integration extends to order twelve, with each order converging at its nominal rate on distributionally sourced wave equations.

desk verdict Useful extension in intent, but the printed 12th-order formula is internally inconsistent and the PDE convergence evidence stops at 6th order. read the letter →

arxiv 2411.14399 v1 pith:A4F6NPGX submitted 2024-11-21 math.NA astro-ph.IMcs.NAgr-qcmath-phmath.MP

Lidia J. Gomes Da Silva This is my paper · ORCID
classification math.NAastro-ph.IMcs.NAgr-qcmath-phmath.MP MSC 65M7065L06
keywords DiscoTEXdistributionallysourcedwaveequationDiracdeltasourcediscontinuouscollocationHermitetimeintegrationhigher-orderconvergencejumpconditionshyperboloidalcoordinates
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper extends the DiscoTEX algorithm—a discontinuous collocation plus implicit-turned-explicit (IMTEX) Hermite time-integration method—to solve wave equations whose sources contain Dirac deltas and their derivatives, up to twelfth order in time. The author's aim is to supply the full higher-order jump conditions that correct the time integrator when a point-particle trajectory crosses a grid point, and to verify the scheme by computing numerical weak-form solutions to the distributionally sourced wave equation and comparing them with exact solutions. The reported convergence tests show that orders 2 through 12 achieve their nominal convergence rates, and that the 12th-order scheme reaches the same error accuracy as lower orders for the tested settings. The practical conclusion is that accuracy stagnates past sixth order while computational cost rises steeply, so the fourth-order variant is the best cost-accuracy compromise under the current control factors. If correct, this gives a concrete recipe for long evolutions, such as extreme-mass-ratio inspirals, that do not require a second interpolation step.

What carries the argument

The mechanism is the discontinuous Hermite time integrator: a smooth Hermite rule whose kernel is replaced by two polynomials $p_-$ (before the discontinuity) and $p_+$ (after), meeting at $t_\times$ with prescribed jumps in value and derivatives. The jump-corrected rule at order 12 is written as the smooth trapezoidal/Hermite combination plus the correction $J_{H12}(\Delta t_\times, \Delta t)$, which packages all $J_0,\dots,J_{11}$ terms produced by the 24-condition linear solve. In the full DiscoTEX scheme the same philosophy is applied in space via a discontinuous Lagrange basis, and the time jumps are combined with differential-operator corrections through the $g$-vectors and the $\Upsilon$ source term; the paper explicitly gives the $g^{(3)},g^{(4)},g^{(5)}$ time derivatives needed by the H8, H10, and H12 variants.

What would settle it

Recompute the 12th-order jump correction $J_{H12}$ symbolically from the 24 collocation-and-jump conditions stated in Section 2.3 and verify that Eq. (47) is reproduced exactly, then run the DH12 scheme on the Legendre test integral (48): if the error does not decay at the 12th-order rate, the central claim is falsified.

Watch

Extended reading notes

Core claim

The central claim is that the DiscoTEX construction generalizes from the second- and fourth-order schemes of the companion paper to a twelfth-order scheme. On each time step the integrand is represented as a pair of degree-11 polynomials matched at the crossing time $t_\times$; the 24 unknowns are fixed by 12 endpoint collocation conditions and 12 jump conditions $J_m = p_+^{(m)}(t_\times)-p_-^{(m)}(t_\times)$, $m=0,\dots,11$. Integrating the piecewise polynomial yields the corrected rule with jump term $J_{H12}(\Delta t_\times, \Delta t)$; the paper also lists the full set of time jumps for the wave-equation application in a hyperboloidal chart. Applied to $\Box\Psi = F(\tau)\delta'(\sigma-\xi_p)+G(\tau)\delta(\sigma-\xi_p)$, the scheme reproduces the exact weak-form solution and the reported convergence rates match the nominal order of each integrator. The paper's own numerical summary is that, for equal control factors, higher-order integrators past sixth order do not improve accuracy but increase wall-clock time substantially.

Load-bearing premise

The extension stands on the correctness of the printed high-order jump formulas, especially the 12th-order correction term and the time jumps in the appendix, which the paper does not fully derive and one of which appears to contain typographical inconsistencies; if those formulas are wrong, the claimed orders of convergence do not follow.

Editorial extensions

If this is right

  • Orders 2, 4, 6, 8, 10, and 12 of the jump-corrected Hermite integrator converge at their nominal rates on the distributionally sourced wave equation, verified against exact solutions.
  • The explicit jump formulas in Appendix A.4 give a template for applying DiscoTEX to other time-dependent point-particle trajectories, including eccentric or generic orbits.
  • For the tested setup, running past sixth order does not lower the error but multiplies the wall-clock time, so fourth order is the recommended setting among the available variants.
  • Because the higher-order schemes reach a given accuracy with larger timesteps, they are candidates for long evolutions where a second interpolation step would be undesirable.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, the same 24-condition construction could be checked automatically by symbolic computation, which would settle whether the apparent typographical inconsistencies in Eq. (47) are typographical or substantive.
  • If the jump machinery is correct, the procedure should extend to sources with higher-order delta derivatives (for example, $\delta''$), where exact solutions are available and the 12th-order convergence could be tested directly.
  • The accuracy stagnation observed past sixth order may be an artifact of the time-dependent trajectory jumps dominating the error; for fixed-position sources, where only the time-stepper order matters, higher-order variants might show a clearer benefit.
  • A natural testable extension is to apply the H12 scheme to the same wave equation with a self-consistently coupled trajectory, where no exact solution exists, and compare against lower-order variants at fixed accuracy.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript extends the DiscoTEX algorithm to time-integration orders 8, 10, and 12. The method uses piecewise Hermite interpolants in time, corrected by jump conditions, and combines them with the discontinuous collocation spatial discretization introduced in the companion paper [1]. The numerical evidence consists of a discontinuous Legendre integral test (Figure 1) and weak-form solutions of a distributionally sourced wave equation compared against exact solutions (Figures 2-5, Table 2). The paper also reports wall-clock times and concludes that orders beyond six are not cost-effective for the tested configuration.

Significance. If the printed formulas are correct, the manuscript would provide a systematic route to high-order weak-form solvers for distributionally sourced PDEs, with potential applications in extreme-mass-ratio inspiral modelling and black-hole spectroscopy. The numerical comparison against independent exact solutions, with no fitted constants, is a genuine strength, as is the explicit reporting of computational cost. However, the central new objects -- the order-8, order-10, and order-12 discontinuous Hermite rules -- are not verified at the PDE level, and the printed 12th-order jump correction contains apparent transcription errors. The significance of the claimed extension therefore cannot be assessed from the manuscript in its present form.

major comments (4)
  1. [Sec. 2.3, Eq. (47)] The printed JH12 expression is internally inconsistent. The second term, (1/44)(5 Delta t^2 - 22 Delta t Delta t_x + 22 Delta t_x^2), is not multiplied by J1, and no other term in JH12 carries J1; since J1 = 15/8 in the paper's own Legendre test (Eq. (49)), this either drops a nonzero jump contribution or introduces a dimensionally inconsistent term. In addition, two different polynomials are labelled J5, with denominators 665280 and 30240, although each jump condition (34)-(45) should contribute exactly one term. As printed, Eq. (47) cannot be the correct 12th-order jump correction, so the claimed 12th-order convergence in Figure 1 and Table 2 is not supported by the manuscript.
  2. [Sec. A.1, Eq. (A.10)] The JH10 expression has the same missing-factor pattern: the term (1/18)(Delta t - 3 Delta t_x)(2 Delta t - 3 Delta t_x) appears without a J1 factor, while all neighbouring terms carry the expected jumps. This affects the printed 10th-order recipe and, together with the JH12 issue, indicates that the higher-order jump formulas need a systematic check rather than isolated typographical fixes.
  3. [Sec. 3.2, Fig. 5 and Table 2] PDE-level convergence is demonstrated only for orders 2, 4, and 6. Figure 5 contains DH2, DH4, and DH6 only, and Table 2 shows stagnating errors (7.7e-11 for H4, 7.2e-11 for H6, 7.2e-9 for H8, 7.2e-11 for H10 and H12) with no convergence-order study. The abstract's claim of an extension to twelve orders therefore rests on the Legendre integral test alone, which uses the very formulas questioned above. A revision must either add H8-H12 convergence studies for the wave equation or restrict the claims.
  4. [Appendix A.4, after Eq. (A.55)] The higher-order time-jump formulas are asserted without derivation. The text states that the remaining terms are determined by 'an analogous trivial process which we omit for brevity,' but these jumps are the central new object of the paper. Without a derivation, a symbolic notebook, or a fully specified recurrence, the reader cannot verify the printed formulas; this gap is especially serious in view of the apparent errors in Eqs. (47) and (A.10).
minor comments (5)
  1. [Eq. (A.8)] The JH4 term printed as 'Delta t_x^2/12 (Delta t^2 - 6 Delta t Delta t_x + 6 Delta t_x^2) J1' appears to have the wrong prefactor; dimensionally it should be Delta t^2/12 (Delta t^2 - 6 Delta t Delta t_x + 6 Delta t_x^2) J1, so please check.
  2. [Eq. (73)] The second component of the second time-derivative source term is printed as ~s^{(2)}_{\Psi,(2)} twice; one occurrence should presumably be ~s^{(2)}_{\Psi,(1)}.
  3. [Table 1 caption] The caption says 'all the 12th order discontinuous Hermite interpolation schemes' but the table lists orders 2 through 12; the wording should be corrected.
  4. [Reference [5]] The exact-solution formulas are cited to an unpublished work; for reproducibility they should either be included in an appendix or linked to a public preprint.
  5. [Eq. (46)] The notation f(t)_{DH12} for the quadrature result is not defined and could be confused with the integrand f(t); a different symbol for the integrated quantity would improve clarity.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the 12th-order extension is validated against independent exact solutions and no fitted parameter is renamed as a prediction; the main risk is a probable typo in Eq. (47), which is a correctness issue, not a circular one.

full rationale

The central claim is that the DiscoTEX time integrator can be extended to 12th order via the discontinuous Hermite rules. The derivation of the new rules is presented in this paper (Eqs. (46)-(47)), and the convergence test in Fig. 1 compares the quadrature of a Legendre function with analytically known jumps (Eqs. (48)-(51)) to its exact integral. The jumps and endpoint derivatives are computed from the test function itself; they are not fitted to force the integral value, so the resulting convergence rates are an independent check of the construction. For the PDE benchmark, the DiscoTEX solutions are compared to exact solutions (Eq. (3)) taken from external works [2-5]. The control factors N, J, and Delta_tau are standard discretization parameters selected through convergence studies; no empirical constant is fitted to the exact solution and then renamed as a prediction. The heavy use of ref. [1] supplies the previously established jump recurrence (Eq. (15)) and spatial discontinuous-collocation machinery; these are methodological antecedents, not premises that assume the order-12 result. There is no imported uniqueness theorem and no ansatz introduced solely by citation. Two qualifications belong in the verdict but are not circularity: Eq. (47) appears to contain typographical inconsistencies (the second J1 term is missing its J1 factor and two different J5 polynomials appear), and Appendix A.4 omits the derivation of the higher-order time jumps with a pointer to [1]. Those omissions and the absence of PDE convergence tests above H6 are correctness/verifiability risks for the printed formulas; they do not make the derivation equivalent to its inputs. The score of 2 reflects only a minor self-citation burden, not any circular reduction.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new physical entities. The free parameters are standard numerical discretization controls chosen by convergence tests. The central method relies on the author's prior work [1] for the jump machinery and on standard hyperboloidal slicing techniques, with no independent verification or machine-checked proofs provided.

free parameters (3)
  • Number of collocation nodes N = 45
    User-specified control factor [CTRL F1] chosen via convergence tests in Figure 4; not derived from theory.
  • Number of jumps J = 19
    User-specified control factor [CTRL F2] chosen via convergence tests in Figure 4; not derived from theory.
  • Timestep size Δτ = 0.00666667
    User-specified control factor [CTRL F3]; set to the same value used in [1] for fair comparison, not derived from theory.
assumptions (4)
  • domain assumption The exact solution for the distributionally-sourced wave equation (Eq. 3) is correct.
    Taken from Field et al. [3] and the unpublished work [5]; used as the reference solution for error computation.
  • domain assumption The hyperboloidal coordinate map (Eqs. 62-63) provides a well-posed outflow boundary setting for the wave equation.
    Standard 'scri-fixing' technique from [9-14]; assumed without analysis in this paper.
  • domain assumption The jump recurrence relation (Eq. 15) correctly determines all higher-order jumps J_m for a moving particle.
    Derived in the author's previous paper [1]; not re-derived or independently verified here.
  • domain assumption The IMTEX replacement, expressing U^(k) in terms of powers of L, is algebraically valid and stable for these orders.
    Used implicitly in Eqs. (70), (A.16), (A.20), and (A.24); no stability analysis is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of DiscoTEX 1.0: Discontinuous collocation and implicit-turned-explicit (IMTEX) integration symplectic, symmetric numerical algorithms with high order jumps for differential equations II: extension to higher-orders of numerical convergence." pith.science (2026). https://pith.science/paper/A4F6NPGX

@misc{pith2026241114399,
  author       = {Pith},
  title        = {Pith review of: DiscoTEX 1.0: Discontinuous collocation and implicit-turned-explicit (IMTEX) integration symplectic, symmetric numerical algorithms with high order jumps for differential equations II: extension to higher-orders of numerical convergence},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/A4F6NPGX}},
  note         = {Machine review of arXiv:2411.14399}
}
read the original abstract

\texttt{DiscoTEX} is a highly accurate numerical algorithm for computing numerical weak-form solutions to distributionally sourced partial differential equations (PDE)s. The aim of this second paper, succeeding \cite{da2024discotex}, is to present its extension up to twelve orders. This will be demonstrated by computing numerical weak-form solutions to the distributionally sourced wave equation and comparing it to its exact solutions. The full details of the numerical scheme at higher orders will be presented.

Figures

Figures reproduced from arXiv: 2411.14399 by the authors.

Figure 1
Figure 1. Numerical error associated with the numerical evaluation of the integral in equation (48) with the numerical scheme of equations ((A.1) - (A.10)) corrected [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Numerical weak-form solution to Ψ(τ, σ) obtained via the DiscoTEX H6 6th- order algorithm. Left: Numerical field Ψ(τ, σ) for a point-particle in time-dependent linear motion ξp(τc), where v is the particle’s velocity. Specifically, here ξp ≈ 0.608 at the coordinate-time τc ≈ 0.876 and v = 1/4. Right: Waveform for the point-particle computed on the numerical domain σ ∈ [0, 1] and τ ∈ [−1.52, 4.50]. As expected the re… view at source ↗
Figure 3
Figure 3. Phase portrait for the numerical weak-form solution obtained via [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Numerical convergence studies assessing the optimal user-specifiable control factors: number of [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: Numerical error associated with computation of the numerical weak-form solution [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 16 canonical work pages

  1. [1]

    L. J. G. Da Silva, Discotex: Discontinuous collocation and implicit- turned-explicit (imtex) integration symplectic, symmetric numerical al- gorithms with higher order jumps for di fferential equations with nu- merical black hole perturbation theory applications, arXiv preprint arXiv:2401.08758 (2024)

    arXiv 2024

  2. [5]

    S. E. Field, Exact solutions to distributionally sourced wave equations ex- tending [2] Appendix B to higher order deltas, unpublished work (2023)

    work page 2023

  3. [2]

    S. E. Field, J. S. Hesthaven, S. R. Lau, Discontinuous galerkin method for computing gravitational waveforms from extreme mass ratio binaries, Classical and Quantum Gravity 26 (16) (2009) 165010

    work page 2009

  4. [3]

    S. E. Field, J. S. Hesthaven, S. R. Lau, Persistent junk solutions in time- domain modeling of extreme mass ratio binaries, Physical Review D 81 (12) (2010) 124030

    work page 2010

  5. [4]

    S. E. Field, S. Gottlieb, G. Khanna, E. McClain, Discontinuous galerkin method for linear wave equations involving derivatives of the dirac delta distribution, in: Spectral and High Order Methods for Partial Differential Equations ICOSAHOM 2020+ 1: Selected Papers from the ICOSAHOM Conference, Vienna, Austria, July 12-16, 2021, Springer, 2022, pp. 307– 321

    work page 2020

  6. [6]

    Afshordi, S

    N. Afshordi, S. Akçay, P. A. Seoane, A. Antonelli, J. C. Aurrekoetxea, L. Barack, E. Barausse, R. Benkel, L. Bernard, S. Bernuzzi, et al., Wave- form modelling for the laser interferometer space antenna, arXiv preprint arXiv:2311.01300 (2023)

    arXiv 2023

  7. [7]

    L. J. Gomes Da Silva, R. P. Macedo, J. E. Thompson, J. A. V . Kroon, L. Durkan, O. Long, Hyperboloidal discontinuous time-symmetric nu- merical algorithm with higher order jumps for gravitational self-force computations in the time domain, arXiv preprint arXiv:2306.13153 (2023)

    arXiv 2023

  8. [8]

    L. J. Gomes Da Silva, Advanced numerical algorithms for EMRIs mod- elling in the time domain via DiscoTEX and relatives - ground-zero of the workpackage [E/(X)MRIs@DiscoTEX], PhD thesis, QMUL Univer- sity, School of Mathematical Sciences, https://qmro.qmul.ac.uk/ xmlui/handle/123456789/99539 (2024)

    arXiv 2024

Show all 20 references
  1. [9]

    Zengino ˘glu, Hyperboloidal foliations and scri-fixing, Classical and Quantum Gravity 25 (14) (2008) 145002

    A. Zengino ˘glu, Hyperboloidal foliations and scri-fixing, Classical and Quantum Gravity 25 (14) (2008) 145002

    2008

  2. [10]

    Zengino ˘glu, D

    A. Zengino ˘glu, D. Nunez, S. Husa, Gravitational perturbations of schwarzschild spacetime at null infinity and the hyperboloidal initial value problem, Classical and Quantum Gravity 26 (3) (2009) 035009

    2009

  3. [11]

    Ansorg, R

    M. Ansorg, R. P. Macedo, Spectral decomposition of black-hole perturba- tions on hyperboloidal slices, Physical Review D 93 (12) (2016) 124016

    2016

  4. [12]

    J. L. Jaramillo, R. P. Macedo, L. Al Sheikh, Pseudospectrum and black hole quasinormal mode instability, Physical Review X 11 (3) (2021) 031003

    2021

  5. [13]

    R. P. Macedo, J. L. Jaramillo, M. Ansorg, Hyperboloidal slicing approach to quasinormal mode expansions: The reissner-nordström case, Physical Review D 98 (12) (2018) 124005

    2018

  6. [14]

    R. P. Macedo, Hyperboloidal framework for the kerr spacetime, Classical and Quantum Gravity 37 (6) (2020) 065019

    2020

  7. [15]

    Amaro-Seoane, Extremely large mass-ratio inspirals, Physical Review D 99 (12) (2019) 123025

    P. Amaro-Seoane, Extremely large mass-ratio inspirals, Physical Review D 99 (12) (2019) 123025

    2019

  8. [16]

    Gourgoulhon, A

    E. Gourgoulhon, A. Le Tiec, F. H. Vincent, N. Warburton, Gravitational waves from bodies orbiting the galactic center black hole and their de- tectability by lisa, Astronomy & Astrophysics 627 (2019) A92

    2019

  9. [17]

    Barsanti, V

    S. Barsanti, V . De Luca, A. Maselli, P. Pani, Detecting subsolar-mass pri- mordial black holes in extreme mass-ratio inspirals with lisa and einstein telescope, Physical Review Letters 128 (11) (2022) 111104

    2022

  10. [18]

    Vázquez-Aceves, L

    V . Vázquez-Aceves, L. Zwick, E. Bortolas, P. R. Capelo, P. Amaro Seoane, L. Mayer, X. Chen, Revised event rates for ex- treme and extremely large mass-ratio inspirals, Monthly Notices of the Royal Astronomical Society 510 (2) (2022) 2379–2390

    2022

  11. [19]

    Vázquez-Aceves, Y

    V . Vázquez-Aceves, Y . Lin, A. Torres-Orjuela, Sgr a* spin and mass esti- mates through the detection of an extremely large mass-ratio inspiral, The Astrophysical Journal 952 (2) (2023) 139

    2023

  12. [20]

    A. Sekhar, Subrahmanyan Chandrasekhar: A man of stellar physics and stellar grace, https://thewire.in/the-sciences/ subrahmanyan-chandrasekhar-a-man-of-stellar-physics-and-stellar-grace , accessed: 31/10/2023. 15

    2023

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.