A new iterative method for solving complex symmetric linear systems

Based on implementation of the quasi-minimal residual (QMR) and biconjugate A-orthogonal residual (BiCOR) method, a new Krylov subspace method is presented for solving complex symmetric linear systems. The new method can be combined with arbitrary symmetric preconditioners. The preconditioned modified Hermitian and Skew-Hermitian splitting (PMHSS) preconditioner is used to accelerate the convergence rate of this method. Numerical experiments indicate that the proposed method and its preconditioned version are efficient and robust, in comparison with other Krylov subspace methods.