[f2py] I: double free or corruption after calling a Slicot routine wrapped with f2py

Hello,

i thought maybe this mailing list would be better for finding some clue about my problem.

regards,

/Enrico

--- Mar 23/2/10, enrico avventi <[email protected]> ha scritto:

Da: enrico avventi <[email protected]>
Oggetto: double free or corruption after calling a Slicot routine wrapped with f2py
A: [email protected]
Data: Martedì 23 febbraio 2010, 17:32

hello,

first of all, as i'm new here, i would like to greet everyone in the list and thank the developers of numpy/scipy. i'm transitioning my work from matlab to python and this software is very helpfull indeed.

the reason i'm writing is that i got to a stumbling block last week. i tried wrote some wrappers of Slicot routines aided by f2py that i keep on github (http://github.com/avventi/Slycot).
the latest routine i tried to wrap, SB02OD, make python to crash with the glibc error double free or corruption. precisely i get this error if the wrapper slycot.sb02od is called within a method, i.e

Python 2.5.2 (r252:60911, Jan 24 2010, 14:53:14) 
[GCC 4.3.2] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> import slycot
>>> slycot.examples.sb02od_example()
--- Example for
 sb01od ...
The solution X is
[[ 1.  0.]
 [ 0.  2.]]
rcond = 0.632455532034
*** glibc detected *** python: double free or corruption (!prev): 0x082ec3b8 ***

as you can see the routine does indeed return the correct output but make python crash afterwards. On the other hand if i type ieach step interactively

Python 2.5.2 (r252:60911, Jan 24 2010, 14:53:14) 
[GCC 4.3.2] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> from numpy import *
>>> import slycot
>>> A = array([[0,1],[0,0]])
>>> B = array([[0],[1]])
>>> C = array([[1,0],[0,1],[0,0]])
>>> Q = dot(C.T,C)
>>> R = zeros((1,1))
>>> L = zeros((2,1))
>>> out = slycot.sb02od('D',2,1,3,A,B,Q,R,L)
>>> out[1]
array([[ 1.,  0.],
       [ 0., 
 2.]])
>>> out[0]
0.63245553203367577
>>> out = slycot.sb02od('D',2,1,3,A,B,Q,R,L)
*** glibc detected *** python: double free or corruption (!prev): 0x0832b428 ***

it works the first time but not the second. i tried Debian lenny 32/64bit and a virtualized Fedora 12 32bit and the error persists.

do you have any ideas why this happens or where should i look to start solving it?

thanks in advance.
regards,

/Enrico

_______________________________________________
f2py-users mailing list
f2py-users-Y4l6ocDipWCuvFJfX82//[email protected]
http://cens.ioc.ee/mailman/listinfo/f2py-users