ENH: output U matrix
ENH: u matrix output