vix.ing · top · new · best · stats · spec

Rake-Compress Riccati Recursions for Parallel Scenario-Tree Model Predictive Control

2026/08/02 by João Sousa-Pinto
Mathematics · Computer Science · #math.OC #cs.RO #msc:49N10 #msc:65Y05 #msc:68W10 #msc:49M27 #acm:49N10 #acm:65Y05 #acm:68W10 #acm:49M27

paper · pdf

16 pages

arxiv created 2026/08/02 · arxiv updated 2026/08/04

Abstract

Scenario-tree model predictive control (MPC) represents future information by a rooted tree and optimizes a nonanticipative policy over that tree. Numerical methods for solving the resulting nonlinear program typically compute their search directions through a sequence of branched linear-quadratic regulator (LQR) subproblems. The standard tree Riccati recursion requires linear work but has a dependency chain proportional to tree height. We present an algebraically exact parallel solver based on rake-compress tree contraction. After independent local control condensation, its two operations act on node and edge data that represent conditional quadratic functions. A rake eliminates a leaf and its parent edge, adding their reduced contribution to the parent-node data. A compress eliminates a unary node and replaces its two adjacent edges by one edge, using the same conditional-value composition as parallel Riccati methods on a chain. Together they contract an arbitrary rooted tree to its root; reversing the contraction recovers every Riccati coefficient, state, control, and multiplier. Given a reusable topology plan, a solve with N nodes and fixed state and control dimensions has O(N) arithmetic work and storage and O(log N) span, independently of tree height, balance, and maximum out-degree. The formulation allows positive-semidefinite dual regularization, including the unregularized case, and an exact linear-size lifting covers the standard scenario-MPC convention of one control per information node. We prove the contraction identities and equivalence to the Karush-Kuhn-Tucker (KKT) system. Three MIT-licensed JAX packages implement the bidirectional contraction, the dual-regularized LQR solver, and a user-facing primal-dual interior-point solver for tree-structured optimal control.

Citations