Improve the efficiency of the Cartesian tensor based fast multipole method for Coulomb interaction using the traces
- Old Dominion Univ., Norfolk, VA (United States)
- Thomas Jefferson National Accelerator Facility (TJNAF), Newport News, VA (United States)
To compute the non-oscillating mutual interaction for a systems with N points, the fast multipole method (FMM) has an efficiency that scales linearly with the number of points. Specifically, for Coulomb interaction, FMM can be constructed using either the spherical harmonic functions or the totally symmetric Cartesian tensors. In this paper, we will present that the effciency of the Cartesian tensor-based FMM for the Coulomb interaction can be significantly improved by implementing the traces of the Cartesian tensors in calculation to reduce the independent elements of the n-th rank totally symmetric Cartesian tensor from (n + 1)(n + 2)=2 to 2n + 1. The computation complexity for the operations in FMM are analyzed and expressed as polynomials of the highest rank of the Cartesian tensors. For most operations, the complexity is reduced by one order. Numerical examples regarding the convergence and the effciency of the new algorithm are demonstrated. As a result, a reduction of computation time up to 50% has been observed for a moderate number of points and rank of tensors.
- Research Organization:
- Thomas Jefferson National Accelerator Facility (TJNAF), Newport News, VA (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Nuclear Physics (NP)
- Grant/Contract Number:
- AC05-06OR23177
- OSTI ID:
- 1439411
- Alternate ID(s):
- OSTI ID: 1582863
- Report Number(s):
- JLAB-ACP-17-2560; DOE/OR/23177-4223; R&D Project: 2016-LDRD-7; PII: S0021999118303280; TRN: US1900604
- Journal Information:
- Journal of Computational Physics, Vol. 371, Issue C; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
A New Relatively Simple Approach to Multipole Interactions in Either Spherical Harmonics or Cartesians, Suitable for Implementation into Ewald Sums
|
journal | December 2019 |
Similar Records
Traceless Totally Symmetric Tensor Based Fast Multipole Method for Space Charge Field Calculation
Fast multipole method using Cartesian tensor in beam dynamic simulation