Multiplying transposed matrices in mir

[jmh530]

On Sunday, 19 April 2020 at 17:55:06 UTC, p.shkadzko wrote:
> snip
>
> So, lubeck mtimes is equivalent to NumPy "a.dot(a.transpose())".

There are elementwise operation on two matrices of the same size 
and then there is matrix multiplication. Two different things. 
You had initially said using an mxn matrix to do the calculation. 
Elementwise multiplication only works for matrices of the same 
size, which is only true in your transpose case when they are 
square. The mtimes function is like dot or @ in python and does 
real matrix multiplication, which works for generic mxn matrices. 
If you want elementwise multiplication of a square matrix and 
it’s transpose in mir, then I believe you need to call 
assumeContiguous after transposed.