二维函数的三角剖分

我有一个函数 f(x,y) = z，我想对它进行三角剖分，以便在线性程序中处理它。我的输入 x、y 和 z 定义如下：

x = arange(0,5,0.25)
y = arange(0,0.25)
z = curve(x,y)

我现在用 tri = delaunay(np.array([x,y]).T) 尝试了 scipy.spatial.Delaunay，但出现错误：

QH6154 Qhull 精度误差：初始单纯形是平坦的（facet 1 与内点共面）

这似乎是因为我的 x 和 y 是一样的。但它们应该是，因为我想在这个“网格”上进行三角测量。

我也尝试了 triang = mtri.Triangulation(x,y) 中的 matplotlib.tri 但我得到了同样的错误（因为 delaunay 在内部使用）

所以我的问题是如何在 f(x,y) = z 网格上对函数 x,y 进行三角测量。提前致谢。