Skip to content
DraconiticResearch
Search Sources API
197,787 studies31 sourcesUpdated 45 min ago
PreprintarXiv licencearXiv

Exploiting Micro-Structure in the Generalized Riccati Recursion for Optimal Control

Abstract

Optimal Control Problems (OCPs) have been applied to a variety of applications. Such problems are often transcribed into a Nonlinear Program (NLP) for which solution methods typically solve a Linear Quadratic (LQ) subproblem. To efficiently solve this LQ subproblem, the Riccati recursion exploits the block-diagonal structure that arises from discretizing an OCP using a Multiple Shooting method. However, reformulating general dynamics constraints to fit the required problem structure to use the Riccati recursion introduces additional micro-structure within these block matrices that remains unexploited. This work demonstrates how structure-exploiting solvers such as Fatrop can benefit from exploiting not only the block-diagonal structure but also the micro-structure within those blocks. We propose a structure-exploiting LU decomposition algorithm and integrate it within the Riccati recursion. Numerical results on randomized linear systems and two OCP examples demonstrate a reduction in computation time of the search direction of up to 20%, with the largest speedup observed for relative zero block sizes around 0.25 and large stage-wise equality constraint jacobians.

Related

  1. Near-Optimal Quantum Algorithm and Complexity Analysis for Riccati Problems
    Jingyao Wang, Yanqiao Wang, Bowen Li +2
    arXivPreprintarXiv licenceFull textMathematical physicsQuantum physics & quantum informationarXiv
  2. Sub-Infinite Horizon Stochastic Linear-Quadratic Optimal Control Problems and Delayed Backward Riccati Equations
    Yutao Chen, Hongwei Lou, Hanxiao Wang
    arXivPreprintCC BY 4.0Full textOptimisation & controlarXiv
  3. RapidDOC: Riccati-Preconditioned Implicit Differentiation for Differentiable Optimal Control under Operator Mismatch
    Xiangda Zhu, Yuankun Chen, Zifei Nie +1
    arXivPreprintarXiv licenceFull textOptimisation & controlarXiv
  4. Exploiting Mutual Coupling Structure for Channel Estimation of Active RIS-Assisted Links
    Simon Tarboush, Pinjun Zheng, Hadi Sarieddeen +2
    arXivPreprintarXiv licenceFull textGeneral mathematicsarXiv
  5. Optimal Two‐Phase Sampling Designs for Generalized Raking Estimators With Multiple Parameters of Interest
    Jasper B. Yang, Bryan E. Shepherd, Thomas Lumley +1
    Statistics in medicinePeer-reviewedCC BYFull textStatisticsEurope PMC
BibTeX
@misc{callens2026exploiting,
  title = {{Exploiting Micro-Structure in the Generalized Riccati Recursion for Optimal Control}},
  author = {Louis Callens and Wilm Decré},
  howpublished = {arXiv},
  year = {2026},
  month = aug,
  doi = {10.48550/arxiv.2608.29238},
  eprint = {2608.29238},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2608.29238},
}
RIS
TY  - UNPB
TI  - Exploiting Micro-Structure in the Generalized Riccati Recursion for Optimal Control
AU  - Louis Callens
AU  - Wilm Decré
PY  - 2026
DA  - 2026/08/29
PB  - arXiv
DO  - 10.48550/arxiv.2608.29238
UR  - https://arxiv.org/abs/2608.29238
AB  - Optimal Control Problems (OCPs) have been applied to a variety of applications. Such problems are often transcribed into a Nonlinear Program (NLP) for which solution methods typically solve a Linear Quadratic (LQ) subproblem. To efficiently solve this LQ subproblem, the Riccati recursion exploits the block-diagonal structure that arises from discretizing an OCP using a Multiple Shooting method. However, reformulating general dynamics constraints to fit the required problem structure to use the Riccati recursion introduces additional micro-structure within these block matrices that remains unexploited. This work demonstrates how structure-exploiting solvers such as Fatrop can benefit from exploiting not only the block-diagonal structure but also the micro-structure within those blocks. We propose a structure-exploiting LU decomposition algorithm and integrate it within the Riccati recursion. Numerical results on randomized linear systems and two OCP examples demonstrate a reduction in computation time of the search direction of up to 20%, with the largest speedup observed for relative zero block sizes around 0.25 and large stage-wise equality constraint jacobians.
ER  -