33
7

Faster Kernel Matrix Algebra via Density Estimation

Abstract

We study fast algorithms for computing fundamental properties of a positive semidefinite kernel matrix KRn×nK \in \mathbb{R}^{n \times n} corresponding to nn points x1,,xnRdx_1,\ldots,x_n \in \mathbb{R}^d. In particular, we consider estimating the sum of kernel matrix entries, along with its top eigenvalue and eigenvector. We show that the sum of matrix entries can be estimated to 1+ϵ1+\epsilon relative error in time sublinearsublinear in nn and linear in dd for many popular kernels, including the Gaussian, exponential, and rational quadratic kernels. For these kernels, we also show that the top eigenvalue (and an approximate eigenvector) can be approximated to 1+ϵ1+\epsilon relative error in time subquadraticsubquadratic in nn and linear in dd. Our algorithms represent significant advances in the best known runtimes for these problems. They leverage the positive definiteness of the kernel matrix, along with a recent line of work on efficient kernel density estimation.

View on arXiv
Comments on this paper