English

Everything is Vecchia: Unifying low-rank and sparse inverse Cholesky approximations

Numerical Analysis 2026-03-09 v1 Numerical Analysis

Abstract

The partial pivoted Cholesky approximation accurately represents matrices that are close to being low-rank. Meanwhile, the Vecchia approximation accurately represents matrices with inverse Cholesky factors that are close to being sparse. What happens if a partial Cholesky approximation is combined with a Vecchia approximation of the residual? This paper shows how the sum is exactly a Vecchia approximation of the original matrix with an augmented sparsity pattern. Thus, Vecchia approximations subsume a class of existing matrix approximations and have broad applicability.

Keywords

Cite

@article{arxiv.2603.05709,
  title  = {Everything is Vecchia: Unifying low-rank and sparse inverse Cholesky approximations},
  author = {Eagan Kaminetz and Robert J. Webber},
  journal= {arXiv preprint arXiv:2603.05709},
  year   = {2026}
}

Comments

32 pages, 9 figures