I wonder if a whole set of infix matrix operations could be added to reduce the need to load numpy for simple tasks.
A @* B A @. B A @+ B