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.
