lagrange - signal interpolation python

Scipy griddata no funciona dentro de un bucle/pérdida de memoria

Question

Tengo un problema al utilizar la cuadrícula de Scipy dentro de un ciclo. Básicamente lo que sucede es que la memoria crece sin límites mientras el ciclo se está ejecutando.

Para reproducir el problema simplemente pon el ejemplo en

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

dentro de un bucle:

for i in range(100000):

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

La versión de My Python es 2.7.3, mi versión numpy es 1.7.0 y mi versión scipy es 0.12.0b1. Lo estoy ejecutando en Windows 7.

¿Es esto un error? ¿Cómo puedo repetir la interpolación muchas veces sin incurrir en un problema de pérdida de memoria?

Con el resto del código:

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'')

Gracias por adelantado.

Answer 1

Es un error en Cython , debería solucionarse en la versión final de Scipy 0.12.0.