Parallel Solution for 3D Volume Integral Equation using PETSc Library

The Method of Moments (MoM) is a popular method for solving volume integral equations arising in electromagnetic scattering from a Three-Dimensional (3D) in-homogeneous body problem. The MoM results in dense matrice, which makes it difficult in the application to realistic and complex large scale problems due to the high solution time. The solution time limitation can be overcome by going for distributed system computation. In this work we use the open source package ‘PETSc’ library which stands for ‘Portable, Extensible Toolkit for Scientific Computation’. PETSc is a suite of data structures and (iterative) solvers based on MPI to solve the matrix/linear system in parallel. The Krylov subspace method, GMRES (Generalized minimum residual) iterative scheme is used to solve the linear system of equations. In this work we show the implementation of PETSc for dense and complex valued MoM matrices. The algorithm is shown to scale well.