Finds the unique commutation matrix K: K vec(A) = vec(A^t)
commutation(m, n = m)
m
rows
n
columns
Klaus K. Holst