pith. sign in

arxiv: 2604.13785 · v1 · submitted 2026-04-15 · 🧮 math.NA · cs.NA

Pathwise convergence of a linearization scheme for stochastic differential-algebraic equations under the local Lipschitz coefficients

Guy Tsafack , Antoine Tambue This is my paper

Pith reviewed 2026-05-10 12:36 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords stochastic differential-algebraic equationsindex-1 SDAElocal linearizationpathwise convergencenumerical methodslocal Lipschitz conditionKhasminskii growth
0
0 comments X

The pith

A local linearization scheme for index-1 stochastic differential-algebraic equations converges pathwise at rate 1/2 minus any positive epsilon.

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

The paper introduces a numerical method for index-1 SDAEs whose drift and diffusion coefficients satisfy local Lipschitz and Khasminskii-type growth conditions. It decomposes the drift locally by Taylor expansion into a linear component solved implicitly to handle the singular matrix at each step and a nonlinear component treated explicitly. Under these conditions the scheme converges in the pathwise sense with rate 1/2 minus epsilon for arbitrary epsilon greater than zero. The approach is presented as efficient for high-dimensional problems. A reader would care because it supplies a practical way to simulate constrained stochastic systems without requiring global Lipschitz assumptions on the coefficients.

Core claim

The novel local linearization scheme for index-1 stochastic differential-algebraic equations converges pathwise with rate 1/2 minus epsilon for any epsilon greater than zero when the coefficients obey the local Lipschitz and Khasminskii growth conditions and the equation remains index-1 at every time despite a singular leading matrix.

What carries the argument

The local linearization technique that splits the drift via Taylor expansion into an implicitly treated linear part resolving the singularity and an explicitly treated nonlinear part.

If this is right

  • The scheme can be implemented efficiently in high dimensions.
  • It resolves the singularity issue at each time step without requiring a global Lipschitz condition.
  • Pathwise error bounds of nearly order 1/2 hold for arbitrary positive epsilon.
  • Numerical experiments can be used to verify the predicted convergence rate in practice.

Where Pith is reading between the lines

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

  • The splitting idea might extend to index-2 or higher SDAEs if the index reduction step is adapted accordingly.
  • The method could be paired with adaptive step-size control to improve practical performance on stiff trajectories.
  • Similar linearization steps may apply to other singular stochastic systems arising in constrained optimization or circuit simulation.

Load-bearing premise

The SDAE must remain index-1 at each time even though the leading matrix is singular, and the coefficients must satisfy local Lipschitz and Khasminskii-type growth conditions.

What would settle it

A counterexample computation or numerical test in which the pathwise convergence rate falls below 1/2 minus epsilon while the index-1 property holds but the local Lipschitz or Khasminskii condition is violated.

Figures

Figures reproduced from arXiv: 2604.13785 by Antoine Tambue, Guy Tsafack.

Figure 1
Figure 1. Figure 1: An example of an electrical circuit with [PITH_FULL_IMAGE:figures/full_fig_p030_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Pathwise convergence with our semi-implicit linearization scheme with three [PITH_FULL_IMAGE:figures/full_fig_p033_2.png] view at source ↗
read the original abstract

The paper deals with the numerical treatment of index-1 stochastic differential-algebraic equations (SDAEs) with nonlinear coefficients that satisfy the local Lipschitz and the Khasminskii conditions. The key challenge here is the presence of a singular and non-autonomous matrix in the equation, which makes the numerical method challenging to analyze. To tackle this challenge, we develop a more general numerical method using a local linearization technique. More precisely, we use the Taylor expansion to decompose locally the drift component of the SDAEs in linear and nonlinear parts. The linear part is approximated implicitly and must resolve the singularity issue of each time step, while the nonlinear part is approximated explicitly. This method is fascinating due to the fact that it is efficient in high dimension. We prove that this novel numerical method converges in the pathwise sense with rate $\frac{1}{2}-\epsilon$, for arbitrary $\epsilon >0$. The implementation of this novel numerical method is also carried out to verify our theoretical result.

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

0 major / 3 minor

Summary. The manuscript develops a local-linearization numerical scheme for index-1 SDAEs whose drift and diffusion coefficients satisfy only local Lipschitz continuity together with a Khasminskii-type growth bound. The drift is split via Taylor expansion into an implicitly treated linear piece (chosen to respect the singular, time-dependent leading matrix) and an explicitly treated nonlinear remainder; the resulting scheme is shown to be well-defined at each step and to converge pathwise with rate 1/2−ε for every ε>0. Numerical experiments are included to illustrate the theoretical rate.

Significance. If the proof is complete, the result extends the known pathwise convergence theory for SDEs under local Lipschitz/Khasminskii hypotheses to the SDAE setting, where the singular non-autonomous matrix introduces additional analytic difficulties. The local-linearization approach is a natural and efficient way to retain solvability of the implicit step while avoiding global Lipschitz assumptions, and the paper supplies both the convergence analysis and reproducible numerical verification.

minor comments (3)
  1. The abstract asserts the pathwise rate 1/2−ε but supplies no outline of the key estimates (control of the linearization remainder, preservation of the index-1 property, or application of the Khasminskii condition). A one-sentence sketch of the main argument would improve accessibility.
  2. Notation for the local linearization operator and the implicit algebraic solve at each time step should be introduced with an explicit equation reference (e.g., the definition of the linear coefficient matrix and the resulting linear SDAE) so that the reader can follow the well-posedness argument without ambiguity.
  3. The numerical section would benefit from a brief statement of the exact discretization parameters (step-size sequence, number of paths, and how the implicit solve is implemented) to allow direct reproduction of the reported convergence plots.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our work on the pathwise convergence of the local-linearization scheme for index-1 SDAEs under local Lipschitz and Khasminskii conditions. The recommendation for minor revision is appreciated. No specific major comments were listed in the report, so we have performed a careful re-reading and incorporated minor editorial improvements for clarity, notation consistency, and reproducibility of the numerical section.

Circularity Check

0 steps flagged

No significant circularity; proof is self-contained under stated assumptions

full rationale

The paper establishes pathwise convergence of a local-linearization scheme for index-1 SDAEs by Taylor-splitting the drift into an implicitly treated linear piece (respecting the singular leading matrix) and an explicit nonlinear remainder, then proving the rate 1/2−ε under local Lipschitz continuity plus Khasminskii growth. This chain uses standard stochastic-analysis estimates and does not reduce the claimed rate or well-posedness to any fitted quantity, self-definition, or load-bearing self-citation; the index-1 property is preserved by construction of the scheme and the hypotheses, with no renaming of known results or ansatz smuggling. The result is therefore independent of its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on two domain assumptions standard in the SDAE literature plus the index-1 regularity of the leading matrix; no free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption Drift and diffusion coefficients satisfy local Lipschitz continuity and Khasminskii-type growth
    Explicitly required in the abstract for the equations under study.
  • domain assumption The SDAE is of index 1 with singular non-autonomous leading matrix
    Stated as the key structural property that the linearization must resolve at each step.

pith-pipeline@v0.9.0 · 5477 in / 1286 out tokens · 28256 ms · 2026-05-10T12:36:40.057445+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

36 extracted references · 36 canonical work pages

  1. [1]

    Oana-Silvia Serea, Antoine Tambue, Guy Tsafack Existence and unique- ness for the solutions of non-autonomous stochastic differential al- gebraic equations with locally Lipschitz coefficients , Nonlinear Dif- ferential Equations and Applications NoDEA; vol 32, num 4, p 69, https://doi.org/10.1007/s00030-025-01079-x springer (2025)

    work page doi:10.1007/s00030-025-01079-x 2025

  2. [2]

    Tsafack Guy and Tambue Antoine , Pathwise convergence of a novel nu- merical scheme based on semi-implicit method for stochastic differential- algebraic equations with non-global Lipschiz coefficients , arXiv preprint arXiv:2506.15627 (2025)

    work page arXiv 2025

  3. [3]

    Tsafack Guy and Tambue Antoine , Strong convergence of a semi tamed scheme for stochastic differential algebraic equation under non-global Lipschitz coefficients , arXiv preprint arXiv:2509.09032 (2025)

    work page arXiv 2025

  4. [4]

    R. Z. Khasminskii , Khasminskii R Z. Stochastic Stability of Differential Equations. vol Maryland: Sijthoff and Noordhoff, springer, (1980). 33

    work page 1980

  5. [5]

    Journal of Computational and Applied Mathe- matics, vol 290, p 370-384, Elsevier (2015)

    Mao Xuerong, The truncated Euler– Maruyama method for stochastic differential equations . Journal of Computational and Applied Mathe- matics, vol 290, p 370-384, Elsevier (2015)

    work page 2015

  6. [6]

    Journal of computational and applied mathematics, vol 157, num 2, p 477-505, Elsevier (2003)

    Winkler Renate, Stochastic differential-algebraic equations of index 1 and applications in circuit simulation . Journal of computational and applied mathematics, vol 157, num 2, p 477-505, Elsevier (2003)

    work page 2003

  7. [7]

    Potential Analysis, vol 8, p 205-216, Springer(1998)

    Gyöngy István, A note on Euler’s approximations . Potential Analysis, vol 8, p 205-216, Springer(1998)

    work page 1998

  8. [8]

    Probability theory and related fields, vol 105, num 2, p 143-158, Springer (1996)

    Gyöngy István and Krylov Nicolai, Existence of strong solutions for Itô’s stochastic equations via approximations . Probability theory and related fields, vol 105, num 2, p 143-158, Springer (1996)

    work page 1996

  9. [9]

    Commu- nications in Nonlinear Science and Numerical Simulation, vol 94, p 105- 574, Elsevier(2021)

    İzgi Burhaneddin and Çetin Coşkun,Strong convergence of semi-implicit split-step methods for SDE with locally Lipschitz coefficients . Commu- nications in Nonlinear Science and Numerical Simulation, vol 94, p 105- 574, Elsevier(2021)

    work page 2021

  10. [10]

    Journal of Computational and Applied Mathematics, vol 313, p 202-217, Elsevier(2017)

    Jimenez Juan Carlos and Mora C and Selva Mónica, A weak Local Lin- earization scheme for stochastic differential equations with multiplicative noise. Journal of Computational and Applied Mathematics, vol 313, p 202-217, Elsevier(2017)

    work page 2017

  11. [11]

    Stochas- tic Processes and their Applications, vol 130, num 8, p 4968-5005, 2020, Elsevier(2020)

    Mukam Jean Daniel and Tambue Antoine, Strong convergence of a stochastic Rosenbrock-type scheme for the finite element discretization of semilinear SPDEs driven by multiplicative and additive noise . Stochas- tic Processes and their Applications, vol 130, num 8, p 4968-5005, 2020, Elsevier(2020)

    work page 2020

  12. [12]

    Higham Desmond J and Mao Xuerong and Stuart Andrew M, Strong convergence of Euler-type methods for nonlinear stochastic differential equations, SIAM journal on numerical analysis, vol 40, num 3, p 1041- 1063, SIAM(2002)

    work page 2002

  13. [13]

    Else- vier(2007) 34

    Mao, X., Stochastic differential equations and applications . Else- vier(2007) 34

    work page 2007

  14. [14]

    Numerical Algorithms, vol 89 (2), p 721-747, Springer(2022)

    Kelly Cónall and Lord Gabriel J, Adaptive Euler methods for stochastic systems with non-globally Lipschitz coefficients . Numerical Algorithms, vol 89 (2), p 721-747, Springer(2022)

    work page 2022

  15. [15]

    BIT Numerical Mathematics, vol 2, p 437-455, Springer(2012)

    Küpper Dominique and Kværnø Anne and Rößler Andreas, A Runge- Kutta method for index 1 stochastic differential-algebraic equations with scalar noise . BIT Numerical Mathematics, vol 2, p 437-455, Springer(2012)

    work page 2012

  16. [16]

    Springer, vol 1409 (2006)

    Hairer Ernst and Lubich Christian and Roche Michel, The numeri- cal solution of differential-algebraic systems by Runge-Kutta methods . Springer, vol 1409 (2006)

    work page 2006

  17. [17]

    Springer Berlin Heidelberg New York, vol 375 , (1996)

    Wanner Gerhard and Hairer Ernst, Solving ordinary differential equa- tions II . Springer Berlin Heidelberg New York, vol 375 , (1996)

    work page 1996

  18. [18]

    Chapman and Hall/CRC (2018)

    Lawler Gregory F, Introduction to stochastic processes . Chapman and Hall/CRC (2018)

    work page 2018

  19. [19]

    World Scientific Publishing Company (2012)

    Klebaner Fima C, Introduction to stochastic calculus with applications . World Scientific Publishing Company (2012)

    work page 2012

  20. [20]

    Schein Oliver and Denk Georg, Numerical solution of stochastic differential-algebraic equations with applications to transient noise simu- lation of microelectronic circuits, Journal of Computational and Applied Mathematics, vol 100(1), p 77-92, Elsevier (1998)

    work page 1998

  21. [21]

    Preprint TUM-M9907, Technische Universität, München, Cite- seer(1999)

    Penski Christian, Analysis and numerical integration of stochastic differential-algebraic equations and applications in electronic circuit sim- ulation. Preprint TUM-M9907, Technische Universität, München, Cite- seer(1999)

    work page 1999

  22. [22]

    Math., vol 38, num 1, p 117-131, (2010)

    Cong Nguyen Dinh and The Nguyen Thi, Stochastic differential algebraic equation of index 1 , Vietnam J. Math., vol 38, num 1, p 117-131, (2010)

    work page 2010

  23. [23]

    Texts in applied mathematics, vol 60, Springer(2014)

    Pavliotis Grigorios A, Stochastic processes and applications . Texts in applied mathematics, vol 60, Springer(2014)

    work page 2014

  24. [24]

    Chapman and Hall/CRC (2018)

    Guttorp Peter, Stochastic modeling of scientific data . Chapman and Hall/CRC (2018). 35

    work page 2018

  25. [25]

    Springer Science & Business Media (2013)

    Grigoriu Mircea, Stochastic calculus: applications in science and engi- neering. Springer Science & Business Media (2013)

    work page 2013

  26. [26]

    World Scientific (2006)

    Henderson Douglas and Plaschko Peter, Stochastic Differential Equa- tions In Science And Engineering (With Cd-rom) . World Scientific (2006)

    work page 2006

  27. [27]

    Springer (2019)

    Campbell Stephen and Ilchmann Achim and Mehrmann Volker and Reis Timo and others, Applications of differential-algebraic equations: exam- ples and benchmarks . Springer (2019)

    work page 2019

  28. [28]

    Shaker (1999)

    Schein Oliver, Stochastic differential-algebraic equations in circuit sim- ulation. Shaker (1999)

    work page 1999

  29. [29]

    Springer (2008)

    Denk Georg and Römisch Werner and Sickenberger Thorsten and Win- kler Renate, Efficient transient noise analysis in circuit simulation . Springer (2008)

    work page 2008

  30. [30]

    Part I: Mathematical structure and index of network equations

    Gœnther M and Feldmann U, CAD based electric circuit modeling in industry. Part I: Mathematical structure and index of network equations . Surveys on Mathematics for Industry, vol 8, (2008)

    work page 2008

  31. [31]

    Springer Science & Business Media, Vol 9 (2007)

    Anile Angelo Marcello and Alì Giuseppe and Mascali Giovanni, Scien- tific computing in electrical engineering . Springer Science & Business Media, Vol 9 (2007)

    work page 2007

  32. [32]

    Springer Science & Business Media vol 425 (2012)

    Demir Alper and Sangiovanni-Vincentelli Alberto, Analysis and simu- lation of noise in nonlinear electronic circuits and systems . Springer Science & Business Media vol 425 (2012)

    work page 2012

  33. [33]

    vol 10, p 235-253, Cam- bridge University Press (2007)

    Kloeden Peter E and Neuenkirch Andreas, The pathwise convergence of approximation schemes for stochastic differential equations , LMS journal of Computation and Mathematics. vol 10, p 235-253, Cam- bridge University Press (2007)

    work page 2007

  34. [34]

    London Mathematical Society Lecture Note Series, vol 363, p 140-161, Cambridge University Press (2009)

    Jentzen Arnulf Kloeden and Peter E Neuenkirch Andreas, Pathwise con- vergence of numerical schemes for random and stochastic differential equations. London Mathematical Society Lecture Note Series, vol 363, p 140-161, Cambridge University Press (2009). 36

    work page 2009

  35. [35]

    IMA Journal of Numerical Analysis, vol 32, num 2, p 701-723, Oxford University Press (2012)

    Nguyen Son Luu and Yin G, Pathwise convergence rate for numerical solutions of stochastic differential equations . IMA Journal of Numerical Analysis, vol 32, num 2, p 701-723, Oxford University Press (2012)

    work page 2012

  36. [36]

    Römisch, Werner and Winkler, Renate, title=Stochastic DAEs in cir- cuit simulation, Modeling, Simulation, and Optimization of Integrated Circuits: Proceedings of a Conference held at the Mathematisches Forschungsinstitut, Oberwolfach, November 25-December 1, 2001, p 303- 318, 2003, Springer(2003) 37

    work page 2001