2022/02/24 by Huang, Baihe, Song, Zhao, Weinstein, Omri +3
#Data Structures and Algorithms (cs.DS) #FOS: Computer and information sciences
paper · doi:10.48550/arxiv.2202.12329
The Fast Gaussian Transform (FGT) enables subquadratic-time multiplication of an n× n Gaussian kernel matrix Ki,j= exp ( - ‖ xi - xj ‖22 ) with an arbitrary vector h ∈ ℝn, where x1,…, xn ∈ ℝd are a set of fixed source points. This kernel plays a central role in machine learning and random feature maps. Nevertheless, in most modern data analysis applications, datasets are dynamically changing (yet often have low rank), and recomputing the FGT from scratch in (kernel-based) algorithms incurs a major computational overhead (\gtrsim n time for a single source update ∈ ℝd). These applications motivate a dynamic FGT algorithm, which maintains a dynamic set of sources under kernel-density estimation (KDE) queries in sublinear time while retaining Mat-Vec multiplication accuracy and speed. Assuming the dynamic data-points xi lie in a (possibly changing) k-dimensional subspace (k≤ d), our main result is an efficient dynamic FGT algorithm, supporting the following operations in logO(k)(n/ε) time: (1) Adding or deleting a source point, and (2) Estimating the ``kernel-density'' of a query point with respect to sources with ε additive accuracy. The core of the algorithm is a dynamic data structure for maintaining the projected ``interaction rank'' between source and target boxes, decoupled into finite truncation of Taylor and Hermite expansions.