×

Adaptive recompression of \(\mathcal H\)-matrices for BEM. (English) Zbl 1070.65028

Summary: The efficient treatment of dense matrices arising, e.g., from the finite element discretisation of integral operators requires special compression techniques. In this article, we use a hierarchical low-rank approximation, the so-called \(\mathcal H\)-matrix, that approximates the dense stiffness matrix in admissible blocks (corresponding to subdomains where the underlying kernel function is smooth) by low rank matrices. The low rank matrices are assembled by the \(ACA+\) algorithm, an improved variant of the well-known ACA method. We present an algorithm that can determine a coarser block structure that minimises the storage requirements (enhanced compression) and speeds up the arithmetic (e.g., inversion) in the \(\mathcal H\)-matrix format. This coarse approximation is done adaptively and on-the-fly to a given accuracy such that the matrix is assembled with minimal storage requirements while keeping the desired approximation quality. The benefits of this new recompression technique are demonstrated by numerical examples.

MSC:

65F30 Other matrix algorithms (MSC2010)
65F05 Direct numerical methods for linear systems and matrix inversion
65F50 Computational methods for sparse matrices
35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation
65N38 Boundary element methods for boundary value problems involving PDEs

Software:

hlib
Full Text: DOI

References:

[1] Bebendorf, M.: Approximation of boundary element matrices. Numer. Math. 86(4), 565–589 (2000). · Zbl 0966.65094
[2] Bebendorf, M.: Hierarchical LU decomposition based preconditioners for BEM. Technical Report 28, Max-Planck-Institute for Mathematics in the Sciences 2004. · Zbl 1071.65031
[3] Bebendorf, M., Rjasanov, S.: Adaptive low-rank approximation of collocation matrices. Computing 70, 1–24 (2003). · Zbl 1068.41052
[4] Börm, S., Grasedyck, L.: HLib – a library for - and -matrices, 1999. Available at http://www.hlib.org.
[5] Börm, S., Grasedyck, L.: Low-rank approximation of integral operators by interpolation. Computing 72, 325–332 (2004). · Zbl 1059.65121
[6] Börm, S., Grasedyck, L., Hackbusch, W.: Introduction to hierarchical matrices with applications. Engineering Analysis with Boundary Elements 27, 405–422 (2003). · Zbl 1035.65042
[7] Dahmen, W., Schneider, R.: Wavelets on manifolds I: Construction and domain decomposition. SIAM J. Math Anal. 31, 184–230 (1999). · Zbl 0955.42025
[8] Sauter, S. A., Erichsen, S.: Efficient automatic quadrature in 3-D Galerkin BEM. Comput. Meth. Appl. Mech. Eng. 157, 215–224 (1998). · Zbl 0943.65139
[9] Grasedyck, L.: Theorie und Anwendungen Hierarchischer Matrizen. PhD thesis, Universität Kiel 2001.
[10] Grasedyck, L., Hackbusch, W.: Construction and arithmetics of -matrices. Computing 70, 295–334 (2003). · Zbl 1030.65033
[11] Greengard, L., Rokhlin, V.: A new version of the fast multipole method for the Laplace in three dimensions. In: Acta Numerica 1997, pp. 229–269. Cambridge University Press 1997. · Zbl 0889.65115
[12] Hackbusch, W.: A sparse matrix arithmetic based on -matrices. Part I: Introduction to -matrices. Computing, 62, 89–108 (1999). · Zbl 0927.65063
[13] Hackbusch, W., Khoromskij, B.: A sparse matrix arithmetic based on -matrices. Part II: Application to multi-dimensional problems. Computing 64, 21–47 (2000). · Zbl 0962.65029
[14] Hackbusch, W., Khoromskij, B., Kriemann, R.: Hierarchical matrices based on a weak admissibility criterion. Technical Report2, Max-Planck Institute for Mathematics in the Sciences, 2003. Computing 73, 207–243 (2004). · Zbl 1063.65035
[15] Hackbusch, W., Nowak, Z. P.: On the fast matrix multiplication in the boundary element method by panel clustering. Numer. Math. 54, 463–491 (1989). · Zbl 0641.65038
[16] Lintner, M.: The eigenvalue problem for the 2d Laplacian in -matrix arithmetic and application to the heat and wave equation. Computing 72, 293–323 (2004). · Zbl 1055.65049
[17] Rokhlin, V.: Rapid solution of integral equations of classical potential theory. J. Comput. Phys. 60, 187–207 (1985). · Zbl 0629.65122
[18] Saad, Y., Schultz, M.: GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Statist. Comput. 7, 856–869 (1986). · Zbl 0599.65018
[19] Sauter, S. A.: Variable order panel clustering (extended version). Technical Report52, Max-Planck-Institute for Mathematik, Leipzig, Germany, 1999.
[20] Tyrtyshnikov, E.: Incomplete cross approximation in the mosaic-skeleton method. Computing 64, 367–380 (2000). · Zbl 0964.65048
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.