[SciPy-User] Fix for griddata unsorted data

Wolfgang Kerzendorf wkerzendorf@googlemail....
Fri Mar 18 08:55:06 CDT 2011

Dear all,

I realized that for the LinearNDInterpolator to work the data needs to 
be sorted. This is explained in griddata, but unfortunatley is not 
included in the LinearNDInterpolator description. It would be great if 
that would be updated. I also think that (as most things in scipy) it 
should just work.

Here's my proposed fix. Before giving the points and values to qhull we 
do this:

sortIDx = np.lexsort(points.transpose()[::-1])
LinearNDInterpolator(points[sortIDx], newValues[sortIDx])

I think that this should work.

I really think that Qhull is an amazing addition to scipy. It basically 
sped up my computations by orders of magnitude.

While playing around with Qhull, I stumbled across another gem for 
numerics: PyMinuit. It is an extremely good optimizer, very easy to use 
and I was wondering if that would be interesting to include in scipy.


More information about the SciPy-User mailing list