On Mon, Sep 2, 2013 at 7:31 PM, Pauli Virtanen <pav@iki.fi> wrote:
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Hi,

02.09.2013 21:12, David Cournapeau kirjoitti:
[clip]
> The pb is specific to MKL, yes, but I've found a workaround: in the
> case of MKL, you just need to workaround the g77 ABI of the MKL vs
> gfortran ABI, but the LAPACK interface is the actual LAPACK, not
> CLAPACK.
>
> So for MKL, you need non dummy wrappers iff the function returns a
> complex. A dirty patch seems to confirm that this fixes the issue,
> I will prepare an actual patch tomorrow,

Thanks, that sounds like an OK solution. I don't see a problem having
separate wrappers for MKL; it's sort of a maintenance burden, but this
probably can't be helped as the ABIs don't sound like they are compatible.

The new ABI wrappers were obviously not tested against MKL, so we have
to do this now...

See  https://github.com/scipy/scipy/pull/2824

I decided to split our wrappers into the G77 ABI-specific part (returning complex) and the CBLAS/CLAPACK hacks specific to Accelerate. Not sure the complexity is worth it compared to simply have three cases (MKL on OS X, Accelerate, rest), but it may help integration with other libraries.

David

Best,
Pauli
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.12 (GNU/Linux)

iEYEARECAAYFAlIk2ZoACgkQ6BQxb7O0pWDUKwCfR35y4AS2zECzr2y7enwaPEbH
JUoAoKEG1ptACEiuShbhr3Tg/QtrpFSz
=irId
-----END PGP SIGNATURE-----

_______________________________________________
SciPy-Dev mailing list
SciPy-Dev@scipy.org
http://mail.scipy.org/mailman/listinfo/scipy-dev