INLINE Vector<max_n> solve_using_qr( Int m, Int n, Matrix<max_m, max_n> a, Vector<max_m> b) { auto qr = factorize_qr_householder(m, n, a); auto qtb = implicit_q_trans_b(m, n, qr.v, b); auto x = solve_upper_triangular(n, qr.r, qtb); return x; }
void eqnsys<nr_type_t>::solve_qr_ls (void) { A->transpose (); factorize_qr_householder (); substitute_qr_householder_ls (); }
void eqnsys<nr_type_t>::solve_qr (void) { factorize_qr_householder (); substitute_qr_householder (); }