2015/10/23 by Steffen Börm, Börm, Steffen · 1 citation
Engineering · Physics and Astronomy · #35J05 #65F30 #65N38 #Electromagnetic Scattering and Analysis #FOS: Mathematics #Numerical Analysis (math.NA) #Sparse and Compressive Sensing Techniques #Synthetic Aperture Radar (SAR) Applications and Techniques
paper · pdf · doi:10.48550/arxiv.1510.07087
openalex publication_date 2015/10/23 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/01
Standard numerical algorithms like the fast multipole method or H-matrix schemes rely on low-rank approximations of the underlying kernel function. For high-frequency problems, the ranks grow rapidly as the mesh is refined, and standard techniques are no longer attractive. Directional compression techniques solve this problem by using decompositions based on plane waves. Taking advantage of hierarchical relations between these waves' directions, an efficient approximation is obtained. This paper is dedicated to directional H2-matrices that employ local low-rank approximations to handle directional representations efficiently. The key result is an algorithm that takes an arbitrary matrix and finds a quasi-optimal approximation of this matrix as a directional H2-matrix using a prescribed block tree. The algorithm can reach any given accuracy, and the approximation requires only O(n k + κ2 k2 log n) units of storage, where n is the matrix dimension, κ is the wave number, and k is the local rank. In particular, we have a complexity of O(n k) if κ is constant and O(n k2 log n) for high-frequency problems characterized by κ2 ∼ n. Since the algorithm can be applied to arbitrary matrices, it can serve as the foundation of fast techniques for constructing preconditioners.