Safety-critical control of autonomous systems requires formal safety certificates, such as controlled invariant sets, that must be computed online as conditions change. Although standard synthesis algorithms scale poorly with state dimension, monotone dynamical systems with lower-closed safety specifications allow for accelerated computation of controlled invariant sets. In particular, lazy fixed-point algorithms exploit monotonicity and track only the antichain basis of the set. However, membership tests and redundancy checks against an evolving basis remain major bottlenecks. We introduce a threshold-function reformulation in which a lower-closed set on a d-dimensional grid is represented by its column heights along a designated axis. This reformulates the greatest-fixed-point iteration as independent one-dimensional binary searches, one per grid column, yielding an embarrassingly parallel iteration with asymptotically lower computational complexity than the lazy fixed-point algorithm. Experiments synthesize invariant sets on 3D grids with 10^9 cells in under 50 ms and 10^14 cells in under two minutes. We further demonstrate online re-synthesis in a safety-informed model predictive controller example.
Figures & tables
Fig. 1: Threshold representation and update. (a) A lower set K is stored by the column height τK(k) , so membership is the O(1) test (k,v)∈K⟺v≤τK(k) . (b) For each column, Φ(K) is a prefix whose height f(k) is found by binary search; all columns are updated independently.
Lazy [ 2 ]
Threshold (ours)
Predecessor tests
O(CN2(d−1))
O(CNd−1logN)
sequential over evolving B
parallel over columns
Basis update
O(N2(d−1))
—
sequential
—
TABLE I: Per-iteration complexity comparison.
Case
∣X∣
tpre
Lazy
Threshold
Rlazy
tlazy
tlazy+τ
Rthr
tthr
Speedup
ACC
106
.569ms
182
63.1s
106ms
26
3.18ms
33.4×
ACC
109
439ms
1595
Timeout
126s
24
16.9ms
7,458×
ACC
1010
Timeout
24
318ms
-
ACC
1012
Timeout
24
750ms
-
ACC
1014
Timeout
24
7.94s
-
TABLE II: Lazy vs. Threshold Comparison. tpre is the successor-table generation time used in Lazy methods and is excluded from the reported speedups, which equal tlazy+τ/tthr ; R = iterations; timeout = 20 minutes.
Fig. 2: Real-time controller for an unprotected left turn. An MPPI controller samples rollouts and queries the threshold table once per rollout step to evaluate the set-membership penalty. The resulting optimized control input is applied. When oncoming velocities change, the invariant set is resynthesized.
Quantity
Value
Grid cells ∣X∣
1010=104×102×104
Sampling period Ts
0.1 s
Prediction and control horizon Hp,Hc
20 steps
MPPI iterations / step Nit
10
Rollouts / iteration Nr
512
Perturbation standard deviation σu
1000 Nm
TABLE III: Real-time unprotected left turn: configuration and timings.