コード例 #1
0
ファイル: qr.hpp プロジェクト: ibaned/omega_h
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;
}
コード例 #2
0
ファイル: eqnsys.cpp プロジェクト: FoxMarts/qucs
void eqnsys<nr_type_t>::solve_qr_ls (void) {
  A->transpose ();
  factorize_qr_householder ();
  substitute_qr_householder_ls ();
}
コード例 #3
0
ファイル: eqnsys.cpp プロジェクト: FoxMarts/qucs
void eqnsys<nr_type_t>::solve_qr (void) {
  factorize_qr_householder ();
  substitute_qr_householder ();
}