• 研究課題をさがす
  • 研究者をさがす
  • KAKENの使い方
  1. 課題ページに戻る

2011 年度 実績報告書

大規模線形方程式系に対する新数値解法,GBi-CGSTAB法の研究

研究課題

研究課題/領域番号 22560060
研究機関東京大学

研究代表者

杉原 正顯  東京大学, 情報理工学(系)研究科, 教授 (80154483)

キーワード線形方程式系の数値解法 / クリロフ部分空間法 / IDR(s)法 / 偽収束 / シフト方程式
研究概要

谷尾,杉原によって開発された大規模線形方程式系の数値解法GBi-CGSTAB(s,L) 法に関して以下のような研究成果を得た.
[1] 偽収束に関する新しい知見の発見とそれを用いた偽収束を改善するための手法のGBi-CGSTAB法への組込み: 昨年,偽収束現象を改善する手法をIDR(s) 法の場合に倣って開発した.その手法を多くの例に適用する中,残差と真の残差の差のノルムは反復初期からほとんど変わらず一定となっていることを発見した.そこでこれを利用した偽収束を改善する自動残差修正機構を構築し,GBi-CGSTAB法に組込んだ.そして,昨年度開発した自動残差修正機構の場合と比較し,同等以上の有効性をもつことを示した.
[2] GBi-CGSTAB(s,L) 法に基づく多数のシフト方程式を効率よく解くため解法の確立: 多数のシフト方程式効率よく解くためのクリロフ部分空間法として,Shifted GMRES法,Shifted BiCGstab(L)法,Shifted IDR(s)法など,様々なものが提案されている.ここではまず手法の基礎になる残差の共線性をGBiCG(s) 法の場合に確立し,GBiCG(s) 法によるシフト線形方程式の解法を与えた.そして,これにshifted BiCGSTAB(L) 法に倣いながら,安定化多項式を付加し,GBi-CGSTAB(s,L) 法を基にしたシフト線形方程式の解法を構築した.数値実験により,既存手法と比較し,同等以上の性能を持つことを確かめた.特に,歪対称に近い係数行列を持つ,IDR(s) 法をベースとした手法では解けなかった問題に対して,有効であることを確かめた.

現在までの達成度 (区分)
現在までの達成度 (区分)

2: おおむね順調に進展している

理由

研究の目標である:偽収束を改善する手法の開発,多数のシフト方程式を効率よく解くため解法の確立が達成され,研究は順調に進展している.

今後の研究の推進方策

偽収束に関する新しい知見(残差と真の残差の差のノルムは反復初期からほとんど変わらず一定となっていること)は,GBi-CGSTAB法のみならず,クリロフ部分空間法に共通する性質のようであり,これを解明することはこの分野におけるブレークスルーとなることが期待される.是非,この方向での研究を進めたい.

  • 研究成果

    (2件)

すべて その他

すべて 学会発表 (2件)

  • [学会発表] GBiCGSTAB (s,L) 法の偽収束性について

    • 著者名/発表者名
      深堀 康紀,杉原 正顯
    • 学会等名
      日本応用数理学会2011年度年会
    • 発表場所
      同志社大学
  • [学会発表] シフト方程式に対するGBiCGSTAB(s,L)法

    • 著者名/発表者名
      深堀 康紀,杉原 正顯
    • 学会等名
      日本応用数理学会 2012年 研究部会 連合発表会
    • 発表場所
      九州大学

URL: 

公開日: 2014-07-24  

サービス概要 検索マニュアル よくある質問 お知らせ 利用規程 科研費による研究の帰属

Powered by NII kakenhi