社団法人 物理探査学会
第117回(平成19年度秋季)学術講演会


電気・電磁探査法数値解析における連立方程式数値解法の基礎研究

講演要旨(和文)
電気電磁探査法数値計算問題における線形連立方程式求解について,より簡便にこの問題を解くための数値計算法について検討した.具体的には,主に大規模疎行列係数を持つ連立方程式を解くための,汎用計算システムを MATLAB (GNU Octave) および C++ 言語で試作し,数値計算例として差分法による線形方程式解法に適用した.複数の前処理と解法を組み合わせて収束性を検証し,前処理や解法によって収束性が変化することを確認した.現時点では,定量的に最適な解法をあらかじめ求めたり,どんな問題でも最速に解くことができる解法というものが存在しない以上,今回試作したような,複数の前処理と解法を選択して利用可能なシステムが,より効率良い数値計算を行う助けとなる可能性がある.今後は,試作したシステムの改良を続け,より物理探査数値計算の線形方程式求解に特化して利便性も高め,将来的な公開を目指したい.

講演要旨(英文)
An important problem in numerical analysis of electric/electromagnetic methods is developing a fast and stable solver for the system of linear equations to solve the partial differential equation (PDE). In order to overcome this problem, I have compiled several iterative solvers with various preconditioners for the large sparse linear system as a test version of geophysics-oriented iterative solver system. The compiled system has been evaluated with two different types of PDEs using various iterative solvers and preconditioners. The results show the efficiency of the new system and importance of developing faster solvers for numerical problems in geophysics.