Skip to main content
SpringerLink
Log in

A prewavelet-based algorithm for the solution of second-order elliptic differential equations with variable coefficients on sparse grids

  • Original Paper
  • Published:
Numerical Algorithms Aims and scope Submit manuscript

Abstract

We present a Ritz-Galerkin discretization on sparse grids using prewavelets, which allows us to solve elliptic differential equations with variable coefficients for dimensions d ≥ 2. The method applies multilinear finite elements. We introduce an efficient algorithm for matrix vector multiplication using a Ritz-Galerkin discretization and semi-orthogonality. This algorithm is based on standard 1-dimensional restrictions and prolongations, a simple prewavelet stencil, and the classical operator-dependent stencil for multilinear finite elements. Numerical simulation results are presented for a three-dimensional problem on a curvilinear bounded domain and for a six-dimensional problem with variable coefficients. Simulation results show a convergence of the discretization according to the approximation properties of the finite element space. The condition number of the stiffness matrix can be bounded below 10 using a standard diagonal preconditioner.

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

Access this article

Log in via an institution

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. Achatz, S.: Higher order sparse grid methods for elliptic partial differential equations with variable coefficients. Computing 71(1), 1–15 (2003)

    Article  MathSciNet  MATH  Google Scholar 

  2. Balder, R., Zenger, C.: The solution of multidimensional real Helmholtz equations of sparse grids. SIAM J. Sci. Comp. 17(3), 631–646 (1996)

    Article  MathSciNet  MATH  Google Scholar 

  3. Barinka, A., Barsch, T., Charton, P., Cohen, A., Dahlke, S., Dahmen, W., Urban, K.: Adaptive wavelet schemes for elliptic problems—implementation and numerical experiments. SIAM J. Sci. Comput. 23(3), 910–939 (2002)

    Article  MathSciNet  MATH  Google Scholar 

  4. Bungartz, H.-J.: Du̇nne Gitter und deren Anwendung bei der adaptiven Lȯsung der dreidimensionalen Poisson-Gleichung. Dissertation, Fakultȧt fu̇r Informatik, Technische Universitȧt Mu̇nchen (1992)

  5. Dornseifer, T., Pflaum, C.: Discretization of elliptic differential equations on curvilinear bounded domains with sparse grids. Computing 56(3), 197–213 (1996)

    Article  MathSciNet  MATH  Google Scholar 

  6. Feuersȧnger, C.: Sparse grid methods for higher dimensional approximation. PhD thesis, University of Bonn (2010)

  7. Gordon, W.J., Hall, C.A.: Construction of curvilinear co-ordinate systems and applications to mesh generation. Int. J. Numer. Methods Eng. 7(4), 461–477 (1973)

    Article  MathSciNet  MATH  Google Scholar 

  8. Griebel, M.: Adaptive sparse grid multilevel methods for elliptic PDEs based on finite differences. Computing 61, 151–179 (1998)

    Article  MathSciNet  MATH  Google Scholar 

  9. Griebel, M., Hamaekers, J.: A wavelet based sparse grid method for electronic Schrödinger equation. In: Sanz-Solë, M, Soria, J., Varona, J.L., Verdera, J. (eds.) Proceedings of the International Congress of Mathematicians Madrid, August 22–30, 2006, vol. 1, pp. 1473–1506. European Mathematical Society, Publishing House (2007)

  10. Griebel, M., Hullmann, A.: On a multilevel preconditioner and its condition numbers for the discretized laplacian on full and sparse grids in higher dimensions. In: Singular Phenomena and Scaling in Mathematical Models, pp. 263–296. Springer International Publishing (2014)

  11. Kestler, S., Stevenson, R.: Fast evaluation of system matrices w.r.t. multi-tree collections of tensor product refinable basis functions. J. Comput. Appl Math. 260, 103–116 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  12. Niedermeier, A., Zimmer, S.: Implementational aspects of prewavelet sparse grid methods. In: Lai, C.-H., Bjoerstad, P., Cross, M., Widlund, O. (eds.) Eleventh International Conference of Domain Decomposition Methods (1999)

  13. Oswald, P.: Multilevel finite element approximation. Teubner Skripten zur Numerik. Teubner, Stuttgart (1994)

    Book  Google Scholar 

  14. Peherstorfer, B., Zimmer, S., Zenger, C., Bungartz, H.-J.: A multigrid method for adaptive sparse grids. SIAM J. Sci. Comput. 37(5), 51–70 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  15. Pflaum, C.: Diskretisierung elliptischer Differentialgleichungen mit dünnen Gittern. Dissertation, Technische Universität München (1995)

  16. Pflaum, C.: A multilevel algorithm for the solution of second order elliptic differential equations on sparse grids. Numer. Math. 79, 141–155 (1998)

    Article  MathSciNet  MATH  Google Scholar 

  17. Pflaum, C.: Robust convergence of multilevel algorithms for convection-diffusion equations. SIAM J. Numer. Anal. 37(2), 443–469 (2000)

    Article  MathSciNet  MATH  Google Scholar 

  18. Pflaum, C., Hartmann, R.: A sparse grid discretization of the helmholtz equation with variable coefficients in high dimensions. SIAM J. Numer. Anal. 54(4), 2707–2727 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  19. Vassilevski, P.S., Wang, J.: Stabilizing the hierarchical basis by approximate wavelets II implementation and numerical results. SIAM J. Sci. Comput. 20(2), 490–514 (1998)

    Article  MathSciNet  MATH  Google Scholar 

  20. Zeiser, A.: Fast matrix-vector multiplication in the sparse-grid Galerkin method. J. Sci. Comput. 47(3), 328–346 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  21. Zenger, C.: Sparse grids. In: Hackbusch, W. (ed.) Proceedings of the Sixth GAMM-Seminar Parallel Algorithms for Partial Differential Equations, Kiel, 1990, vol. 31 of Notes on Numerical Fluid Mechanics, pp. 240–251. Kiel, Vieweg (1991)

Download references

Acknowledgments

The authors gratefully acknowledge funding of the Erlangen Graduate School in Advanced Optical Technologies (SAOT) by the German Research Foundation (DFG) in the framework of the German excellence initiative.

Author information

Authors and Affiliations

  1. Department of Computer Science, System Simulation, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany

    Rainer Hartmann & Christoph Pflaum

Authors
  1. Rainer Hartmann
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Christoph Pflaum
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Rainer Hartmann.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hartmann, R., Pflaum, C. A prewavelet-based algorithm for the solution of second-order elliptic differential equations with variable coefficients on sparse grids. Numer Algor 78, 929–956 (2018). https://doi.org/10.1007/s11075-017-0407-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11075-017-0407-9

Keywords

Search

Navigation