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

Para reproducir el problema solo 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')

Mi versión de 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.