Dgelqf Dgeqp3 Dgeqrf Dorglq Dorgqr Dormlq Dormqr

Dgelqf() Sub Dgelqf ( M As  Long, N As  Long, A() As  Double, Tau() As  Double, Info As  Long  ) LQ分解 目的 本ルーチンはm×n行列AのLQ分解を計算する. A = L * Q 引数 [in] M 行列 A の行数. (M >= 0) (M = 0 の場合, 処理を行わずに戻る) [in] N 行列 A の列数. (N >= 0) (N = 0 の場合, 処理を行わずに戻る) [in,out] A() 配列 A(LA1 - 1, LA2 - 1) (LA1 >= M, LA2 >= N) [in] M×N行列 A. [out] 対角およびその下の要素にM×min(M, N)下台形行列Lが入る(M <= Nであれば, Lは下三角行列である). 対角より上の要素は, 配列Tau()と合わせて, 基本鏡映変換の積として直交行列Qを表す (詳細を参照のこと). [out] Tau() 配列 Tau(LTau - 1) (LTau >= min(M, N)) 基本鏡映変換のスカラー因子 (詳細を参照のこと). [out] Info = 0: 正常終了. = -1: パラメータ M の誤り. (M < 0) = -2: パラメータ N の誤り. (N < 0) = -3: パラメータ A() の誤り. = -4: パラメータ Tau() の誤り. 詳細 行列Qは基本鏡映変換の積で表される. Q = H(k) . . . H(2) H(1), ただし k = min(M, N). 各H(i)は次のように表される. H(i) = I - tau * v * v^T ただし, tauは実数スカラー, また, vは実数ベクトルで, v(1〜i-1) = 0, v(i) = 1 である. v(i+1〜N)はA(i-1, i〜N-1)に, tauはTau(i-1)に格納される. 出典 LAPACK 使用例 行列Aの行列AのLQ分解を求める. ただし, ( 0.20 -0.11 -0.93 ) A = ( -0.32 0.81 0.37 ) ( -0.80 -0.92 -0.29 ) である. Sub Ex_Dgelqf() Const M = 3, N = 3, K = N Dim A(M - 1, N - 1) As Double, Tau(N - 1) As Double, Info As Long A(0, 0) = 0.2: A(0, 1) = -0.11: A(0, 2) = -0.93 A(1, 0) = -0.32: A(1, 1) = 0.81: A(1, 2) = 0.37 A(2, 0) = -0.8: A(2, 1) = -0.92: A(2, 2) = -0.29 '-- Compute LQ factorization of A Call Dgelqf(M, N, A(), Tau(), Info) Debug.Print "L =" Debug.Print A(0, 0) Debug.Print A(1, 0), A(1, 1) Debug.Print A(2, 0), A(2, 1), A(2, 2) Debug.Print "Info =", Info '-- Compute Q Call Dorglq(M, N, K, A(), Tau(), Info) Debug.Print "Q =" Debug.Print A(0, 0), A(0, 1), A(0, 2) Debug.Print A(1, 0), A(1, 1), A(1, 2) Debug.Print A(2, 0), A(2, 1), A(2, 2) Debug.Print "Info =", Info End Sub 実行結果 L = -0.957601169589929 0.51921406926948 -0.791085804620857 -0.220237826244838 0.609478035393939 -1.07262846515618 Info = 0 Q = -0.208855216922558 0.114870369307407 0.971176758689895 0.267429186186539 -0.948516183365263 0.169701739147835 0.940670573973608 0.295164103761576 0.167382863850482 Info = 0