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 ();
}
