[Numpy-discussion] How to Force Storage Order

Sturla Molden sturla.molden at gmail.com
Wed Apr 1 17:12:21 EDT 2015


"Klemm, Michael" <michael.klemm at 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




More information about the NumPy-Discussion mailing list