English

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

R2 v1 2026-07-01T11:16:46.827Z