def METHOD_NAME(p: npt.NDArray[np.floating], q: npt.NDArray[np.floating]) -> npt.NDArray[np.floating]:
    """Compute the distance between two convex bodies p and q, each defined by a set of points.
    Uses the Gilbert–Johnson–Keerthi (GJK) distance algorithm.
    Args:
        p: Body 1 list of points (``shape=(num_points, gdim)``).
        q: Body 2 list of points (``shape=(num_points, gdim)``).
    Returns:
        Shortest vector between the two bodies.
    """
    return _cpp.geometry.METHOD_NAME(p, q)