PDEを解くためのマルチグリッド法


15

マルチグリッド法の簡単な説明またはこれに関するいくつかの文献が必要です。

私は、BiCGStab、CG、GS、Jacobiおよび事前調整などの反復法に精通していますが、マルチグリッド法の初心者です。

非常に初心者向けの優れた文献であっても、誰かがこれを詳細に説明したり、少なくとも明確に擬似コードまたはソースコードを提供したりできますか?ありがとう!

回答:


15

マルチグリッドの背後にある主なアイデアは、投影です。私はそれについて次のように考えようとします:

かなりの精度でPDEを解きたいと思うので、たくさんの点がたくさんある非常に細かいグリッドでドメインを離散化することにしましょう(たとえば、有限差分法を使用して)。最後に、方程式系をセットアップし、それを解決する準備ができました。私はお気に入りの反復ソルバー(jacobi、gauss seidel、共役勾配など)を使用してみます。私は1日以上待って、コンピューターがまだ答えを計算しようとしていることに気づきました!!!

これらの反復法がすぐに動作しない理由は、(通常)このような方程式の大規模なシステムをセットアップすると、行列自体の固有値が1に非常に近いためです。多くの反復法の収束速度は最大の固有値に反比例するためです(27ページのブリッグのマルチグリッドチュートリアルスライド、パート1へのクリスチャンクラソンのリンクを参照)。そのため、最大固有値が1に近いほど、反復法は遅くなります。(注:これは少し単純化しすぎていますが、マルチグリッドの必要性を動機付けるのに役立ちます)。

明らかに、未知数が少ない場合(つまり、グリッドポイントの少ない粗いグリッド上)に問題を解決する方が常に高速です。しかし、もっと重要なのは、より粗いグリッドの解(または近似解)が、より細かいグリッドの問題を解決するための良い出発点です。 これは、ほとんどの(すべてではないにしても)マルチグリッドメソッドの背後にある重要な考え方です。これはなぜですか?直感的には理にかなっていますが、これを正当化する数学的に厳密な方法があります。

元の細かいグリッドの問題に適用された反復法(引数のために、ジャコビまたはガウスサイデル)のエラーのフーリエモードを見てみましょう。最初の数回の反復で、ほとんどの高周波(高度に振動する)エラーが除去されることがわかります。これは素晴らしいことですが、低周波(振動が少ない)エラーが残っており、すぐには消えません。実際、標準の反復法が迅速に収束するのを妨げるのは低周波誤差です。

粗いグリッド(たとえば、jacobiやgauss-seidelなどの反復法)で問題を解決するとき、基本的には、細かいグリッドよりもはるかに迅速に(つまり、より少ない反復で)低周波数エラーを除去できます。したがって、粗いグリッドの問題を解決する場合、より低い周波数誤差が大幅に減少したソリューションがあります。したがって、より細かいグリッドでの反復法の開始点として役立ちます。

さまざまなマルチグリッド方式がありますが、それらのほとんどは以下のいくつかのバリエーションによって動作します。

  1. 細かいグリッドの問題から始めます
  2. 粗いグリッドに投影する(制限とも呼ばれます
  3. 粗いグリッド上の解を近似します(他のソルバーを使用)
  4. 粗いグリッドソリューションをより細かいグリッドに投影します(延長とも呼ばれます)。
  5. 4.からの投影を初期推測として使用し、反復法により細かいグリッドの問題を解きます。

私にとって、マルチグリッド方式の最も難しい部分は、グリッド間の投影です。@ChristianClasonによって提案されたBriggsのチュートリアルは、このテーマを私よりもはるかにうまく処理します。


詳細な回答ありがとうございます!これで、multigrdiメソッドに関する基本的な知識が得られました。今、私は制限と延長プロセスについて具体的な質問があります。いくつかの制限行列Rと補間行列M、およびA2 = RAIのような式がグリッド間の投影を行うために使用されたことを読みました。しかし、行列RとIを構築する方法を理解していませんでした。これに関するアイデアはありますか。
ヌルラン

ChristianClasonが投稿したブリッグスマルチグリッドチュートリアル(パート1)のスライド45〜57をご覧ください。そこで、ブリッグスは幾何学的マルチグリッド法のプロセスを説明しています。それは私が見つけることができる最も簡単な説明です。さらに質問がある場合は、新しい質問を投稿してください。:)
ポール

15

このサイトは、擬似コードで詳細な説明を求めるのに適した場所ではない可能性があります(FAQに記載されているように、「質問に答える本全体を想像できるなら、あまりにも多くを求めています。」)このトピックに関する古典的な本の1つ(以下にリスト)から始めて、問題がある具体的な詳細についての特定の質問に戻ってください。


2
ブリッグスは本当に良いです!
vanCompute

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.