Massively parallel computation of globally optimal shortest paths with curvature penalization

Jean-Marie Mirebeau (1) , Lionel Gayraud (2) , Rémi Barrère (2) , Da Chen (3) , François Desquilbet (4)

1 CB - CB - Centre Borelli - UMR 9010
2 Thales Research and Technology [Palaiseau]
3 Shandong Artificial Intelligence Institute
4 EDP - Equations aux Dérivées Partielles

We address the computation of paths globally minimizing an energy involving their curvature, with given endpoints and tangents at these endpoints, according to models known as the Reeds-Shepp car (reversible or forward only), the Euler-Mumford elasticae, and the Dubins car. For that purpose, we numerically solve degenerate variants of the eikonal equation, on a three dimensional domain, in a massively manner on a graphical processing unit. Due to the high anisotropy and non-linearity of the addressed PDE, the discretization stencil is rather wide, has numerous elements, and is costly to generate, which leads to subtle compromises between computational cost, memory usage, and cache coherency. Accelerations by a factor 30 to 120 are obtained w.r.t a sequential implementation. The efficiency and robustness of the method is illustrated in various contexts, ranging from motion planning to vessel segmentation and radar configuration.

Mots clés: Eikonal equation; curvature penalization; GPU computing

Domaines: Analyse numérique [math.NA]; Equations aux dérivées partielles [math.AP]

Jean-Marie Mirebeau, Lionel Gayraud, Rémi Barrère, Da Chen, François Desquilbet. Massively parallel computation of globally optimal shortest paths with curvature penalization. Concurrency and Computation: Practice and Experience, 2022, 35 (2), pp.e7472. ⟨10.1002/cpe.7472⟩. ⟨hal-03171069v2⟩