Fast direct solver for surface integral equation methods for electromagnetic wave scattering from PEC targets

A parallel fast direct solver for surface integral equations based on interpolative decomposition (ID) is proposed. The method of moments (MoM) matrix can be approximated by using ID skeletonization. The error of this decomposition can be well-controlled. By introducing an artificial cube, the cost of the skeletonization can be reduced greatly. OpenMP parallel computing technology is employed to accelerate the calculation. The parallelization saves a lot of CPU time in comparison with the traditional direct method. As a result, the proposed parallel algorithm can achieve accurate and efficient decomposition where the traditional direct method can hardly handle.