Pith. sign in

REVIEW 3 major objections 4 minor 95 references

Flux-Corrected Diagonal Frog: second order and positivity at all time steps

T0 review · 3 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash

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

Pith's one-line read A nonlinear flux limiter lets one-dimensional Fokker-Planck solvers stay positive and second-order at every time step.

desk verdict A genuine new construction with clean small-step theory, but the 'all time steps' claim leans on an unproven active-set nonsingularity; worth refereeing. read the letter →

arxiv 2607.20415 v1 pith:H2T5EWAA submitted 2026-07-22 math.NA cs.CEcs.NAphysics.comp-phq-fin.CP

Andrey Itkin This is my paper · ORCID
classification math.NAcs.CEcs.NAphysics.comp-phq-fin.CP MSC 65M0665M1235Q84
keywords Fokker-Planckequationpositivitypreservationflux-correctedtransportM-matrixsplitsecond-orderaccuracymassconservationdefectcorrectionactive-setNewton
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

The paper's central claim is that the Diagonal Frog finite-difference framework, augmented by a nonlinear flux limiter acting inside the implicit solve, can preserve positivity of the computed density and exact discrete mass at every time step while retaining second-order accuracy wherever the density is resolved. The scheme splits the second-order directional operator into a monotone M-matrix core and an antidiffusive flux correction, then clamps each interface flux so no node can donate more than half its budget. Because limiting is driven by the sign of the right-hand side rather than the cell Péclet number, the limiter stays inactive in smooth regions and the global L1 error stays O(h^2) even at large Péclet; the price of unconditional positivity is paid only at unresolved fronts, and even there it is confined to a few cells. A defect-corrected two-stage variant restores second-order accuracy in time, and an active-set solver covers arbitrarily large steps. If the claims hold, the method removes the small-time-step limitation of the earlier Diagonal Frog solvers.

What carries the argument

The load-bearing object is the split A2 = A1 + C: A1 is the tridiagonal M-matrix core (central diffusion plus first-order upwind convection), and C is the antidiffusive correction written as a flux difference. A per-interface limiter with two-point caps limits each donor node's withdrawal to half its current budget, so the clamped flux is the nonexpansive clamp of the unlimited flux onto a fixed interval. The core resolvent (I - γA1)^{-1} has ℓ1 norm exactly one, which gives the contraction constant q = 2γ‖μ‖∞/h and the threshold γ_pic = h/(2‖μ‖∞). The fixed-point equation is read row by row against the unlimited second-order system; the defect identity Δt²/2 A2² Y0 provides the second-order

What would settle it

Take a smooth, strictly positive advection-diffusion density (e.g. a Gaussian) with cell Péclet number around 4 on meshes n = 101..801, and measure the L1 order of the FCDF-B fixed point against the semi-discrete exact solution; the claim fails if the order drops below about 1.9 while the unlimited second-order operator measures 1.99. A second check looks at the limited region: if any active limiter interface carries density scale above the truncation level on such smooth data, the localization assumption is violated.

Watch

Extended reading notes

Core claim

The paper's central claim is that the Diagonal Frog second-order directional discretization, split into a monotone M-matrix core and an antidiffusive flux correction with per-interface cap limiting inside the implicit solve, is unconditionally positive and exactly mass-conservative for every time step. The fixed point of the limited iteration solves the full second-order system at every node where the flux caps are slack; because the limiter is driven by the sign of the right-hand side rather than the cell Péclet number, it stays inactive in resolved regions, yielding a global L1 error O(h^2) uniformly in Péclet except at order-one unresolved fronts, where it degrades locally to O(h). The de

Load-bearing premise

The uniform second-order statement rests on the assumption that in every resolved region the limiter caps remain slack, i.e. the fixed-point density is comparable to the budget in the sense of Eq. (A.6); the paper verifies this a posteriori on examples but leaves the a priori Harnack-type estimate open.

Editorial extensions

If this is right

  • Positivity and exact mass conservation become unconditional: no minimum time step, no threshold in γ, and no limiter-value dependence.
  • On smooth, strictly positive densities the scheme converges at O(h^2) uniformly in the cell Péclet number; at an unresolved order-one front it falls to O(h) only in a few cells, unlike exponential-fitting schemes that degrade globally.
  • Second-order temporal accuracy is achievable from backward Euler solves alone via defect correction, with local error O(Δt^3) wherever the caps are slack.
  • Positivity and conservation cover every step size when the computed Picard threshold exceeds the smaller linear positivity threshold; full second order at every step follows when it also exceeds the Padé threshold, with the active-set solver closing the gap when it does not.

Reading between the lines

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

  • The unresolved-front O(h) floor is stated as an observed floor, not a proved lower bound; a front-tracking or subcell reconstruction extension might beat it, since no lower-bound theorem for nonlinear positivity-preserving schemes is known.
  • The a priori localization claim (budgets comparable to density) is left open; a Harnack-type estimate for the core resolvent would turn the a posteriori check into a theorem, and its absence is the main condition on which the uniform second-order claim depends.
  • Because the limiter works on fluxes and the solver is one-dimensional, the scheme should compose cleanly with directional splitting for multidimensional Fokker-Planck problems; the mass-conservation property would then transfer to each sweep.
  • In option pricing, a density that stays nonnegative at every time step removes butterfly-arbitrage violations in local-volatility calibration without restricting the time step; this is a direct but unstated practical consequence.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper proposes a nonlinear extension of the Diagonal Frog finite-difference framework for one-dimensional Fokker–Planck equations. The spatial operator is split into a monotone M-matrix core and an antidiffusive flux correction; a Zalesak-type limiter is applied per interface inside implicit banded solves. Two variants are analyzed: FCDF-A (global stopping rule) and FCDF-B (per-interface limiting), plus a defect-corrected two-stage time stepper FCDF-DC. The main claimed results are: unconditional positivity and exact mass conservation for every time step and limiter value; ℓ1 contraction of the Picard iteration under a purely convective step restriction γ < h/(2μ̄); Péclet-uniform second-order L1 accuracy modulo a layer term; and coverage of all step sizes by handing over to the linear resolvent/Padé windows of a companion paper, with an active-set semismooth Newton solver filling the remaining gap. Numerical experiments on Ornstein–Uhlenbeck and advection-dominated benchmarks are reported.

Significance. If the main claims held as stated, this would be a significant advance: it would provide a positivity-preserving, mass-conservative, second-order scheme that avoids the global first-order degradation of Chang–Cooper and extends to arbitrary time steps without violating Godunov-type barriers, by making the nonlinearity local and adaptive. The paper has real strengths: the proof of unconditional positivity and conservation in Proposition 1(i)–(ii) is clean and self-contained; the ℓ1 contraction estimate is explicit and correctly isolates the convective constant; the defect-correction identity in Lemma 1 and the resulting stability function are elegant; and the two-term error decomposition in Proposition 2 is a sensible way to separate smooth and layer contributions. The paper is also unusually candid in stating its own open questions, which is commendable. However, the two headline claims — 'every time step' and 'second order uniformly in Péclet' — are not fully established: one rests on an invalid block-diagonalization argument for the active-set Newton matrices, and the other on an a posteriori budget–density comparability that is explicitly left open.

major comments (3)
  1. [Appendix D, Lemma 2 / Proposition 5(b)] The block-diagonalization is invalid. A clamped interface i+1/2 removes only the correction C_S at that interface; the monotone core A1 is tridiagonal and still couples nodes i and i+1 across every interface, clamped or free. After the permutation in Eq. (D.1), the off-diagonal blocks between free and clamped blocks contain these A1 entries and do not vanish, so V_S is not block diagonal as claimed in Eq. (D.2). Consequently the Neumann-series bound on the diagonal blocks does not imply nonsingularity, and Proposition 5(b) is unproved. This is load-bearing: because Eq. (23) fails on the tested meshes, the active-set solver is the only large-step second-order mechanism, and the 'every time step' claim rests on an empirical no-failure statement rather than on a proof.
  2. [Proposition 1(iv), Eq. (A.6), Section 7] The Péclet-uniform L1 second-order claim is conditional on budget–density comparability b_j ~ p*_j (Eq. (A.6)), which the paper explicitly leaves open and verifies only a posteriori. The abstract's statement that the method 'ensures the global L1 convergence remains second-order uniformly in the cell Péclet number' overstates what is proved: Proposition 2 assumes the layer-structure and layer-strength hypotheses, and no Harnack-type estimate is supplied to show that caps remain slack in resolved regions. The localization claim should either be proved or stated in the abstract/introduction with the same conditionality as in Section 7.
  3. [Section 5.2, Eqs. (24)–(25)] For γ > γ_pic, no existence or uniqueness of a zero of the piecewise-linear residual F is proved; the contraction argument stops at γ_pic, and Proposition 5(c) covers only the all-free and all-clamped patterns. Since Lemma 2 fails as a proof of nonsingularity for mixed patterns, the semismooth Newton iteration is not proved to be well-defined in the regime where the active-set solver is needed. The numerical pattern-update counts in Tables 8–9 are valuable evidence but do not replace a theorem. This gap should be closed or explicitly quarantined from the unconditional claims in the abstract and Section 7.
minor comments (4)
  1. [Section 5.2] The sentence 'Because F is piecewise linear, the iteration terminates as soon as the pattern at the solution is identified; one further solve then lands exactly on the zero.' is repeated verbatim two sentences apart.
  2. [Table 3 and surrounding text] The caption says 'the last two columns' but the columns showing the temporal-floor decay are FCDF-B, unlimited, and CC; the last two numeric columns are CC and core. Please clarify which columns are meant.
  3. [Section 6.1] Typo: 'Pad’e' should be 'Padé'. Also, the claim that the thresholds are 'exact for the OU benchmark by eigendecomposition' is imprecise; the table reports thresholds located by bisection on the assembled matrix, not closed-form expressions.
  4. [Section 3.1 / Proposition 2] The layer-interface set A and the layer-strength assumption S_L are introduced after the informal three-regime taxonomy. Moving the formal assumptions before the taxonomy would make it easier to see exactly which cases Proposition 2 covers and which it excludes.

Circularity Check

0 steps flagged · score 2.0 of 10

No material circularity: the FCDF positivity/conservation proofs are self-contained; the second-order statement is honestly conditional, and the main self-citation supplies a large-step dependency without reducing the derivation to its inputs.

full rationale

The central derivation is not circular. Proposition 1 proves unconditional positivity and exact conservation directly from the flux-form clamp (Eqs. 13-15), the M-matrix core resolvent (Lemma 3), and the nonexpansiveness of the clamp; no accuracy target is used to fit the limiter. The thresholds gamma_pic, gamma_0, and gamma_r are computed from the assembled matrix and mesh, not calibrated to the reported errors, so the coverage claims are not fitted inputs renamed as predictions. The second-order statement is explicitly conditional: Proposition 1(iv) says the fixed point coincides with the unlimited second-order system only where caps are slack, which is definitional in the slack regime, and the paper verifies limiter localization a posteriori (Section 6.3 and Eq. A.6) while leaving a Harnack estimate open (Section 7). That is a conditionality caveat, not a circular reduction. The main self-citation is the companion paper [Itkin and Kazbek, 2026], which supplies the large-step linear positivity windows used in Corollary 4 and Proposition 5(c); this is a dependency for the all-step coverage corollary, but it is not equivalent to the present paper's own inputs, and the thresholds are also measured numerically here. The admitted gap in Proposition 5(c) (mixed-pattern nonsingularity for gamma > gamma_pic remains open, supported only numerically) is a proof gap rather than an input-output circularity.

Assumptions & free parameters 1 free parameters · 6 assumptions · 0 invented entities

The ledger is light: no fitted physics constants and no invented entities. The main burden is the a posteriori budget-density comparability used for localization, plus imported companion-paper windows.

free parameters (1)
  • even budget split kappa_j = 2
    Chosen by hand so each of a node's two adjacent interfaces can withdraw at most half the node's budget. This design choice directly powers the unconditional positivity proof (Prop 1(i)); it is not fitted to data.
assumptions (6)
  • standard math Godunov's theorem and Bolley-Crouzeix temporal barrier: linear second-order positivity-preserving schemes are impossible.
    Used to motivate the need for nonlinearity; not needed as a step in the proofs.
  • standard math The core A1 is an irreducible M-matrix generator with 1^T A1 = 0 under the zero-flux closure, so M = I - gamma A1 is inverse-positive with unit column sums (Lemma 3).
    This is the backbone of the positivity and conservation proofs.
  • standard math Nonexpansiveness of clamping onto a fixed interval containing zero.
    Used in the contraction bound of Prop 1(iii) and Prop 3(iii).
  • domain assumption Existence of a layer set A with Péclet-uniform smoothness bounds M_S and layer strength S_L (Prop 2, assumptions 2 and 3).
    Needed for the O(h^2) Péclet-uniform error bound; explicitly excludes sub-cell boundary accumulation.
  • ad hoc to paper Log-Lipschitz resolution condition and budget-density comparability (Eq. A.6) so that caps are slack in resolved regions.
    The paper states this as an a posteriori check, not a proved theorem; the localization claim rests on it, and Sec. 7 leaves the Harnack estimate open.
  • domain assumption Linear positivity windows for the resolvent and Pade(0,2) maps from the companion papers [Itkin 2026; Itkin and Kazbek 2026].
    Used for the large-step hand-over; the companion paper is self-cited and 'in preparation', not reproduced here.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Flux-Corrected Diagonal Frog: second order and positivity at all time steps." pith.science (2026). https://pith.science/paper/H2T5EWAA

@misc{pith2026260720415,
  author       = {Pith},
  title        = {Pith review of: Flux-Corrected Diagonal Frog: second order and positivity at all time steps},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/H2T5EWAA}},
  note         = {Machine review of arXiv:2607.20415}
}
abstract

By Godunov's theorem, linear second-order finite-difference schemes for the Fokker-Planck equation cannot preserve positivity. The Diagonal Frog (DF) framework previously bypassed this barrier using eventual positivity, but required a strict minimum time step. This paper resolves the small-step limitation using a nonlinear extension of the DF solvers. We split the second-order directional operator into a monotone M-matrix core and an antidiffusive flux correction. A Zalesak-type limiter is then applied iteratively within the implicit banded solve. The resulting Flux-Corrected DF (FCDF) schemes (variants A and B) are unconditionally positive across all time steps. Because the limiter acts on fluxes rather than point values, these schemes conserve discrete mass exactly and maintain second-order accuracy. Crucially, the limiter activates only within unresolved layers. This ensures the global $L_1$ convergence remains second-order uniformly in the cell P\'eclet number, avoiding the first-order degradation seen in the Chang-Cooper scheme. The method's Picard iteration is contractive under a purely convective step restriction. To support arbitrary step sizes, we develop an active-set reformulation. This solves the system using a semismooth Newton iteration, where computational cost scales only with the number of nodes where positivity binds. Finally, we introduce a defect-corrected time stepping approach that restores second-order time accuracy. Numerical experiments on Ornstein-Uhlenbeck and advection-dominated benchmarks confirm our claims.

Figures

Figures reproduced from arXiv: 2607.20415 by the authors.

Figure 1
Figure 1. Coverage thresholds against the mesh on the OU benchmark. The Picard threshold falls linearly in h. The resolvent threshold falls faster and crosses below it, so condition Eq. (22) holds from n = 101 onward. The Padé threshold stays three orders of magnitude above and escapes upward, so condition Eq. (23) fails on every mesh. Open symbols mark meshes on which no Padé window was found below the search bound 106 . h F… view at source ↗
Figure 2
Figure 2. Limiter activity at an unresolved front, D = 10−4 and n = 801, on a logarithmic density scale. The shaded interfaces are those where the clamp is active. They begin only where the density has fallen several orders of magnitude below the plateau, and the advected packet itself is untouched. The clamp is active on 38% of the interfaces, which carry 1.7 × 10−5 of the total mass. 6.4 Long-time front smearing Proposition… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

95 extracted references · 1 canonical work pages

  1. [1]

    and Wu, J

    Gao, Z. and Wu, J. , title =. SIAM Journal on Scientific Computing , volume =

  2. [2]

    McFarland, D. M. and Ye, F. and Zong, C. and Zhu, R. and Han, T. and Fu, H. and others , title =. Mathematics , volume =

  3. [3]

    Park, B. T. and Petrosian, V. , title =. The Astrophysical Journal Supplement Series , year =

  4. [4]

    Communications in Computational Physics , year =

    Qi, Chen and Su, Zhenpeng and Wu, Zhaoguo and Zheng, Huinan and Wang, Yuming , title =. Communications in Computational Physics , year =

  5. [5]

    Journal of Geophysical Research: Space Physics , year =

    Camporeale, Enrico and Delzanno, Gian Luca and Zaharia, Sorin and Koller, Josef , title =. Journal of Geophysical Research: Space Physics , year =

  6. [6]

    , title =

    Albert, Jay M. , title =. Journal of Geophysical Research: Space Physics , year =

  7. [7]

    Chang, J. S. and Cooper, G. , title =. Journal of Computational Physics , year =

  8. [8]

    Journal of Computational Physics , year =

    Harten, Ami , title =. Journal of Computational Physics , year =

Show all 95 references
  1. [9]

    Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences , year =

    Pareschi, Lorenzo and Russo, Giovanni , title =. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences , year =

  2. [10]

    Communications in Computational Physics , year =

    Yuan, Li and Shu, Chi-Wang , title =. Communications in Computational Physics , year =

  3. [11]

    SIAM Journal on Numerical Analysis , year =

    Osher, Stanley , title =. SIAM Journal on Numerical Analysis , year =

  4. [12]

    SIAM Journal on Numerical Analysis , year =

    Osher, Stanley and Chakravarthy, Sukumar , title =. SIAM Journal on Numerical Analysis , year =

  5. [13]

    Mathematics of Computation , year =

    Osher, Stanley and Solomon, Fred , title =. Mathematics of Computation , year =

  6. [14]

    and Lubuma, J

    Chapwanya, M. and Lubuma, J. M. S. , title =. Computers & Mathematics with Applications , year =

  7. [15]

    IMA Journal of Numerical Analysis , year =

    Hu, Jingwei and Zhang, Xiangxiong , title =. IMA Journal of Numerical Analysis , year =

  8. [16]

    Journal of Scientific Computing , year =

    Liu, Chen and Gao, Yuan and Zhang, Xiangxiong , title =. Journal of Scientific Computing , year =

  9. [17]

    and Cates, M

    Tailleur, J. and Cates, M. E. , title =. Physical Review Letters , year =

  10. [18]

    Journal of Geophysical Research: Space Physics , year =

    Gong, Zheng and Fu, Suiyan and Ma, Qiugang and Xiang, Zheng and Ni, Binbin and Gan, Leping , title =. Journal of Geophysical Research: Space Physics , year =

  11. [19]

    A Deep Smoothness

    Kossaczk. A Deep Smoothness. Progress in Industrial Mathematics at. 2021 , pages =

    2021

  12. [20]

    Deep smoothness

    Kossaczk. Deep smoothness. Physics of Fluids , year =

  13. [21]

    Applied Mathematics and Computation , year =

    Hoang, Manh Tuan and Ehrhardt, Matthias , title =. Applied Mathematics and Computation , year =

  14. [22]

    International Journal of Numerical Analysis and Modeling , Year =

    Jumps without Tears: A New Splitting Technology for Barrier Options , Author =. International Journal of Numerical Analysis and Modeling , Year =

  15. [23]

    Itkin , Journal =

    A. Itkin , Journal =. Efficient Solution of Backward Jump-Diffusion. 2016 , Pages =

    2016

  16. [24]

    International Journal of Computer Mathematics , Year =

    LSV models with stochastic interest rates and correlated jumps , Author =. International Journal of Computer Mathematics , Year =

  17. [25]

    Itkin , Journal =

    A. Itkin , Journal =. 2015 , Number =

    2015

  18. [26]

    2018 , address =

    Kwok, Sau Fa , title =. 2018 , address =

    2018

  19. [27]

    1996 , series =

    Risken, Hannes , title =. 1996 , series =

    1996

  20. [28]

    Thomas, J. W. , title =. 1995 , series =

    1995

  21. [29]

    , title =

    Elhashash, Amna and Szyld, Daniel B. , title =. Linear Algebra and its Applications , year =

  22. [30]

    The Electronic Journal of Linear Algebra , year =

    Olesky, Dale and Tsatsomeros, Michael and Van den Driessche, Pauline , title =. The Electronic Journal of Linear Algebra , year =

  23. [31]

    Linear Algebra and its Applications , year =

    Noutsos, Dimitrios , title =. Linear Algebra and its Applications , year =

  24. [32]

    , title =

    Noutsos, Dimitrios and Tsatsomeros, Michael J. , title =. SIAM Journal on Matrix Analysis and Applications , year =

  25. [33]

    , title =

    Berman, Abraham and Plemmons, Robert J. , title =. 1994 , doi =

    1994

  26. [34]

    2009 , publisher=

    Stochastic Methods: A Handbook for the Natural and Social Sciences , author=. 2009 , publisher=

    2009

  27. [35]

    On the construction and comparison of difference schemes , Author =. SIAM J. Numerical Analysis , Year =

  28. [36]

    BIT Numerical Mathematics , volume =

    Jahnke, Tobias and Lubich, Christian , title =. BIT Numerical Mathematics , volume =. 2000 , pages =

    2000

  29. [37]

    Communications on Pure and Applied Mathematics , volume =

    Magnus, Wilhelm , title =. Communications on Pure and Applied Mathematics , volume =. 1954 , pages =

    1954

  30. [38]

    Blanes, Sergio and Casas, Fernando and Oteo, Jos. The. Physics Reports , volume =. 2009 , pages =

    2009

  31. [39]

    An Exact Local Error Representation of Exponential Operator Splitting Methods for Evolutionary Problems and Applications to Linear

    Descombes, St. An Exact Local Error Representation of Exponential Operator Splitting Methods for Evolutionary Problems and Applications to Linear. BIT Numerical Mathematics , volume =. 2010 , pages =

    2010

  32. [40]

    , title =

    Hundsdorfer, Willem and Verwer, Jan G. , title =

  33. [41]

    SIAM Review , volume =

    Moler, Cleve and Van Loan, Charles , title =. SIAM Review , volume =. 2003 , pages =

    2003

  34. [42]

    , title =

    Higham, Nicholas J. , title =. SIAM Review , volume =. 2009 , pages =

    2009

  35. [43]

    Saad, Yousef , title =

  36. [44]

    , title =

    Varah, James M. , title =. Linear Algebra and its Applications , volume =

  37. [45]

    in 't Hout, K. J. and Welfert, B. D. , title =. Applied Numerical Mathematics , volume =

  38. [46]

    Numerical Methods for Partial Differential Equations , volume =

    Ikonen, Samuli and Toivanen, Jari , title =. Numerical Methods for Partial Differential Equations , volume =

  39. [47]

    Computational Methods in Applied Sciences , volume =

    Toivanen, Jari , title =. Computational Methods in Applied Sciences , volume =

  40. [48]

    and Ziogas, Andrew , title =

    Chiarella, Carl and Kang, Boda and Meyer, Gunter H. and Ziogas, Andrew , title =. International Journal of Theoretical and Applied Finance , volume =

  41. [49]

    Flux-Corrected Transport: Principles, Algorithms, and Applications , edition =

  42. [50]

    Bellman, Richard , title =

  43. [51]

    International Journal of Computer Mathematics , Year =

    Efficient solution of structural default models with correlated jumps and mutual obligations , Author =. International Journal of Computer Mathematics , Year =

  44. [52]

    Ikonen and J

    S. Ikonen and J. Toivanen , Journal =. Componentwise splitting methods for pricing. 2007 , Pages =

    2007

  45. [53]

    RAIRO Analyse Num

    Bolley, Catherine and Crouzeix, Michel , title =. RAIRO Analyse Num. 1978 , pages =

    1978

  46. [54]

    and Coughran, William M

    Bank, Randolph E. and Coughran, William M. and Fichtner, Wolfgang and Grosse, Eric H. and Rose, Donald J. and Smith, R. Kent , title =. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems , volume =. 1985 , pages =

    1985

  47. [55]

    Rational

    G\". Rational. GAMM-Mitteilungen , year =

  48. [56]

    Hull , publisher =

    J.C. Hull , publisher =. Options, Futures, and Other Derivatives , year =

  49. [57]

    , title =

    LeVeque, Randall J. , title =

  50. [58]

    Solving Ordinary Differential Equations

    Hairer, Ernst and N. Solving Ordinary Differential Equations

  51. [59]

    Hairer, Ernst and Lubich, Christian and Wanner, Gerhard , title =

  52. [60]

    Hosea, M. E. and Shampine, Lawrence F. , title =. Applied Numerical Mathematics , volume =. 1996 , pages =

    1996

  53. [61]

    2026 , eprint =

    Itkin, Andrey , title =. 2026 , eprint =. doi:10.48550/arXiv.2606.23980 , url =

    2026 doi arXiv

  54. [62]

    K. J. International journal of numerical analysis and modeling , Year =

  55. [63]

    Scharfetter, D. L. and Gummel, H. K. , title =. IEEE Transactions on Electron Devices , volume =

  56. [64]

    Larsen, E. W. and Levermore, C. D. and Pomraning, G. C. and Sanderson, J. G. , title =. Journal of Computational Physics , volume =

  57. [65]

    and Dellacherie, S

    Buet, C. and Dellacherie, S. , title =. Communications in Mathematical Sciences , volume =

  58. [66]

    and Borz

    Mohammadi, M. and Borz. Analysis of the. Journal of Numerical Mathematics , volume =

  59. [67]

    , title =

    Bessemoulin-Chatard, M. , title =. Numerische Mathematik , volume =

  60. [68]

    and Zanella, M

    Pareschi, L. and Zanella, M. , title =. Journal of Scientific Computing , volume =

  61. [69]

    Sweby, P. K. , title =. SIAM Journal on Numerical Analysis , volume =

  62. [70]

    and Turek, S

    Kuzmin, D. and Turek, S. , title =. Journal of Computational Physics , volume =

  63. [71]

    , title =

    Kuzmin, D. , title =. Journal of Computational Physics , volume =

  64. [72]

    Barrenechea, G. R. and John, V. and Knobloch, P. , title =. SIAM Journal on Numerical Analysis , volume =

  65. [73]

    Barrenechea, G. R. and John, V. and Knobloch, P. , title =. SIAM Review , volume =

  66. [74]

    and Shu, C.-W

    Zhang, X. and Shu, C.-W. , title =. Journal of Computational Physics , volume =

  67. [75]

    and Shu, C.-W

    Gottlieb, S. and Shu, C.-W. and Tadmor, E. , title =. SIAM Review , volume =

  68. [76]

    and Meister, A

    Kopecz, S. and Meister, A. , title =. Applied Numerical Mathematics , volume =

  69. [77]

    Daners, D. and Gl. Eventually positive semigroups of linear operators , journal =

  70. [78]

    and Sun, J

    Qi, L. and Sun, J. , title =. Mathematical Programming , volume =

  71. [79]

    and Deleersnijder, E

    Burchard, H. and Deleersnijder, E. and Meister, A. , title =. Applied Numerical Mathematics , volume =

  72. [80]

    The primal-dual active set strategy as a semismooth

    Hinterm. The primal-dual active set strategy as a semismooth. SIAM Journal on Optimization , volume =

  73. [81]

    Journal of Computational Physics , year =

    Zalesak, Steven T , title =. Journal of Computational Physics , year =

  74. [82]

    2026 , note=

    Diagonal Frog meets ADI: trading matrix exponentials for resolvents in the Fokker-Planck equation , author=. 2026 , note=

    2026

  75. [83]

    Applied Numerical Mathematics , volume=

    A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations , author=. Applied Numerical Mathematics , volume=. 2003 , publisher=

    2003

  76. [84]

    SIAM Journal on Optimization , volume=

    The primal-dual active set strategy as a semismooth Newton method , author=. SIAM Journal on Optimization , volume=. 2002 , publisher=

    2002

  77. [85]

    Journal of Computational Physics , volume=

    Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms , author=. Journal of Computational Physics , volume=. 2011 , publisher=

    2011

  78. [86]

    Journal of Computational and Applied Mathematics , volume=

    Linearity-preserving flux correction and convergence acceleration for constrained Galerkin schemes , author=. Journal of Computational and Applied Mathematics , volume=. 2012 , publisher=

    2012

  79. [87]

    Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , author=. Journal of Computational Physics , volume=. 1973 , publisher=

    1973

  80. [88]

    A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics , author=. Mat. Sb. (N.S.) , volume=

  81. [89]

    2026 , eprint=

    Diagonal Frog: High-order positivity-preserving FD schemes for anisotropic Fokker-Planck equations , author=. 2026 , eprint=

    2026

  82. [90]

    Itkin , title =

    A. Itkin , title =. 2020 , publisher =

    2020

  83. [91]

    and Chakravarthy, S

    Osher, S. and Chakravarthy, S. , title =. SIAM Journal on Numerical Analysis , volume =

  84. [92]

    Kuznetsov, N. N. , title =. USSR Computational Mathematics and Mathematical Physics , volume =

  85. [93]

    The optimal convergence rate of monotone finite difference methods for hyperbolic conservation laws , journal =

  86. [94]

    and Shu, C.-W

    Zhang, X. and Shu, C.-W. , title =. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences , volume =

  87. [95]

    A Harnack inequality for parabolic differential equations , journal =

    J. A Harnack inequality for parabolic differential equations , journal =. 1964 , publisher =

    1964

Pith tools

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