pith. sign in

arxiv: 2604.25279 · v1 · submitted 2026-04-28 · 🧮 math.OC

Sakawa-Shindo algorithm for optimal control of time-delay systems, with applications to epidemiology

Rami Katz , Francesca Cal\`a Campana , Giulia Giordano This is my paper

Pith reviewed 2026-05-07 15:56 UTC · model grok-4.3

classification 🧮 math.OC
keywords Sakawa-Shindo algorithmoptimal controltime-delay systemsdelay differential equationsepidemic modelsnon-pharmaceutical interventionsvaccination strategies
0
0 comments X

The pith

The Sakawa-Shindo algorithm extends to optimal control problems with any number of discrete state delays and converges to first-order optimal solutions.

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

The paper adapts the Sakawa-Shindo algorithm to optimal control settings where the system state evolves according to delay differential equations with an arbitrary number of fixed delays. It proves that the adapted method terminates after finitely many steps, that the sequence of controls it produces becomes asymptotically first-order optimal, and that a subsequence converges to a control satisfying the first-order necessary conditions. The authors apply the method to epidemic models that incorporate delays for incubation periods and the time lag before vaccination becomes effective, using it to optimize the timing and intensity of non-pharmaceutical interventions and vaccination schedules.

Core claim

We extend the Sakawa-Shindo algorithm to solve optimal control problems where the system dynamics involve an arbitrary number of discrete state delays. We prove that the algorithm guarantees termination in a finite number of steps, asymptotic first-order optimality of the generated control sequence and convergence of a subsequence to a control satisfying first-order optimality, and we apply it to the optimal design of non-pharmaceutical interventions and vaccination plans for epidemic models with delays associated with incubation period and vaccination.

What carries the argument

The extended Sakawa-Shindo algorithm, which iteratively updates candidate controls using adjoint equations modified to handle multiple discrete state delays.

If this is right

  • The algorithm terminates after a finite number of iterations for any finite number of discrete delays.
  • The generated control sequence achieves asymptotic first-order optimality.
  • A subsequence of controls converges to one that satisfies the first-order necessary optimality conditions.
  • The method can be used to compute optimal non-pharmaceutical intervention schedules and vaccination plans in epidemic models that include incubation and vaccination delays.

Where Pith is reading between the lines

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

  • The same delay-handling technique could be tested on models that replace discrete delays with distributed delays to check whether convergence properties survive.
  • Optimal controls produced by the algorithm could be compared against real-world timing data from past outbreaks to assess whether they improve predicted outcomes.
  • The approach might transfer to other biological control problems, such as optimizing treatment schedules in models of chronic disease progression with maturation delays.

Load-bearing premise

The underlying delay differential equations and cost functionals must satisfy regularity conditions such as Lipschitz continuity and differentiability.

What would settle it

Numerical execution of the algorithm on a concrete epidemic model with known incubation and vaccination delays that either fails to terminate in finite steps or produces a final control violating the first-order optimality condition to machine precision.

Figures

Figures reproduced from arXiv: 2604.25279 by Francesca Cal\`a Campana, Giulia Giordano, Rami Katz.

Figure 1
Figure 1. Figure 1: Leftmost two panels: Extended SIRV OCP with delays h1 = 5 and h2 = 7, control bounds umax = 0.4 and vmax = 0.8, parameters T = 350 days, Λ = 2.91 · 10−5, µi = 2.90 · 10−5 ∀i, β = 1 and γ = 1/6, σV = σR = 1/500, θV = 0.0013, θR = 0.0021, initial conditions I0 = 10−6, S0 = 1 − I0, R0 = 0, V0 = 0, and cost functional weights wI = 104, wu = 1, wv = 10. Rightmost two panels: Extended SIDARTHE-V OCP with delays … view at source ↗
read the original abstract

We extend the Sakawa-Shindo algorithm to solve optimal control problems where the system dynamics involve an arbitrary number of discrete state delays. We prove that the algorithm guarantees termination in a finite number of steps, asymptotic first-order optimality of the generated control sequence and convergence of a subsequence to a control satisfying first-order optimality, and we apply it to the optimal design of non-pharmaceutical interventions and vaccination plans for epidemic models with delays associated with incubation period and vaccination.

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 / 2 minor

Summary. The manuscript extends the Sakawa-Shindo algorithm to optimal control problems governed by delay differential equations with an arbitrary number of discrete state delays. It proves finite termination of the algorithm, asymptotic first-order optimality of the generated control sequence, and convergence of a subsequence to a first-order stationary control. The method is applied to the optimal design of non-pharmaceutical interventions and vaccination plans in epidemic models that incorporate delays associated with incubation periods and vaccination.

Significance. If the stated regularity conditions hold, the work supplies a convergent, implementable algorithm for a practically important class of delayed optimal control problems. The finite-termination argument, the embedding of the multi-delay adjoint system, and the concrete epidemiological application constitute a useful contribution to the literature on numerical methods for DDE control.

minor comments (2)
  1. [Abstract] Abstract: the phrase 'asymptotic first-order optimality of the generated control sequence' would benefit from a parenthetical clarification of the topology (e.g., L^2 or L^∞) in which the limit is taken.
  2. [Adjoint derivation] §3 (or wherever the adjoint equations appear): the notation for the delayed adjoint variables could be made more uniform across the multiple-delay case to ease verification of the first-order condition.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our work, the accurate summary of the contributions, and the recommendation for minor revision. We are pleased that the extension of the Sakawa-Shindo algorithm to multi-delay systems, the finite-termination and convergence results, and the epidemiological application are viewed as a useful addition to the literature.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper extends the Sakawa-Shindo algorithm to DDEs with an arbitrary number of discrete delays, derives the corresponding adjoint system, and proves finite termination plus convergence to first-order stationarity by showing that each successful step yields a strict cost decrease bounded away from zero. These arguments follow the standard structure of the original algorithm once the delayed state and adjoint are embedded; they do not reduce to the inputs by definition, rely on fitted parameters renamed as predictions, or depend on load-bearing self-citations. The epidemiological application is a straightforward instantiation of the extended method under the stated regularity assumptions. The derivation chain is therefore self-contained and independent.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Review based on abstract only; no explicit free parameters, axioms, or invented entities are stated. Standard regularity assumptions for delay systems are implicitly required for the algorithm to be well-defined.

pith-pipeline@v0.9.0 · 5368 in / 1062 out tokens · 38587 ms · 2026-05-07T15:56:42.244343+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

20 extracted references · 20 canonical work pages

  1. [1]

    V . B. Kolmanovskii and L. E. Shaikhet,Control of systems with aftereffect. American Mathematical Soc., 1996, vol. 157

    work page 1996

  2. [2]

    The maximum principle for optimal control problems with time delays,

    A. Boccia and R. B. Vinter, “The maximum principle for optimal control problems with time delays,”SIAM J. Control Optim., vol. 55, no. 5, pp. 2905–2935, 2017

    work page 2017

  3. [3]

    Numerical optimal control for delay differential equations: A simultaneous approach based on linearization of the delayed state,

    T. K. S. Ritschel and S. Stange, “Numerical optimal control for delay differential equations: A simultaneous approach based on linearization of the delayed state,” inEuropean Control Conf., 2025, pp. 777–782

    work page 2025

  4. [4]

    Optimal control problems with delays in state and control variables subject to mixed control–state constraints,

    L. Göllmann, D. Kern, and H. Maurer, “Optimal control problems with delays in state and control variables subject to mixed control–state constraints,”Optim. Control Appl. Meth., vol. 30, no. 4, pp. 341–365, 2008

    work page 2008

  5. [5]

    On global convergence of an algorithm for optimal control,

    Y . Sakawa and Y . Shindo, “On global convergence of an algorithm for optimal control,”IEEE Trans. Autom. Control, vol. 25, no. 6, pp. 1149– 1153, 1980

    work page 1980

  6. [6]

    A contribution to the mathematical theory of epidemics,

    W. Kermack and A. McKendrick, “A contribution to the mathematical theory of epidemics,”Proc. Royal Society London, vol. 115:772, 1927

    work page 1927

  7. [7]

    Mathematical models in population biology and epidemiology,

    F. Brauer and C. Castillo-Chavez, “Mathematical models in population biology and epidemiology,”Springer, vol. 2nd ed., 2012

    work page 2012

  8. [8]

    Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy,

    G. Giordano, F. Blanchini, R. Bruno, P. Colaneri, A. Di Filippo, A. Di Matteo, and M. Colaneri, “Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy,”Nature Medicine, vol. 26, pp. 855–860, 2020

    work page 2020

  9. [9]

    Modeling vaccination rollouts, SARS-CoV-2 variants and the requirement for non- pharmaceutical interventions in Italy,

    G. Giordano, M. Colaneri, A. Di Filippo, F. Blanchini, P. Bolzern, G. De Nicolao, P. Sacchi, P. Colaneri, and R. Bruno, “Modeling vaccination rollouts, SARS-CoV-2 variants and the requirement for non- pharmaceutical interventions in Italy,”Nature Medicine, 2021

    work page 2021

  10. [10]

    Modelling and control of epidemics across scales,

    E. Hernandez-Vargas, A. González, C. Beck, X. Bi, F. Calà Campana, and G. Giordano, “Modelling and control of epidemics across scales,” Proc. IEEE Conf. Dec. Control (CDC), pp. 4963–4980, 2022

    work page 2022

  11. [11]

    Optimal control applied to biological models,

    S. Lenhart and J. T. Workman, “Optimal control applied to biological models,”Chapman and Hall, CRC, 2007

    work page 2007

  12. [12]

    Optimal control in epidemiology,

    O. Sharomi and T. Malik, “Optimal control in epidemiology,”Ann. Oper . Res., vol. 251, pp. 55–71, 2017

    work page 2017

  13. [13]

    Optimal time-profiles of public health intervention to shape voluntary vaccination for childhood diseases,

    B. Buonomo, P. Manfredi, and A. d’Onofrio, “Optimal time-profiles of public health intervention to shape voluntary vaccination for childhood diseases,”J. Math. Biol., vol. 78, p. 1089–1113, 2019

    work page 2019

  14. [14]

    Sequential-quadratic- hamiltonian optimal control of epidemic models with an arbitrary number of infected and non-infected compartments,

    F. Calà Campana, R. Katz, and G. Giordano, “Sequential-quadratic- hamiltonian optimal control of epidemic models with an arbitrary number of infected and non-infected compartments,”IEEE Control Systems Letters, vol. 8, pp. 1805–1810, 2024

    work page 2024

  15. [15]

    F. A. Rihan,Delay Differential Equations and Applications to Biology, 2nd ed. Cham: Springer, 2026

    work page 2026

  16. [16]

    Optimal control of an epidemiological model with multiple time delays,

    E. B. Bashier and K. C. Patidar, “Optimal control of an epidemiological model with multiple time delays,”Appl Math Comput, vol. 292, pp. 47–56, 2017

    work page 2017

  17. [17]

    Application of a time-delay SIR model with vaccination in COVID-19 prediction and its optimal control strategy,

    S. Dong, L. Xu, Y . A, Z.-Z. Lan, D. Xiao, and B. Gao, “Application of a time-delay SIR model with vaccination in COVID-19 prediction and its optimal control strategy,”Nonlinear Dynamics, vol. 111, pp. 10 677– 10 692, 2023

    work page 2023

  18. [18]

    J. K. Hale,Ordinary differential equations. Courier Corporation, 2009

    work page 2009

  19. [19]

    R. E. Megginson,An introduction to Banach space theory. Springer Science & Business Media, 2012, vol. 183

    work page 2012

  20. [20]

    Sakawa-Shindo algorithm for optimal control of time-delay systems, with applications to epidemiology,

    R. Katz, F. Calà Campana, and G. Giordano, “Sakawa-Shindo algorithm for optimal control of time-delay systems, with applications to epidemiology,” 2026. [Online]. Available: https://giuliagiordano.dii. unitn.it/docs/Katz2026.pdf

    work page 2026