Blender qhull stuff

test.py

import bpy
import numpy as np
from scipy.spatial import Delaunay

points = np.random.rand(20,3)
tri = Delaunay(points)

mesh = bpy.data.meshes.new("myUnrulyMess")  # add the new mesh
obj = bpy.data.objects.new(mesh.name, mesh)
col = bpy.data.collections.get("Collection")
col.objects.link(obj)
bpy.context.view_layer.objects.active = obj

verts = [list(v) for v in points]
edges = []
faces = np.ndarray.tolist(tri.simplices)

mesh.from_pydata(verts, edges, faces)