Function bidiag_unpack (o2scl_linalg)¶
-
template<class mat_t, class vec_t, class mat2_t, class vec2_t, class mat3_t, class vec3_t, class vec4_t>
int o2scl_linalg::bidiag_unpack(size_t M, size_t N, const mat_t &A, const vec_t &tau_U, mat2_t &U, const vec2_t &tau_V, mat3_t &V, vec3_t &diag, vec4_t &superdiag)¶ Unpack a matrix
Awith the bidiagonal decomposition and create matricesU,V, diagonaldiagand superdiagonalsuperdiag.Given a matrix
Aof size(M,N)with \( M \geq N \) created by bidiag_decomp(), this function creates the matrixUof size(M,N), the matrixVof size(N,N), the diagonaldiagof sizeNand the super-diagonalsuperdiagof sizeN-1.