Mixed precision thin SVD algorithms based on the Gram matrix
Numerical Analysis
2026-03-13 v1 Numerical Analysis
Abstract
In this work, we present a mixed precision algorithm that leverages the Gram matrix and Jacobi methods to compute the singular value decomposition (SVD) of tall-and-skinny matrices. By constructing the Gram matrix in higher precision and coupling it with a Jacobi algorithm, our theoretical analysis and numerical experiments both indicate that the singular values computed by this mixed precision thin SVD algorithm attain high relative accuracy. In practice, our mixed precision thin SVD algorithm yields speedups of over 10x on a single CPU and about 2x on distributed memory systems when compared with traditional thin SVD methods.
Cite
@article{arxiv.2603.11953,
title = {Mixed precision thin SVD algorithms based on the Gram matrix},
author = {Erin Carson and Yuxin Ma and Meiyue Shao},
journal= {arXiv preprint arXiv:2603.11953},
year = {2026}
}
Comments
20 pages, 4 figures