Hi everyone, I've come across a memory issue when trying to assign data to slices of a Numpy memory-mapped array. The short story is that if I create a memory mapped array and try to add data to subsets of the array many times in a loop, the memory usage of my code grows over time, suggesting there is some kind of memory leak. More specifically, if I run the following script: import random import numpy as np image = np.memmap('image.np', mode='w+', dtype=np.float32, shape=(10000, 10000)) print("Before assignment") for i in range(1000): x = random.uniform(1000, 9000) y = random.uniform(1000, 9000) imin = int(x) - 128 imax = int(x) + 128 jmin = int(y) - 128 jmax = int(y) + 128 data = np.random.random((256,256)) image[imin:imax, jmin:jmax] = image[imin:imax, jmin:jmax] + data del x, y, imin, imax, jmin, jmax, data the memory usage goes up to ~300Mb after 1000 iterations (and proportionally more if I increase the number of iterations). I've written up a more detailed overview of the issue on stackoverflow (with memory profiling): http://stackoverflow.com/questions/15473377/memory-issue-with-numpy-memory-m... Does anyone have any idea what is going on, and how I can avoid this issue? Thanks! Tom