int gsl_linalg_hessenberg_unpack(gsl_matrix * H, gsl_vector * tau, gsl_matrix * U) { int s; gsl_matrix_set_identity(U); s = gsl_linalg_hessenberg_unpack_accum(H, tau, U); return s; } /* gsl_linalg_hessenberg_unpack() */
/** * C++ version of gsl_linalg_hessenberg_unpack_accum(). * @param H An upper Hessenberg matrix * @param tau A vector * @param U An orthogonal matrix * @return Error code on failure */ inline int hessenberg_unpack_accum( matrix& H, vector& tau, matrix& U ){ return gsl_linalg_hessenberg_unpack_accum( H.get(), tau.get(), U.get() ); }