概要
本サンプルはFortran言語によりLAPACKルーチンDGESVを利用するサンプルプログラムです。
:6x4の行列の特異値と左および右特異ベクトルを求めます。
計算された特異値と特異ベクトルの誤差限界近似値も合わせて求めます。
DGESDDの例題プログラムは の場合の特異値分解を示します。
入力データ
(本ルーチンの詳細はDGESV のマニュアルページを参照)1 2 3 4 5 6 7 8 9 10
このデータをダウンロード |
DGESV Example Program Data 4 :Value of N 1.80 2.88 2.05 -0.89 5.25 -2.95 -0.95 -3.80 1.58 -2.69 -2.90 -1.04 -1.11 -0.66 -0.59 0.80 :End of matrix A 9.52 24.35 0.77 -6.22 :End of vector b
出力結果
(本ルーチンの詳細はDGESV のマニュアルページを参照)1 2 3 4 5 6 7 8 9 10 11 12 13 14
この出力例をダウンロード |
DGESV Example Program Results Solution 1.0000 -1.0000 3.0000 -5.0000 Details of factorization 1 2 3 4 1 5.2500 -2.9500 -0.9500 -3.8000 2 0.3429 3.8914 2.3757 0.4129 3 0.3010 -0.4631 -1.5139 0.2948 4 -0.2114 -0.3299 0.0047 0.1314 Pivot indices 2 2 3 4
ソースコード
(本ルーチンの詳細はDGESV のマニュアルページを参照)※本サンプルソースコードのご利用手順は「サンプルのコンパイル及び実行方法」をご参照下さい。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71
このソースコードをダウンロード |
Program dgesv_example ! DGESV Example Program Text ! Copyright 2017, Numerical Algorithms Group Ltd. http://www.nag.com ! .. Use Statements .. Use lapack_example_aux, Only: nagf_file_print_matrix_real_gen Use lapack_interfaces, Only: dgesv Use lapack_precision, Only: dp ! .. Implicit None Statement .. Implicit None ! .. Parameters .. Integer, Parameter :: nin = 5, nout = 6 ! .. Local Scalars .. Integer :: i, ifail, info, lda, ldb, n ! .. Local Arrays .. Real (Kind=dp), Allocatable :: a(:, :), b(:) Integer, Allocatable :: ipiv(:) ! .. Executable Statements .. Write (nout, *) 'DGESV Example Program Results' Write (nout, *) ! Skip heading in data file Read (nin, *) Read (nin, *) n lda = n ldb = n Allocate (a(lda,n), b(ldb), ipiv(n)) ! Read A and B from data file Read (nin, *)(a(i,1:n), i=1, n) Read (nin, *) b(1:n) ! Solve the equations Ax = b for x Call dgesv(n, 1, a, lda, ipiv, b, ldb, info) If (info==0) Then ! Print solution Write (nout, *) 'Solution' Write (nout, 100) b(1:n) ! Print details of factorization Write (nout, *) Flush (nout) ! ifail: behaviour on error exit ! =0 for hard exit, =1 for quiet-soft, =-1 for noisy-soft ifail = 0 Call nagf_file_print_matrix_real_gen('General', ' ', n, n, a, lda, & 'Details of factorization', ifail) ! Print pivot indices Write (nout, *) Write (nout, *) 'Pivot indices' Write (nout, 110) ipiv(1:n) Else Write (nout, 120) 'The (', info, ',', info, ')', & ' element of the factor U is zero' End If 100 Format ((3X,7F11.4)) 110 Format ((3X,7I11)) 120 Format (1X, A, I3, A, I3, A, A) End Program