A domain decomposition finite-difference method for parallel numerical implementation of time-dependent Maxwell´s equations

A domain decomposition technique together with an implicit finite-difference scheme is used to design a parallel algorithm to solve for electromagnetic scattering in the time-domain by an infinite square metallic cylinder. The implicit difference scheme yields second-order accuracy, unconditional stability, and at each time step, a large system of linear equations. The domain decomposition technique reduces the solution of this large system to that of many independent smaller subsystems. The present algorithm can be easily implemented on coarse-grain parallel vector supercomputers to obtain a speedup close to the number of available central processing units (CPUs)