Compressed Fast Multipole Representations for Homogeneous 3-D Kernels
DOI:
https://doi.org/10.13052/2024.ACES.J.390201Keywords:
fast multipole method, integral equationAbstract
For homogeneous kernels, the memory requirements associated with H2 representations of integral equation matrices can be reduced by incorporating translational invariance. Starting with a non-translationally invariant H2 representation, this can be accomplished using a left/right iterative algorithm. In this paper, it is shown that a similar algorithm can also be used to compress an existing fast multipole method (FMM). It is observed that the iterative compression converges faster when used to compress an FMM than when it is applied to an H2 representation. Resulting savings in floating-point operations are indicated, and extensions of the reported method are discussed.
Downloads
References
L. Greengard and V. Rokhlin, “A fast algorithm for particle simulations,” J. Comput. Phys., vol. 73, no. 2, pp. 325-348, 1987.
W. Fong and E. Darve, “The black-box fast multipole method,” J. Comput. Phys., vol. 228, no. 23, pp. 8712-8725, Dec. 2009.
W. Hackbusch, Hierarchical Matrices: Algorithms and Analysis (Springer Series in Computational Mathematics). Berlin: Springer-Verlag, p. 511, 2015.
X. Xu and R. J. Adams, “Sparse matrix factorization using overlapped localizing LOGOS modes on a shifted grid,” IEEE T. Antenn. Propag., vol. 60, no. 3, pp. 1414-1424, 2012.
R. J. Adams, J. C. Young, and S. D. Gedney, “Compressing H2 matrices for translationally invariant kernels,” Applied Computational Electromagnetics Society (ACES) Journal, vol. 35, no. 11, pp. 1392-1393, Nov. 2020.
R. J. Adams, J. C. Young, and S. D. Gedney,“Efficiency improvements in compressing H2 matrices for translationally invariant kernels,” presented at the 2022 IEEE International Symposium on Antennas and Propagation and USNC-URSI Radio Science Meeting, Denver, CO, July 10-15, 2022.
R. J. Adams, J. C. Young, and S. D. Gedney, “Compressing a fast multipole method representation of an integral equation matrix,” presented at the 2023 International Applied Computational Electromagnetics Society Symposium, Monterey/Seaside, CA, 2023.