A Mixed Precision Direct Electromagnetic Finite Element Solver on GPUs

In this paper, we present a mixed-precision direct electromagnetic finite element solver on GPU machines. In our solver, we propose a mixed-precision factorization algorithm aimed at obtaining results of the same quality as the double-precision algorithm at a much lower cost. In addition, to accelerate the direct solution of finite element sparse linear systems using GPU machines more efficiently, we propose a novel pipelined offloading technique that utilizes CUDA streams to achieve pipeline scheduling with overlapping computations and data transfers. Our numerical experiments show that by combining our proposed mixed-precision algorithm with the pipelined offloading technique, our proposed solver produces exact solutions with tremendous speedup.