[論文レビュー] An Integral Equation Formulation of the $N$-body Dielectric Spheres Problem. Part II: Complexity Analysis

Bérenger Bramas, Muhammad Hassan+1 more|arXiv (Cornell University)|Nov 17, 2019

Electromagnetic Scattering and Analysis参考文献 83被引用数 4

ひとこと要約

本稿では、一様な媒質内に存在するN体の誘電体球の問題を解くために、ガレルキン境界積分方程式法の複雑性解析を提示している。その結果、線形系を解くために必要なGMRESおよびCG反復回数がNに依存せずに有界であることを示した。FMMと組み合わせることで、計算コストが線形スケーリングすることを保証し、大規模な極性誘電体粒子の電気的状態シミュレーションにおいて、本手法が精度面でも線形スケーリングであることを確立した。

ABSTRACT

This article is the second in a series of two papers concerning the mathematical study of a boundary integral equation of the second kind that describes the interaction of $N$ dielectric spherical particles undergoing mutual polarisation. The first article presented the numerical analysis of the Galerkin method used to solve this boundary integral equation and derived $N$-independent convergence rates for the induced surface charges and total electrostatic energy. The current article will focus on computational aspects of the algorithm. We provide a convergence analysis of the iterative method used to solve the underlying linear system and show that the number of liner solver iterations required to obtain a solution is independent of $N$. Additionally, we present two linear scaling solution strategies for the computation of the approximate induced surface charges. Finally, we consider a series of numerical experiments designed to validate our theoretical results and explore the dependence of the numerical errors and computational cost of solving the underlying linear system on different system parameters.

研究の動機と目的

• N個の誘電体球に対するガレルキン境界積分方程式法が計算コスト面で線形スケーリングであることを確立すること。
• 線形系を解くために必要なGMRESおよびCG反復回数が粒子数Nに依存しないことを証明すること。
• 変動するシステムパラメータの下で理論的に導出された収束境界を数値実験により検証すること。
• O(N)の計算複雑性を達成しつつ、N-誤差安定性を維持できることを示すこと。
• 本手法をラプラス方程式を超える他の演算子、特に静電的力の計算に拡張する基盤を築くこと。

提案手法

• 本手法は、均一な媒質内に存在するN個の誘電体球に対して定式化された第二種境界積分方程式をガレルキン離散化することで構成される。
• 得られる線形系は、反復解法（GMRESおよび、対称化後に適用可能な共役勾配法（CG））を用いて解かれる。
• 行列-ベクトル積の計算には、O(N)の計算量で実行可能な高速多重極法（FMM）が用いられ、スケーラブルな線形系の解法を可能にする。
• GMRESおよびCGの収束解析が実施され、与えられた許容誤差に対して反復回数にN依存性がない上界が存在することを証明した。
• 線形系を対称化することで、CG法の適用が可能となり、計算効率が向上するとともに、並列的な収束解析が可能になった。
• 数値実験は、ラップトップ上で実施され、Nとℓ_max（5および10）を変化させ、解法の収束許容誤差をそれぞれ10^−6および10^−9に設定し、スケーラビリティの検証が行われた。

実験結果

リサーチクエスチョン

• RQ1 線形系を解くために必要なGMRES反復回数は、誘電体球の数Nに伴い増加するか？
• RQ2 ガレルキン離散化から生じる対称化線形系に対して、共役勾配法（CG）を効果的に適用できるか？その反復回数はNに依存しないか？
• RQ3 FMMとNに依存しない収束性を有する反復解法を組み合わせた場合、アルゴリズム全体の計算コストはNに対して線形スケーリングするか？
• RQ4 実際の計算では、離散化誤差とFMM誤差はどのように相互作用するか？FMM誤差が解の精度を支配するか？
• RQ5 本手法は、類似したスケーラビリティと精度特性を有する静電的力の計算に拡張可能か？

主な発見

• 与えられた相対誤差許容誤差を達成するために必要なGMRES反復回数は、Nに依存せずに有界であり、反復解法のNに依存しない収束性を証明した。
• 対称化された線形系を解くために必要なCG反復回数も、Nに依存せずに有界であり、GMRESと同程度の性質を示す収束境界が存在する。
• 数値実験により、GMRESおよびCGの両アプローチにおいて、計算時間の線形スケーリングが確認された。Nおよびℓ_maxの変化に関わらず、ほぼ同一の性能が得られた。
• テストされた幾何的配置において、FMM誤差は離散化誤差を支配しないことが示され、漸近的複雑性解析の妥当性が裏付けられた。
• N-誤差安定性を有するガレルキン離散化と、FMMを用いたO(N)の反復解法の組み合わせにより、本手法が精度面でも線形スケーリングであることが確立された。
• 結果は、今後の研究計画として想定されている静電的力の計算へのフレームワークの拡張が現実的であることを支持している。