دورية أكاديمية

Combining Fast Multipole Techniques and an Approximate Inverse Preconditioner for Large Electromagnetism Calculations.

التفاصيل البيبلوغرافية
العنوان: Combining Fast Multipole Techniques and an Approximate Inverse Preconditioner for Large Electromagnetism Calculations.
المؤلفون: Carpentieri, B., Duff, I. S., Giraud, L., Sylvand, G.
المصدر: SIAM Journal on Scientific Computing; 2006, Vol. 27 Issue 3, p774, 19p, 3 Diagrams, 9 Charts, 3 Graphs
مصطلحات موضوعية: SPARSE matrices, MAXWELL equations, PARTIAL differential equations, BOUNDARY element methods, ELECTRIC fields, ELECTROMAGNETISM, INTEGRAL equations, LINEAR systems
مستخلص: The boundary element method has become a popular tool for the solution of Maxwell's equations in electromagnetism. From a linear algebra point of view, this leads to the solution of large dense complex linear systems, where the unknowns are associated with the edges of the mesh defined on the surface of the illuminated object. In this paper, we address the iterative solution of these linear systems via preconditioned Krylov solvers. Our primary focus is on the design of an efficient parallelizable preconditioner. In that respect, we consider an approximate inverse method based on the Frobenius-norm minimization. The preconditioner is constructed from a sparse approximation of the dense coefficient matrix, and the patterns both for the preconditioner and for the coefficient matrix are computed a priori using geometric information from the mesh. We describe how such a preconditioner can be naturally implemented in a parallel code that implements the multipole technique for the matrix-vector product calculation. We investigate the numerical scalability of our preconditioner on realistic industrial test problems and show that it exhibits some limitations on very large problems of size close to one million unknowns. To improve its robustness on those large problems we propose an embedded iterative scheme that combines nested GMRES solvers with different fast multipole computations. We show through extensive numerical experiments that this new scheme is extremely robust at affordable memory and CPU costs for the solution of very large and challenging problems. [ABSTRACT FROM AUTHOR]
Copyright of SIAM Journal on Scientific Computing is the property of Society for Industrial & Applied Mathematics and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
قاعدة البيانات: Complementary Index
الوصف
تدمد:10648275
DOI:10.1137/040603917