pith. machine review for the scientific record. sign in

arxiv: 2604.07644 · v1 · submitted 2026-04-08 · 💻 cs.RO · cs.AI· cs.SY· eess.SY· math.OC

Recognition: unknown

Safe Large-Scale Robust Nonlinear MPC in Milliseconds via Reachability-Constrained System Level Synthesis on the GPU

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:56 UTC · model grok-4.3

classification 💻 cs.RO cs.AIcs.SYeess.SYmath.OC
keywords robust nonlinear MPCsystem level synthesisGPU accelerationreachability analysissafe robotic controlreal-time optimizationADMM quadratic programming
0
0 comments X

The pith

GPU-parallelized system level synthesis computes safe robust nonlinear MPC for 75D humanoids in 20 milliseconds on average.

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

The paper develops a GPU framework that simultaneously optimizes a feasible nominal trajectory, a tracking controller, and an over-approximated reachable set under bounded disturbances. This joint computation is performed via a custom ADMM-based quadratic program solver that exploits parallel associative scans and adaptive caching. The approach scales to problems with hundreds of thousands of decision variables while delivering 100 percent empirical safety on whole-body quadruped and humanoid models. A sympathetic reader would care because it turns high-dimensional robust control from an offline or slow process into a real-time capability for uncertain robotic systems.

Core claim

By embedding reachability constraints directly into system level synthesis and executing the resulting large-scale quadratic programs on the GPU, the method produces dynamically feasible nominal trajectories, robust tracking controllers, and closed-loop reachable sets in a single online optimization loop. This yields average solve times of 20 milliseconds for 61- to 75-dimensional nonlinear systems while satisfying up to 80,000 constraints and achieving 100 percent empirical safety under disturbance.

What carries the argument

Reachability-constrained system level synthesis (SLS) executed through a GPU-accelerated ADMM quadratic program solver that uses parallel associative scans and adaptive caching to optimize closed-loop reachable sets alongside nominal trajectories and controllers.

If this is right

  • Real-time robust control becomes feasible for whole-body humanoid and quadruped robots with 60-75 state dimensions.
  • Solve times for nominal trajectory optimization drop by 97.7 percent relative to state-of-the-art CPU solvers.
  • SLS-based robust controller and reachable-set synthesis accelerate by a factor of 237 compared with prior GPU implementations.
  • The method maintains 100 percent empirical safety on high-dimensional nonlinear systems where learning-based reachability baselines fail.
  • Problem sizes up to 200,000 decision variables and 80,000 constraints can be handled in fixed- or receding-horizon settings.

Where Pith is reading between the lines

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

  • The same GPU backend could be reused for other large-scale robust optimization problems that admit SLS formulations.
  • Combining the reachability constraints with learned disturbance models might further reduce conservatism without losing formal guarantees.
  • Longer planning horizons may become practical in real-time if the current scaling continues with future GPU hardware.

Load-bearing premise

The system dynamics must be known to sufficient accuracy and disturbances must lie in compact sets whose images under the closed-loop maps remain tractable to over-approximate.

What would settle it

Run the synthesized controller on a physical or simulated robot subjected to disturbances drawn from a distribution larger than the assumed compact sets and observe whether safety violations occur despite the computed reachable sets remaining invariant.

Figures

Figures reproduced from arXiv: 2604.07644 by Glen Chou, Jeffrey Fang.

Figure 1
Figure 1. Figure 1: (a): GPU-accelerated nonlinear constrained whole-body control (61 states, 12 controls) executed on hardware with a Unitree Go2 EDU quadruped, navigating through an obstacle field in real-time at 50 Hz. (b): Overhead view of simulated positions reconstructed from hardware-experiment data showing successful navigation around the obstacles. (c): Graph illustrating minimum distance to any obstacle, demonstrati… view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of our method, GPU-SLS. At each SQP iteration (s), we compute a nominal trajectory update (δx, δu) using a GPU-parallelized ADMM-based QP solver that exploits caching and parallel associative scans (PASs) for acceleration. The resulting dual variables are post-processed into SLS￾compatible duals τ (App. B), which inform the objective in the SLS controller optimization problem (18a). This optimiza… view at source ↗
Figure 3
Figure 3. Figure 3: (a): Comparison of average solve-time scaling with horizon length for various NMPC solvers on a torque-constrained 10-link pendulum stabilization task under random external disturbances. Our method consistently outperforms all baselines for long-horizon problems. (b): Solve-time scaling of GPU-SLS with a family of n-link torque-constrained pendulums. constraints. We warm-start the controller update using t… view at source ↗
Figure 1
Figure 1. Figure 1: In RTI, only a single SQP iteration is performed per [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 6
Figure 6. Figure 6: Dubins car experiment using GPU-SLS to calculate a robust [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: Runtime comparison of FastSLS vs GPU-SLS for solving the [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 10
Figure 10. Figure 10: Runtime comparisons between our GPU solver and the CPU-based [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 8
Figure 8. Figure 8: (a): Whole-body humanoid (75D, 19C) navigating around an obstacle using (14) calculated from GPU-SLS. (b): A visualization of the robust forward tubes (blue squares) generated by (18) and the humanoid’s rolled out trajectory remaining within the tubes. 0.0 0.3 0.6 0.9 p x [m] −0.75 −0.50 −0.25 0.00 p y [m] 0.8 0.9 1.0 1.1 p z [m] 0 1 Time [s] 0.0 0.4 0.8 1.2 v x [m / s] 0 1 Time [s] −0.8 −0.4 0.0 v y [m / … view at source ↗
Figure 11
Figure 11. Figure 11: (a): GPU-accelerated robust whole-body control (61 states, 12 controls) using RTI GPU-SLS executed on hardware with a Unitree Go2 EDU quadruped. The robot successfully navigates robustly around an obstacle in real-time at 50 Hz. (b): Graph illustrating the robot’s distance to the obstacle and the maximum tube size for each time step along the trajectory. The tubes are used to tighten the constraints in th… view at source ↗
Figure 12
Figure 12. Figure 12: Upsweep phase of the parallel associative scan. Leaf nodes [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Downsweep phase of the parallel associative scan. Using the subtree [PITH_FULL_IMAGE:figures/full_fig_p013_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Whole-body humanoid (75D, 19C) robustly navigating through two [PITH_FULL_IMAGE:figures/full_fig_p015_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Visualization of (px, py)-slices of the robust reachable tubes at certain time steps for the whole-body humanoid rollout in [PITH_FULL_IMAGE:figures/full_fig_p015_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Simulated rollouts of Fig. 8, showing disturbed states and robust tubes under adversarial disturbances. [PITH_FULL_IMAGE:figures/full_fig_p017_16.png] view at source ↗
read the original abstract

We present GPU-SLS, a GPU-parallelized framework for safe, robust nonlinear model predictive control (MPC) that scales to high-dimensional uncertain robotic systems and long planning horizons. Our method jointly optimizes an inequality-constrained, dynamically-feasible nominal trajectory, a tracking controller, and a closed-loop reachable set under disturbance, all in real-time. To efficiently compute nominal trajectories, we develop a sequential quadratic programming procedure with a novel GPU-accelerated quadratic program (QP) solver that uses parallel associative scans and adaptive caching within an alternating direction method of multipliers (ADMM) framework. The same GPU QP backend is used to optimize robust tracking controllers and closed-loop reachable sets via system level synthesis (SLS), enabling reachability-constrained control in both fixed- and receding-horizon settings. We achieve substantial performance gains, reducing nominal trajectory solve times by 97.7% relative to state-of-the-art CPU solvers and 71.8% compared to GPU solvers, while accelerating SLS-based control and reachability by 237x. Despite large problem scales, our method achieves 100% empirical safety, unlike high-dimensional learning-based reachability baselines. We validate our approach on complex nonlinear systems, including whole-body quadrupeds (61D) and humanoids (75D), synthesizing robust control policies online on the GPU in 20 milliseconds on average and scaling to problems with 2 x 10^5 decision variables and 8 x 10^4 constraints. The implementation of our method is available at https://github.com/Jeff300fang/gpu_sls.

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

2 major / 3 minor

Summary. The manuscript presents GPU-SLS, a GPU-parallelized framework for safe robust nonlinear MPC that jointly optimizes an inequality-constrained nominal trajectory via ADMM-based SQP, a tracking controller, and closed-loop reachable-set over-approximations via system level synthesis (SLS). It reports 97.7% reduction in nominal solve time versus CPU solvers, 71.8% versus other GPU solvers, 237x acceleration of SLS components, average 20 ms synthesis time, scaling to 2×10^5 decision variables and 8×10^4 constraints, and 100% empirical safety on 61D quadruped and 75D humanoid systems, with open-source code provided.

Significance. If the central claims hold, the work enables real-time formally constrained robust control for high-dimensional robotic systems that were previously intractable, bridging nonlinear MPC with SLS reachability in a scalable GPU implementation. The open-source release at the cited GitHub repository is a clear strength that supports reproducibility and further development.

major comments (2)
  1. [Experiments] Experiments section: the 100% empirical safety result is central to the safety claim, yet the manuscript provides no breakdown of the number of Monte Carlo trials, the specific disturbance sets used for reachable-set validation, or quantitative comparison against the over-approximation tightness; without these, it is difficult to judge whether the result generalizes beyond the reported scenarios.
  2. [GPU-accelerated QP solver] Section on the GPU ADMM QP solver: the adaptive caching and parallel associative scan technique is load-bearing for the reported speedups, but the manuscript does not include a convergence analysis or iteration-count comparison under the same termination tolerances as the CPU baselines; this leaves open whether the 97.7% speedup preserves solution quality.
minor comments (3)
  1. [Abstract] Abstract: the phrase 'accelerating SLS-based control and reachability by 237x' should specify the exact baseline implementation and hardware to avoid ambiguity in the speedup claim.
  2. [Preliminaries] Notation: the distinction between open-loop and closed-loop reachable sets is introduced without an explicit equation reference in the early sections; adding a forward reference to the SLS parameterization would improve readability.
  3. [Figures] Figure captions: several timing plots lack error bars or mention of the number of repeated runs, which would help readers assess variability of the millisecond-scale results.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments and positive overall assessment of the work. We address each major comment below and will revise the manuscript to incorporate the suggested clarifications.

read point-by-point responses
  1. Referee: Experiments section: the 100% empirical safety result is central to the safety claim, yet the manuscript provides no breakdown of the number of Monte Carlo trials, the specific disturbance sets used for reachable-set validation, or quantitative comparison against the over-approximation tightness; without these, it is difficult to judge whether the result generalizes beyond the reported scenarios.

    Authors: We agree that additional details on the Monte Carlo validation procedure will strengthen the safety claims. In the revised manuscript we will expand the Experiments section to report the number of trials performed, the precise disturbance sets used (including their bounds and structure), and quantitative metrics comparing the closed-loop reachable-set over-approximations to the empirical sets (e.g., volume ratios or maximum observed deviations). These additions will make the generalization of the 100% safety result clearer. revision: yes

  2. Referee: Section on the GPU ADMM QP solver: the adaptive caching and parallel associative scan technique is load-bearing for the reported speedups, but the manuscript does not include a convergence analysis or iteration-count comparison under the same termination tolerances as the CPU baselines; this leaves open whether the 97.7% speedup preserves solution quality.

    Authors: We acknowledge that an explicit convergence comparison would be helpful. Although the GPU ADMM solver is terminated at identical residual tolerances to the CPU baselines, we will add to the revised manuscript a table and brief analysis of iteration counts across representative problem sizes together with verification that final objective values and constraint violations match those of the CPU solvers to within numerical tolerance. This will confirm that the reported speedups preserve solution quality. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper presents an algorithmic construction for GPU-parallelized SLS-based robust nonlinear MPC, including a novel ADMM QP solver with associative scans and caching, jointly optimizing nominal trajectories, tracking controllers, and closed-loop reachable sets. No equations or claims reduce by construction to fitted parameters, self-definitions, or renamed known results. Self-citations to SLS literature are standard external references and not load-bearing for the core GPU scaling and empirical safety claims, which rest on benchmarks (97.7% faster solves, 237x acceleration, 100% safety on 61D/75D systems) rather than internal circular steps. The derivation chain is self-contained against external validation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard robust-control assumptions (bounded disturbances, accurate nominal model) rather than new free parameters or invented entities; no ad-hoc constants are mentioned.

axioms (1)
  • domain assumption Disturbances are bounded and the closed-loop maps admit tractable reachable-set over-approximations via SLS
    Invoked to guarantee safety while optimizing reachable sets

pith-pipeline@v0.9.0 · 5599 in / 1216 out tokens · 53332 ms · 2026-05-10T16:56:20.299818+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Local Linearity of LLMs Enables Activation Steering via Model-Based Linear Optimal Control

    cs.LG 2026-04 conditional novelty 7.0

    Local linearity of LLM layers enables LQR-based closed-loop activation steering with theoretical tracking guarantees.

  2. Over-Approximating Minimizer Sets of Constrained Convex Programs with Parametric Uncertainty via Reachability Analysis

    math.OC 2026-04 unverdicted novelty 6.0

    A reachability-analysis method on projected gradient descent dynamics produces certified outer approximations to the minimizer sets of strongly convex programs whose costs depend on bounded uncertain parameters.

  3. VISION-SLS: Safe Perception-Based Control from Learned Visual Representations via System Level Synthesis

    cs.RO 2026-04 conditional novelty 6.0

    VISION-SLS learns visual features with state-dependent error bounds and optimizes causal affine output-feedback policies via system level synthesis to achieve safe nonlinear control from RGB images.

Reference graph

Works this paper leans on

57 extracted references · 14 canonical work pages · cited by 3 Pith papers · 1 internal anchor

  1. [1]

    In: 2025 IEEE 64th Conference on Decision and Control (CDC)

    Arshiya Taj Abdul, Augustinos D. Saravanos, and Evan- gelos A. Theodorou. Nonlinear robust optimization for planning and control. In2025 IEEE 64th Conference on Decision and Control (CDC), pages 3383–3390, 2025. doi: 10.1109/CDC57313.2025.11312909

  2. [2]

    Mpcgpu: Real-time nonlinear model predictive control through preconditioned conjugate gradient on the gpu

    Emre Adabag, Miloni Atal, William Gerard, and Brian Plancher. Mpcgpu: Real-time nonlinear model predictive control through preconditioned conjugate gradient on the gpu. In2024 IEEE International Conference on Robotics and Automation (ICRA), pages 9787–9794. IEEE, 2024

  3. [3]

    Differentiable model predictive control on the gpu

    Emre Adabag, Marcus Greiff, John Subosits, and Thomas Lew. Differentiable model predictive control on the gpu. arXiv preprint arXiv:2510.06179, 2025

  4. [4]

    Set propagation techniques for reachability analysis.Annual Review of Control, Robotics, and Autonomous Systems, 4(1):369–395, 2021

    Matthias Althoff, Goran Frehse, and Antoine Girard. Set propagation techniques for reachability analysis.Annual Review of Control, Robotics, and Autonomous Systems, 4(1):369–395, 2021

  5. [5]

    Primal-dual ilqr for gpu-accelerated learning and control in legged robots.arXiv preprint arXiv:2506.07823, 2025

    Lorenzo Amatucci, Jo ˜ao Sousa-Pinto, Giulio Turrisi, Dominique Orban, Victor Barasuol, and Claudio Semini. Primal-dual ilqr for gpu-accelerated learning and control in legged robots.arXiv preprint arXiv:2506.07823, 2025

  6. [6]

    Control barrier functions: Theory and applications

    Aaron D Ames, Samuel Coogan, Magnus Egerstedt, Gen- naro Notomista, Koushil Sreenath, and Paulo Tabuada. Control barrier functions: Theory and applications. In 2019 18th European control conference (ECC), pages 3420–3431. Ieee, 2019

  7. [7]

    Doyle, Steven H

    James Anderson, John C. Doyle, Steven H. Low, and Nikolai Matni. System level synthesis.Annu. Rev. Control., 47:364–393, 2019

  8. [8]

    System level synthesis.Annual Reviews in Control, 47:364–393, 2019

    James Anderson, John C Doyle, Steven H Low, and Nikolai Matni. System level synthesis.Annual Reviews in Control, 47:364–393, 2019

  9. [9]

    Deepreach: A deep learning approach to high-dimensional reachability

    Somil Bansal and Claire J Tomlin. Deepreach: A deep learning approach to high-dimensional reachability. In 2021 IEEE International Conference on Robotics and Automation (ICRA), pages 1817–1824. IEEE, 2021

  10. [10]

    Hamilton-jacobi reachability: A brief overview and recent advances

    Somil Bansal, Mo Chen, Sylvia Herbert, and Claire J Tomlin. Hamilton-jacobi reachability: A brief overview and recent advances. In2017 IEEE 56th Annual Confer- ence on Decision and Control (CDC), pages 2242–2253. IEEE, 2017

  11. [11]

    Stochastic mpc with online-optimized policies and closed-loop guaran- tees.arXiv preprint arXiv:2502.06469, 2025

    Marcell Bartos, Alexandre Didier, Jerome Sieber, Jo- hannes K ¨ohler, and Melanie N Zeilinger. Stochastic mpc with online-optimized policies and closed-loop guaran- tees.arXiv preprint arXiv:2502.06469, 2025

  12. [12]

    Relu- qp: A gpu-accelerated quadratic programming solver for model-predictive control

    Arun L Bishop, John Z Zhang, Swaminathan Guru- murthy, Kevin Tracy, and Zachary Manchester. Relu- qp: A gpu-accelerated quadratic programming solver for model-predictive control. In2024 IEEE International Conference on Robotics and Automation (ICRA), pages 13285–13292. IEEE, 2024

  13. [13]

    Scans as primitive parallel operations

    Guy E Blelloch. Scans as primitive parallel operations. IEEE Transactions on computers, 38(11):1526–1538, 2002

  14. [14]

    Distributed optimization and statistical learning via the alternating direction method of multipliers.Foundations and Trends® in Machine learning, 3(1):1–122, 2011

    Stephen Boyd, Neal Parikh, Eric Chu, Borja Peleato, Jonathan Eckstein, et al. Distributed optimization and statistical learning via the alternating direction method of multipliers.Foundations and Trends® in Machine learning, 3(1):1–122, 2011

  15. [15]

    Robust model predictive control with polytopic model uncertainty through system level syn- thesis.Automatica, 162:111431, 2024

    Shaoru Chen, Victor M Preciado, Manfred Morari, and Nikolai Matni. Robust model predictive control with polytopic model uncertainty through system level syn- thesis.Automatica, 162:111431, 2024

  16. [16]

    Model error propagation via learned contraction metrics for safe feedback motion planning of unknown systems

    Glen Chou, Necmiye Ozay, and Dmitry Berenson. Model error propagation via learned contraction metrics for safe feedback motion planning of unknown systems. In2021 60th IEEE Conference on Decision and Control (CDC), pages 3576–3583. IEEE, 2021

  17. [17]

    Safe output feedback motion planning from images via learned perception modules and contraction theory

    Glen Chou, Necmiye Ozay, and Dmitry Berenson. Safe output feedback motion planning from images via learned perception modules and contraction theory. In International Workshop on the Algorithmic Foundations of Robotics, pages 349–367. Springer, 2022

  18. [18]

    Goal-reaching trajectory design near danger with piecewise affine reach-avoid computation

    Long Kiu Chung, Wonsuhk Jung, Chuizheng Kong, and Shreyas Kousik. Goal-reaching trajectory design near danger with piecewise affine reach-avoid computation. arXiv preprint arXiv:2402.15604, 2024

  19. [19]

    Fast parallel matrix inversion algorithms

    Laszlo Csanky. Fast parallel matrix inversion algorithms. In16th Annual Symposium on Foundations of Computer Science (sfcs 1975), pages 11–12. IEEE, 1975

  20. [20]

    Convex synthesis and verification of control-lyapunov and barrier functions with input constraints.arXiv preprint arXiv:2210.00629, 2022

    Hongkai Dai and Frank Permenter. Convex synthesis and verification of control-lyapunov and barrier functions with input constraints.arXiv preprint arXiv:2210.00629, 2022

  21. [21]

    Safe nonlinear control using robust neural lyapunov- barrier functions

    Charles Dawson, Zengyi Qin, Sicun Gao, and Chuchu Fan. Safe nonlinear control using robust neural lyapunov- barrier functions. InConference on Robot Learning, pages 1724–1735. PMLR, 2022

  22. [22]

    Augmented lagrangian and alternating direction methods for convex optimiza- tion: A tutorial and some illustrative computational re- sults.RUTCOR Research Reports, 32(3):44, 2012

    Jonathan Eckstein and Wang Yao. Augmented lagrangian and alternating direction methods for convex optimiza- tion: A tutorial and some illustrative computational re- sults.RUTCOR Research Reports, 32(3):44, 2012

  23. [23]

    Two innovations in inexact augmented lagrangian methods for convex optimization.arXiv preprint arXiv:2503.11809, 2025

    Jonathan Eckstein and Chang Yu. Two innovations in inexact augmented lagrangian methods for convex optimization.arXiv preprint arXiv:2503.11809, 2025

  24. [24]

    Hpipm: a high- performance quadratic programming framework for model predictive control.IFAC-PapersOnLine, 53(2): 6563–6569, 2020

    Gianluca Frison and Moritz Diehl. Hpipm: a high- performance quadratic programming framework for model predictive control.IFAC-PapersOnLine, 53(2): 6563–6569, 2020

  25. [25]

    Optimization over state feedback policies for robust con- trol with constraints.Automatica, 42(4):523–533, 2006

    Paul J Goulart, Eric C Kerrigan, and Jan M Maciejowski. Optimization over state feedback policies for robust con- trol with constraints.Automatica, 42(4):523–533, 2006

  26. [26]

    Fastrack: A modular framework for fast and guaranteed safe motion planning

    Sylvia L Herbert, Mo Chen, SooJean Han, Somil Bansal, Jaime F Fisac, and Claire J Tomlin. Fastrack: A modular framework for fast and guaranteed safe motion planning. In2017 IEEE 56th Annual Conference on Decision and Control (CDC), pages 1517–1522. IEEE, 2017

  27. [27]

    A parallel-in-time newton’s method for nonlinear model predictive control.IEEE Transactions on Control Systems Technology, 2025

    Casian Iacob, Hany Abdulsamad, and Simo S ¨arkk¨a. A parallel-in-time newton’s method for nonlinear model predictive control.IEEE Transactions on Control Systems Technology, 2025

  28. [28]

    Parallel and proximal con- strained linear-quadratic methods for real-time nonlinear mpc.arXiv preprint arXiv:2405.09197, 2024

    Wilson Jallet, Ewen Dantec, Etienne Arlaud, Justin Car- pentier, and Nicolas Mansard. Parallel and proximal con- strained linear-quadratic methods for real-time nonlinear mpc.arXiv preprint arXiv:2405.09197, 2024

  29. [29]

    Proxddp: Proximal constrained trajectory optimization

    Wilson Jallet, Antoine Bambade, Etienne Arlaud, Sarah El-Kazdadi, Nicolas Mansard, and Justin Carpentier. Proxddp: Proximal constrained trajectory optimization. IEEE Transactions on Robotics, 2025

  30. [30]

    Cusadi: A gpu paralleliza- tion framework for symbolic expressions and optimal control.IEEE Robotics and Automation Letters, 2024

    Se Hwan Jeon, Seungwoo Hong, Ho Jae Lee, Charles Khazoom, and Sangbae Kim. Cusadi: A gpu paralleliza- tion framework for symbolic expressions and optimal control.IEEE Robotics and Automation Letters, 2024

  31. [31]

    Using neural networks to compute approximate and guar- anteed feasible hamilton-jacobi-bellman pde solutions

    Frank Jiang, Glen Chou, Mo Chen, and Claire J Tomlin. Using neural networks to compute approximate and guar- anteed feasible hamilton-jacobi-bellman pde solutions. arXiv preprint arXiv:1611.03158, 2016

  32. [32]

    Planning with learned dynamics: Probabilis- tic guarantees on safety and reachability via lipschitz constants.IEEE Robotics and Automation Letters, 6(3): 5129–5136, 2021

    Craig Knuth, Glen Chou, Necmiye Ozay, and Dmitry Berenson. Planning with learned dynamics: Probabilis- tic guarantees on safety and reachability via lipschitz constants.IEEE Robotics and Automation Letters, 6(3): 5129–5136, 2021

  33. [33]

    Statistical safety and robustness guarantees for feedback motion planning of unknown underactuated stochastic systems.arXiv preprint arXiv:2212.06874, 2022

    Craig Knuth, Glen Chou, Jamie Reese, and Joe Moore. Statistical safety and robustness guarantees for feedback motion planning of unknown underactuated stochastic systems.arXiv preprint arXiv:2212.06874, 2022

  34. [34]

    A computationally efficient robust model predictive control framework for uncertain nonlin- ear systems.IEEE Transactions on Automatic Control, 66(2):794–801, 2020

    Johannes K ¨ohler, Raffaele Soloperto, Matthias A M ¨uller, and Frank Allg ¨ower. A computationally efficient robust model predictive control framework for uncertain nonlin- ear systems.IEEE Transactions on Automatic Control, 66(2):794–801, 2020

  35. [35]

    Shreyas Kousik, Sean Vaskov, Fan Bu, Matthew Johnson- Roberson, and Ram Vasudevan. Bridging the gap be- tween safety and real-time performance in receding- horizon trajectory design for mobile robots.The Interna- tional Journal of Robotics Research, 39(12):1419–1469, 2020

  36. [36]

    Reach-avoid problems via sum-or-squares opti- mization and dynamic programming

    Benoit Landry, Mo Chen, Scott Hemley, and Marco Pavone. Reach-avoid problems via sum-or-squares opti- mization and dynamic programming. In2018 IEEE/RSJ International Conference on Intelligent Robots and Sys- tems (IROS), pages 4325–4332. IEEE, 2018

  37. [37]

    Robust model predictive control using tubes.Automatica, 40(1):125–133, 2004

    Wilbur Langson, Ioannis Chryssochoos, SV Rakovi ´c, and David Q Mayne. Robust model predictive control using tubes.Automatica, 40(1):125–133, 2004

  38. [38]

    Fast system level synthesis: Robust model predictive control using riccati recursions.IFAC-PapersOnLine, 58(18): 173–180, 2024

    Antoine P Leeman, Johannes Kohler, Florian Messerer, Amon Lahr, Moritz Diehl, and Melanie N Zeilinger. Fast system level synthesis: Robust model predictive control using riccati recursions.IFAC-PapersOnLine, 58(18): 173–180, 2024

  39. [39]

    Robust non- linear optimal control via system level synthesis.IEEE Transactions on Automatic Control, 2025

    Antoine P Leeman, Johannes K ¨ohler, Andrea Zanelli, Samir Bennani, and Melanie N Zeilinger. Robust non- linear optimal control via system level synthesis.IEEE Transactions on Automatic Control, 2025

  40. [40]

    Funnel libraries for real-time robust feedback motion planning.The International Journal of Robotics Research, 36(8):947– 982, 2017

    Anirudha Majumdar and Russ Tedrake. Funnel libraries for real-time robust feedback motion planning.The International Journal of Robotics Research, 36(8):947– 982, 2017

  41. [41]

    Mayne, Maria M

    David Q. Mayne, Maria M. Seron, and Sa ˇsa V . Rakovi´c. Robust output feedback model predictive control of constrained linear systems.Automatica, 41(2):219–224, 2005

  42. [42]

    Tube-based robust nonlinear model pre- dictive control.International journal of robust and nonlinear control, 21(11):1341–1353, 2011

    David Q Mayne, Erric C Kerrigan, EJ Van Wyk, and Paola Falugi. Tube-based robust nonlinear model pre- dictive control.International journal of robust and nonlinear control, 21(11):1341–1353, 2011

  43. [43]

    Formal safety verification and refinement for generative motion planners via certified local stabilization.arXiv preprint arXiv:2509.19688, 2025

    Devesh Nath, Haoran Yin, and Glen Chou. Formal safety verification and refinement for generative motion planners via certified local stabilization.arXiv preprint arXiv:2509.19688, 2025

  44. [44]

    Tinympc: Model- predictive control on resource-constrained microcon- trollers

    Khai Nguyen, Sam Schoedel, Anoushka Alavilli, Brian Plancher, and Zachary Manchester. Tinympc: Model- predictive control on resource-constrained microcon- trollers. In2024 IEEE International Conference on Robotics and Automation (ICRA), pages 1–7. IEEE, 2024

  45. [45]

    Wright.Numerical Opti- mization

    Jorge Nocedal and Stephen J. Wright.Numerical Opti- mization. Springer, New York, NY , 2nd edition, 2006

  46. [46]

    Cyqlone: A parallel, high-performance linear solver for optimal control.arXiv preprint arXiv:2512.09058, 2025

    Pieter Pas and Panagiotis Patrinos. Cyqlone: A parallel, high-performance linear solver for optimal control.arXiv preprint arXiv:2512.09058, 2025

  47. [47]

    Rawlings, David Q

    James B. Rawlings, David Q. Mayne, and Moritz Diehl. Model Predictive Control: Theory, Computation, and Design. Nob Hill Publishing, Madison, WI, 2nd edition, 2017

  48. [48]

    Learning control barrier functions from expert demonstrations

    Alexander Robey, Haimin Hu, Lars Lindemann, Hanwen Zhang, Dimos V Dimarogonas, Stephen Tu, and Nikolai Matni. Learning control barrier functions from expert demonstrations. In2020 59th IEEE Conference on Decision and Control (CDC), pages 3717–3724. Ieee, 2020

  49. [49]

    Tempo- ral parallelization of dynamic programming and linear quadratic control.IEEE Transactions on Automatic Control, 68(2):851–866, 2022

    Simo S ¨arkk¨a and ´Angel F Garc ´ıa-Fern´andez. Tempo- ral parallelization of dynamic programming and linear quadratic control.IEEE Transactions on Automatic Control, 68(2):851–866, 2022

  50. [50]

    Learning barrier functions for constrained motion planning with dynami- cal systems

    Matteo Saveriano and Dongheui Lee. Learning barrier functions for constrained motion planning with dynami- cal systems. In2019 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), pages 112–

  51. [51]

    Robust tracking with model mismatch for fast and safe planning: an sos optimization approach

    Sumeet Singh, Mo Chen, Sylvia L Herbert, Claire J Tomlin, and Marco Pavone. Robust tracking with model mismatch for fast and safe planning: an sos optimization approach. InInternational Workshop on the Algorithmic Foundations of Robotics, pages 545–564. Springer, 2018

  52. [52]

    Robust feedback motion planning via contraction theory.The International Journal of Robotics Research, 42(9):655–688, 2023

    Sumeet Singh, Benoit Landry, Anirudha Majumdar, Jean- Jacques Slotine, and Marco Pavone. Robust feedback motion planning via contraction theory.The International Journal of Robotics Research, 42(9):655–688, 2023

  53. [53]

    How to train your neural control barrier function: Learning safety filters for complex input-constrained systems

    Oswin So, Zachary Serlin, Makai Mann, Jake Gonzales, Kwesi Rutledge, Nicholas Roy, and Chuchu Fan. How to train your neural control barrier function: Learning safety filters for complex input-constrained systems. In 2024 IEEE International Conference on Robotics and Automation (ICRA), pages 11532–11539. IEEE, 2024

  54. [54]

    Safety beyond the training data: Robust out-of- distribution mpc via conformalized system level synthe- sis.arXiv preprint arXiv:2602.12047, 2026

    Anutam Srinivasan, Antoine Leeman, and Glen Chou. Safety beyond the training data: Robust out-of- distribution mpc via conformalized system level synthe- sis.arXiv preprint arXiv:2602.12047, 2026

  55. [55]

    Osqp: An operator splitting solver for quadratic programs.Mathematical Programming Computation, 12(4):637–672, 2020

    Bartolomeo Stellato, Goran Banjac, Paul Goulart, Al- berto Bemporad, and Stephen Boyd. Osqp: An operator splitting solver for quadratic programs.Mathematical Programming Computation, 12(4):637–672, 2020

  56. [56]

    A generalized itera- tive lqg method for locally-optimal feedback control of constrained nonlinear stochastic systems

    Emanuel Todorov and Weiwei Li. A generalized itera- tive lqg method for locally-optimal feedback control of constrained nonlinear stochastic systems. InProceedings of the 2005, American Control Conference, 2005., pages 300–306. IEEE, 2005

  57. [57]

    Robustly Constrained Dynamic Games for Uncertain Nonlinear Dynamics

    Shuyu Zhan, Chih-Yuan Chiu, Antoine P Leeman, and Glen Chou. Robustly constrained dynamic games for uncertain nonlinear dynamics.arXiv preprint arXiv:2509.16826, 2025. APPENDIX In this section, we provide supplemental material to sup- port the information provided in our methodology. App. A discusses the mechanism of parallel associative scans, which are ...