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的。