English

Faster Kernel Interpolation for Gaussian Processes

Machine Learning 2021-08-16 v2 Artificial Intelligence

Abstract

A key challenge in scaling Gaussian Process (GP) regression to massive datasets is that exact inference requires computation with a dense n x n kernel matrix, where n is the number of data points. Significant work focuses on approximating the kernel matrix via interpolation using a smaller set of m inducing points. Structured kernel interpolation (SKI) is among the most scalable methods: by placing inducing points on a dense grid and using structured matrix algebra, SKI achieves per-iteration time of O(n + m log m) for approximate inference. This linear scaling in n enables inference for very large data sets; however the cost is per-iteration, which remains a limitation for extremely large n. We show that the SKI per-iteration time can be reduced to O(m log m) after a single O(n) time precomputation step by reframing SKI as solving a natural Bayesian linear regression problem with a fixed set of m compact basis functions. With per-iteration complexity independent of the dataset size n for a fixed grid, our method scales to truly massive data sets. We demonstrate speedups in practice for a wide range of m and n and apply the method to GP inference on a three-dimensional weather radar dataset with over 100 million points.

Keywords

Cite

@article{arxiv.2101.11751,
  title  = {Faster Kernel Interpolation for Gaussian Processes},
  author = {Mohit Yadav and Daniel Sheldon and Cameron Musco},
  journal= {arXiv preprint arXiv:2101.11751},
  year   = {2021}
}

Comments

To appear, Artificial Intelligence and Statistics (AISTATS) 2021