Segfault in Python 2.5

[Nikhil Padmanabhan]

Hi,

I recently tried to upgrade to Python 2.5 (MacPython), and I am  
getting segfaults on some of my f2py wrapped fortran codes. I'm  
running OS X 10.4.8 (ppc), and this occurs for both numpy 1.0 and  
1.0rc2 (the two versions I tried). Running with Python 2.4 does not  
give me this error. This is using the IBM xlf compiler....

A simple test case is below.....

Not sure where to go from here -- any help would be great!
Thanks,

Here is a stripped down version that still crashes :

subroutine foo(xx, nxx, bj)
   implicit none
   integer nxx
   real*8 xx(nxx),bj(nxx)
!f2py intent(out) bj

   bj = xx

end subroutine foo

Compiling this to foo.so, I get :
Python 2.5 (r25:51918, Sep 19 2006, 08:49:13)
[GCC 4.0.1 (Apple Computer, Inc. build 5341)] on darwin
Type "help", "copyright", "credits" or "license" for more information.
 >>> from numpy import *
 >>> import foo as p
 >>> a = arange(10)
 >>> p.foo(a)
Segmentation fault

Intriguingly, swapping the import orders fixes this :
Python 2.5 (r25:51918, Sep 19 2006, 08:49:13)
[GCC 4.0.1 (Apple Computer, Inc. build 5341)] on darwin
Type "help", "copyright", "credits" or "license" for more information.
 >>> import foo as p
 >>> from numpy import *
 >>> a = arange(10)
 >>> p.foo(a)
array([ 0.,  1.,  2.,  3.,  4.,  5.,  6.,  7.,  8.,  9.])
 >>>



------------------------------------
Nikhil Padmanabhan
NPadmanabhan at lbl.gov
nikhil@(510) 495-2943




-------------------------------------------------------------------------
Using Tomcat but need to do more? Need to support web services, security?
Get stuff done quickly with pre-integrated technology to make your job easier
Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642