Testing the output

> evalm(htranspose(S)&*F&*S); The columns of S are F-ortogonal