what is the fastest way to do this
(a.transpose() @ b.transpose()).transpose()
# a and b are 4x4 numpy matrices
dot gives a different result
import numpy as np
a = np.array((
(1, 2, 3, 4),
(0, 1, 0, 0),
(0, 0, 1, 0),
(0, 0, 0, 1),
))
b = np.array((
(2, 0, 7, 0),
(0, 2, 0, 0),
(0, 2, 2, 0),
(0, 2, 0, 2),
))
print(a @ b)
print(np.dot(a, b))
print((a.transpose() @ b.transpose()).transpose())