線形システム(Aはls_arnoldi2_projection.eq1.gif)

ls_arnoldi2_projection.eq2.gif

ls_arnoldi2_projection.eq3.gifステップ目の近似解ls_arnoldi2_projection.eq4.gifを初期値ls_arnoldi2_projection.eq5.gifから求める. いま,2つの部分空間ls_arnoldi2_projection.eq6.gifを考える. そして,

  • ls_arnoldi2_projection.eq7.gif
  • ls_arnoldi2_projection.eq8.gif

として,解を探索するのがProjection法である. Kは解候補が含まれるsearch subspaceで, Lは近似解を得るための制約空間(subspace of constraints)である. Projection法では残差ベクトルls_arnoldi2_projection.eq9.gifがLに対して垂直になるようにする. この条件をPetrov-Galerkin conditionという.

ここで,ls_arnoldi2_projection.eq10.gifとすると,

ls_arnoldi2_projection.eq11.gif

よって,

ls_arnoldi2_projection.eq12.gif

2つめの条件により,ls_arnoldi2_projection.eq13.gifはLに垂直なベクトルとなる. これが基本的なProjection法のステップである. また,K,Lはステップごとに適切なものを選ぶ必要がある.

Krylov部分空間によるProjection法

共役勾配法を代表的なものとして,多くの方法が Projection法において,K,LにKrylov部分空間を用いたものになっている. Krylov部分空間をls_arnoldi2_projection.eq14.gifとすると, これらの方法は大まかに以下のように分類される.

  • ls_arnoldi2_projection.eq15.gif : FOMなど
  • ls_arnoldi2_projection.eq16.gif : Lanczos, BiCG, CGS, QMRなど
  • ls_arnoldi2_projection.eq17.gif : GMRESなど

添付ファイル: filels_arnoldi2_projection.eq9.gif 687件 [詳細] filels_arnoldi2_projection.eq11.gif 689件 [詳細] filels_arnoldi2_projection.eq12.gif 713件 [詳細] filels_arnoldi2_projection.eq13.gif 725件 [詳細] filels_arnoldi2_projection.eq14.gif 806件 [詳細] filels_arnoldi2_projection.eq15.gif 699件 [詳細] filels_arnoldi2_projection.eq16.gif 684件 [詳細] filels_arnoldi2_projection.eq17.gif 819件 [詳細] filels_arnoldi2_projection.eq2.gif 656件 [詳細] filels_arnoldi2_projection.eq3.gif 709件 [詳細] filels_arnoldi2_projection.eq4.gif 702件 [詳細] filels_arnoldi2_projection.eq5.gif 692件 [詳細] filels_arnoldi2_projection.eq6.gif 679件 [詳細] filels_arnoldi2_projection.eq7.gif 694件 [詳細] filels_arnoldi2_projection.eq8.gif 843件 [詳細] filels_arnoldi2_projection.eq1.gif 827件 [詳細] filels_arnoldi2_projection.eq10.gif 808件 [詳細]

トップ   編集 凍結 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2024-03-08 (金) 18:06:05