Parallelising a simulator for the analysis of electromagnetic radiation using MUMPS library

The practical experience of parallelising a simulator of general scattering and radiation electromagnetic problems is presented. The simulator stems from an existing sequential simulator in the frequency domain and can be fruitfully used in applications such as the test of coverage of a wireless network, analysis of complex structures, and so on. After the analysis of a test case, two steps were carried out: firstly, a "hand-crafted" code parallelisation was developed within the kernel of the simulator. Secondly, the sequential library, used in the existing simulator, was replaced by the parallel MUMPS library in order to solve the associated linear algebra problem in parallel. For factorising the matrix with MUMPS, two different ordering methods have been considered.

[1]  Tapan K. Sarkar,et al.  Fully coupled hybrid‐method FEM/high‐frequency technique for the analysis of 3D scattering and radiation problems , 2005 .

[2]  Vipin Kumar,et al.  Analysis of Multilevel Graph Partitioning , 1995, Proceedings of the IEEE/ACM SC95 Conference.

[3]  Patrick Amestoy,et al.  Hybrid scheduling for the parallel solution of linear systems , 2006, Parallel Comput..

[4]  Jack Dongarra,et al.  Basic Linear Algebra Subprograms (BLAS) , 2011, Encyclopedia of Parallel Computing.

[5]  Jennifer A. Scott.  On ordering elements for a frontal solver , 1999 .

[6]  R. Fernandez-Recio,et al.  Fully Coupled Hybrid FEM-UTD Method Using NURBS for the Analysis of Radiation Problems , 2008, IEEE Transactions on Antennas and Propagation.

[7]  Joseph W. H. Liu,et al.  The Multifrontal Method for Sparse Matrix Solution: Theory and Practice , 1992, SIAM Rev..

[8]  B. McDonald,et al.  Finite-Element Solution of Unbounded Field Problems , 1972 .

[9]  Jürgen Schulze.  Towards a Tighter Coupling of Bottom-Up and Top-Down Sparse Matrix Ordering Methods , 2001 .

[10]  Jianming Jin,et al.  Application of hybrid finite element method to electromagnetic scattering from coated cylinders , 1988 .

[11]  S. Sloan.  An algorithm for profile and wavefront reduction of sparse matrices , 1986 .

[12]  M. Salazar-Palma,et al.  A finite element method for the analysis of radiation and scattering of electromagnetic waves on complex environments , 2005 .

[13]  Tapan K. Sarkar,et al.  Iterative and Self-Adaptive Finite-Elements in Electromagnetic Modeling , 1998 .

[14]  O. Zienkiewicz,et al.  The coupling of the finite element method and boundary solution procedures , 1977 .

[15]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..

[16]  Patrick R. Amestoy,et al.  Multifrontal parallel distributed symmetric and unsymmetric solvers , 2000 .

[17]  Magdalena Salazar-Palma,et al.  Iterative and self-adaptive finite-elements in electromagnetic modeling , 1998 .