本研究では、不規則な構造をもつ非構造格子系を離散化して得られる大規模連立1次方程式の並列解法について研究した。対称正定値な方程式に対しては前処理付き共役勾配法(PCG法、preconditioned CG method)が有効な解法であるが、PCG法の有効性は、前処理による収束の加速と、前処理の並列演算可能性に依存する。本研究では、楕円型偏微分方程式を不規則な有限要素法で離散化した方程式に対し、幾何学的構造に基づき、Maximum Independent Set法により自動的にマルチグリッドを生成するアルゴリズムを考案し、これにGalerkin法を用いてより粗いグリッド上の方程式を構成する。この方法は等方的問題については有効であるが、異方性のある問題では必ずしも収束の早いマルチグリッド法にはならない。そこで代数的な情報を用いながら異方的問題のsemicourseningを自動的に行う方法を提案した。数値実験により、不完全Cholesky前処理法(ICCG)などと比較した。多くの問題でICCGに近い性能を示し、ICCGが並列化に困難があることを考えると十分実用的になると考えられる。等方的でないもんだいでは、異方性が100以上ではICCGにかなわないが、より現実的な問題では高速に収束できることを見いだした。今後は、代数的にマルチグリッドを構成するAMGやそれを前処理とする解法について研究する必要がある。
|