IOMにおいて問題となるのは,ls_arnoldi4_diom.eq1.gifの計算ですべてのls_arnoldi4_diom.eq2.gifが必要になることである. そこで,ls_arnoldi4_diom.eq3.gifls_arnoldi4_diom.eq4.gifからでなく,ls_arnoldi4_diom.eq5.gifから計算するように修正したのがDIOM(Direct IOM)である.

IOMでi=j-k+1〜jとしたことで,ls_arnoldi4_diom.eq6.gifは値を持つ部分の幅がk+1な行列になる.以下にm=5, k=3の場合の例を示す.

ls_arnoldi4_diom.eq7.gif

まずls_arnoldi4_diom.eq8.gifとLU分解する. ピボッティングがなければ,分解した行列は以下のようになる.

ls_arnoldi4_diom.eq9.gif

DIOMではls_arnoldi4_diom.eq5.gifからls_arnoldi4_diom.eq10.gifを求める形にする. FOMのls_arnoldi4_diom.eq3.gifに関する式は,

ls_arnoldi4_diom.eq11.gif

であり,ls_arnoldi4_diom.eq6.gifをLU分解した行列で置き換える.

ls_arnoldi4_diom.eq12.gif

いま,

ls_arnoldi4_diom.eq13.gif

とおくと以下のように変形できる.

ls_arnoldi4_diom.eq14.gif

ls_arnoldi4_diom.eq15.gifls_arnoldi4_diom.eq16.gifそれぞれの中身について考えてみる.

  • ls_arnoldi4_diom.eq15.gifについて
    ls_arnoldi4_diom.eq15.gifの式はls_arnoldi4_diom.eq17.gifとなり,ls_arnoldi4_diom.eq15.gifの各列をls_arnoldi4_diom.eq18.gifとすると,ls_arnoldi4_diom.eq19.gifが上三角行列であることから, ls_arnoldi4_diom.eq20.gifはガウスの消去法の後退代入を使って以下のように計算される.
    ls_arnoldi4_diom.eq21.gif
    よって,ls_arnoldi4_diom.eq20.gifはm-1までのls_arnoldi4_diom.eq18.gifを使って算出される.ここでls_arnoldi4_diom.eq15.gifを以下のように書くことにする.
    ls_arnoldi4_diom.eq22.gif
  • ls_arnoldi4_diom.eq16.gifについて
    ls_arnoldi4_diom.eq16.gifについてもls_arnoldi4_diom.eq23.gifであり,例えば,m=5の場合,
    ls_arnoldi4_diom.eq24.gif
    よって,
    ls_arnoldi4_diom.eq25.gif
    ただし,ls_arnoldi4_diom.eq26.gifである. このようにls_arnoldi4_diom.eq27.gifもm-1以下のls_arnoldi4_diom.eq28.gifから計算される.こちらもls_arnoldi4_diom.eq15.gifの時と同様に以下のように書く.
    ls_arnoldi4_diom.eq29.gif

これらのことを使って式を書き換えると,

ls_arnoldi4_diom.eq30.gif

ここで,ls_arnoldi4_diom.eq31.gifls_arnoldi4_diom.eq5.gifなので,

ls_arnoldi4_diom.eq32.gif

添付ファイル: filels_arnoldi4_diom.eq23.gif 545件 [詳細] filels_arnoldi4_diom.eq13.gif 577件 [詳細] filels_arnoldi4_diom.eq12.gif 526件 [詳細] filels_arnoldi4_diom.eq29.gif 537件 [詳細] filels_arnoldi4_diom.eq14.gif 541件 [詳細] filels_arnoldi4_diom.eq26.gif 524件 [詳細] filels_arnoldi4_diom.eq1.gif 483件 [詳細] filels_arnoldi4_diom.eq18.gif 548件 [詳細] filels_arnoldi4_diom.eq7.gif 512件 [詳細] filels_arnoldi4_diom.eq16.gif 526件 [詳細] filels_arnoldi4_diom.eq24.gif 531件 [詳細] filels_arnoldi4_diom.eq3.gif 449件 [詳細] filels_arnoldi4_diom.eq15.gif 523件 [詳細] filels_arnoldi4_diom.eq11.gif 495件 [詳細] filels_arnoldi4_diom.eq19.gif 521件 [詳細] filels_arnoldi4_diom.eq22.gif 484件 [詳細] filels_arnoldi4_diom.eq30.gif 517件 [詳細] filels_arnoldi4_diom.eq31.gif 473件 [詳細] filels_arnoldi4_diom.eq32.gif 487件 [詳細] filels_arnoldi4_diom.eq9.gif 526件 [詳細] filels_arnoldi4_diom.eq28.gif 512件 [詳細] filels_arnoldi4_diom.eq10.gif 455件 [詳細] filels_arnoldi4_diom.eq2.gif 510件 [詳細] filels_arnoldi4_diom.eq25.gif 486件 [詳細] filels_arnoldi4_diom.eq20.gif 489件 [詳細] filels_arnoldi4_diom.eq27.gif 514件 [詳細] filels_arnoldi4_diom.eq21.gif 457件 [詳細] filels_arnoldi4_diom.eq4.gif 439件 [詳細] filels_arnoldi4_diom.eq8.gif 477件 [詳細] filels_arnoldi4_diom.eq5.gif 536件 [詳細] filels_arnoldi4_diom.eq17.gif 493件 [詳細] filels_arnoldi4_diom.eq6.gif 499件 [詳細]

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2012-07-13 (金) 14:54:11 (2994d)