Right now we store G, B in `TLMat`, but storing orthogonal bases separately would allow for more convenient operations.
Right now we store G, B in
TLMat, but storing orthogonal bases separately would allow for more convenient operations.