Skip to main content
Log in

Inexact GMRES iterations and relaxation strategies with fast-multipole boundary element method

  • Published:
Advances in Computational Mathematics Aims and scope Submit manuscript

Abstract

Boundary element methods produce dense linear systems that can be accelerated via multipole expansions. Solved with Krylov methods, this implies computing the matrix-vector products within each iteration with some error, at an accuracy controlled by the order of the expansion, p. We take advantage of a unique property of Krylov iterations that allows lower accuracy of the matrix-vector products as convergence proceeds, and propose a relaxation strategy based on progressively decreasing p. In extensive numerical tests of the relaxed Krylov iterations, we obtained speed-ups of between 1.5 × and 2.3 × for Laplace problems and between 2.7 × and 3.3 × for Stokes problems. We include an application to Stokes flow around red blood cells, computing with up to 64 cells and problem size up to 131k boundary elements and nearly 400k unknowns. The study was done with an in-house multi-threaded C++ code, on a hexa-core CPU. The code is available on its version-control repository, https://github.com/barbagroup/fmm-bem-relaxed, and we share reproducibility packages for all results in https://github.com/barbagroup/inexact-gmres/.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Similar content being viewed by others

References

  1. Appel, A.W.: An efficient program for many-body simulation. SIAM J. Sci. Stat. Comput. 6(1), 85–103 (1985). https://doi.org/10.1137/0906008

    Article  MathSciNet  Google Scholar 

  2. Barnes, J., Hut, P.: A hierarchical \(O(N \log n)\) force-calculation algorithm. Nature 324, 446–449 (1986). https://doi.org/10.1038/324446a0

    Article  Google Scholar 

  3. Bouras, A., Frayssé, V.: A relaxation strategy for inexact matrix-vector products for Krylov methods. CERFACS TR0PA000015 European Centre for Research and Advanced Training in Scientific Computation (2000)

  4. Bouras, A., Frayssé, V.: Inexact matrix-vector products in Krylov methods for solving linear systems: A relaxation strategy. SIAM J. Matrix Analysis Applications 26(3), 660–678 (2005)

    Article  MathSciNet  Google Scholar 

  5. Brebbia, C., Dominguez, I.: Boundary Elements An Introductory Course, 2nd edn. WIT Press, USA (1992)

    MATH  Google Scholar 

  6. Cheng, H., Greengard, L., Rokhlin, V.: A fast adaptive multipole algorithm in three dimensions. J. Comput. Phys. 155, 468–498 (1999)

    Article  MathSciNet  Google Scholar 

  7. Van den Eshof, J., Sleijpen, G.: Inexact Krylov subspace methods for linear systems. SIAM J. Matrix Analysis Applications 26(1), 125–153 (2004)

    Article  MathSciNet  Google Scholar 

  8. Evans, E., Fung, Y.C.: Improved measurements of the erythrocyte geometry. Microvascular Research 4(4), 335–347 (1972)

    Article  Google Scholar 

  9. Fata, S.N.: Explicit expressions for 3D boundary integrals in potential theory. Int. J. Num. Meth. Engineering 78, 32–47 (2009)

    Article  MathSciNet  Google Scholar 

  10. Fata, S.N.: Explicit expressions for three-dimensional boundary integrals in linear elasticity. J. Comp. App. Mathematics 235, 4480–4495 (2011)

    Article  MathSciNet  Google Scholar 

  11. Fedosov, D.A., Caswell, B., Suresh, S., Karniadakis, G.E.: Quantifying the biophysical characteristics of Plasmodium-falciparum-parasitized red blood cells in microcirculation. Proc. Natl. Acad. Sci. 108(1), 35–39 (2011). https://doi.org/10.1073/pnas.1009492108. http://www.pnas.org/content/108/1/35.abstract

    Article  Google Scholar 

  12. Freund, J.B.: Numerical simulation of flowing blood cells. Ann. Rev. Fluid Mech. 46, 67–95 (2014)

    Article  MathSciNet  Google Scholar 

  13. Greengard, L.: The rapid evaluation of potential fields in particle systems. The MIT Press, Cambridge (1987)

    MATH  Google Scholar 

  14. Greengard, L., Rokhlin, V.: A fast algorithm for particle simulations. J. Comput. Phys. 73(2), 325–348 (1987). https://doi.org/10.1016/0021-9991

    Article  MathSciNet  Google Scholar 

  15. Hess, J.L., Smith, A.: Calculation of potential flow about arbitrary bodies. Progress in Aerospace Sciences 8, 1–138 (1967)

    Article  Google Scholar 

  16. Liu, Y., Nishimura, N.: The fast multipole boundary element method for potential problems: a tutorial. Engineering Analysis with Boundary Elements 30(5), 371–381 (2006). https://doi.org/10.1016/j.enganabound.2005.11.006

    Article  Google Scholar 

  17. Nishimura, N.: Fast multipole accelerated boundary integral equation methods. Appl. Mech. Rev. 55(4), 299–324 (2002). https://doi.org/10.1115/1.1482087. http://link.aip.org/link/?AMR/55/299/1

    Article  Google Scholar 

  18. Roache, P.J.: Verification and validation in computational science and engineering. Hermosa Albuquerque (1998)

  19. Saad, Y., Schultz, M.: GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput. 7, 856–869 (1986)

    Article  MathSciNet  Google Scholar 

  20. Simoncini, V., Szyld, D.: Theory of inexact Krylov subspace methods and applications to scientific computing. SIAM J. Sci. Comput. 25(2), 454–477 (2003)

    Article  MathSciNet  Google Scholar 

  21. Wang, T., Layton, S.K., Barba, L.A.: Inexact GMRES: FMM scaling figshare. https://doi.org/10.6084/m9.figshare.3100816.v1 (2016)

  22. Wang, T., Layton, S.K., Barba, L.A.: Inexact GMRES: Laplace equation figshare. https://doi.org/10.6084/m9.figshare.3100882.v2 (2016)

  23. Wang, T., Layton, S.K., Barba, L.A.: Inexact GMRES: Red blood cells figshare. https://doi.org/10.6084/m9.figshare.3100891.v1 (2016)

  24. Wang, T., Layton, S.K., Barba, L.A.: Inexact GMRES: Stokes equation figshare. https://doi.org/10.6084/m9.figshare.3100885.v3 (2016)

  25. White, C.A., Head-Gordon, M.: Rotating around the quartic angular momentum barrier in fast multipole method calculations. J. Chem. Phys. 105(12), 5061–5067 (1996)

    Article  Google Scholar 

  26. Yokota, R.: An FMM based on dual tree traversal for many-core architectures. Journal of Algorithms & Computational Technology 7(3), 301–324 (2013). https://doi.org/10.1260/1748-3018.7.3.301. Preprint arXiv:1209.3516v3

    Article  Google Scholar 

  27. Yokota, R., Barba, L.A.: A tuned and scalable fast multipole method as a preeminent algorithm for exascale systems. Int. J. High-perf. Comput. Applic. Published online 24 Jan. 2012 http://doi.org/fzgqgm preprint on arXiv:1106.2176 (2012)

  28. Yokota, R., Bardhan, J.P., Knepley, M.G., Barba, L.A., Hamada, T.: Biomolecular electrostatics using a fast multipole BEM on up to 512 GPUs and a billion unknowns. Comput. Phys. Comm. 182(6), 1271–1283 (2011). https://doi.org/10.1016/j.cpc.2011.02.013

    Article  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Science Foundation via NSF CAREER award OCI- 1149784. LAB acknowledges support from NVIDIA Corp. via the CUDA Fellows Program. Dr. Cris Cecka (previously at Harvard University, currently at Nvidia Corp.) contributed to the development of the code, particularly writing the octree and base evaluator. He later continued developing his fast-multipole framework, which evolved into his FMMTL project (see https://github.com/ccecka/fmmtl). The authors also wish to acknowledge valuable interactions with Dr. Christopher Cooper (previously at Boston University, currently at Universidad Técnica Federico Santa María) that helped with the implementation of the boundary element method.

Funding

This work was supported by the National Science Foundation via NSF CAREER award OCI-1149784.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Lorena A. Barba.

Ethics declarations

Conflict of interest

The authors declare no competing interests.

Additional information

Communicated by: Michael O’Neil

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This article belongs to the Topical Collection: Advances in Computational Integral Equations Guest Editors: Stephanie Chaillat, Adrianna Gillman, Per-Gunnar Martinsson, Michael O’Neil, Mary-Catherine Kropinski, Timo Betcke, Alex Barnett

Appendix: Algorithm listings

Appendix: Algorithm listings

figure r

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, T., Layton, S.K. & Barba, L.A. Inexact GMRES iterations and relaxation strategies with fast-multipole boundary element method. Adv Comput Math 48, 32 (2022). https://doi.org/10.1007/s10444-022-09932-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1007/s10444-022-09932-8

Keywords

Mathematics Subject Classification (2010)

Navigation