[Numpy-discussion] inplace dot products

[David Warde-Farley]

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