[Numpy-discussion] Optimisation of matrix multiplication

[Andrew Nelson]

Hi all,
I have a matrix multiplication that I'd like to optimize.

I have a matrix `a` (dtype=complex) with shape (N, M, 2, 2). I'd like to do
the following multiplication:

a[:, 0] @ a[:, 1] @ ... @ a[:, M-1]

where the first dimension, N, is element wise (and hopefully vectorisable)
and M>=2. So for each row do M-1 matrix multiplications of 2x2 matrices.
The output should have shape (N, 2, 2).

What would be the best way of going about this?

-- 
_____________________________________
Dr. Andrew Nelson


_____________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20180130/5336e536/attachment.html>