Matrix multiply ignores strides

[mateogianolio]

Reproducible example:

const x = range(0, 9).reshape(3, 3);

x.multiply(x)
/* array([
  [ 15, 18, 21 ],
  [ 42, 54, 66 ],
  [ 69, 90, 111 ]
], dtype=float64) */

x.T.multiply(x)
/* array([
  [ 15, 18, 21 ],
  [ 42, 54, 66 ],
  [ 69, 90, 111 ]
], dtype=float64)
... should be ...
array([
  [ 45, 54, 63 ],
  [ 54, 66, 78 ],
  [ 63, 78, 93 ]
], dtype=float64) */

Temporary workaround:

x.T.copy().multiply(x)

Proposed fix is to call copy() in multiply to normalise strides before performing the operation.