SciPy

scipy.linalg.find_best_blas_type

scipy.linalg.find_best_blas_type(arrays=(), dtype=None)[source]

Find best-matching BLAS/LAPACK type.

Arrays are used to determine the optimal prefix of BLAS routines.

Parameters:

arrays : sequence of ndarrays, optional

Arrays can be given to determine optimal prefix of BLAS routines. If not given, double-precision routines will be used, otherwise the most generic type in arrays will be used.

dtype : str or dtype, optional

Data-type specifier. Not used if arrays is non-empty.

Returns:

prefix : str

BLAS/LAPACK prefix character.

dtype : dtype

Inferred Numpy data type.

prefer_fortran : bool

Whether to prefer Fortran order routines over C order.

Examples

>>> import scipy.linalg.blas as bla
>>> a = np.random.rand(10,15)
>>> b = np.asfortranarray(a)  # Change the memory layout order
>>> bla.find_best_blas_type((a,))
('d', dtype('float64'), False)
>>> bla.find_best_blas_type((a*1j,))
('z', dtype('complex128'), False)
>>> bla.find_best_blas_type((b,))
('d', dtype('float64'), True)