On Thu, May 21, 2015 at 7:06 PM, Alexander Belopolsky <ndarray@mac.com> wrote:
1. Is there a simple expression using existing numpy functions that implements PEP 465 semantics for @?

2. Suppose I have a function that takes two vectors x and y, and a matrix M and returns x.dot(M.dot(y)).  I would like to "vectorize" this function so that it works with x and y of any ndim >= 1 and M of any ndim >= 2 treating multi-dimensional x and y as arrays of vectors and M as an array of matrices (broadcasting as necessary).  The result should be an array of xMy products.  How would I achieve that using  PEP 465's @?


If you are willing to run Python 3.5 (use 3.6.0a3, a4 crawls with the bugs), you can use gh-5878. The override mechanisms are still in process in Nathaniel's PR, so that may change. I'd welcome any feedback.

Chuck