This paper deals with a finite element analysis (FEA) of an electromagnetic problem with a quadruple precision floating point number. A Conjugate Gradient (CG) method for real symmetric matrices and a Conjugate Orthogonal Conjugate Gradient (COCG) method for complex symmetric matrices are implemented in the quadruple precision. Then, an FEA solver with the quadruple precision arithmetic for the electromagnetic problem is implemented by the C language. As a result, the convergence of the iterative method is improved, and the minimum residual norm is reduced by approximately three orders of magnitude.
抄録全体を表示