pith. sign in

arxiv: 2606.25908 · v1 · pith:MLYAFUFInew · submitted 2026-06-24 · 🧮 math.NA · cs.NA

Spectral deferred corrections parallelized across the method for differential-algebraic equations

Matthias Bolten , Lisa Wimmer This is my paper

Pith reviewed 2026-06-25 19:34 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords spectral deferred correctionsdifferential-algebraic equationsparallel methodsindex-one DAEscollocation methodsnumerical integration
0
0 comments X

The pith

Parallel spectral deferred correction methods solve index-one DAEs faster than sequential versions while keeping high accuracy.

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

The paper tests spectral deferred correction methods on initial value problems for differential-algebraic equations of index one. It applies a parallelization across the method that uses a number of processes equal to the collocation nodes and employs analytically derived optimal diagonal coefficients. The work shows that this parallel version runs faster than the matching sequential SDC scheme and still reaches high accuracy in the small-scale parallel regime. A reader would care because many simulation tasks involve DAEs and modest parallel speedups without loss of accuracy or reformulation could reduce run times on modest hardware.

Core claim

The parallel SDC method with optimal diagonal coefficients solves DAE problems of index one with high accuracy and obtains speedup compared to sequential SDC schemes in a small-scale parallel regime where the number of processes matches the number of collocation nodes.

What carries the argument

Spectral deferred corrections parallelized across the method, using optimal diagonal coefficients in the correction equations applied to index-one DAEs.

If this is right

  • Parallel SDC achieves measurable speedup over sequential SDC for index-one DAE problems.
  • High accuracy is retained when the optimal coefficients are used without reformulation.
  • Small-scale parallelism limited to the number of collocation nodes is sufficient to obtain the reported gains.

Where Pith is reading between the lines

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

  • The same coefficients might allow parallel SDC on certain higher-index DAEs if stability carries over.
  • Combining method-level parallelism with spatial or other parallel techniques could produce larger overall speedups.
  • The approach could be checked on application-derived DAE systems from multibody or circuit models to measure practical gains.

Load-bearing premise

The optimal diagonal coefficients derived analytically for ordinary differential equations remain effective and stable when the same SDC iteration is applied to index-one differential-algebraic equations without additional stabilization or reformulation.

What would settle it

A test on a standard index-one DAE problem where the parallel SDC method loses stability or accuracy relative to the sequential version.

read the original abstract

In this work, the performance of different spectral deferred corrections (SDC) methods applied to initial value problems for differential-algebraic equations (DAEs) of index one is analyzed. The SDC method solves a series of correction equations, and after each iteration, the numerical solution is corrected by adding the obtained approximation of the error. The formulation of the SDC method allows for a parallelization across the method to obtain small-scale parallelism, for which a number of processes equal to the number of collocation nodes can be used. Recently, an analytical approach to finding optimal diagonal coefficients for parallel SDC was proposed. So far, no analysis with the new coefficients for parallel SDC applied to DAEs was done. We demonstrate that parallel SDC methods solve DAE problems with high accuracy faster than the associated sequential SDC schemes, obtaining speedup in a small-scale parallel regime.

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

2 major / 0 minor

Summary. The manuscript analyzes the performance of spectral deferred correction (SDC) methods for initial-value problems in index-one differential-algebraic equations. It applies recently derived optimal diagonal coefficients (originally obtained for ODEs) inside a parallel-across-the-method formulation that uses one process per collocation node, and reports that the resulting parallel SDC schemes achieve high accuracy with measurable speedup over the corresponding sequential SDC integrators in a small-scale parallel regime.

Significance. If the numerical evidence and stability properties hold, the work would supply a concrete route to small-scale parallelism for high-order DAE integrators without requiring additional stabilization or reformulation. Such methods are relevant to constrained mechanical systems and circuit simulation, where index-one DAEs appear routinely and where the algebraic constraints can degrade the contraction rate of standard deferred-correction iterations.

major comments (2)
  1. [Abstract] The central construction re-uses the analytically optimal diagonal coefficients derived for ODEs inside the parallel SDC iteration for semi-explicit index-one DAEs. Because the algebraic constraint couples the error equation at each collocation node, the contraction mapping and the expected convergence order are not guaranteed to survive without a separate DAE-specific stability or order analysis; no such re-derivation or modified-coefficient calculation is supplied.
  2. [Abstract] The abstract asserts that the parallel schemes solve DAE problems “with high accuracy” and obtain speedup, yet the manuscript contains no convergence proof, no DAE error-propagation estimate, and no statement that the algebraic variables were accounted for when the diagonal entries were chosen. This omission is load-bearing for the claim that the ODE-optimal coefficients remain effective.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for highlighting the distinction between the ODE-derived coefficients and their application to index-1 DAEs. Our manuscript is an empirical study of parallel performance; we address the two major comments below.

read point-by-point responses

  1. Referee: [Abstract] The central construction re-uses the analytically optimal diagonal coefficients derived for ODEs inside the parallel SDC iteration for semi-explicit index-one DAEs. Because the algebraic constraint couples the error equation at each collocation node, the contraction mapping and the expected convergence order are not guaranteed to survive without a separate DAE-specific stability or order analysis; no such re-derivation or modified-coefficient calculation is supplied.

    Authors: We agree that the algebraic constraints introduce coupling not present in the ODE case and that a dedicated stability analysis for DAEs would be desirable. The present work re-uses the ODE-optimal diagonal coefficients without modification and evaluates the resulting schemes through numerical experiments on standard index-1 test problems. No DAE-specific re-derivation is performed because the manuscript scope is limited to demonstrating practical accuracy and small-scale parallel speedup; such an analysis would constitute a separate theoretical contribution. revision: no

  2. Referee: [Abstract] The abstract asserts that the parallel schemes solve DAE problems “with high accuracy” and obtain speedup, yet the manuscript contains no convergence proof, no DAE error-propagation estimate, and no statement that the algebraic variables were accounted for when the diagonal entries were chosen. This omission is load-bearing for the claim that the ODE-optimal coefficients remain effective.

    Authors: The abstract claims are grounded in the numerical results shown in the paper, which document both accuracy and wall-clock speedup relative to sequential SDC on the chosen test set. The manuscript does not contain a convergence proof or DAE-specific error estimate. We can revise the abstract and introduction to state explicitly that effectiveness is demonstrated numerically rather than established by a DAE-tailored analysis. revision: partial

Circularity Check

0 steps flagged

No circularity; empirical speedup claims rest on independent timing measurements

full rationale

The paper applies analytically derived diagonal coefficients from prior ODE work to index-1 DAE problems and reports measured wall-clock speedups versus sequential SDC. No equation or claim reduces a reported result to a fit on the same data, a self-definition, or a self-citation chain whose validity is presupposed. The central demonstration consists of numerical experiments whose outcomes are falsifiable by re-running the timings; the coefficient derivation itself is external to this manuscript and does not appear inside the reported performance numbers. Self-citation of the coefficient method is present but does not carry the speedup claim.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; the work relies on the prior derivation of optimal diagonal coefficients whose validity for DAEs is the central untested premise.

pith-pipeline@v0.9.1-grok · 5676 in / 1013 out tokens · 16106 ms · 2026-06-25T19:34:58.979901+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

17 extracted references · 15 canonical work pages · 1 internal anchor

  1. [1]

    and Petzold, Linda R

    Ascher, Uri M. and Petzold, Linda R. , title =. 1998 , isbn =. doi:10.1137/1.9781611971392 , language =

    work page doi:10.1137/1.9781611971392 1998

  2. [2]

    Solving Ordinary Differential Equations

    Hairer, Ernst and Wanner, Gerhard , title =. 2010 , isbn =. doi:10.1007/978-3-642-05221-7 , language =

    work page doi:10.1007/978-3-642-05221-7 2010

  3. [3]

    Andrews, G. C. and Ormrod, M. K. , title =. Proceedings of the ASME Design Engineering Technical Conference (DETC) , year =

  4. [4]

    2025 , howpublished =

    Benabdallah, Seyyid Ali and Souilah, Messoud , title =. 2025 , howpublished =. doi:10.48550/arXiv.2411.15658 , keywords =

    work page doi:10.48550/arxiv.2411.15658 2025

  5. [5]

    Gear, C. W. , title =. Calcolo , fjournal =. 1988 , issn =. doi:10.1007/BF02575744 , language =

    work page doi:10.1007/bf02575744 1988

  6. [6]

    Bit Numerical Mathematics , author =

    Dutt, Alok and Greengard, Leslie and Rokhlin, Vladimir , title =. 2000 , issn =. doi:10.1023/A:1022338906936 , language =

    work page doi:10.1023/a:1022338906936 2000

  7. [7]

    Huang, Jingfang and Jia, Jun and Minion, Michael , title =. J. Comput. Phys. , volume =. 2006 , issn =. doi:10.1016/j.jcp.2005.10.004 , language =

    work page doi:10.1016/j.jcp.2005.10.004 2006

  8. [8]

    Huang, Jingfang and Jia, Jun and Minion, Michael , title =. J. Comput. Phys. , volume =. 2007 , issn =. doi:10.1016/j.jcp.2006.06.040 , language =

    work page doi:10.1016/j.jcp.2006.06.040 2007

  9. [9]

    Discrete Contin

    Xia, Yinhua and Xu, Yan and Shu, Chi-Wang , title =. Discrete Contin. Dyn. Syst., Ser. B , issn =. 2007 , language =. doi:10.3934/dcdsb.2007.8.677 , keywords =

    work page doi:10.3934/dcdsb.2007.8.677 2007

  10. [10]

    2015 , issn =

    Weiser, Martin , title =. 2015 , issn =. doi:10.1007/s10543-014-0540-y , language =

    work page doi:10.1007/s10543-014-0540-y 2015

  11. [11]

    Speck, Robert , title =. Comput. Vis. Sci. , volume =. 2018 , issn =. doi:10.1007/s00791-018-0298-x , zbl =

    work page doi:10.1007/s00791-018-0298-x 2018

  12. [12]

    SIAM J. Sci. Comput. , volume =. 2025 , issn =. doi:10.1137/24M1649800 , language =

    work page doi:10.1137/24m1649800 2025

  13. [13]

    Parallel Performance of Shared Memory Parallel Spectral Deferred Corrections , journal =

    Freese, Philip and G. Parallel Performance of Shared Memory Parallel Spectral Deferred Corrections , journal =. 2025 , issn =. doi:10.1177/10943420251400406 , publisher =

    work page doi:10.1177/10943420251400406 2025

  14. [14]

    On the analysis of spectral deferred corrections for differential-algebraic equations of index one

    Bolten, Matthias and Wimmer, Lisa , title =. 2026 , howpublished =. doi:10.48550/arXiv.2601.16744 , note =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2601.16744 2026

  15. [15]

    Dalc. Adv. Water Resour. , fjournal =. 2011 , issn =. doi:10.1016/j.advwatres.2011.04.013 , language =

    work page doi:10.1016/j.advwatres.2011.04.013 2011

  16. [16]

    2025 , month = apr, version =

    Speck, Robert and Lunet, Thibaut and Baumann, Thomas and Wimmer, Lisa and Akramov, Ikrom and Rosilho de Souza, Giacomo and Fritz, Jakob and Shipton, Jemma , title =. 2025 , month = apr, version =. doi:10.5281/zenodo.15196003 , url =

    work page doi:10.5281/zenodo.15196003 2025

  17. [17]

    Pleiades High-Performance Computing Cluster , year =