[論文レビュー] Semi-analytic integration for a parallel space-time boundary element method modeling the heat equation

Jan Zapletal, Raphael Watschinger+2 more|arXiv (Cornell University)|Feb 19, 2021

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

ひとこと要約

本稿では、3次元熱方程式を解くための空間時間境界要素法における半解析的積分技法を提示する。時間成分を解析的に統合することで、行列構築と表現式評価を効率的に行える。この手法は、オープンソースのC++ライブラリbestheaに実装されており、1ソケットあたり最大18スレッドまで強スケーラビリティを示すマルチコアシステムで高い並列効率を達成しているが、NUMA効果によりソケット間のスケーリングが制限されている。

ABSTRACT

The presented paper concentrates on the boundary element method (BEM) for the heat equation in three spatial dimensions. In particular, we deal with tensor product space-time meshes allowing for quadrature schemes analytic in time and numerical in space. The spatial integrals can be treated by standard BEM techniques known from three dimensional stationary problems. The contribution of the paper is twofold. First, we provide temporal antiderivatives of the heat kernel necessary for the assembly of BEM matrices and the evaluation of the representation formula. Secondly, the presented approach has been implemented in a publicly available library besthea allowing researchers to reuse the formulae and BEM routines straightaway. The results are validated by numerical experiments in an HPC environment.

研究の動機と目的

• 3次元熱方程式を解くための、高効率な空間時間境界要素法のための堅牢で再利用可能なフレームワークの開発。
• 正確で高速なガレルキン行列構築に必要な、熱核の正確な時間積分の導出。
• 公開可能なC++ライブラリ（besthea）への実装により、将来のBEMプロジェクトにおける採用と再利用を促進。
• HPC環境における数値実験を通じた手法の妥当性とスケーラビリティの検証。
• 特に並列処理およびベクトル化実行を想定した、空間時間BEMの詳細で実装可能なガイダンスの不足に対処。

提案手法

• テンソル積による空間時間メッシュを採用し、時間方向には解析的積分、空間方向には数値積分を適用する。
• ガレルキン行列の構築および表現式評価に用いるため、熱核およびその法線微分の正確な逆微分（原始関数）を導出する。
• 単層、双層およびそれらの双対作用素を、標準的な空間BEM技術を用いたガレルキンBEMにより離散化する。
• 共有メモリシステムにおけるパフォーマンス向上のため、OpenMPによるスレッド処理とAVX512によるベクトル化を活用する。
• bestheaライブラリは再利用可能なC++ルーチンと数式を提供し、将来の高速多重極法（FMM）との統合に向けたニアフィールド計算に重点を置いている。
• 区分的定数および線形空間基底関数の両方をサポートし、スレッド数とメモリアクセスパターンの変化に伴うパフォーマンスを評価する。

実験結果

リサーチクエスチョン

• RQ1 3次元熱方程式の空間時間境界要素行列の時間積分をどのように解析的に実行すれば、計算コストを低減できるか？
• RQ2 行列構築および表現式評価に必要な、熱核およびその法線微分の正確な半解析的原始関数は何か？
• RQ3 現代のマルチコアアーキテクチャ上で、得られたBEMシステムはどの程度並列に効率的に構築・解法できるか？
• RQ4 NUMAノード（ソケット）をまたいでスケーリングする際の実装の限界は何か？
• RQ5 詳細なドキュメンテーション付きのオープンソースライブラリを通じて、この手法を研究コミュニティがどのように容易に再利用できるか？

主な発見

• 導出された熱核およびその法線微分の半解析的原始関数は正しいものであり、空間時間BEMにおける正確な時間積分を可能にする。
• bestheaライブラリはこの手法を効果的に実装しており、即座に再利用可能な数式とBEMルーチンを提供している。
• 1ソケット内では18スレッドまで、行列構築の効率が90％以上に達し、区分的定数および線形基底関数の両方でソケット内でのスケーリングが良好に実現されている。
• ポテンシャル評価についても同様のスケーリングが得られ、ソケット内では90％以上の効率を示すが、非NUMA対応のメモリ割り当てのため、両ソケットへのアクセス時に効率が低下している。
• std::vectorと標準アロケータを用いた実装のため、ソケット間でのスケーリングが最適でないが、これは将来のFMMベースで分散メモリを想定したバージョンで解消される見込みである。
• 数値実験により、HPC環境下での手法の正しさとパフォーマンスが確認され、理論的導出の妥当性が検証された。