Python scipy find_simplex 卡住

Question

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

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

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

我已经分解了整个域，因此每个点都应该在一个三角形/四边形内。

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

 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 版本。