Research map / Probabilistic and scientific modeling / Uncertainty and Gaussian processes
Conditioning Gaussian Processes on Almost Anything
Read the paperCheck your idea against 108,967 papers
AI-extracted summary. Written by AI from the paper's full text. Numbers and claims can be wrong or missing context, so check the paper before you cite or build on it.
Shows GPs are an exact closed-form linear diffusion/flow-matching model and uses a Monte Carlo guided probability-flow ODE (FlowGP, with whitening) to condition GPs on non-linear, non-Gaussian constraints, PDEs and LLM-derived text.
Problem
Exact GP inference is limited to linear-Gaussian conditioning; non-Gaussian likelihoods and non-linear constraints (physics, text) require approximate inference schemes that are bespoke and hard to generalise.
Approach
Recast GP predictive sampling as a linear diffusion/probability-flow ODE with closed-form Gaussian dynamics, add a guidance term estimated by self-normalised importance sampling using only point-wise likelihood evaluations, and whiten the covariance to reduce stiffness and Wasserstein-2 transport cost.
Key results
- On the Allen-Cahn benchmark, FlowGP with whitening reaches RMSE 0.13 and NLPD -0.83 versus PHYSS 0.17/1.69 and AUTOIP 0.17/-0.29, at roughly 4 s versus 1e4 s runtime.
- On the non-linear pendulum benchmark, FlowGP-W matches PHYSS/AUTOIP RMSE of 0.05 with NLPD -1.05 (best) and ~4 s runtime (about 25x lower cost).
- FlowGP generates 100 samples for an LLM text prompt in under 2 seconds, with the paper stating runtimes from milliseconds to under 5 seconds.
Surprising findings
- The paper reports its simple Monte Carlo guidance estimator outperforms DPS and MPGD even at S=1 samples.
- Without whitening, FlowGP on Allen-Cahn has high RMSE variance (0.38 with std 1.39), whereas whitening stabilises it.
Limitations
- Cubic scaling in the number of grid points makes large-scale problems intractable.
- Theoretical guarantees are strong only for log-concave likelihoods; non-log-concave likelihoods admit only worst-case exponential bounds.
- Kernel hyperparameters are fitted only on the linear-Gaussian component, which may bias results when the conditioning information is highly informative; possible importance weight collapse in high dimensions.
Experimental setup
Benchmarks: monotone bounded regression benchmark non-linear pendulum Allen-Cahn equation Burgers equation
Metrics: RMSE NLPD runtime
Baselines: PHYSS AUTOIP DPS MPGD
Builds on: probability flow ODE diffusion guidance (DPS) stochastic interpolants pathwise conditioning probit likelihood
Closest papers by meaning
- Computation-Aware Gaussian Processes: Model Selection And Linear-Time Inference (NeurIPS 2024 · similarity 0.88)
- Turbocharging Gaussian Process Inference with Approximate Sketch-and-Project (NeurIPS 2025 · similarity 0.88)
- Robust and Conjugate Gaussian Process Regression (ICML 2024 · similarity 0.87)
- Conditioning non-linear and infinite-dimensional diffusion processes (NeurIPS 2024 · similarity 0.87)
- Exact Inference for Continuous-Time Gaussian Process Dynamics (AAAI 2024 · similarity 0.87)
- Scalable Derivative Gaussian Processes via Exact Gradient Reduction (NeurIPS 2026 · similarity 0.87)
Writing a paper in this area? scime draft is a free LaTeX editor with AI included, real-time co-authors and a check of your references and venue rules before you submit. It also imports OpenAI Prism projects.
