English

Block-encoding dense and full-rank kernels using hierarchical matrices: applications in quantum numerical linear algebra

Quantum Physics 2022-12-14 v4 Numerical Analysis Numerical Analysis

Abstract

Many quantum algorithms for numerical linear algebra assume black-box access to a block-encoding of the matrix of interest, which is a strong assumption when the matrix is not sparse. Kernel matrices, which arise from discretizing a kernel function k(x,x)k(x,x'), have a variety of applications in mathematics and engineering. They are generally dense and full-rank. Classically, the celebrated fast multipole method performs matrix multiplication on kernel matrices of dimension NN in time almost linear in NN by using the linear algebraic framework of hierarchical matrices. In light of this success, we propose a block-encoding scheme of the hierarchical matrix structure on a quantum computer. When applied to many physical kernel matrices, our method can improve the runtime of solving quantum linear systems of dimension NN to O(κpolylog(Nε))O(\kappa \operatorname{polylog}(\frac{N}{\varepsilon})), where κ\kappa and ε\varepsilon are the condition number and error bound of the matrix operation. This runtime is near-optimal and, in terms of NN, exponentially improves over prior quantum linear systems algorithms in the case of dense and full-rank kernel matrices. We discuss possible applications of our methodology in solving integral equations and accelerating computations in N-body problems.

Keywords

Cite

@article{arxiv.2201.11329,
  title  = {Block-encoding dense and full-rank kernels using hierarchical matrices: applications in quantum numerical linear algebra},
  author = {Quynh T. Nguyen and Bobak T. Kiani and Seth Lloyd},
  journal= {arXiv preprint arXiv:2201.11329},
  year   = {2022}
}

Comments

Added affiliations and acknowledgments

R2 v1 2026-06-24T09:04:54.580Z