Randomly pivoted Cholesky: Practical approximation of a kernel matrix with few entry evaluations
Not on the map. Why?
Abstract (arXiv)
The randomly pivoted partial Cholesky algorithm (RPCholesky) computes a factorized rank-k approximation of an N x N positive-semidefinite (psd) matrix. RPCholesky requires only (k + 1) N entry evaluations and O(k^2 N) additional arithmetic operations, and it can be implemented with just a few lines of code. The method is particularly useful for approximating a kernel matrix. This paper offers a thorough new investigation of the empirical and theoretical behavior of this fundamental algorithm. For matrix approximation problems that arise in scientific machine learning, experiments show that RPCholesky matches or beats the performance of alternative algorithms. Moreover, RPCholesky provably returns low-rank approximations that are nearly optimal. The simplicity, effectiveness, and robustness of RPCholesky strongly support its use in scientific computing and machine learning applications.
Key facts. Cited by 31 papers in the map (data through 7 Oct 2026). Builds on 19 of its 89 references that are in the map. First posted on arXiv 13 Jul 2022.
Builds on
19 of 89 references.
References · 19
- 2014OpenML: networked science in machine learningcited by 523
- 2015An Introduction to Matrix Concentration Inequalitiescited by 391
- 2020Finite Versus Infinite Neural Networks: an Empirical Studycited by 156
- 2020Kernel methods through the roof: handling billions of points efficientlycited by 88
- 2019Harnessing the Power of Infinitely Wide Deep Nets on Small-data Taskscited by 87
- 2017FALKON: An Optimal Large Scale Kernel Methodcited by 82
- 2017Time-lagged autoencoders: Deep learning of slow collective variables for molecular kineticscited by 64
- 2022A Framework and Benchmark for Deep Batch Active Learning for Regressioncited by 47
- 2020Determinantal Point Processes in Randomized Numerical Linear Algebracited by 41
- 2018On Fast Leverage Score Sampling and Optimal Learningcited by 40
- 2016Monte Carlo Markov Chain Algorithms for Sampling Strongly Rayleigh Distributions and Determinantal Point Processescited by 24
- 2017Fixed-Rank Approximation of a Positive-Semidefinite Matrix from Streaming Datacited by 22
7 more
- 2013Improving CUR Matrix Decomposition and the Nystr"om Approximation via Adaptive Samplingcited by 21
- 2019Exact sampling of determinantal point processes with sublinear time preprocessingcited by 19
- 2017Sublinear Time Low-Rank Approximation of Positive Semidefinite Matricescited by 18
- 2021Simple, Fast, and Flexible Framework for Matrix Completion with Infinite Width Neural Networkscited by 17
- 2018Distributed Adaptive Sampling for Kernel Matrix Approximationcited by 9
- 2020Diversity sampling is an implicit regularization for kernel methodscited by 4
- 2021Kernel approximation on algebraic varietiescited by 4
Built on by
Citing papers per year: 2022: 1; 2023: 5; 2024: 7; 2025: 11; 2026: 7.
31 papers.
Sort:
- 2023Sparse Cholesky Factorization for Solving Nonlinear PDEs via Gaussian Processescited by 31
- 2025xRFM: Accurate, scalable, and interpretable feature learning models for tabular datacited by 14
- 2023Toward Large Kernel Modelscited by 12
- 2023Kernel quadrature with randomly pivoted Choleskycited by 11
- 2023Robust, randomized preconditioning for kernel ridge regressioncited by 8
- 2022Multiresolution kernel matrix algebracited by 6
- 2025Data-Efficient Kernel Methods for Learning Differential Equations and Their Solution Operators: Algorithms and Error Analysiscited by 6
- 2025Operator Learning at Machine Precisioncited by 6
- 2026WildCat: Near-Linear Attention in Theory and Practicecited by 6
- 2024The fast committor machine: Interpretable prediction with kernelscited by 5
- 2024Sketchy Moment Matching: Toward Fast and Provable Data Selection for Finetuningcited by 5
- 2024Embrace rejection: Kernel matrix approximation by accelerated randomly pivoted Choleskycited by 5
19 more
- 2023A Preconditioned Interior Point Method for Support Vector Machines Using an ANOVA-Decomposition and NFFT-Based Matrix-Vector Productscited by 3
- 2024Supervised Kernel Thinningcited by 3
- 2025Kernel Density Machinescited by 2
- 2026Linearization Explains Fine-Tuning in Large Language Modelscited by 2
- 2024Randomly Pivoted Partial Cholesky: Random How?cited by 1
- 2024Joint Estimation of Conditional Mean and Covariance for Unbalanced Panelscited by 1
- 2025Accelerating Natural Gradient Descent for PINNs with Randomized Numerical Linear Algebracited by 1
- 2025Learning where to learn: Training data distribution optimization for scientific machine learningcited by 1
- 2025Diagonally-Weighted Generalized Method of Moments Estimation for Gaussian Mixture Modelingcited by 1
- 2025Inference on Variable Importance for Treatment Effect Heterogeneity: Shapley Values and Beyondcited by 1
- 2026GPTQ-intrinsic LoRA: A Near-optimal Algorithm for Low-precision Quantization with Low-rank Adaptationcited by 1
- 2024Variational Markov chain mixtures with automatic component selectioncited by 0
- 2025Metric Learning in an RKHScited by 0
- 2025Adaptive randomized pivoting and volume samplingcited by 0
- 2025A Streaming Sparse Cholesky Method for Derivative-Informed Gaussian Process Surrogates Within Digital Twin Applicationscited by 0
- 2026Communication-Avoiding Linear Algebraic Kernel K-Means on GPUscited by 0
- 2026Deep Convolutional Large-Margin ℓp-SVDD for Visual Anomaly Detectioncited by 0
- 2026Fast Length-Squared Sampling for Positive-Semidefinite Matricescited by 0
- 2026SuperPCA: subspace analysis and an efficient algorithm for high-dimensional PCAcited by 0
- 2026SuperPCA: subspace analysis and an efficient algorithm for high-dimensional PCA22 Sep 2026
- 2026Fast Length-Squared Sampling for Positive-Semidefinite Matrices12 Aug 2026
- 2026Deep Convolutional Large-Margin ℓp-SVDD for Visual Anomaly Detection24 Jul 2026
- 2026GPTQ-intrinsic LoRA: A Near-optimal Algorithm for Low-precision Quantization with Low-rank Adaptation31 May 2026
- 2026WildCat: Near-Linear Attention in Theory and Practice10 Feb 2026
- 2026Linearization Explains Fine-Tuning in Large Language Models9 Feb 2026
- 2026Communication-Avoiding Linear Algebraic Kernel K-Means on GPUs23 Jan 2026
- 2025Operator Learning at Machine Precision25 Nov 2025
- 2025A Streaming Sparse Cholesky Method for Derivative-Informed Gaussian Process Surrogates Within Digital Twin Applications1 Nov 2025
- 2025Inference on Variable Importance for Treatment Effect Heterogeneity: Shapley Values and Beyond21 Oct 2025
- 2025Adaptive randomized pivoting and volume sampling2 Oct 2025
- 2025xRFM: Accurate, scalable, and interpretable feature learning models for tabular data12 Aug 2025
19 more
- 2025Metric Learning in an RKHS6 Aug 2025
- 2025Diagonally-Weighted Generalized Method of Moments Estimation for Gaussian Mixture Modeling28 Jul 2025
- 2025Learning where to learn: Training data distribution optimization for scientific machine learning27 May 2025
- 2025Accelerating Natural Gradient Descent for PINNs with Randomized Numerical Linear Algebra16 May 2025
- 2025Kernel Density Machines30 Apr 2025
- 2025Data-Efficient Kernel Methods for Learning Differential Equations and Their Solution Operators: Algorithms and Error Analysis2 Mar 2025
- 2024Joint Estimation of Conditional Mean and Covariance for Unbalanced Panels29 Oct 2024
- 2024Supervised Kernel Thinning17 Oct 2024
- 2024Embrace rejection: Kernel matrix approximation by accelerated randomly pivoted Cholesky4 Oct 2024
- 2024Sketchy Moment Matching: Toward Fast and Provable Data Selection for Finetuning8 Jul 2024
- 2024Variational Markov chain mixtures with automatic component selection7 Jun 2024
- 2024The fast committor machine: Interpretable prediction with kernels16 May 2024
- 2024Randomly Pivoted Partial Cholesky: Random How?17 Apr 2024
- 2023A Preconditioned Interior Point Method for Support Vector Machines Using an ANOVA-Decomposition and NFFT-Based Matrix-Vector Products1 Dec 2023
- 2023Kernel quadrature with randomly pivoted Cholesky6 Jun 2023
- 2023Robust, randomized preconditioning for kernel ridge regression24 Apr 2023
- 2023Sparse Cholesky Factorization for Solving Nonlinear PDEs via Gaussian Processes3 Apr 2023
- 2023Toward Large Kernel Models6 Feb 2023
- 2022Multiresolution kernel matrix algebra21 Nov 2022
- 2023Sparse Cholesky Factorization for Solving Nonlinear PDEs via Gaussian Processescited by 31
- 2025xRFM: Accurate, scalable, and interpretable feature learning models for tabular datacited by 14
- 2023Toward Large Kernel Modelscited by 12
- 2023Kernel quadrature with randomly pivoted Choleskycited by 11
- 2023Robust, randomized preconditioning for kernel ridge regressioncited by 8
- 2022Multiresolution kernel matrix algebracited by 6
- 2025Data-Efficient Kernel Methods for Learning Differential Equations and Their Solution Operators: Algorithms and Error Analysiscited by 6
- 2025Operator Learning at Machine Precisioncited by 6
- 2026WildCat: Near-Linear Attention in Theory and Practicecited by 6
- 2024The fast committor machine: Interpretable prediction with kernelscited by 5
- 2024Sketchy Moment Matching: Toward Fast and Provable Data Selection for Finetuningcited by 5
- 2024Embrace rejection: Kernel matrix approximation by accelerated randomly pivoted Choleskycited by 5
12 most cited of 31 All 31 (JSON)
- 2026SuperPCA: subspace analysis and an efficient algorithm for high-dimensional PCA22 Sep 2026
- 2026Fast Length-Squared Sampling for Positive-Semidefinite Matrices12 Aug 2026
- 2026Deep Convolutional Large-Margin ℓp-SVDD for Visual Anomaly Detection24 Jul 2026
- 2026GPTQ-intrinsic LoRA: A Near-optimal Algorithm for Low-precision Quantization with Low-rank Adaptation31 May 2026
- 2026WildCat: Near-Linear Attention in Theory and Practice10 Feb 2026
- 2026Linearization Explains Fine-Tuning in Large Language Models9 Feb 2026
- 2026Communication-Avoiding Linear Algebraic Kernel K-Means on GPUs23 Jan 2026
- 2025Operator Learning at Machine Precision25 Nov 2025
- 2025A Streaming Sparse Cholesky Method for Derivative-Informed Gaussian Process Surrogates Within Digital Twin Applications1 Nov 2025
- 2025Inference on Variable Importance for Treatment Effect Heterogeneity: Shapley Values and Beyond21 Oct 2025
- 2025Adaptive randomized pivoting and volume sampling2 Oct 2025
- 2025xRFM: Accurate, scalable, and interpretable feature learning models for tabular data12 Aug 2025
12 newest of 31 All 31 (JSON)