Orthogonal SVD — orthonormalize named axes via QR, fill the rest from SVD
Source:R/RcppExports.R
orthogonal_svd.RdMirrors rENA's orthogonal_svd() in ena.rotate.by.mean.R:
the named axes in the output are the orthonormalized Q columns, not the
original weights columns. Use complete_rotation to keep
the named axes verbatim.