[SciPy-user] amd64 specific error in fitpack spline interpolation

[George Nurser]

On 30 Aug 2006, at 11:00, Angus McMorland wrote:

> Hi all,
>
> I posted about this problem a few months ago, but we didn't get very
> far, but the annoyance level of the problem is enough now for me to  
> have
> another go.
>
> I'm trying to run the spline example from the cookbook (my exact code,
> called spl.py, is inlined below). On an ix86 computer this works fine,
> and returns the expected graphs, but on an amd64 machine I get the  
> error
> below. Both computers are running debian etch, with recent svn  
> numpy and
> scipy.
>
> TypeError: array cannot be safely cast to required type
>> /usr/lib/python2.4/site-packages/scipy/interpolate/fitpack.py(217) 
>> splprep()
>     216
> t,c,o=_fitpack._parcur(ravel(transpose(x)),w,u,ub,ue,k,task,ipar,s,t,
> --> 217                              nest,wrk,iwrk,per)
>     218     _parcur_cache['u']=o['u']
>
> I believe the error is generated by the _parcur call -> once in pdb,
> calling it again re-raises the same error.
>
> ipdb>_fitpack._parcur(ravel(transpose(x)),w,u,ub,ue,k,task,ipar,s,t,
> nest,wrk,iwrk,per)

Same problem on our Opterons.
numpy v 2631
scipy  v 1614
built & linked to acml

[8]nohow@/noc/users/agn/python> python amdtest.py
Traceback (most recent call last):
   File "amdtest.py", line 23, in ?
     tckp,u = splprep([x,y,z],s=s,k=k,nest=-1)
   File "/data/jrd/mod1/agn/ext/Linux/lib64/python2.3/site-packages/ 
scipy/interpolate/fitpack.py", line 217, in splprep
TypeError: array cannot be safely cast to required type


-George.