Python 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和Python2.6.5

我知道这是有限的信息，因为我没有提供数据集，但有谁知道是什么原因造成这一点，或可能建议一个解决方案？

提前感谢

Answer 1

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