組み込みのマルチフロンタル直接法ソルバーがほしい
課題
外部ライブラリに依存しない直接法ソルバーは組み込みの DIRECT のみで、逐次はスカラー実装にとどまり、MPI 並列版はあるもののプロセス数が 2^n+1 個に限られる制約がある。スレッド並列・MPI 並列・大規模問題・接触解析で実用になる直接法は MUMPS または MKL を有効にしたビルドが前提であり、ライブラリを用意できない環境では選択肢がない。
また、ソルバー内部の改良(ピボット戦略、低ランク圧縮、分散方式など)を FrontISTR 側で完結して行える直接法の土台がない。
目標
hecmw にマルチフロンタル法の直接法ソルバーを組み込みで追加し、!SOLVER の METHOD から選択できるようにする。満たしたい要件:
- 外部の直接法ライブラリなしでビルド・実行できる
- スレッド並列・MPI 並列で実用的な性能が出る
- 不定値を含む対称行列と数値非対称な行列(摩擦接触・非関連塑性など)の双方を扱える
- 接触解析や固有値解析など、頑健性の面で直接法が実用上選ばれる解析種別で、既存の MUMPS/MKL と同様に選択できる
- SA-AMG 前処理の粗グリッド用ソルバーとしても使える
- MUMPS 利用時と同等の解が得られることをテストで確認できる