Function Reference
— Function File: commutation_matrix (m, n)

Return the commutation matrix which is the unique matrix such that for all matrices

If only one argument m is given, is returned.

See Magnus and Neudecker (1988), Matrix differential calculus with applications in statistics and econometrics.