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