P = randint(0,2, (20,13))

for i in range(50):

    for ai in ndindex((2,)*13):

        j = P.dot(ai)

If you replace the last line with something like j = ai, the memory leak goes away. I’m not exactly sure what’s going on but the .dot seems to be causing the memory taken by the tuple ai to be held.

This devours RAM in python 2.7.5 (OS X Mavericks default I believe), numpy version 1.8.0.dev-3084618. I’m upgrading to the latest Superpack (numpy 1.9) right now but I somehow doubt this behavior will change.

Any thoughts?

Best, Chris