[SciPy-Dev] flapack and the GIL

[Pauli Virtanen]

Hi,

04.02.2013 13:24, Sturla Molden kirjoitti:
> It seems arbitrary if the flapack f2py interfaces have a 'threadsafe' 
> statement or not. Why is that?
> 
> https://github.com/scipy/scipy/blob/master/scipy/linalg/flapack.pyf.src
> 
> I'd assume the whole of BLAS/LAPACK would be reeantrant or not. And this 
> would possibly depend on the libraries itself, e.g. if they are compiled 
> to assume all local variables are 'STATIC' by default.
> Which means that the right place to determine threadsafety of flapack 
> routines would be in the build process.

I think those have been added randomly, so there's no logic behind it

AFAIK, *all* of Lapack is threadsafe --- there are no SAVE or COMMON
statements. The same is probably true for any modified LAPACK routines
in MKL, Accelerate and the ilk.

-- 
Pauli Virtanen