Matti Picus email@example.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.