Year: 2026
Type: preprint
Abstract: Many dense matrices in scientific computing have hierarchically low-rank off-diagonal structure but are accessible only through routines that apply the matrix and, in general, its transpose or adjoint to vectors, that is, through matrix--vector products (matvecs). We present a fully matrix-free algorithm that compresses such a matrix into Hierarchically Semi-Separable (HSS) form and factorizes it simultaneously with the off-diagonal ranks chosen adaptively from a prescribed tolerance. The method uses $\mathcal{O}(k)$ matvecs with each operator, where $k$ is the largest off-diagonal rank found. Building on fixed-rank block nullification, it samples the matrix in small random batches and enlarges each local basis until a randomized error estimate falls below the tolerance. Unitary nullspace updates give $\mathcal{O}(k^3)$ local dense work, the same asymptotic cost as when the ranks are known in advance. Because compression and elimination proceed together, no separate factorization step is required. Excluding the matvecs, the factorization uses $\mathcal{O}(N)$ arithmetic and storage when the off-diagonal ranks are bounded, and its inverse is applied in $\mathcal{O}(N)$ time. Experiments on Poisson and Helmholtz Schur complements, Matérn covariance matrices, and a complex Lippmann--Schwinger equation demonstrate the factorization as both a direct solver and a preconditioner on problems with up to millions of unknowns.
Source: arXiv (Cornell University)
Authors: Do, Elly, Chen, Chao
Institution: North Carolina State University
Language: English
Topic: Matrix Theory and Algorithms
Subfield: Computational Theory and Mathematics
Field: Computer Science
Domain: Physical Sciences
Open Access status: green
Awards: G1453181110, Track 4: Advanced CI Coordination Ecosystem: Monitoring and Measurement Services, G3562386289, Track 1: ACCESS Resource Allocations Marketplace and Platform Services (RAMPS), G4010157413, ACO: An Open CI Ecosystem to Advance Scientific Discovery (OpenCI), Track 3: COre National Ecosystem for CyberinfrasTructure (CONECT), G5367127336, G5464919783, Track 2: Customized Multi-tier Assistance, Training, and Computational Help (MATCH) for End User ACCESS to CI
OpenAlex ID: W7221770426