[Numpy-discussion] inplace dot products
Mon Mar 16 03:11:20 CDT 2009
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
> 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.
More information about the Numpy-discussion