Function rgsl::linear_algebra::bidiag_unpack

pub fn bidiag_unpack(a: &MatrixF64, tau_u: &VectorF64, u: &MatrixF64, tau_v: &VectorF64, v: &MatrixF64, diag: &VectorF64, superdiag: &VectorF64) -> Value

This function unpacks the bidiagonal decomposition of A produced by gsl_linalg_bidiag_decomp, (A, tau_U, tau_V) into the separate orthogonal matrices U, V and the diagonal vector diag and superdiagonal superdiag. Note that U is stored as a compact M-by-N orthogonal matrix satisfying U^T U = I for efficiency.