Matrizes de SWAP

Bom, estou trabalhando com um zilhão dessas, então preciso de uma função para criá-las.

qiRowSwap[dimension_, i_, j_] :=
IdentityMatrix[dimension] -
qiOrthoMatrix[dimension, i] -
qiOrthoMatrix[dimension, j] +
KroneckerProduct[UnitVector[dimension, i],
UnitVector[dimension, j]] +
KroneckerProduct[UnitVector[dimension, j], UnitVector[dimension, i]]