"Klemm, Michael" <michael.kl...@intel.com> wrote: > I have found that the numpy.linalg.svd algorithm creates the resulting U, > sigma, and V matrixes with Fortran storage. Is there any way to force > these kind of algorithms to not change the storage order? That would > make passing the matrixes to the native dgemm operation much easier.
NumPy's dot function will call cblas_dgemm in the most efficient way regardless of storage. It knows what to do with C and Fortran arrays. Sturla _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion