Charles R Harris <charlesr.har...@gmail.com> wrote: > In any case, we support the `@` operator in 1.10, but not the `@=` > operator. The `@=` operator is tricky to have true inplace matrix > multiplication, as not only are elements on the left overwritten, but the > dimensions need to be preserved.
As long as we use BLAS, we can never have true inplace matrix multiplication because Fortran prohibits pointer aliasing. We therefore need to make a temporary copy before we call BLAS. But as for preserving dimensions: This should be allowed if the array is square. Sturla _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org https://mail.scipy.org/mailman/listinfo/numpy-discussion
