Andrew M. Bradley
Abstract:We describe a new method to map the requested error tolerance on an H-matrix approximation to the block error tolerances. Numerical experiments show that the method produces more efficient approximations than the standard method for kernels having singularity order greater than one, often by factors of 1.5 to 5 and at a lower computational cost.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: how to effectively map the error tolerance of the H - matrix (H - matrix) requested by the user to the error tolerances of each sub - block, thereby improving computational efficiency while ensuring accuracy. Specifically, the author proposes a new method to achieve this mapping, so that for kernel functions with a higher singularity order, more efficient approximations can be generated, and better compression effects can be achieved at a lower computational cost.
### Problem Background
The H - matrix is a technique for efficiently approximating fully dense matrices, and is especially suitable for singular kernel functions that occur in problems such as integral equations or interactions between particles. The H - matrix divides the matrix through a hierarchical structure and uses low - rank approximation (LRA) to represent each sub - block, thereby reducing storage and computational complexity.
### Core Challenges
1. **Error Control**: Ensure that the error of the overall approximation does not exceed the user - specified tolerance.
2. **Computational Efficiency**: Minimize the consumption of computational resources, including time and space complexity, on the premise of meeting the error requirements.
### Paper Contributions
The author proposes an inequality - based method (matrix - wise relative - error method, MREM) to map the global error tolerance to the error tolerance of each sub - block. Compared with the traditional block - wise relative - error method (BREM), MREM shows significant advantages when dealing with kernel functions with a high singularity order:
- **Higher Compression Rate**: MREM can achieve a higher compression rate under the same error tolerance, and usually can achieve an improvement of 1.5 to 5 times.
- **Better Robustness**: MREM hardly produces a lower compression rate than BREM, and in most cases, the actual error is close to the user - requested tolerance.
### Numerical Experiments
The paper verifies the effectiveness of MREM through a series of numerical experiments, including particle systems with different geometric distributions and the fault - slip problem in geophysics. The experimental results show that when dealing with problems with a high singularity order, MREM can not only provide a higher compression rate, but also maintain good error control.
### Conclusion
The MREM method provides a more effective error tolerance mapping strategy for H - matrix approximation, especially when dealing with kernel functions with a high singularity order, it can significantly improve computational efficiency and compression rate. Therefore, for research fields involving high - singularity problems, such as geophysics and materials science, MREM is a recommended choice.
### Related Formulas
- **Frobenius Norm Relative Error**:
\[
\left\| E \right\|_F \leq \varepsilon \left\| B \right\|_F
\]
where \(E = B-\bar{B}\), and \(\bar{B}\) is the approximation matrix of \(B\).
- **Maximum Element Error**:
\[
\left\| E \right\|_{\max} \leq \varepsilon N^{-1} \left\| B \right\|_1
\]
- **MREM Error Condition**:
\[
\left\| E_i \right\|_F \leq \varepsilon \sqrt{\frac{m_i n_i}{N}} \left\| B \right\|_F
\]
These formulas ensure that while meeting the global error tolerance, the error tolerance of each sub - block can be effectively allocated, thereby improving the efficiency of the overall approximation.