A Stabilized Multilevel B-Spline-Based Fast Integral Method for the Solution of the Electric Field Integral Equation

技術要約：EFIEのための安定化した多レベルBスプラインに基づく高速積分法

問題提起
電気界積分方程式（EFIE）は、表面離散化のみを必要とする電磁散乱および放射問題を解くための標準的なツールである。しかし、標準的な境界要素法（BEM）の離散化では、未知数の数を$N$とすると、行列のアセンブリ、格納、および行列ベクトル積（MVP）において$O(N^2)$の計算量が必要となる。計算量を$O(N)$またはそれに近いオーダーに低減するために、多レベル高速多重極展開法（MLFMM）や、FFT加速型のアダプティブ積分法（AIM）のような高速積分法が開発されてきた。

等間隔のラグランジュ補間を用いるFFT加速型の手法には、特定の制限が存在する。これらの手法は、FFT構造を利用するために等間隔の補間点を必要とする。しかし、等間隔のラグランジュ補間は、高次の多項式においてルンゲ現象（Runge-type instabilities）による不安定性を引き起こし、精度の崩壊を招き、堅牢な高次圧縮を妨げる。Bスプラインは、等間隔格子におけるルンゲ現象に対する解決策となるが、多レベル高速積分法への適用は非自明である。ラグランジュ多項式とは異なり、Bスプラインの係数は補間点における関数値と一致しないため、サンプリング行列の逆行列計算を必要とする。これらの行列は高次において極めて悪条件（ill-conditioned）になりやすく、数値的安定性を損なう。さらに、多レベルの効率性に不可欠な「正確なレベル間転送」を維持しつつ、安定性を確保することは、これまで大きな課題であった。

手法
著者らは、FFT適合型のカーネル補間と堅牢な高次補間を組み合わせた、安定化した多レベルBスプラインに基づく高速積分法を提案している。核となる手法は、主に以下の3つのコンポーネントで構成される。

1. Bスプライン・カーネル補間：
   グローバルなラグランジュ多項類の代わりに、等間隔のノットベクトル上に定義されたBスプラインを採用する。この選択により、FFT加速に必要な等間隔サンプリングとの互換性を維持しつつ、ルンゲ現象を緩和できる。グリーン関数はテンソル積多項式を用いて因数分解される。Bスプラインの場合、補間重みはグリーン関数の直接的なサンプルではなく、サンプリング行列$B$を介して導出される。Bスプラインのテンソル積構造により、サンプリング行列の逆行列計算を、ナイーブな$O(s^6)$ではなく、効率的な$O(s^3)$の演算（ここで$s$はスプラインの数）で行うことが可能となり、MVPの全体的な$O(s^3 \log s)$の計算量を維持できる。

2. ノット除去による安定化：
   高次におけるBスプライン・サンプリング行列の悪条件に対処するため、著者らは「ノット除去（knot-removal）」による安定化戦略を導入している。経験的な分析により、等間隔ノットベクトルの境界付近の特定のノットを除去することで、多レベルスキームに必要な等間隔構造を破壊することなく、安定性定数（逆サンプリング行列のノルム）を大幅に減少させられることが示されている。

• バッファ・ノット： 基底関数が関連するボックスから突出する問題（BEMにおける一般的な問題）に対処するため、ノットベクトルに「バッファ・ノット」を追加して拡張する。これにより、すべての基底関数のサポートが補間領域内に含まれることが保証される。
• 正確なレベル間転送： 多レベル手法における重要な要件は、親ボックスと子ボックスの間でモーメントを正確に転送できることである。提案された安定化戦略は、ノットベクトルの入れ子構造を保持する。「ノット挿入（knot insertion）」（標準的なスプライン細分化技術）を用いることで、レベル間の正確な転送行列を構築する。これにより、安定化（ノット除去）が適用されたグリッド上で計算されたモーメントを、カーネル補間に使用される元の等間隔グリッドへ正確にマッピングでき、レベル間でのマシン精度での正確な転送が維持される。

3. 因数分解戦略：
   本論文では、離散化されたEFIEに対して、既存の直接因数分解およびダイアド（dyadic）因数分解に加え、新しい「スターベース（star-based）」因数分解を導入している。

• スターベース因数分解： ループ・スター分解の概念を用いて、スカラーポテンシャル演算子の表現を利用する。これは、スカラーポテンシャル演算子を、疎なスター行列と単層（single-layer）演算子へと因数分解するものである。この戦略は、スカラーポテンシャル演算子の零空間（null space）をマシン精度まで保持するため、低周波プリコンディショニング（例：カルデロン・プリコンディショニング）の効果にとって極めて重要である。これには、RWGおよび定数（piecewise-constant）モーメントの両方の計算が必要となる。
• 直接およびダイアド因数分解： これらは代替案としてレビューされている。直接因数分解はカーネルを直接補間し、ダイアド形式はダイアド・カーネルを使用する。スターベースのアプローチは、零空間に関する優れた安定性特性を持つことから強調されている。

主な貢献

• 安定化したBスプライン補間： 等間隔のBスプライン補間におけるルンゲ現象とサンプリング行列の悪条件を排除し、堅牢な高次圧縮を可能にするノット除去戦略の導入。
• 正確な多レベル転送： ノット挿入を用いることで、安定化（ノット除去）が適用されている場合でも、正確なレベル間転送を実現するメカニズムの開発。これは、安定した補間と正確な階層的転送の間の矛盾を解決する。
• スターベースの因数分解： スカラーポテンシャル演算子の零空間をマシン精度まで保持する、EFIEのための新しい因数分解戦略の提案。これにより、低周波安定化およびプリコンディショニングとの適合性が向上する。
• 計算量の維持： Bスプラインのコンパクトなサポートとサンプリング行列のテンソル積構造を活用することで、MVP、アセンブリ、およびストレージにおいて$O(N)$の計算量を維持することを実証。

結果
標準的な形状（球体）および現実的な形状（航空機）の幾何学形状に対し、様々な電気的サイズ（$0.2\lambda$から$4\lambda$）にわたって数値結果が示された。

• 安定性と精度： Bスプラインベースの手法は、ラグランジュベースのアプローチで見られた崩壊を起こすことなく、安定した高次補間を示した。誤差解析により、多項式の次数$p$およびスプライン数$s$が増加するにつれて、理論的限界と一致する収束が示された。安定化戦略は、高次（$p=5, 6$）においても誤差を制御することに成功したが、ラグランジュ法は不安定性を示した。
• 因数分解の性能： スターベースの因数分解は、テストされた3つの戦略の中で最も堅牢な誤差挙動を示した。
• スケーラビリティ： 性能分析により、球体および航空機モデルの両方において、ストレージ、アセンブリ時間、およびMVP時間の$O(N)$スケーリングが確認された。未知数が約300万の球体の場合、メモリ消費量は$p$に応じて91 GBから182 GBの範囲であり、MVP時間は7.6秒から54秒であった。最大約400万の未知数を持つ航空機モデルについても、同様の線形スケーリングが観察された。

意義
本論文は、提案された安定化した多レベルBスプライン法が、高速積分方程式ソルバーのための、実用的かつ高精度なラグランジュ補間の代替手段を提供すると主張している。ルンゲ現象とサンプリング行列の悪条件を克服することで、本手法は、従来の等間隔FFT加速型スキームでは到達不可能であった堅牢な高次圧縮を可能にする。スカラーポテンシャルの零空間をスターベースの因数分解によって保持することは、低周波問題への有用性をさらに高める。著者らは、予備的な結果が会議で発表されていること、および本手法が広帯域定式化のためのMLFMMのような高周波技術と組み合わせることが可能であることを述べているが、主な貢献は低・中周波カーネル補間フレームワークの安定化にある。