[Scipy-tickets] [SciPy] #1727: LinearNDInterpolator hangs on 4D data set

SciPy Trac scipy-tickets@scipy....
Thu Sep 13 18:12:08 CDT 2012


#1727: LinearNDInterpolator hangs on 4D data set
-----------------------------------------------+----------------------------
 Reporter:  miked983                           |       Owner:  somebody
     Type:  defect                             |      Status:  new     
 Priority:  high                               |   Milestone:          
Component:  scipy.interpolate                  |     Version:  0.10.0  
 Keywords:  LinearNDInterpolator, interpolate  |  
-----------------------------------------------+----------------------------
 The set is fairly large with 250,000 entries but it obviously is getting
 stuck somewhere. I've left it running 30 minutes and the memory usage
 steadily grows with no result, while MATLAB's interpn solves the system
 immediately. The code is as follows:

 {{{
 import numpy as np
 import scipy.io
 from scipy.interpolate import LinearNDInterpolator

 data = scipy.io.loadmat('data.mat')['data']
 pts =
 np.array((data['var1'][0][0].ravel(),data['var2'][0][0].ravel(),data['var3'][0][0].ravel(),data['var4'][0][0].ravel())).T
 vals = data['var5'][0][0].ravel()
 interpFun = LinearNDInterpolator(pts,vals)
 }}}

 A similar issue is mentioned here:[http://mail.scipy.org/pipermail/scipy-
 user/2011-March/028841.html]
 and in the ticket that that post points to. However, it is stated that the
 problem has been fixed; perhaps this is something new...

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1727>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list