def METHOD_NAME(self):
    np.random.seed(0)
    pts = np.random.rand(50, 3)
    tri = op.network.Delaunay(points=pts, shape=[1, 1, 1])
    assert tri.coords.shape == (50, 3)
    assert np.all(tri.coords == pts)