Publications
Publications.
Journal Articles
- PinTRKConvergence of ParaOpt for general Runge-Kutta time discretizationsF. Kwok, J. Salomon, and D. TognonTo appear in SIAM Journal on Numerical Analysis, 2026
ParaOpt is a time parallel method based on Parareal for solving optimality systems arising in optimal control problems. The method was presented in [M.J. Gander, F. Kwok and J. Salomon, SIAM J. Sci. Comput., 42 (2020), A2773–A2802] together with a convergence analysis in the case where implicit Euler is used to discretize the differential equations governing the system dynamics. However, its convergence behaviour for higher order time discretizations has not been considered. In this paper, we use an operator norm analysis to prove that the convergence rate of ParaOpt applied to a linear-quadratic optimal control problem has the same order as the Runge-Kutta time integration method used, provided that a few auxiliary order conditions are satisfied. We illustrate our theoretical results with numerical examples, before showing an additional test case not covered by our analysis, namely, a nonlinear optimal control problem involving a Schrödinger type system.
@article{PinTRK, author = {Kwok, F. and Salomon, J. and Tognon, D.}, journal = {To appear in SIAM Journal on Numerical Analysis}, title = {Convergence of ParaOpt for general Runge-Kutta time discretizations}, url = {https://hal.science/hal-05014759}, volume = {accepted, February}, year = {2026} } - PinT-USParaOpt for unstable systemsDjahou N. TognonSubmitted, Feb 2025
We present a convergence analysis of the ParaOpt algorithm applied to linear optimal control problems where the linear operator involved in the constraint – namely, a linear ordinary differential equation (ODE) – has at least one positive eigenvalue, which usually corresponds to unstable systems. When the control problem is discretized using the explicit Euler scheme, the resulting formulation can be expressed as a stationary iterative method, in which a coarse Jacobian approximation acts as a preconditioner. An explicit bound on the spectral radius of the corresponding iteration matrix is proposed, depending on the discretization parameters as well as on the smallest negative and largest positive eigenvalues of the considered linear operator. Our theoretical results are completed by numerical experiments performed on both linear and nonlinear examples.
@article{PinT-US, author = {Tognon, Djahou N.}, journal = {Submitted}, month = feb, title = {ParaOpt for unstable systems}, url = {https://inria.hal.science/hal-04966081}, year = {2025} }
Proceedings
- QuantumQuantum approach for electrical system maintenance scheduling problemQ. Zhou, D. Tognon, S. Chalal, and 4 more authorsIn CEMRACS 2025, Feb 2026
Maintenance scheduling is essential to ensure the reliability of power transmission systems while limiting operational risk and respecting technical and resource constraints. Increasing system complexity due to asset aging and renewable integration challenges classical optimization approaches traditionally used by transmission system operators. This work focuses on the optimization stage of a risk-based maintenance planning framework, where predefined risk indicators guide the construction of feasible schedules. After discussing the limitations of classical mixed-integer programming formulations in terms of scalability and computational effort, we investigate the applicability of quantum optimization methods to the maintenance scheduling problem. In particular, we explore the use of the Quantum Approximate Optimization Algorithm (QAOA) and Grover Adaptive Search to minimize operational risk under realistic constraints. While current quantum hardware limits practical deployment, this study provides insights into the potential of quantum approaches as complementary tools for future large-scale power system maintenance planning.
@inproceedings{Quantum, author = {Zhou, Q. and Tognon, D. and Chalal, S. and Henka, N. and Tazi, S. and Hafid, F. and Titouan, C.}, booktitle = {CEMRACS 2025}, title = {Quantum approach for electrical system maintenance scheduling problem}, url = {https://hal.science/hal-05568055v1}, year = {2026} } - PinTOCPA Parallel in Time Algorithm Based on ParaExp for Optimal Control ProblemsFelix Kwok, and Djahou N. TognonIn 2024 IEEE 63rd Conference on Decision and Control (CDC), Feb 2024
We propose a new parallel-in-time algorithm for solving optimal control problems constrained by discretized partial differential equations. Our approach, which is based on a deeper understanding of ParaExp, considers an overlapping time-domain decomposition in which we combine the solution of homogeneous problems using exponential propagation with the local solutions of inhomogeneous problems. The algorithm yields a linear system whose matrix-vector product can be fully performed in parallel. We then propose a preconditioner to speed up the convergence of GMRES in the special cases of the heat and wave equations. Numerical experiments are provided to illustrate the efficiency of our preconditioners.
@inproceedings{PinTOCP, author = {Kwok, Felix and Tognon, Djahou N.}, booktitle = {2024 IEEE 63rd Conference on Decision and Control (CDC)}, doi = {10.1109/CDC56724.2024.10886085}, pages = {6685-6690}, title = {A Parallel in Time Algorithm Based on ParaExp for Optimal Control Problems}, url = {https://ieeexplore.ieee.org/document/10886085/metrics#metrics}, year = {2024} } - NGSFPA Dynamical Neural Galerkin Scheme for Filtering ProblemsJoubine Aghili, Joy Atokple Zialesi, Marie Billaud-Friess, and 3 more authorsIn ESAIM: Proceedings and Surveys, Feb 2025
This paper considers the filtering problem which consists in reconstructing the state of a dynamical system with partial observations coming from sensor measurements, and the knowledge that the dynamics are governed by a physical PDE model with unknown parameters. We present a filtering algorithm where the reconstruction of the dynamics is done with neural network approximations whose weights are dynamically updated using observational data. In addition to the estimate of the state, we also obtain time-dependent parameter estimations of the PDE parameters governing the observed evolution. We illustrate the behavior of the method in a one-dimensional KdV equation involving the transport of solutions with local support. Our numerical investigation reveals the importance of the location and number of the observations. In particular, it suggests to consider dynamical sensor placement.
@inproceedings{NGSFP, author = {Aghili, Joubine and Atokple Zialesi, Joy and Billaud-Friess, Marie and Garnier, Guillaume and Mula, Olga and Tognon, Djahou N.}, booktitle = {ESAIM: Proceedings and Surveys}, doi = {10.1051/proc/202581002}, pages = {2-15}, title = {A Dynamical Neural Galerkin Scheme for Filtering Problems}, url = {https://www.esaim-proc.org/articles/proc/abs/2024/01/proc240008/proc240008.html}, volume = {81}, year = {2025} }
Theses
- PhD ThesisTime parallelization and machine learning for optimal control and inverse problemsDjahou N. TognonSorbonne Université, Feb 2025
Time-parallel algorithms (PinT) such as Parareal, ParaExp, etc., are well known in the literature for their ability to exploit the parallel architecture of today’s computers to solve initial value porblem. Recently, ParaOpt, an algorithm based on the principles of the Parareal algorithm, has made it possible to parallelize the solution of optimality systems. A first analysis of the convergence of this algorithm was presented in [M.J. Gander, F. Kwok and J. Salomon, SIAM J. Sci. Comput.,42 (2020), A2773-A2802] in the restricted case of the implicit Euler method for linear quadratic optimal control problems (LQOCP) involving dissipative systems. In this thesis, we first present a convergence analysis when the system under study is unstable. Secondly, we focus on the influence of the solver used for time resolution on convergence, considering the more general case of LQOCP discretization by Runge-Kutta methods. We show that the convergence rate of ParaOpt has the same order as the Runge-Kutta time integration method used, provided that the Runge-Kutta method satisfies some additional order conditions. We then consider a preconditioning problem, for which we introduce a new PinT algorithm, this time based on the ParaExp algorithm, for LQOCP solving. Our approach is based on an overlapping time-interval decomposition in which we combine the solution of homogeneous subproblems of optimality systems using exponential propagation with local solutions of inhomogeneous subproblems. The formulation leads to a linear system whose matrix-vector product can be fully computed in parallel. We then propose two preconditioners to accelerate the convergence of GMRES in the special cases of heat and wave equations. In a final independent chapter, we analyze a learning-based model correction method. The approach studied follows the principles of Aphynity [Yuan Yin et al J. Stat. Mech. (2021) 124012], a method consisting of introducing a corrective term in the form of a neural network into the PDE under consideration. Learning is then performed in the outer loop of the solver used, so that training takes place indirectly through a time scheme. We study the influence of the solver on the resulting network, and show in particular that the order of approximation of the corrected model is equal to the order of the scheme used. Finally, we show how a Richardson-type acceleration strategy can speed up convergence by considering two smaller, independently trained networks with different time steps.
@phdthesis{PhD Thesis, author = {Tognon, Djahou N.}, school = {Sorbonne Université}, title = {Time parallelization and machine learning for optimal control and inverse problems}, url = {https://theses.hal.science/tel-05262960}, year = {2025} }