並列ODEメソッドの最新技術とは何ですか?


39

現在、ODE統合のための並列メソッドを検討しています。さまざまなアプローチを説明する新しい文献や古い文献がたくさんありますが、最近の調査や一般的なトピックを説明する概要記事は見つかりませんでした。

Burrageによる本[1]がありますが、それはほぼ20年前であり、したがって、パラリアルアルゴリズムのようなより現代的なアイデアの多くをカバーしていません。

[1] K.バラージ、常微分方程式の並列法および逐次法、クラレンドンプレス、オックスフォード、1995

回答:


35

私は最近の概要記事を知りませんが、PFASSTアルゴリズムの開発に積極的に関与しているので、いくつかの考えを共有できます。

私が知っている時系列技術には3つの広範なクラスがあります。

  • メソッド全体— RKまたは外挿積分器の独立した段階を並行して評価できます。RIDC(修正主義の積分遅延補正アルゴリズム)も参照してください。
  • 問題全体—波形緩和
  • 時間領域全体—パラリアル。PITA(時間並列アルゴリズム); およびPFASST(空間と時間の並列完全近似スキーム)。

メソッド全体で並列化するメソッドは、通常、仕様に非常に近いパフォーマンスを発揮しますが、少数の(時間)プロセッサーを超えてスケ​​ーリングすることはありません。通常、それらは他の方法よりも実装が比較的簡単であり、いくつかのコアを追加して、予測可能で控えめな高速化を探している場合に適しています。

時間領域で並列化するメソッドには、Parareal、PITA、PFASSTが含まれます。これらの方法はすべて反復的で、安価な(ただし不正確な)「粗い」プロパゲーターと高価な(しかし正確な)「細かい」プロパゲーターで構成されています。粗いプロパゲーターを使用して得られたシリアルソリューションを改善するために、細かいプロパゲーターを並列で繰り返し評価することにより、並列効率を達成します。

EE<1/KK

これらのすべての方法で多くのゲームをプレイして、それらを試して速度を上げることができます。これらのドメイン全体のテクニックのパフォーマンスは、解決している問題と粗いものをスピードアップするために利用可能なテクニックに依存しているようです伝播者(粗いグリッド、粗い演算子、粗い物理など)。

いくつかの参考文献(論文に記載されている参考文献も参照):

ネットで利用可能なPFASSTの2つの実装を作成しました:PyPFASSTlibpfasst


1
私は現在、パラリアルを学んでいます。そして、それは私にとって大きな助けになると思います。
eccstartup

これは素晴らしい概要です。ただし、ODEはPDEの空間的離散化後に解決されることが多いことを明示的に言及する必要があります。したがって、空間ドメインが十分に大きい場合、メソッド全体の並列処理により、数千のコアに対する優れたスケーラビリティが得られます。これは、計算時間の大半が、たとえばRKステージのRHS評価の計算に費やされるためです。
NoseKnowsAll

15

この投稿は現在2年ですが、誰かがつまずいた場合に備えて、簡単な更新をさせてください。

Martin Ganderは最近、この分野の歴史的な視点を提供し、多くの異なるPINTメソッドについて議論する素晴らしいレビュー記事を書きました:http : //www.unige.ch/~gander/Preprints/50YearsTimeParallel.pdf

また、非常に多くのリファレンスをリストし、さまざまな方法の説明を提供するコミュニティWebサイトもあります。http//www.parallel-in-time.org/

特にPararealの並列時間アルゴリズムについては、https//en.wikipedia.org/wiki/Pararealで説明しています。


1
特に素晴らしいソフトウェア(XBraid)に支えられているように見えるため、GanderがFalgoutなどによるMGRITアプローチについて話していないことに少し驚いていますが、MGRITの論文が最近出てきたのを知っています。
ジェフオックスベリー

1
こんにちはGeoff、MGRITの論文が発表される前にMartin Ganderが論文を書いたと確信しています。レビュー論文は2015年に登場する予定ですが、プレプリントは2013年後半にすでにオンラインになったと思います。
ダニエル

1
一見すると、このレビューでは「メソッド全体で並列」が省略されているように見えます。たとえば、外挿は決して言及されていません。
デビッドケッチャソン16

4

u0u(t)=exp(λt)u0, Reλ>0.


私が言ったように、私はすでに個々のトピックに関する多くの記事を見つけました。私が欠けているのは、アプローチの一般的な概要です。
フロリアンブラッカー

1
FWIW、PFASSTアルゴリズムは、多数(数百)のタイムプロセッサであっても、ハミルトニアンシステムに対して非常に優れた収束(近日公開予定)を示します。そうは言っても、かなりのスピードアップを得るには、もう一度、粗いプロパゲーターを細かいプロパゲーターよりもはるかに安くすることに依存します-粒子系の高速化を得るには、多重極展開または他のマルチフィジックスアプローチが必要と思われます。
マシューエメット
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.