Kernel Operations on the GPU, with Autodiff, without Memory Overflows
Not on the map. Why?
Abstract (arXiv)
The KeOps library provides a fast and memory-efficient GPU support for tensors whose entries are given by a mathematical formula, such as kernel and distance matrices. KeOps alleviates the major bottleneck of tensor-centric libraries for kernel and geometric applications: memory consumption. It also supports automatic differentiation and outperforms standard GPU baselines, including PyTorch CUDA tensors or the Halide and TVM libraries. KeOps combines optimized C++/CUDA schemes with binders for high-level languages: Python (Numpy and PyTorch), Matlab and GNU R. As a result, high-level "quadratic" codes can now scale up to large data sets with millions of samples processed in seconds. KeOps brings graphics-like performances for kernel methods and is freely available on standard repositories (PyPi, CRAN). To showcase its versatility, we provide tutorials in a wide range of settings online at urlwww.kernel-operations.io.
Key facts. Cited by 101 papers in the map (data through 7 Oct 2026). First posted on arXiv 27 Mar 2020.
Builds on
No reference list yet.
Built on by
Citing papers per year: 2021: 20; 2022: 25; 2023: 14; 2024: 13; 2025: 16; 2026: 13.
101 papers.
Sort:
By area:
- 2022FlashAttention: Fast and Memory-Efficient Exact Attention with IO-Awarenesscited by 3,611
- 2022Diagonal State Spaces are as Effective as Structured State Spacescited by 350
- 2021Sinkformers: Transformers with Doubly Stochastic Attentioncited by 136
- 2022Liquid Structural State-Space Modelscited by 134
- 2021Shape registration in the time of transformerscited by 51
- 2021Accurate Point Cloud Registration with Robust Optimal Transportcited by 38
- 2021On UMAP's true loss functioncited by 32
- 2021Towards a mathematical theory of trajectory inferencecited by 31
- 2022From t-SNE to UMAP with contrastive learningcited by 28
- 2022Posterior and Computational Uncertainty in Gaussian Processescited by 23
- 2022Latent Graph Inference using Product Manifoldscited by 23
- 2022Simplifying and Understanding State Space Models with Diagonal Linear RNNscited by 23
38 more
- 2021Preconditioning for Scalable Gaussian Process Hyperparameter Optimizationcited by 22
- 2021Tighter Bounds on the Log Marginal Likelihood of Gaussian Process Regression Using Conjugate Gradientscited by 20
- 2022Elastic shape analysis of surfaces with second-order Sobolev metrics: a comprehensive numerical frameworkcited by 18
- 2021Decoupled Exploration and Exploitation Policies for Sample-Efficient Reinforcement Learningcited by 14
- 2024Memory-Scalable and Simplified Functional Map Learningcited by 14
- 2021ResNet-LDDMM: Advancing the LDDMM Framework using Deep Residual Networkscited by 13
- 2021Fast and Scalable Optimal Transport for Brain Tractogramscited by 11
- 2022Efficient Hyperparameter Tuning for Large Scale Kernel Ridge Regressioncited by 11
- 2022Fast Nonlinear Vector Quantile Regressioncited by 11
- 2023qPOTS: Efficient batch multiobjective Bayesian optimization via Pareto optimal Thompson samplingcited by 11
- 2023Large-Scale Gaussian Processes via Alternating Projectioncited by 11
- 2022Low-Precision Arithmetic for Fast Gaussian Processescited by 10
- 2021SKIing on Simplices: Kernel Interpolation on the Permutohedral Lattice for Scalable Gaussian Processescited by 9
- 2021Towards Efficient Point Cloud Graph Neural Networks Through Architectural Simplificationcited by 9
- 2024Conditional simulation via entropic optimal transport: Toward non-parametric estimation of conditional Brenier mapscited by 9
- 2025Position: Curvature Matrices Should Be Democratized via Linear Operatorscited by 9
- 2025ELECTRA: A Cartesian Network for 3D Charge Density Prediction with Floating Orbitalscited by 9
- 2022Single-Photon Structured Lightcited by 8
- 20223D Shape Sequence of Human Comparison and Classification using Current and Varifoldscited by 8
- 2023CoLA: Exploiting Compositional Structure for Automatic and Efficient Numerical Linear Algebracited by 7
- 2024Advancing Deep Active Learning & Data Subset Selection: Unifying Principles with Information-Theory Intuitionscited by 7
- 2024Fast Kernel Summation in High Dimensions via Slicing and Fourier Transformscited by 7
- 2024Fast Summation of Radial Kernels via QMC Slicingcited by 7
- 2026FlashSinkhorn: IO-Aware Entropic Optimal Transport on GPUcited by 7
- 2022Unsupervised Voice-Face Representation Learning by Cross-Modal Prototype Contrastcited by 6
- 2023Toward Mesh-Invariant 3D Generative Deep Learning with Geometric Measurescited by 6
- 2023Efficient and Accurate Optimal Transport with Mirror Descent and Conjugate Gradientscited by 6
- 2023Neural Snowflakes: Universal Latent Graph Inference via Trainable Latent Geometriescited by 6
- 2022Self-Supervised Correspondence Estimation via Multiview Registrationcited by 5
- 2023Label Shift Quantification with Robustness Guarantees via Distribution Feature Matchingcited by 5
- 2025A Truncated Newton Method for Optimal Transportcited by 5
- 2021Unrolled Variational Bayesian Algorithm for Image Blind Deconvolutioncited by 4
- 2021Partial Matching in the Space of Varifoldscited by 3
- 2021Low-rank statistical finite elements for scalable model-data synthesiscited by 3
- 2021Nonlinearities in Steerable SO(2)-Equivariant CNNscited by 3
- 2022Neural Bregman Divergences for Distance Learningcited by 3
- 2022Memory Safe Computations with XLA Compilercited by 3
- 2022GNPM: Geometric-Aware Neural Parametric Modelscited by 3
50 most cited of 101 All 101 (JSON)
- 2026A Unified Dual Method for Matching Problems30 Sep 2026
- 2026Quasi Linear Kernel Attention with Infinite Capacity28 Sep 2026
- 2026DanLing NestedTensor: Composable Multi-Ragged Tensors for Deep Learning24 Sep 2026
- 2026Fast Gauss Sums via Flash Attention4 Sep 2026
- 2026Gleam: Adaptive Network-Efficient CUDA API Remoting for Cross-Device GPU Sharing over LANs25 Jul 2026
- 2026Sobek: Streaming Equivariant Tensor Product Convolutions20 Jul 2026
- 2026SIKA-GP: Accelerating Gaussian Process Inference with Sparse Inducing Kernel Approximations for Bayesian Deep Learning26 May 2026
- 2026Lightweight Gaussian Process Inference in C++ on Metal and CUDA18 May 2026
- 2026Differentiable latent structure discovery for interpretable forecasting in clinical time series30 Apr 2026
- 2026Fast Log-Domain Sinkhorn Optimal Transport with Warp-Level GPU Reductions4 Apr 2026
- 2026k-Maximum Inner Product Attention for Graph Transformers and the Expressive Power of GraphGPS4 Apr 2026
- 2026Flash-SD-KDE: Accelerating SD-KDE with Tensor Cores10 Feb 2026
38 more
- 2026FlashSinkhorn: IO-Aware Entropic Optimal Transport on GPU3 Feb 2026
- 2025From Theory to Throughput: CUDA-Optimized APML for Large-Batch 3D Learning17 Dec 2025
- 2025Differentiable Expectation-Maximisation and Applications to Gaussian Mixture Model Optimal Transport2 Sep 2025
- 20253D Human Mesh Estimation from Single View RGBD11 Aug 2025
- 2025Fast Gaussian process inference by exact Mat'ern kernel decomposition3 Aug 2025
- 2025Sinkhorn Normalization of Diffusion Kernels8 Jul 2025
- 2025ClusterUCB: Efficient Gradient-Based Data Selection for Targeted Fine-Tuning of LLMs12 Jun 2025
- 2025A new approach for image segmentation based on diffeomorphic registration and gradient fields11 Jun 2025
- 2025Unsupervised Learning for Optimal Transport plan prediction between unbalanced graphs21 May 2025
- 2025Turbocharging Gaussian Process Inference with Approximate Sketch-and-Project19 May 2025
- 2025A Truncated Newton Method for Optimal Transport2 Apr 2025
- 2025ELECTRA: A Cartesian Network for 3D Charge Density Prediction with Floating Orbitals11 Mar 2025
- 2025Dimensionality reduction for homological stability and global structure preservation5 Mar 2025
- 2025ARC-Flow : Articulated, Resolution-Agnostic, Correspondence-Free Matching and Interpolation of 3D Shapes Under Flow Fields4 Mar 2025
- 2025HiGP: A high-performance Python package for Gaussian Process4 Mar 2025
- 2025Position: Curvature Matrices Should Be Democratized via Linear Operators31 Jan 2025
- 2025The Artificial Scientist -- in-transit Machine Learning of Plasma Simulations6 Jan 2025
- 2024Conditional simulation via entropic optimal transport: Toward non-parametric estimation of conditional Brenier maps11 Nov 2024
- 2024High-Dimensional Gaussian Process Regression with Soft Kernel Interpolation28 Oct 2024
- 2024Batched Energy-Entropy acquisition for Bayesian Optimization11 Oct 2024
- 2024Stochastic Kernel Regularisation Improves Generalisation in Deep Kernel Machines8 Oct 2024
- 2024Fast Summation of Radial Kernels via QMC Slicing2 Oct 2024
- 2024Mat'ern Kernels for Tunable Implicit Surface Reconstruction23 Sep 2024
- 2024Learning to Embed Distributions via Maximum Kernel Entropy1 Aug 2024
- 2024Neural varifolds: an aggregate representation for quantifying the geometry of point clouds5 Jul 2024
- 2024Gradients of Functions of Large Matrices27 May 2024
- 2024Fast Evaluation of Additive Kernels: Feature Arrangement, Fourier Methods, and Kernel Derivatives26 Apr 2024
- 2024Memory-Scalable and Simplified Functional Map Learning30 Mar 2024
- 2024Fast Kernel Summation in High Dimensions via Slicing and Fourier Transforms16 Jan 2024
- 2024Advancing Deep Active Learning & Data Subset Selection: Unifying Principles with Information-Theory Intuitions9 Jan 2024
- 2023On the Nystrom Approximation for Preconditioning in Kernel Machines6 Dec 2023
- 2023Persistent Homology for High-dimensional Data Based on Spectral Methods6 Nov 2023
- 2023Accelerating Non-Conjugate Gaussian Processes By Trading Off Computation For Uncertainty31 Oct 2023
- 2023Large-Scale Gaussian Processes via Alternating Projection26 Oct 2023
- 2023qPOTS: Efficient batch multiobjective Bayesian optimization via Pareto optimal Thompson sampling24 Oct 2023
- 2023Neural Snowflakes: Universal Latent Graph Inference via Trainable Latent Geometries23 Oct 2023
- 2023Compressed online Sinkhorn8 Oct 2023
- 2023Tracking Control for a Spherical Pendulum via Curriculum Reinforcement Learning25 Sep 2023
50 newest of 101 All 101 (JSON)
- 2022FlashAttention: Fast and Memory-Efficient Exact Attention with IO-Awarenesscited by 3,611
- 2021Sinkformers: Transformers with Doubly Stochastic Attentioncited by 136
- 2022Liquid Structural State-Space Modelscited by 134
- 2021Shape registration in the time of transformerscited by 51
- 2021Accurate Point Cloud Registration with Robust Optimal Transportcited by 38
- 2021On UMAP's true loss functioncited by 32
- 2021Towards a mathematical theory of trajectory inferencecited by 31
- 2022From t-SNE to UMAP with contrastive learningcited by 28
- 2022Posterior and Computational Uncertainty in Gaussian Processescited by 23
- 2022Latent Graph Inference using Product Manifoldscited by 23
- 2022Simplifying and Understanding State Space Models with Diagonal Linear RNNscited by 23
- 2021Preconditioning for Scalable Gaussian Process Hyperparameter Optimizationcited by 22
12 most cited of 100 All 100 (JSON)
- 2026A Unified Dual Method for Matching Problems30 Sep 2026
- 2026Quasi Linear Kernel Attention with Infinite Capacity28 Sep 2026
- 2026DanLing NestedTensor: Composable Multi-Ragged Tensors for Deep Learning24 Sep 2026
- 2026Fast Gauss Sums via Flash Attention4 Sep 2026
- 2026Gleam: Adaptive Network-Efficient CUDA API Remoting for Cross-Device GPU Sharing over LANs25 Jul 2026
- 2026Sobek: Streaming Equivariant Tensor Product Convolutions20 Jul 2026
- 2026SIKA-GP: Accelerating Gaussian Process Inference with Sparse Inducing Kernel Approximations for Bayesian Deep Learning26 May 2026
- 2026Lightweight Gaussian Process Inference in C++ on Metal and CUDA18 May 2026
- 2026Differentiable latent structure discovery for interpretable forecasting in clinical time series30 Apr 2026
- 2026Fast Log-Domain Sinkhorn Optimal Transport with Warp-Level GPU Reductions4 Apr 2026
- 2026k-Maximum Inner Product Attention for Graph Transformers and the Expressive Power of GraphGPS4 Apr 2026
- 2026Flash-SD-KDE: Accelerating SD-KDE with Tensor Cores10 Feb 2026
12 newest of 100 All 100 (JSON)