Thanks for the quick reply. A was not aware of the fact that this kind of "fuse" function is not compatible with Numpy API. I understand the point. FYI, numba is not able to simplify this kind of calculus: C += A @ B. Nor numexpr which is not compatible with dot product. I did not test pythran.
Thus, the only solution is to use the Blas functions through scipy, as recalled by Kevin. I'll play a bit with transposition and alignment issues ... Regards, -- Sent from: http://numpy-discussion.10968.n7.nabble.com/ _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@python.org https://mail.python.org/mailman/listinfo/numpy-discussion