Testing the output
>
evalm(transpose(S)&*F&*S);