[Numpy-discussion] NumPy-Discussion OpenBLAS and dotblas

Sturla Molden sturla.molden at gmail.com
Tue Aug 12 10:06:15 EDT 2014


Matti Picus <matti.picus at gmail.com> wrote:

> Thanks for your prompt reply. I think numpy is a wonderful project, and 
> you all do a great job moving it forward.
> If you ask what would my vision for maturing numpy, I would like to see 
> a grouping of linalg matrix-operation functionality into a python level 
> package, exactly the opposite of more tightly tying linalg into the core 
> of numpy. 

But with the @ operator in Python 3.5 it makes sence to have both matrix
multiplication and linear algebra solvers in the core of NumP. Just
consider:

A @ B
A.LazyInverse @ B 
A @ B.LazyInverse

(Lazy matrix inversion with O(1) complexity does not exist right now, but
could be added in the future.)

To implement this efficiently, we need BLAS and LAPACK in the core of
NumPy. It does not mean there would not be a linalg namespace for LU, SVD,
et al.

Sturla




More information about the NumPy-Discussion mailing list