Testing the output

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