[SciPy-dev] possible bug in interp1d

[Darren Dale]

Running this short script yields an error message:

from scipy import *
i = interpolate.interp1d(arange(100), arange(100))
i(arange(10, 90, 100))

exceptions.AttributeError                            Traceback (most recent 
call last)

/home/darren/<console>

/usr/lib64/python2.4/site-packages/scipy/interpolate/interpolate.py in 
__call__(self, x_new)
    196         out_of_bounds.shape = sec_shape
    197         new_out = ones(new_shape)*out_of_bounds
--> 198         putmask(y_new, new_out.flat, self.fill_value)
    199         y_new.shape = yshape
    200         # Rotate the values of y_new back so that they coorespond to 
the

/usr/lib64/python2.4/site-packages/scipy/base/oldnumeric.py in putmask(a, 
mask, v)
    186     """
    187     print dir(a)
--> 188     return a.putmask(v, mask)
    189
    190 def swapaxes(a, axis1, axis2):

AttributeError: 'scipy.flatiter' object has no attribute 'putmask'

I checked the namespace of a in oldnumeric.putmask:

['__array__', '__class__', '__delattr__', '__delitem__', '__doc__', 
'__getattribute__', '__getitem__', '__hash__', '__init__', '__iter__', 
'__len__', '__new__', '__reduce__', '__reduce_ex__', '__repr__', 
'__setattr__', '__setitem__', '__str__', 'base', 'copy', 'next']


Darren