Papers › Stochastic Gradient Descent for Gaussian Processes Done Right

Stochastic Gradient Descent for Gaussian Processes Done Right

31 Oct 2023arXiv:2310.20581archive 2025-07-28

Jihao Andreas Lin, Shreyas Padhy, Javier Antorán, Austin Tripp, Alexander Terenin, Csaba Szepesvári, José Miguel Hernández-Lobato, David Janz

As is well known, both sampling from the posterior and computing the mean of the posterior in Gaussian process regression reduces to solving a large linear system of equations. We study the use of stochastic gradient descent for solving this linear system, and show that when \emph{done right} -- by which we mean using specific insights from the optimisation and kernel communities -- stochastic gradient descent is highly effective. To that end, we introduce a particularly simple \emph{stochastic dual descent} algorithm, explain its design in an intuitive manner and illustrate the design choices through a series of ablation studies. Further experiments demonstrate that our new method is highly competitive. In particular, our evaluations on the UCI regression tasks and on Bayesian optimisation set our approach apart from preconditioned conjugate gradients and variational Gaussian process approximations. Moreover, our method places Gaussian process regression on par with state-of-the-art graph neural networks for molecular binding affinity prediction.

PaperPDFCodeCode Syntology ran

In Syntology Open this paper in Syntology's Atlas, the map of the papers in Syntology's graph and their citations.

For agents, Syntology's MCP tool lists every function and class Syntology harvested from this paper and whether it ran (how to connect): get_harvested_code_for_paper(arxiv_id="2310.20581")

Code

Syntology Ran 15 of 19 code samples harvested from 1 repository linked to this paper; 4 have no recorded run. Of those that ran: 15 ran with no contract checked.

By repository: official repository: 19 samples from 1 repository, 15 ran. The run record, sample by sample. “Ran” means executed on a synthesized input, not that the code is correct or reproduces the paper.

cambridge-mlg/sgd-gp officialmentioned in papermentioned on GitHubjaxMIT report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

19 samples harvested; 15 ran; 0 honoured the contract we drafted; 4 have no recorded run. Read from Syntology's graph 2026-09-24; that is when this build read the record, not when the samples ran.

15ran
4unverified

Licence: 0 of the 19 samples are pointer only, meaning Syntology does not serve that copy's text. This page shows no code text for any sample; each one links to its file in the repository.

Harvested from cambridge-mlg/sgd-gp. “Ran” means the sample executed on a synthesized input. It does not mean the output is correct, and nothing here reproduces the paper's results. “Honoured” and “violated” refer to a contract Syntology drafted from the code itself; “our draft was wrong” and “fixture could not drive it” are failures of Syntology's instrument, not of the code.

Each sample ends with its code_sha256, Syntology's identity for that exact code. An agent fetches the stored sample with Syntology's MCP tool get_code(code_sha256="…") (how to connect); click an identity to copy that call.

Repository labels, per sample. official repository: The archive marks this repository official for the paper. named in the paper: The archive records that the paper mentions this repository; it is not marked official. community (archive-listed): In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper. found in paper text by Syntology: Syntology found this repository in the paper's own text; whether it is the authors' implementation is not asserted. community: Not in the archive's code links for this paper; a community repository Syntology harvested. Samples from a repository marked official are listed first. Licence labels name the repository's licence as recorded at harvest. “Pointer only” means Syntology does not serve that copy's text, for one of four reasons: no licence file was found; the licence was not identified; the licence is recorded as permissive but that copy's record is not marked cleared; or the licence is outside the permissive list Syntology serves text under (MIT, Apache-2.0, BSD and similar). Some licences outside that list permit redistribution, such as WTFPL, and GPL-3.0 under its conditions; they are simply not on the list. Hover a licence label for the reason. File links open the file on GitHub at the default branch, which may have changed since the harvest.

add_centroid cambridge-mlg/sgd-gp/scalable_gps/kmeans.py official repository ran MIT (permissive) · 256df157025464e0 · report
custom_cg cambridge-mlg/sgd-gp/scalable_gps/custom_cg.py official repository ran MIT (permissive) · e516cb910a18d1c4 · report
error cambridge-mlg/sgd-gp/scalable_gps/linear_model.py official repository ran MIT (permissive) · e0776f5f9f8da068 · report
fp_dicts_to_arr cambridge-mlg/sgd-gp/scalable_gps/fingerprint_utils.py official repository ran MIT (permissive) · 902854c4bdd7e801 · report
freeze_non_variational_params cambridge-mlg/sgd-gp/scalable_gps/SVGP.py official repository ran MIT (permissive) · 182f12ef71a39631 · report
generate_close_pair_dict cambridge-mlg/sgd-gp/scalable_gps/knn.py official repository ran MIT (permissive) · 5aff9dfc1a529f98 · report
get_concentrating_toy_sin_dataset cambridge-mlg/sgd-gp/scalable_gps/data.py official repository ran MIT (permissive) · c95de35c70b9ab83 · report
get_pruned_indices cambridge-mlg/sgd-gp/scalable_gps/knn.py official repository ran MIT (permissive) · bc4d5ad47891f330 · report
get_split_toy_sin_dataset cambridge-mlg/sgd-gp/scalable_gps/data.py official repository ran MIT (permissive) · 4eee6abddb7dac5c · report
hilbert_space_RMSE cambridge-mlg/sgd-gp/scalable_gps/eval_utils.py official repository ran MIT (permissive) · 2bda1745f6d5e206 · report
i_error cambridge-mlg/sgd-gp/scalable_gps/inducing_linear_model.py official repository ran MIT (permissive) · 776edbdfc8352a25 · report
kmeanspp cambridge-mlg/sgd-gp/scalable_gps/kmeans.py official repository ran MIT (permissive) · 8940eaaee4dd71b2 · report
pred_space_RMSE cambridge-mlg/sgd-gp/scalable_gps/eval_utils.py official repository ran MIT (permissive) · 5eca356948077f17 · report
subsample cambridge-mlg/sgd-gp/scalable_gps/data.py official repository ran MIT (permissive) · 5e56bc4f20bb7af5 · report
vector_quantize cambridge-mlg/sgd-gp/scalable_gps/kmeans.py official repository ran MIT (permissive) · f9ffe6a03a54afb0 · report
i_error_grad_sample cambridge-mlg/sgd-gp/scalable_gps/inducing_linear_model.py official repository unverified MIT (permissive) · 0f413768b7c2cd5a · report
i_regularizer cambridge-mlg/sgd-gp/scalable_gps/inducing_linear_model.py official repository unverified MIT (permissive) · c92ac444fbedc857 · report
regularizer cambridge-mlg/sgd-gp/scalable_gps/linear_model.py official repository unverified MIT (permissive) · 05b0c6b467ba7da1 · report
sample_from_qu cambridge-mlg/sgd-gp/scalable_gps/SVGP.py official repository unverified MIT (permissive) · e79e6d84397801b3 · report

Tasks

Bayesian OptimisationGaussian Processesregression

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Methods

Gaussian ProcessSET

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections