Python scipy find_simplex 冻结

Question

我正在尝试使用scipy.spatial.Delaunay来实现“基于Delaunay图映射的快速动态网格变形”：

http://www.sciencedirect.com/science/article/pii/S0021999105002809

我正在尝试使用 find_simplex 来查找哪些点在哪个三角形/四边形中，但由于某种原因，该过程冻结了一个特定点。

我已经分解了整个域，所以每个点都应该在一个三角形/tet 内。

这些点是围绕机翼段的网格的一大组坐标。翼段表面上的点和远场构成顶点。导致我出现问题的代码段如下：

 for j in xrange(len(points[:,0])):
      simp = tri.find_simplex(points[j,:3]) <<freezes here at j=213877
      vert = tri.vertices[simp]

我正在使用： SciPy 0.9.0 和 Python 2.6.5

我知道这是有限的信息，因为我没有提供数据集，但有人知道是什么原因造成的，或者可能建议解决方法吗？

提前致谢

Answer 1

该解决方案很可能只是升级到更新的 Scipy 版本。