arXiv research

Bayesian quadrature uses probabilistic models for estimating intractable integrals.

Bayesian quadrature improves integration efficiency with invariant priors.

A new type of quadrature is developed. The Gaussian quadrature, for a given measure, finds optimal values of a function's argument (nodes) and the corresponding weights. In contrast, the Lebesgue quadrature developed in this paper, finds optimal values of function (value-nodes) and the corresponding weights. The Gaussi…

The aim of this article is to design a moment transformation for Student- t distributed random variables, which is able to account for the error in the numerically computed mean. We employ Student-t process quadrature, an instance of Bayesian quadrature, which allows us to treat the integral itself as a random variable…

Computation of moments of transformed random variables is a problem appearing in many engineering applications. The current methods for moment transformation are mostly based on the classical quadrature rules which cannot account for the approximation errors. Our aim is to design a method for moment transformation for …

New quadrature method using randomly pivoted Cholesky outperforms existing techniques.

Adaptive Bayesian quadrature (ABQ) is a powerful approach to numerical integration that empirically compares favorably with Monte Carlo integration on problems of medium dimensionality (where non-adaptive quadrature is not competitive). Its key ingredient is an acquisition function that changes as a function of previou…

Unified quadrature framework for large-scale kernel machines.

Bayesian quadrature improves integration on Riemannian manifolds.

Herding and kernel herding are deterministic methods of choosing samples which summarise a probability distribution. A related task is choosing samples for estimating integrals using Bayesian quadrature. We show that the criterion minimised when selecting samples in kernel herding is equivalent to the posterior varianc…

Herding and kernel herding are deterministic methods of choosing samples which summarise a probability distribution. A related task is choosing samples for estimating integrals using Bayesian quadrature. We show that the criterion minimised when selecting samples in kernel herding is equivalent to the posterior varianc…

There is renewed interest in formulating integration as an inference problem, motivated by obtaining a full distribution over numerical error that can be propagated through subsequent computation. Current methods, such as Bayesian Quadrature, demonstrate impressive empirical performance but lack theoretical analysis. A…

Parallelized Bayesian quadrature improves sample efficiency and inference.

This paper introduces repulsive Monte Carlo methods for computing the sliced Wasserstein distance.

In this study, a numerical quadrature for the generalized inverse Gaussian distribution is derived from the Gauss-Hermite quadrature by exploiting its relationship with the normal distribution. The proposed quadrature is not Gaussian, but it exactly integrates the polynomials of both positive and negative orders. Using…

Combines control variates and adaptive importance sampling for Monte Carlo integration.

Integration over non-negative integrands is a central problem in machine learning (e.g. for model averaging, (hyper-)parameter marginalisation, and computing posterior predictive distributions). Bayesian Quadrature is a probabilistic numerical integration technique that performs promisingly when compared to traditional…

In this paper we develope, in a geometric framework, a Hamilton-Jacobi Theory for general dynamical systems. Such a theory contains the classical theory for Hamiltonian systems on a cotangent bundle and recent developments in the framework of general symplectic, Poisson and almost-Poisson manifolds (including some appr…

Adaptive quadrature improves Bayesian inference through active learning.

This article is concerned with Gaussian process quadratures, which are numerical integration methods based on Gaussian process regression methods, and sigma-point methods, which are used in advanced non-linear Kalman filtering and smoothing algorithms. We show that many sigma-point methods can be interpreted as Gaussia…

Improved kernel herding algorithm for faster quadrature rule convergence.

Improved option pricing for SABR model using Gauss-Hermite quadrature.

TQ separates sampling and integration for high-dimensional integrals.

The sigma-point filters, such as the UKF, which exploit numerical quadrature to obtain an additional order of accuracy in the moment transformation step, are popular alternatives to the ubiquitous EKF. The classical quadrature rules used in the sigma-point filters are motivated via polynomial approximation of the integ…

The standard Kernel Quadrature method for numerical integration with random point sets (also called Bayesian Monte Carlo) is known to converge in root mean square error at a rate determined by the ratio $s/d$, where $s$ and $d$ encode the smoothness and dimension of the integrand. However, an empirical investigation re…

Study of Hamilton-Jacobi Theory with symmetries and integrability by quadratures.

Kernel-based quadrature rules are becoming important in machine learning and statistics, as they achieve super-$\sqrt{n}$ convergence rates in numerical integration, and thus provide alternatives to Monte Carlo integration in challenging settings where integrands are expensive to evaluate or where integrands are high d…

A new method calculates accurate SABR model option prices and deltas.

This work introduces a fixed-point optimization for variational inference.

Study geodesic curves on Heisenberg group, classify them, and compute first step of quadrature.

The paper improves probabilistic herding methods using Gibbs distributions.

New Fourier features improve high-precision approximation in large-scale problems.

An important application of Lebesgue integral quadrature arXiv:1807.06007 is developed. Given two random processes, $f(x)$ and $g(x)$, two generalized eigenvalue problems can be formulated and solved. In addition to obtaining two Lebesgue quadratures (for $f$ and $g$) from two eigenproblems, the projections of $f$- and…

The paper introduces new methods for Asian option pricing using Laguerre quadrature.

New theory extends Hamilton-Jacobi for contact systems, ensuring integrability.

Kernel quadrature improves CRPS estimation for probabilistic time-series forecasting.

The paper improves error bounds for Bayesian quadrature in noisy settings.

Bayesian neural networks speed up numerical integration.

We show that kernel-based quadrature rules for computing integrals can be seen as a special case of random feature expansions for positive definite kernels, for a particular decomposition that always exists for such kernels. We provide a theoretical analysis of the number of required samples for a given approximation e…

An infinitely wide model is a weighted integration $\int \varphi(x,v) d μ(v)$ of feature maps. This model excels at handling an infinite number of features, and thus it has been adopted to the theoretical study of deep learning. Kernel quadrature is a kernel-based numerical integration scheme developed for fast approxi…

Faster training of neural ODEs using Gauß-Legendre quadrature.

A new method optimizes Fourier pricing for multi-asset options using adaptive quadrature.

This paper provides a dictionary of closed-form kernel mean embeddings.

We propose and analyze numerical methods for the Heath-Jarrow-Morton (HJM) model. To construct the methods, we first discretize the infinite dimensional HJM equation in maturity time variable using quadrature rules for approximating the arbitrage-free drift. This results in a finite dimensional system of stochastic dif…

The paper addresses numerical integration issues in SV models, proposing a fast regime switching algorithm.

The paper studies symmetries and conservation laws of non-diagonalisable hydrodynamic systems.

New filters improve radar target inference in complex scenarios.

Bayesian Probabilistic Integration uses BART for high-dimensional, non-smooth functions.