Python scipy find_simplex冻结

Question

我试着用scipy.spatial.Delaunay公司实现“基于Delaunay图映射的快速动态网格变形”：

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

我试图使用find_simplex来找出哪些点在哪个三角形/tet中，但由于某些原因，该过程冻结在一个特定点上。在

我已经分解了整个域，所以每个点都应该在一个三角形/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版本。在