5 その他のサンプル

ナビゲーション:前へ   上へ   次へ

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     )



ナビゲーション:前へ   上へ   次へ

関連情報
Privacy Policy  /  Trademarks