[Numpy-discussion] Patch for accelerate cblas_sgemv (Re: Reporting a bug to Apple)

On 10/06/14 14:57, Matthew Brett wrote:

Would you consider doing a PR for that?

Here is a patch you can try before I post a PR. It can also be build independently of NumPy, so you don't need to rebuild NumPy just for testing it. (The change to numpy is in a different folder.) I decided against using cblas_sgemm. Instead it just enforces alignment to 32 byte boundaries. Because cblas_sgemm would require a copy if the vector is strided, it didn't matter. I have tested with Accelerate, OpenBLAS and MKL, clang and icc. From what I can tell it works correctly and does not segfault on misalignment. Sturla