Search papers, labs, and topics across Lattice.
This paper introduces a boundary-row Divide-and-Conquer (D&C) algorithm for eigenvalue-only computation of symmetric tridiagonal matrices, addressing the memory bottleneck of conventional D&C methods. By propagating only boundary rows/columns during the conquer phase, the algorithm achieves linear space complexity, a significant reduction from the quadratic space required by traditional approaches. The authors demonstrate the algorithm's correctness, stability, and performance advantages over QR and standard D&C routines on both CPUs and GPUs.
Solving eigenvalue problems just got a whole lot faster and cheaper: a new divide-and-conquer algorithm slashes memory requirements from quadratic to linear, finally making it competitive with QR-based methods.
Divide and Conquer (D&C) is a widely used algorithmic strategy for symmetric eigenvalue decomposition. Its natural parallelism makes D&C attractive on modern multicore CPUs and GPUs, but existing eigenvalue-only routines often default to QR-based methods because conventional D&C still materializes or replays large transformation matrices during the conquer phase. This paper proposes a boundary-row D&C algorithm for eigenvalue-only computation. The key observation is that the conquer phase only needs selected boundary rows/columns rather than the full accumulated eigenvector matrix. By propagating these boundary rows directly through the recursion, the proposed algorithm reduces the memory requirement from quadratic to linear space while also eliminating unnecessary matrix-vector work in the conventional lazy-replay formulation. We provide the algorithm, its time and space complexity analysis, correctness and stability arguments, optimized CPU and GPU implementations, and an evaluation against QR and D&C routines in standard numerical libraries.
