New integrators for mechanical systems on Lie groups simplify based on group properties.
problem Designing numerical integrators for mechanical systems on Lie groups.
method Leverage retraction maps and Lie group properties to design structure-preserving integrators.
result Simplified design of integrators for Euler-Poincare and Lie-Poisson equations.
i-flow uses normalizing flows for high-dimensional integration and sampling.
problem High-dimensional integration in science and statistics.
method Normalizing flows for bijective mappings between distributions.
result i-flow outperforms other algorithms for high-dimensional correlated integrals.
The paper addresses numerical integration issues in SV models, proposing a fast regime switching algorithm.
problem Numerical integration challenges in SV models, especially with high precision and low computational time.
method Proposes a fast regime switching algorithm to determine when higher precision arithmetic is needed.
result Shows that numerical quadratures need to be carefully chosen based on model parameters and parameter values.
Develops integrators for nonholonomic systems on Lie groups.
problem Nonholonomic constraints on Lie groups.
method Using retraction maps and Hamel formulation.
result Structure-preserving numerical integrators for nonholonomic systems.
Bayesian neural networks speed up numerical integration.
problem Scalability of Bayesian quadrature methods.
method Bayesian Stein networks using neural networks and Laplace approximation.
result Orders of magnitude speed-up on benchmark functions and real-world problems.
Motivated by numerical integration on manifolds, we relate the algebraic properties of invariant connections to their geometric properties. Using this perspective, we generalize some classical results of Cartan and Nomizu to invariant connections on algebroids. This has fundamental consequences for the theory of numeri…
Neural dynamical systems are dynamical systems that are described at least in part by neural networks. The class of continuous-time neural dynamical systems must, however, be numerically integrated for simulation and learning. Here, we present a compact neural circuit for two common numerical integrators: the explicit …
This paper deals with the evaluation of double line integrals of the squared exponential covariance function. We propose a new approach in which the double integral is reduced to a single integral using the error function. This single integral is then computed with efficiently implemented numerical techniques. The perf…
The paper compares numerical schemes for nonholonomic systems using retraction maps.
problem Optimal control of nonholonomic systems with numerical approximations.
method Retraction maps used as seed for geometric integrators of Hamilton equations.
result Performance comparison of symplectic and non-symplectic integrators.
A contour integral method recently proposed by Weideman [IMA J. Numer. Anal., to appear] for integrating semi-discrete advection-diffusion PDEs, is extended for application to some of the important equations of mathematical finance. Using estimates for the numerical range of the spatial operator, optimal contour parame…
A research frontier has emerged in scientific computation, wherein numerical error is regarded as a source of epistemic uncertainty that can be modelled. This raises several statistical challenges, including the design of statistical methods that enable the coherent propagation of probabilities through a (possibly dete…
We discuss two numerical methods, based on a path integral approach described in a previous paper (I), for solving the stochastic equations underlying the financial markets: the Monte Carlo approach, and the Green function deterministic numerical method. Then, we apply the latter to some specific financial problems. In…
Bayesian quadrature uses probabilistic models for estimating intractable integrals.
problem Estimating intractable integrals in complex models.
method Probabilistic, model-based approach using Gaussian processes.
result Comprehensive review and systematic taxonomy of Bayesian quadrature methods.
Bayesian Probabilistic Integration uses BART for high-dimensional, non-smooth functions.
problem Bayesian quadrature's limitations in high-dimensional or non-smooth functions.
method Bayesian Additive Regression Trees (BART) priors for numerical integration.
result Explicit convergence rates can be obtained in various settings.
A new method simulates square-root processes efficiently.
problem Simulating square-root processes accurately and efficiently.
method Simulate the integrated square-root process instead of the square-root process itself.
result High precision with low number of time steps, and exact limiting Inverse Gaussian distributions.
Bayesian quadrature improves integration efficiency with invariant priors.
problem Efficient numerical integration with known structure.
method Invariance priors for bijective transformations in input domain.
result Superior performance in synthetic and real-world applications.
A new method for estimating uncertainties in neural ODEs without numerical integration.
problem Accurate estimation of predictive uncertainties in neural ODEs.
method Distributional Gradient Matching (DGM) algorithm that jointly trains a smoother and a dynamics model.
result Significantly more accurate predictions compared to traditional methods.
Geometric integrator preserves coadjoint orbits in dissipative systems.
problem Preserving coadjoint orbits in dissipative mechanical systems.
method Adapted discrete variational integrators for forced Euler-Poincaré and Lie-Poisson systems.
result Preserves coadjoint orbits exactly, improving over general-purpose methods.
Probabilistic solvers improve stability for stiff systems.
problem Performance penalties for small steps in stiff systems.
method Probabilistic exponential integrators that include fast linear dynamics in the prior.
result Proven L-stability and probabilistic error accounting.
An efficient adaptive direct numerical integration (DNI) algorithm is developed for computing high quantiles and conditional Value at Risk (CVaR) of compound distributions using characteristic functions. A key innovation of the numerical scheme is an effective tail integration approximation that reduces the truncation …
This paper tackles Bayesian system identification with probabilistic numerical methods.
problem Accurately modeling nonlinear dynamic systems from noisy data.
method Probabilistic Sequential Monte Carlo (SMC) combined with probabilistic numerical integration.
result Efficient identification of latent states and system parameters from noisy measurements.
Integration of the form ∫a∞f(x)w(x)dx, where w(x) is either sin(ωx) or cos(ωx), is widely encountered in many engineering and scientific applications, such as those involving Fourier or Laplace transforms. Often such integrals are approximated by a numerical integration…
We consider numerical integrators of ODEs on homogeneous spaces (spheres, affine spaces, hyperbolic spaces). Homogeneous spaces are equipped with a built-in symmetry. A numerical integrator respects this symmetry if it is equivariant. One obtains homogeneous space integrators by combining a Lie group integrator with an…
DNA-SE uses deep learning to solve semiparametric problems efficiently.
problem Solving semiparametric integral equations in high dimensions.
method Formulates semiparametric estimation as a bi-level optimization problem and uses DNN to approximate solutions.
result Demonstrates numerical and statistical advantages over traditional methods.
Gradient flow preserves speed for integral Menger curvature curves.
problem Optimizing curves with integral Menger curvature constraints.
method Projected Sobolev gradient flow in Hilbert space.
result Long-time existence and C1,1-bounds for the flow. The paper efficiently solves a complex option valuation equation for two assets.
problem Valuation of European options under a two-asset Kou jump-diffusion model.
method Extends an efficient algorithm for a one-dimensional integral to a two-dimensional one, using operator splitting schemes for time discretization.
result The method achieves optimal computational cost and stable convergence for various operator splitting schemes.
Efficiently approximates integrals using a subset of samples from a target distribution in RKHS.
problem Approximating integrals with a target distribution using limited pointwise evaluations.
method Proposes a procedure using a small random subset of samples from the target distribution, either uniformly or using approximate leverage scores.
result Upper bound on approximation error for both sampling strategies, achieving optimal rate with reduced evaluations.
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 study examines stability of Hamiltonian Poisson integrators on both integrable and non-integrable systems.
problem Investigating stability properties of Hamiltonian Poisson integrators.
method Examples of Lotka-Volterra dynamics and numerical investigations of a non-integrable system are used.
result The existence of a modified Hamiltonian is crucial for the stability of Hamiltonian Poisson integrators.
Bayesian quadrature improves integration on Riemannian manifolds.
problem Efficiently computing integrals on nonlinear geometric data.
method Probabilistic numerical methods, specifically Bayesian quadrature, on Riemannian manifolds.
result Bayesian quadrature reduces the number of function evaluations compared to Monte Carlo methods.
Paper proposes efficient method to calculate Fisher-Bingham distribution normalizing constant.
problem Efficiently calculating the normalizing constant of Fisher-Bingham distributions.
method Numerical integration with continuous Euler transform to Fourier-type integral representation.
result The method is fast and accurate, applicable to high-dimensional distributions.
Median-of-means sampling outperforms mean-of-means for large sample sizes in numerical integration.
problem Improving numerical integration accuracy in high dimensions.
method Median-of-means sampling compared to mean-of-means using RQMC methods.
result Median-of-means sampling is superior for large sample sizes, while mean-of-means is better for smaller sample sizes.
This paper proposes a new method to learn integration schemes for complex ODEs.
problem Learning efficient integration schemes for non-linear ODEs and their identification.
method A novel framework to learn integration schemes that minimize an integration-related cost function.
result The proposed learning-based approach provides integration schemes close to analytical solutions.
New method for pricing American options in time-dependent models, improving accuracy and efficiency.
problem Pricing American options in time-dependent models with improved accuracy and efficiency.
method Semi-analytical pricing using a nonlinear Volterra integral equation and numerical methods.
result Improved accuracy and efficiency in pricing American options compared to forward finite difference solvers.
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…
In the framework of path integral the evolution operator kernel for the Merton-Garman Hamiltonian is constructed. Based on this kernel option formula is obtained, which generalizes the well-known Black-Scholes result. Possible approximation numerical schemes for path integral calculations are proposed.
Fenrir uses probabilistic numerics to simplify solving initial value problems.
problem Solving initial value problems in ordinary differential equations.
method Probabilistic numerics and Gauss--Markov regression.
result The method simplifies parameter estimation in ODEs, making it easier and more robust.
This paper develops a novel numerical method for pricing American options in a two-asset jump-diffusion model.
problem Pricing American options under correlated two-asset jump-diffusion models using finite difference methods often fails to preserve monotonicity and accurately discretize jump integrals.
method Introduces a novel monotone integration scheme to solve 2-D Partial Integro-Differential Equations (PIDEs) efficiently and accurately.
result The proposed method ensures convergence to the viscosity solution of the variational inequality and is both ℓ∞-stable and consistent. Method calculates function integrals on complex manifolds.
problem Integrating functions on complex manifolds.
method Digital representation and calculation method.
result Integral calculation on compact manifolds.
Research covers geometry, analysis, and integration on infinite-dimensional spaces.
problem Exploring geometric and analytical structures in infinite-dimensional settings.
method Analyzes numerical schemes, Lie groups, connections, and integration theory.
result Developed new methods for integration and analysis on infinite-dimensional manifolds.
The paper finds new metrics for geodesic flows with rational integrals.
problem Finding Riemannian metrics with rational integrals for geodesic flows.
method Explicit construction of metrics and integrals.
result New examples of metrics with rational integrals are provided.
In this note we describe how some objects from generalized geometry appear in the qualitative analysis and numerical simulation of mechanical systems. In particular we discuss double vector bundles and Dirac structures. It turns out that those objects can be naturally associated to systems with constraints -- we recall…
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…
Study discretizes Dirac and port-Hamiltonian systems using manifolds.
problem Discretization of Dirac and port-Hamiltonian systems.
method Retraction and discretization maps on manifolds for Dirac structures, applied to port-Hamiltonian systems.
result Numerical integrators for port-Hamiltonian systems derived from discretization techniques.
We consider an application involving a financial quadratic portfolio of options, when the joint underlying log-returns changes with multivariate elliptic distribution. This motivates the needs for methods for the approximation of multiple integrals over hyperboloids. A transformation is used to reduce the hyperboloid i…
Paper solves stock loan pricing with finite maturity using integral equations.
problem Valuation of margin-call stock loans with finite maturities.
method Fourier Sine transform and Volterra integral equation approach.
result Integral representation of margin-call stock loan value.
Consider a process, stochastic or deterministic, obtained by using a numerical integration scheme, or from Monte-Carlo methods involving an approximation to an integral, or a Newton-Raphson iteration to approximate the root of an equation. We will assume that we can sample from the distribution of the process from time…
Combines normalizing flows and quasi-Monte Carlo for improved numerical integration.
problem Improving the efficiency of numerical integration methods.
method Uses normalizing flows to approximate distributions and quasi-Monte Carlo for sampling.
result Demonstrates an estimator with significantly lower variance.