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 _______________________________________________ Numpy-discussion mailing list [email protected] http://mail.scipy.org/mailman/listinfo/numpy-discussion
