March 16, 2009
1:11 a.m.
On 20-Feb-09, at 6:41 AM, Olivier Grisel wrote:
Alright, thanks for the reply.
Is there a canonical way /sample code to gain low level access to blas / lapack atlas routines using ctypes from numpy / scipy code?
I don't mind fixing the dimensions and the ndtype of my array if it can decrease the memory overhead.
I got some clarification from Pearu Peterson off-list. For gemm the issue is that if the matrix C is not Fortran-ordered, it will be copied, and that copy will be over-written. order='F' when creating the array being overwritten will fix this. DWF