Abstract:The discretization of non-local operators, e.g., solution operators of partial differential equations or integral operators, leads to large densely populated matrices. $\mathcal{H}^2$-matrices take advantage of local low-rank structures in these matrices to provide an efficient data-sparse approximation that allows us to handle large matrices efficiently, e.g., to reduce the storage requirements to $\mathcal{O}(n k)$ for $n$-dimensional matrices with local rank $k$, and to reduce the complexity of the matrix-vector multiplication to $\mathcal{O}(n k)$ operations.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: when dealing with large - scale dense matrices, how to efficiently perform the multiplication operation of H2 matrices and ensure that the sub - matrices of the result have a controllable block - relative - error estimate. Specifically, the author proposes a new algorithm for approximating the product of two H2 matrices and guaranteeing the relative error of each resulting sub - matrix.
### Background and Problem Description
When numerically solving partial differential equations or integral equations, the discretization scheme usually generates large - scale dense matrices. The storage and computational complexity of these matrices are very high and difficult to handle directly. The H2 matrix provides an effective data - sparse approximation by exploiting the local low - rank structure, so that these large - scale matrices can be efficiently processed. For example, the H2 matrix can reduce the storage requirement of an n - dimensional matrix from O(n^2) to O(nk), and reduce the complexity of matrix - vector multiplication from O(n^2) to O(nk).
However, in order to perform more complex operations (such as constructing efficient pre - conditioners or evaluating matrix functions), we need algorithms that can approximate the results of operations between H2 matrices with controllable precision. In particular, for the multiplication of H2 matrices, the existing methods are either insufficient in precision or low in computational efficiency.
### Main Contributions of the Paper
In this paper, the author introduces a new algorithm for approximating the product of two H2 matrices and is able to provide a strict block - relative - error estimate for each sub - matrix in the result. The main features of this algorithm are as follows:
1. **Exact Representation of Intermediate Results**: The algorithm first prepares an exact representation of the product of two H2 matrices, and then applies rank - revealing factorization to obtain an approximate result in the form of an H2 matrix.
2. **Strict Error Control**: Compared with previous algorithms, a major advantage of the new algorithm is that it can provide a strict relative - error estimate for each sub - matrix appearing in the result.
3. **Efficiency**: By using a special tree structure to represent the intermediate steps of the product, the algorithm can complete the calculation within a time complexity of O(nk^2 log n) and can guarantee the given precision requirements.
### Summary of Mathematical Formulas
- Definition of H2 matrix:
\[
G|_{\hat{t} \times \hat{s}} \approx V_t S_{ts} W_s^*,
\]
where \(V_t \in \mathbb{R}^{\hat{t} \times k}\), \(S_{ts} \in \mathbb{R}^{k \times k}\), \(W_s \in \mathbb{R}^{\hat{s} \times k}\).
- Definition of Accumulator:
\[
\text{accu}(A_{tr}) := A_{tr} W_r^* + V_t B_{tr}^* + N_{tr} + \sum_{(s, X, Y) \in P_{tr}} XY.
\]
- Definition of Basis Tree:
\[
B_t =
\begin{cases}
(C, N, M) & \text{if } t \text{ is a leaf in } T_I, \\
(C, M) & \text{if } t \text{ is not a leaf in } T_I, \\
(C, M, \alpha_1, \ldots, \alpha_m) & \text{if } t \text{ has children } t_1, \ldots, t_m.
\end{cases}
\]
Through these methods, the new algorithm proposed in the paper not only improves the computational efficiency of H2 - matrix multiplication, but also ensures the high - precision of the result and strict error control.