dt = Delaunay(pts[ch.vertices])
tets = dt.points[dt.simplices]
return np.sum(tetrahedron_volume(tets[:, (0)], tets[:, (1)], tets[:, (2)],
    tets[:, (3)]))

def convex_hull_volume_bis(pts):
    ch = ConvexHull(pts)
simplices = np.column_stack((np.repeat(ch.vertices[0], ch.nsimplex), ch.
    simplices))
tets = ch.points[simplices]
return np.sum(tetrahedron_volume(tets[:, (0)], tets[:, (1)], tets[:, (2)],
    tets[:, (3)]))
Volume of convex hull with QHull from SciPy