5 その他のサンプル
5.1 Harwell Boeing形式(実対称行列、反復法)のサンプル
このサンプルはHarwell Boeing形式の実対称行列 (MXTYPE=RSA) を読み込みnAGソルバールーチンF11JEFを用いて解くサンプルです。
プリコンディショナーは使用せずにCG法で解いています。
※他のプリコンディショナーや解法を指定することも可能です。詳細はマニュアルページをご参照下さい
※以下のサンプルはnAG数値計算ライブラリ(トライアルのご案内 )を利用します。
サンプルコード: solve_real_symm.f90
係数行列ファイルサンプル: f11jef_matrix.hb
右辺ベクトルファイルサンプル:f11jef_uhen.txt
以下に実行結果例を示します。
Enter Harwell Boeing Matrix File Name f11jef_matrix.hb 7 x 7 (16 non-zeros) mxtype is RSA Enter Vector (txt) File Name f11jef_uhen.txt Converting to CS(SCS)... Sorting CS(SCS)... Solving... Results: number of iterations = 7 1.00000 2.00000 3.00000 4.00000 5.00000 6.00000 7.00000
5.2 Harwell Boeing形式(実非対称行列、反復法)のサンプル
このサンプルはHarwell Boeing形式の実非対称行列 (MXTYPE=RUA) を読み込みnAGソルバールーチンF11DEFを用いて解くサンプルです。
プリコンディショナー利用せずbi-conjugate gradient stabilized法で解いています。
※他のプリコンディショナーや解法を指定することも可能です。詳細はマニュアルページをご参照下さい
※以下のサンプルはnAG数値計算ライブラリ(トライアルのご案内 )を利用します。
サンプルコード: solve_real_nonsymm_itr.f90
係数行列ファイルサンプル: f11def_matrix.hb
右辺ベクトルファイルサンプル:f11def_uhen.txt
以下に実行結果例を示します。
Enter Harwell Boeing Matrix File Name f11def_matrix.hb 5 x 5 (16 non-zeros) mxtype is RUA Enter Vector (txt) File Name f11def_uhen.txt Converting to CS(SCS)... Sorting CS(SCS)... Solving... Results: number of iterations = 6 1.00000 2.00000 3.00000 4.00000 5.00000
5.3 Harwell Boeing形式(複素非エルミート行列、反復法)のサンプル
このサンプルはHarwell Boeing形式の複素非エルミート行列 (MXTYPE=CUA) を読み込みnAGソルバールーチンF11DSFを用いて解くサンプルです。
プリコンディショナーは使用せずにbi-conjugate gradient stabilized法で解いています。
※他のプリコンディショナーや解法を指定することも可能です。詳細はマニュアルページをご参照下さい
※以下のサンプルはnAG数値計算ライブラリ(トライアルのご案内 )を利用します。
サンプルコード: solve_complex_nonhermitian.f90
係数行列ファイルサンプル: f11dsf_matrix.hb
右辺ベクトルファイルサンプル:f11dsf_uhen.txt
以下に実行結果例を示します。
Enter Harwell Boeing Matrix File Name f11dsf_matrix.hb 5 x 5 (16 non-zeros) mxtype is CUA Enter Vector (txt) File Name f11dsf_uhen.txt Converting to CS(SCS)... Sorting CS(SCS)... Solving... Results: number of iterations = 5 ( 1.00000 , 2.00000 ) ( 2.00000 , 3.00000 ) ( 3.00000 , 4.00000 ) ( 4.00000 , 5.00000 ) ( 5.00000 , 6.00000 )
ナビゲーション:前へ 上へ 次へ