Scipy griddata在循环/内存泄漏中不起作用

时间:2013-03-15 15:23:30

标签: python scipy interpolation

我在循环中使用Scipy的griddata时遇到问题。基本上发生的事情是在循环运行时内存不受限制地增长。

要重现问题,请将示例放在

http://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.griddata.html

循环中的

for i in range(100000):

    grid_z1 = griddata(points, values, (grid_x, grid_y), method='linear')

我的Python版本是2.7.3,我的numpy版本是1.7.0,我的scipy版本是0.12.0b1。我在WIndows 7上运行它。

这是一个错误吗?如何在不引发内存泄漏问题的情况下多次重复插值?

使用其余代码:

def func(x, y):
    return x*(1-x)*np.cos(4*np.pi*x) * np.sin(4*np.pi*y**2)**2

grid_x, grid_y = np.mgrid[0:1:100j, 0:1:200j]
points = np.random.rand(1000, 2)
values = func(points[:,0], points[:,1])

for i in range(100000):

    grid_z1 = griddata(points, values, (grid_x, grid_y), method='linear')

提前致谢。

1 个答案:

答案 0 :(得分:0)

这是一个bug in Cython,应该在Scipy最终的0.12.0版本中解决。