[SciPy-Dev] SciPy 0.12.0/Numpy 1.7.1 sgerqf segfault

[Orion Poplawski]

On 04/15/2013 03:32 PM, Pauli Virtanen wrote:
> 16.04.2013 00:13, Orion Poplawski kirjoitti:
> [clip]
>>> (scipy.linalg doesn't have it's own python_xerbla.c, so it ends up
>>> looking for the symbol in numpy.linalg.lapack_lite --- rather unexpected!)
>>
>> So, any suggestions for this?  Should flapack.pyf.src not mark them as thread
>> safe?  Should scipy provide its own xerbla?
>
> Yes, to the latter.
>
> The actual fix to this bug is with high probability a replacement from
>
> integer optional,intent(in),depend(n),check(lwork>=m||lwork==-1) ::
> lwork=3*m
>
> to
>
> integer optional,intent(in),depend(m),check(lwork>=m||lwork==-1) ::
> lwork=3*m
>
> on line 653 of scipy/linalg/flapack.pyf.src. This bug doesn't manifest
> always as it depends on ordering of dict entries inside f2py...

I can confirm that this fixes it.  Thanks!



-- 
Orion Poplawski
Technical Manager                     303-415-9701 x222
NWRA, Boulder/CoRA Office             FAX: 303-415-9702
3380 Mitchell Lane                       orion at nwra.com
Boulder, CO 80301                   http://www.nwra.com
-------------- next part --------------
A non-text attachment was scrubbed...
Name: scipy-gerqf.patch
Type: text/x-patch
Size: 731 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20130416/68b7f425/attachment.bin>