[論文レビュー] lifex: a flexible, high performance library for the numerical solution of complex finite element problems

Pasquale Claudio Africa|arXiv (Cornell University)|Jul 29, 2022

Advanced Numerical Methods in Computational Mathematics被引用数 10

ひとこと要約

lifex は、使いやすさ、スケーラビリティ、パフォーマンスに重点を置き、複雑なマルチフィジックス、マルチスケール、マルチドメインの有限要素問題を解くために、deal.II をベースに構築されたオープンソースで高パフォーマンスな C++ ライブラリです。lifex は、ラップトップからエクサスケール HPC システムに至る多様なアーキテクチャで、効率的かつポータブルなシミュレーションを可能にし、ほぼ理想に近い並列スループットを達成するとともに、任意の有限要素、複雑な結合、実行時パラメータチューニングをサポートし、最小限のコードで実装が可能です。

ABSTRACT

Numerical simulations are ubiquitous in mathematics and computational science. Several industrial and clinical applications entail modeling complex multiphysics systems that evolve over a variety of spatial and temporal scales. This study introduces the design and capabilities of lifex, an open source C++ library for high performance finite element simulations of multiphysics, multiscale, and multidomain problems. lifex meets the emerging need for versatile, efficient computational tools that are easily accessed by users and developers. We showcase its flexibility and effectiveness on a number of illustrative examples and advanced applications of use and demonstrate its parallel performance up to thousands of cores.

研究の動機と目的

• 多様な科学分野において、複雑なマルチフィジックスシステムをシミュレートするための、包括的で効率的かつユーザーフレンドリーな計算フレームワークの増大するニーズに対応する。
• 統合された環境内で、サイクリック実験、感度分析、最適化、不確実性定量化を可能にすることで、ソフトウェア開発を時間のかかるシミュレーションワークフローから分離する。
• 数値的正確性と効率性を維持したまま、複雑な有限要素ソルバーの実装を簡素化する高パフォーマンスで、ポータブルかつ拡張可能なライブラリを提供する。
• 標準ファイルフォーマット（例：VTK）、メッシュタイプ（六面体および四面体）、およびサードパーティライブラリ（PETSc、Trilinos、Boost）とのシームレスな相互運用性をサポートする。
• 良好にドキュメント化され、活発にメンテナンスされているコードベースを通じて、再現可能で再利用可能で透明性のあるサイクリック実験を可能にすることで、持続可能なオープンサイエンスエコシステムを促進する。

提案手法

• 複雑な偏微分方程式（PDE）およびその結合スキームの定義を簡素化するために、現代的な C++17 機能と洗練された抽象化プログラミングインターフェースを活用する。
• アダプティブで高次元、および hp-リファインド離散化のための堅牢なインfra を再利用するため、deal.II 有限要素ライブラリをコアバックエンドとして採用する。
• MPI を用いた分散メモリ並列処理を採用し、数千コアにわたるドメイン分割と負荷分散を実現することで、高いスケーラビリティを確保する。
• 実行時における線形代数バックエンド（PETSc、Trilinos、または deal.II 内蔵）およびソルバーの選択をサポートすることで、コード変更なしにパフォーマンスチューニングが可能となる。
• VTK を統合して可視化およびメッシュ I/O を実現し、標準フォーマットでのメッシュおよび解のデータのインポート／エクスポートをサポートする。
• Spack、Docker、mk、life x-env などの高度なデプロイメントツールを用いることで、Linux、Windows、macOS システム間でのポータビリティと再現性を確保する。

実験結果

リサーチクエスチョン

• RQ1 複雑なマルチフィジックス問題に対して、高パフォーマンス、使いやすさ、拡張性を両立させる有限要素ライブラリは、どのように設計できるか？
• RQ2 現代的な C++ の抽象化とモジュラ設計は、計算効率を損なわせることなく、実装の複雑さをどの程度軽減できるか？
• RQ3 多様な HPC アーキテクチャ上で、柔軟でユーザーフレンドリーな有限要素フレームワークが達成できる並列スケーラビリティの程度はどのようになるか？
• RQ4 実世界の科学的応用において、任意の有限要素、複雑な結合スキーム、実行時パラメータ構成が、ライブラリによってどの程度効果的にサポートできるか？
• RQ5 研究指向のライブラリは、初心者ユーザーとエキスパート開発者の両方にアクセス可能でありながら、生産環境レベルのパフォーマンスとポータビリティを達成できるか？

主な発見

• lifex は、最大数千コアでほぼ理想に近い並列スループットを達成し、高い計算効率と強いスケーラビリティを示している。
• パフォーマンスが重要なコンponentsにおいて、計算オーバーヘッドがほとんどないことから、抽象化レイヤーが実行速度を損なわないことが確認された。
• deal.II バックエンドを介した任意の有限要素のサポートにより、多様な PDE 型および離散化戦略にわたる正確なシミュレーションが可能となった。
• VTK および標準メッシュフォーマットとの相互運用性により、可視化および前処理／後処理ツールとのシームレスなデータ交換が実現された。
• 現代的な C++17 とポータブルなデプロイメントスタック（Docker、Spack、mk）の活用により、Linux、Windows、macOS プラットフォーム間で信頼性の高い再現可能なビルドが可能となった。
• ライブラリのモジュラ設計と実行時構成可能機能により、複雑なソルバーや結合スキームが最小限のコードで実装可能となり、開発生産性とコードの保守性が向上した。