複素対称線形システム向けのCOMINRES-QLP法

抄録

Solving complex symmetric systems efficiently is a significantly important issue in the applications of electromagnetics since the discretization of the numerical model using the finite element method results in complex symmetric systems when computing the electromagnetic field. However, such complex symmetric systems are usually much more difficult to be solved because of the property of non-Hermitian therefore, efficient numerical methods for the solution of complex symmetric systems are demanded. In this study, we proposed an iterative method called COMINRES-QLP which is specially designed for complex symmetric linear systems. Moreover, we proposed a new stopping criterion which could be very useful in practical computations when the condition number of the coefficient matrix is unknown.