Ticket #58: hetland.py

from scikits import delaunay # or wherever your delaunay package lives these days
from numpy.random import rand

xi, yi = rand(2, 1000)
x, y = rand(2, 100)
tri = delaunay.Triangulation(xi, yi)

for n in range(100000):
    interp = tri.nn_interpolator(rand(1000))
    z = interp(x, y)

# Segfaults