Dgelqf Dgeqp3 Dgeqrf Dorglq Dorgqr Dormlq Dormqr

Dgeqp3() Sub Dgeqp3 ( M As  Long, N As  Long, A() As  Double, Jpvt() As  Long, Tau() As  Double, Info As  Long  ) QR分解 (ピボット選択付き) 目的 本ルーチンは行列AのQR分解をピボット選択付きで計算する. A * P = Q * R 計算にはレベル3 BLASを用いる. 引数 [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] 上三角部分にmin(M, N)×N上台形行列Rが入る. 対角下の要素は, 配列Tau()と合わせて, min(M, N)個の基本鏡映変換の積として直交行列Qを表す. [in,out] Jpvt() 配列 Jpvt(LJpvt - 1) (LJpvt >= N) [in] Jpvt(j-1) <> 0 の場合, Aの第j列はA*Pの前(先頭列)に入れ替えられる. Jpvt(j-1) = 0 の場合, Aの第j列はフリーな列である. [out] Jpvt(j-1) = k の場合, A*Pの第j列は第k列であったことを示す. [out] Tau() 配列 Tau(LTau - 1) (LTau >= min(M, N)) 基本鏡映変換のスカラー因子. [out] Info = 0: 正常終了. = -1: パラメータ M の誤り. (M < 0) = -2: パラメータ N の誤り. (N < 0) = -3: パラメータ A() の誤り. = -4: パラメータ Jpvt() の誤り. = -5: パラメータ Tau() の誤り. 出典 LAPACK 使用例 行列Aの行列AのQR分解を求める. ただし, ( 0.20 -0.11 -0.93 ) A = ( -0.32 0.81 0.37 ) ( -0.80 -0.92 -0.29 ) である. Sub Ex_Dgeqp3() Const M = 3, N = 3, K = N Dim A(M - 1, N - 1) As Double, Jpvt(N - 1) As Long, Tau(N - 1) As Double Dim Info As Long, I 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 For I = 0 To N - 1 Jpvt(I) = 1 Next '-- Compute QR factorization of A Call Dgeqp3(M, N, A(), Jpvt(), Tau(), Info) Debug.Print "R =" Debug.Print A(0, 0), A(0, 1), A(0, 2) Debug.Print A(1, 1), A(1, 2) Debug.Print A(2, 2) Debug.Print "Jpvt =", Jpvt(0), Jpvt(1), Jpvt(2) Debug.Print "Info =", Info '-- Compute Q Call Dorgqr(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 実行結果 R = -0.884533775499839 -0.514169173181655 8.18510293279503E-02 -1.11813687057967 -0.635776616242947 -0.821576839520579 Jpvt = 1 2 3 Info = 0 Q = -0.226107815823067 0.202352390539057 0.952852961199602 0.361772505316908 -0.890778487093709 0.275016983707865 0.904431263292269 0.406899492472248 0.128206446816755 Info = 0