Testing the output

> evalm(htranspose(S)&*S);