安全性解析/GPU計算2026/11/3arXiv:2611.04096

ロボット制御ループ上で動くリアルタイムGPU版Hamilton-Jacobi安全性解析ソルバ

Online Hamilton-Jacobi Reachability Onboard the Robot Control Loop: A Real-Time GPU Solver for Safety-Critical Autonomy

WENO5+TVD-RK3+Lax-FriedrichsのHJ-Isaacsレベルセット解法をCUDAで実装し、2〜4次元の安全到達管(BRT)を制御ループ内でリアルタイム計算できるようにした。

詳しい要約

1. どんなもの？

- HJ reachability の BRT をオンラインで計算する GPU ソルバの提案 - 対象: 2D/3D/4D の動力学に対する level-set HJ-Isaacs 計算 - 実装: CUDA による WENO5 + TVD-RK3 + Lax-Friedrichs スキーム - 目的: ロボットの制御ループ内で動かせる決定論的安全層の提供 - 応用: ADAS/AV、humanoid、VLA 等の black-box 方策の補完 - 拡張: in-kernel label propagation による K-agent capture-basin Voronoi 分割

2. 先行研究と比べてどこがすごい？

- 従来 HJ reachability は状態次元に対し指数的に計算量が増え、offline 実行が通例 - 本研究は制御ループ内実行可能な速度を達成 - JAX 参照実装 hj_reachability 比で FP32 は 5.9-17.5x 高速 - 同一スキームの torch.compile 再実装比でも 1.7-16x 高速 - FP64 でも hj_reachability 比 2.0-4.6x 高速かつ一致 - 組込み NVIDIA Jetson AGX Thor 上でも 3.8-8.9x 高速

3. 技術・手法の肝は？

- カノニカルな level-set HJ-Isaacs ソルバを CUDA 実装 - 数値スキーム: WENO5 + TVD-RK3 + Lax-Friedrichs - FP32 と FP64 の両精度に対応 - 単一ワークステーション GPU (NVIDIA RTX 4000 Ada) を想定 - 組込み GPU (Jetson AGX Thor) でも動作 - in-kernel label propagation で単一の統合 HJ 解から K-agent の capture-basin Voronoi 分割を復元し、fleet size に依存しない

4. どうやって有効だと検証した？

- FP32 で 3D Air3D BRT 112 Hz、4D pursuit-evasion 61 Hz、2D wind-affected drifter 74-147 Hz - FP64 で 8 ベンチマークケースにて hj_reachability と一致 - Mitchell の Air3D BRT の参照セル数に対し 0.04% 以内で再現 - Jetson AGX Thor で Air3D 66 Hz (N=41)、2D drifter 54 Hz - 6 機の quadcopter による time-varying・非一様風下の wildfire-response 閉ループシミュレーションで 4.0 ms/tick

5. 議論はある？

- 状態次元の指数的スケーリングへの対処は 2D-4D に限定され、高次元化は要旨からは不明 - FP32 と FP64 の精度差の安全保証への影響は要旨からは不明 - 実機ロボットでの閉ループ検証は wildfire シミュレーションのみで、実機実験は要旨からは不明 - black-box 方策との統合時の安全層設計の詳細は要旨からは不明

6. 次に読むべき論文は？

- hj_reachability (JAX 参照実装) - Mitchell の Air3D BRT 研究 - WENO5 / TVD-RK3 / Lax-Friedrichs の level-set 数値解法 - torch.compile による再実装比較 - capture-basin Voronoi tessellation と multi-agent HJ 関連研究

※ AIが要旨から生成した要約です。正確性は原文をご確認ください。

著者: Kevin Weekly, Alexandre M. Bayen

分類: cs.RO, eess.SY

原文アブストラクト

Hamilton-Jacobi (HJ) reachability provides certified backward reachable tubes (BRTs) for safety-critical robotic systems, but the level-set computation scales exponentially in the state dimension and has historically been performed offline. As autonomy stacks (ADAS/AV, humanoid robotics, and vision-language-action models) increasingly rely on black-box, probabilistic policies, a deterministic safety layer with formal guarantees is an essential complement. We present a CUDA implementation of the canonical WENO5 + TVD-RK3 + Lax-Friedrichs level-set Hamilton-Jacobi-Isaacs solver for 2D, 3D, and 4D dynamics that is fast enough to run inside a robot's control loop. In 32-bit floating point (FP32) on a single workstation GPU (NVIDIA RTX 4000 Ada), it computes the 3D Air3D BRT at 112 Hz, a 4D pursuit-evasion BRT at 61 Hz, and a 2D wind-affected drifter BRT at 74-147 Hz, 5.9-17.5x faster than the JAX reference hj_reachability at matched precision and 1.7-16x faster than a torch.compile reimplementation of the same scheme. In 64-bit floating point (FP64) it agrees with hj_reachability on eight benchmark cases, reproduces Mitchell's Air3D BRT to within 0.04% of the reference cell count, and runs 2.0-4.6x faster. On an embedded NVIDIA Jetson AGX Thor it computes Air3D at 66 Hz (N=41) and the 2D drifter at 54 Hz, 3.8-8.9x faster than hj_reachability on the same device. As a multi-agent application, an in-kernel label-propagation extension recovers the K-agent capture-basin Voronoi tessellation from a single combined HJ solve, independent of fleet size; it runs at 4.0 ms per tick in a closed-loop wildfire-response simulation of six quadcopters under time-varying, non-uniform wind.

arXivで原文を読む →arXiv概要ページ

← 論文一覧