Skip to main content
SpringerLink
Log in
Book cover

Numerical Mathematics and Advanced Applications 2011 pp 653–661Cite as

3D Helmholtz Krylov Solver Preconditioned by a Shifted Laplace Multigrid Method on Multi-GPUs

  • Conference paper
  • First Online:

Abstract

We are focusing on an iterative solver for the three-dimensional Helmholtz equation on multi-GPU using CUDA (Compute Unified Device Architecture). The Helmholtz equation discretized by a second order finite difference scheme is solved with Bi-CGSTAB preconditioned by a shifted Laplace multigrid method. Two multi-GPU approaches are considered: data parallelism and split of the algorithm. Their implementations on multi-GPU architecture are compared to a multi-threaded CPU and single GPU implementation. The results show that the data parallel implementation is suffering from communication between GPUs and CPU, but is still a number of times faster compared to many-cores. The split of the algorithm across GPUs limits communication and delivers speedups comparable to a single GPU implementation.

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

Chapter
USD   29.95
Price excludes VAT (USA)
  • Available as PDF
  • Read on any device
  • Instant download
  • Own it forever
eBook
USD   169.00
Price excludes VAT (USA)
  • Available as EPUB and PDF
  • Read on any device
  • Instant download
  • Own it forever
Softcover Book
USD   219.99
Price excludes VAT (USA)
  • Compact, lightweight edition
  • Dispatched in 3 to 5 business days
  • Free shipping worldwide - see info
Hardcover Book
USD   219.99
Price excludes VAT (USA)
  • Durable hardcover edition
  • Dispatched in 3 to 5 business days
  • Free shipping worldwide - see info

Tax calculation will be finalised at checkout

Purchases are for personal use only

Learn about institutional subscriptions

Notes

  1. 1.

    During the work on this paper, the newer version of CUDA 4.0 has been released. It was not possible to have the newer version installed on all systems for our experiments. That is why for consistency and comparability of experiments, we use the previous version

References

  1. J.J. Dongarra, I.S. Duff, D.C. Sorensen, and H.A. van der Vorst. Solving Linear Systems on Vector and Shared Memory Computers. SIAM, Philadelphia (1991).

    Google Scholar 

  2. B. Engquist and A. Majda. Absorbing boundary conditions for numerical simulation of waves. Math. Comput., 31:629–651 (1977).

    Article  MathSciNet  MATH  Google Scholar 

  3. Y. A. Erlangga, C. W. Oosterlee, and C. Vuik. A novel multigrid based preconditioner for heterogeneous Helmholtz problems. SIAM J. Sci. Comput., 27:1471–1492 (2006).

    Article  MathSciNet  MATH  Google Scholar 

  4. O. Ernst and M. Gander. Why it is difficult to solve Helmholtz problems with classical iterative methods. In Durham Symposium 2010 (2010).

    Google Scholar 

  5. H. Knibbe, C. W. Oosterlee, and C. Vuik. GPU implementation of a Helmholtz Krylov solver preconditioned by a shifted Laplace multigrid method. Journal of Computational and Applied Mathematics, 236:281–293 (2011).

    Article  MathSciNet  MATH  Google Scholar 

  6. http://www.nvidia.com (2011).

  7. E. Zhebel. A Multigrid Method with Matrix-Dependent Transfer Operators for 3D Diffusion Problems with Jump Coefficients. PhD thesis, Technical University Bergakademie Freiberg, Germany (2006).

    Google Scholar 

Download references

Author information

Authors and Affiliations

  1. Delft University of Technology, Delft, Netherlands

    H. Knibbe & C. Vuik

  2. Dutch National Research Centre for Mathematics and Computer Science (CWI), Delft University of Technology, Delft, Netherlands

    C. W. Oosterlee

Authors
  1. H. Knibbe
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. C. W. Oosterlee
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. C. Vuik
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to H. Knibbe .

Editor information

Editors and Affiliations

  1. , Department of Mathematics, University of Leicester, Leicester, LE!1 7RH, United Kingdom

    Andrea Cangiani

  2. , Department of Mathematics, University of Leicester, Leicester, LE1 7RL, United Kingdom

    Ruslan L. Davidchack

  3. , Department of Mathematics, University of Leicester, Leicester, LE1 7RH, United Kingdom

    Emmanuil Georgoulis

  4. , Department of Mathematics, University of Leicester, Leicester, LE1 7RH, United Kingdom

    Alexander N. Gorban

  5. , Department of Mathematics, University of Leicester, Leicester, LE1 7RH, United Kingdom

    Jeremy Levesley

  6. , School of Mathematical Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, United Kingdom

    Michael V. Tretyakov

Rights and permissions

Reprints and permissions

Copyright information

© 2013 Springer-Verlag Berlin Heidelberg

About this paper

Cite this paper

Knibbe, H., Oosterlee, C.W., Vuik, C. (2013). 3D Helmholtz Krylov Solver Preconditioned by a Shifted Laplace Multigrid Method on Multi-GPUs. In: Cangiani, A., Davidchack, R., Georgoulis, E., Gorban, A., Levesley, J., Tretyakov, M. (eds) Numerical Mathematics and Advanced Applications 2011. Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-642-33134-3_69

Download citation

Publish with us

Policies and ethics

Search

Navigation