Study on Quadruple-Precision Arithmetic in Electromagnetic Field Analysis by Hierarchical Domain Decomposition Method

Abstract

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.