Optimized general sparse grid approximation spaces for operator equations

Griebel, M.; Knapek, S.

doi:10.1090/S0025-5718-09-02248-0

Optimized general sparse grid approximation spaces for operator equations
HTML articles powered by AMS MathViewer

by M. Griebel and S. Knapek PDF

Math. Comp. 78 (2009), 2223-2257 Request permission

Abstract:

This paper is concerned with the construction of optimized sparse grid approximation spaces for elliptic pseudodifferential operators of arbitrary order. Based on the framework of tensor-product biorthogonal wavelet bases and stable subspace splittings, we construct operator-adapted subspaces with a dimension smaller than that of the standard full grid spaces but which have the same approximation order as the standard full grid spaces, provided that certain additional regularity assumptions on the solution are fulfilled. Specifically for operators of positive order, their dimension is $O(2^{J})$ independent of the dimension $n$ of the problem, compared to $O(2^{Jn})$ for the full grid space. Also, for operators of negative order the overall cost is significantly in favor of the new approximation spaces. We give cost estimates for the case of continuous linear information. We show these results in a constructive manner by proposing a Galerkin method together with optimal preconditioning. The theory covers elliptic boundary value problems as well as boundary integral equations.

References

Robert A. Adams, Sobolev spaces, Pure and Applied Mathematics, Vol. 65, Academic Press [Harcourt Brace Jovanovich, Publishers], New York-London, 1975. MR 0450957
Pascal Auscher, Wavelets with boundary conditions on the interval, Wavelets, Wavelet Anal. Appl., vol. 2, Academic Press, Boston, MA, 1992, pp. 217–236. MR 1161253
H.-J. Bungartz, Dünne Gitter und deren Anwendung bei der adaptiven Lösung der dreidimensionalen Poisson-Gleichung, Dissertation, TU München, Institut für Informatik, 1992.
H.-J. Bungartz, Finite elements of higher order on sparse grids, Habilitation, TU München, Institut für Informatik, 1998.
Hans-Joachim Bungartz and Michael Griebel, A note on the complexity of solving Poisson’s equation for spaces of bounded mixed derivatives, J. Complexity 15 (1999), no. 2, 167–199. MR 1693892, DOI 10.1006/jcom.1999.0499
Hans-Joachim Bungartz and Michael Griebel, Sparse grids, Acta Numer. 13 (2004), 147–269. MR 2249147, DOI 10.1017/S0962492904000182
Charles K. Chui and Jian-zhong Wang, On compactly supported spline wavelets and a duality principle, Trans. Amer. Math. Soc. 330 (1992), no. 2, 903–915. MR 1076613, DOI 10.1090/S0002-9947-1992-1076613-3
Albert Cohen, Ingrid Daubechies, and Pierre Vial, Wavelets on the interval and fast wavelet transforms, Appl. Comput. Harmon. Anal. 1 (1993), no. 1, 54–81. MR 1256527, DOI 10.1006/acha.1993.1005
Wolfgang Dahmen, Stability of multiscale transformations, J. Fourier Anal. Appl. 2 (1996), no. 4, 341–361. MR 1395769
Wolfgang Dahmen, Wavelet and multiscale methods for operator equations, Acta numerica, 1997, Acta Numer., vol. 6, Cambridge Univ. Press, Cambridge, 1997, pp. 55–228. MR 1489256, DOI 10.1017/S0962492900002713
Wolfgang Dahmen, Helmut Harbrecht, and Reinhold Schneider, Adaptive methods for boundary integral equations: complexity and convergence estimates, Math. Comp. 76 (2007), no. 259, 1243–1274. MR 2299773, DOI 10.1090/S0025-5718-07-01970-9
Wolfgang Dahmen and Angela Kunoth, Multilevel preconditioning, Numer. Math. 63 (1992), no. 3, 315–344. MR 1186345, DOI 10.1007/BF01385864
W. Dahmen, S. Prössdorf, and R. Schneider, Wavelet approximation methods for pseudodifferential equations. II. Matrix compression and fast solution, Adv. Comput. Math. 1 (1993), no. 3-4, 259–335. MR 1242378, DOI 10.1007/BF02072014
Ronald A. DeVore, Nonlinear approximation, Acta numerica, 1998, Acta Numer., vol. 7, Cambridge Univ. Press, Cambridge, 1998, pp. 51–150. MR 1689432, DOI 10.1017/S0962492900002816
Ronald A. DeVore, Björn Jawerth, and Vasil Popov, Compression of wavelet decompositions, Amer. J. Math. 114 (1992), no. 4, 737–785. MR 1175690, DOI 10.2307/2374796
Karin Frank, Stefan Heinrich, and Sergei Pereverzev, Information complexity of multivariate Fredholm integral equations in Sobolev classes, J. Complexity 12 (1996), no. 1, 17–34. MR 1386591, DOI 10.1006/jcom.1996.0004
C. Feuersänger, Dünngitterverfahren für hochdimensionale elliptische partielle Differentialgleichungen, Diplomarbeit, Institut für Numerische Simulation, Universität Bonn, 2005.
M. Griebel, A parallelizable and vectorizable multi-level algorithm on sparse grids, Parallel algorithms for partial differential equations (Kiel, 1990) Notes Numer. Fluid Mech., vol. 31, Friedr. Vieweg, Braunschweig, 1991, pp. 94–100. MR 1167870
Michael Griebel, Multilevelmethoden als Iterationsverfahren über Erzeugendensystemen, Teubner Skripten zur Numerik. [Teubner Scripts on Numerical Mathematics], B. G. Teubner, Stuttgart, 1994 (German, with German summary). MR 1312162, DOI 10.1007/978-3-322-89224-9
M. Griebel, Adaptive sparse grid multilevel methods for elliptic PDEs based on finite differences, Computing 61 (1998), no. 2, 151–179. MR 1650985, DOI 10.1007/BF02684411
Michael Griebel, Sparse grids and related approximation schemes for higher dimensional problems, Foundations of computational mathematics, Santander 2005, London Math. Soc. Lecture Note Ser., vol. 331, Cambridge Univ. Press, Cambridge, 2006, pp. 106–161. MR 2277104, DOI 10.1017/CBO9780511721571.004
Michael Griebel and Jan Hamaekers, Sparse grids for the Schrödinger equation, M2AN Math. Model. Numer. Anal. 41 (2007), no. 2, 215–247. MR 2339626, DOI 10.1051/m2an:2007015
Michael Griebel and Jan Hamaekers, A wavelet based sparse grid method for the electronic Schrödinger equation, International Congress of Mathematicians. Vol. III, Eur. Math. Soc., Zürich, 2006, pp. 1473–1506. MR 2275738
M. Griebel, S. Knapek, Optimized approximation spaces for operator equations, Technical Report 568, SFB 256, Univ. Bonn, 1999.
M. Griebel and S. Knapek, Optimized tensor-product approximation spaces, Constr. Approx. 16 (2000), no. 4, 525–540. MR 1771694, DOI 10.1007/s003650010010
M. Griebel and D. Oeltz, A sparse grid space-time discretization scheme for parabolic problems, Computing 81 (2007), no. 1, 1–34. MR 2369419, DOI 10.1007/s00607-007-0241-3
M. Griebel and P. Oswald, Tensor product type subspace splittings and multilevel iterative methods for anisotropic problems, Adv. Comput. Math. 4 (1995), no. 1-2, 171–206. MR 1338900, DOI 10.1007/BF02123478
M. Griebel, P. Oswald, and T. Schiekofer, Sparse grids for boundary integral equations, Numer. Math. 83 (1999), no. 2, 279–312. MR 1712687, DOI 10.1007/s002110050450
Michael Griebel, Michael Schneider, and Christoph Zenger, A combination technique for the solution of sparse grid problems, Iterative methods in linear algebra (Brussels, 1991) North-Holland, Amsterdam, 1992, pp. 263–281. MR 1159736
R. Hochmuth, S. Knapek, G. Zumbusch, Tensor products of Sobolev spaces and applications, Technical Report 685, SFB 256, Univ. Bonn, 2000.
Stéphane Jaffard, Wavelet methods for fast resolution of elliptic problems, SIAM J. Numer. Anal. 29 (1992), no. 4, 965–986. MR 1173180, DOI 10.1137/0729059
S. Knapek, Hyperbolic cross approximation of integral operators with smooth kernel, Technical Report 665, SFB 256, Univ. Bonn, 2000.
S. Knapek, Compression of anisotropic tensor-product discretizations, Technical Report 0200, Institut für Numerische Simulation, Universität Bonn, 2002.
S. Knapek and F. Koster, Integral operators on sparse grids, SIAM J. Numer. Anal. 39 (2001/02), no. 5, 1794–1809. MR 1885717, DOI 10.1137/S0036142900375426
Lydia Kronsjö and Germund Dahlquist, On the design of nested iterations for elliptic difference equations, Nordisk Tidskr. Informationsbehandling (BIT) 12 (1972), 63–71. MR 312754, DOI 10.1007/bf01932674
Silvano Martello and Paolo Toth, Knapsack problems, Wiley-Interscience Series in Discrete Mathematics and Optimization, John Wiley & Sons, Ltd., Chichester, 1990. Algorithms and computer implementations. MR 1086874
Pál-Andrej Nitsche, Sparse approximation of singularity functions, Constr. Approx. 21 (2005), no. 1, 63–81. MR 2105391, DOI 10.1007/s00365-004-0559-4
P. Oswald, On discrete norm estimates related to multilevel preconditioners in the finite element method, in: Constructive Theory of Functions, K.G. Ivanov, P. Petrushev, B. Sendov (editors), Proc. Int. Conf. Varna, 1991, Bulg. Acad. Sci., Sofia, 203–214, 1992.
P. Oswald, Best N-term capacitance approximation on sparse grids, Proc. 12th Intern. Conf. on Domain Decomposition Methods in Science and Engineering, T. Chan et al. (editors), Chiba, 1999, 437–444, 2001.
S. V. Pereverzev, On the complexity of the problem of finding the solutions of Fredholm equations of the second kind with smooth kernels. II, Ukrain. Mat. Zh. 41 (1989), no. 2, 189–193, 287 (Russian); English transl., Ukrainian Math. J. 41 (1989), no. 2, 169–173. MR 992820, DOI 10.1007/BF01060382
D. Röschke, Über eine Kombinationstechnik zur Lösung partieller Differentialgleichungen, Diplomarbeit, Institut für Informatik, TU München, 1991.
H.-J. Schmeisser and H. Triebel, Topics in Fourier analysis and function spaces, Mathematik und ihre Anwendungen in Physik und Technik [Mathematics and its Applications in Physics and Technology], vol. 42, Akademische Verlagsgesellschaft Geest & Portig K.-G., Leipzig, 1987. MR 900143
Reinhold Schneider, Multiskalen- und Wavelet-Matrixkompression, Advances in Numerical Mathematics, B. G. Teubner, Stuttgart, 1998 (German, with German summary). Analysisbasierte Methoden zur effizienten Lösung großer vollbesetzter Gleichungssysteme. [Analysis-based methods for the efficient solution of large nonsparse systems of equations]. MR 1623209, DOI 10.1007/978-3-663-10851-1
Ian H. Sloan and Henryk Woźniakowski, When are quasi-Monte Carlo algorithms efficient for high-dimensional integrals?, J. Complexity 14 (1998), no. 1, 1–33. MR 1617765, DOI 10.1006/jcom.1997.0463
Ian H. Sloan and Henryk Woźniakowski, Tractability of integration in non-periodic and periodic weighted tensor product Hilbert spaces, J. Complexity 18 (2002), no. 2, 479–499. Algorithms and complexity for continuous problems/Algorithms, computational complexity, and models of computation for nonlinear and multivariate problems (Dagstuhl/South Hadley, MA, 2000). MR 1919445, DOI 10.1006/jcom.2001.0626
S. A. Smoljak, Quadrature and interpolation formulae on tensor products of certain function classes, Dokl. Akad. Nauk SSSR 148 (1963), 1042–1045 (Russian). MR 0147825
Rob Stevenson, On the compressibility operators in wavelet coordinates, SIAM J. Math. Anal. 35 (2004), no. 5, 1110–1132. MR 2050194, DOI 10.1137/S0036141002411520
J. F. Traub, G. W. Wasilkowski, and H. Woźniakowski, Information-based complexity, Computer Science and Scientific Computing, Academic Press, Inc., Boston, MA, 1988. With contributions by A. G. Werschulz and T. Boult. MR 958691
Joachim Weidmann, Linear operators in Hilbert spaces, Graduate Texts in Mathematics, vol. 68, Springer-Verlag, New York-Berlin, 1980. Translated from the German by Joseph Szücs. MR 566954, DOI 10.1007/978-1-4612-6027-1
Arthur G. Werschulz, The computational complexity of differential and integral equations, Oxford Mathematical Monographs, The Clarendon Press, Oxford University Press, New York, 1991. An information-based approach; Oxford Science Publications. MR 1144521
Arthur G. Werschulz, The complexity of the Poisson problem for spaces of bounded mixed derivatives, The mathematics of numerical analysis (Park City, UT, 1995) Lectures in Appl. Math., vol. 32, Amer. Math. Soc., Providence, RI, 1996, pp. 895–914. MR 1421372, DOI 10.1006/jcom.1996.0027
Christoph Zenger, Sparse grids, Parallel algorithms for partial differential equations (Kiel, 1990) Notes Numer. Fluid Mech., vol. 31, Friedr. Vieweg, Braunschweig, 1991, pp. 241–251. MR 1167882