cs.CGSep 14, 2026

Low-Dimensional Embeddings for Gaussian Kernels on Manifolds

Authors: Soumik DuttaKunal Dutta

Abstract

The Gaussian kernel is a widely used similarity measure underlying kernel methods such as kernel PCA and spectral clustering, but computing Gaussian kernel distances for many pairs of points can be expensive. Using Random Fourier Features (RFF), Chen and Phillips [ALT 2017] showed that for points in a dd-dimensional Euclidean ball in RN{\mathbb R}^N, t=Ω((d/ε2)log(dR/ε))t=Ω((d/\varepsilon^2)\log(dR/\varepsilon)) features suffice to preserve all pairwise Gaussian kernel distances within a (1±ε)(1\pm\varepsilon) factor with high probability. We establish a uniform relative-error embedding theorem for the more general setting of an arbitrary positive-reach submanifold MRN\mathcal M\subset{\mathbb R}^N of intrinsic dimension dd. We show that t=O((d/ε2)log(vol(M)2N2d/(vol(B1d(0))2rch(M)2dε2d+1δ)))t=O((d/\varepsilon^2)\log(\operatorname{vol}(\mathcal M)^2N^{2d}/(\operatorname{vol}(B_1^d(0))^2\operatorname{rch}(\mathcal M)^{2d}\varepsilon^{2d+1}δ))), or approximately O((d2/ε2)(logN+log(1/(εδ))))O((d^2/\varepsilon^2)(\log N+\log(1/(\varepsilonδ)))), RFFs suffice, with probability 1δ1-δ, to preserve the Gaussian kernel distance between every pair of manifold points up to relative error ε\varepsilon. Thus the bound depends only logarithmically on the ambient dimension and on manifold parameters such as volume and reach, while retaining the 1/ε21/\varepsilon^2 Euclidean rate. We also prove a topological consequence: under the same RFF embedding, persistent homology is preserved in the sense that weighted Cech and Rips filtrations built from Gaussian kernel power distance are (1±ε)(1\pm\varepsilon_\star)-interleaved, where ε\varepsilon_\star accounts for both distance distortion and kernel-weight approximation.

Explore similar work

CardsList