Testing the output

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