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