Testing the output

> evalm(htranspose(T)&*F&*T):map(evala,%); The columns of T are F-ortonormal