Saved in:
Bibliographic Details
Main Authors: Sommer, David, Gruhlke, Robert, Kirstein, Max, Eigel, Martin, Schillings, Claudia
Format: Preprint
Published: 2024
Subjects:
Online Access:https://arxiv.org/abs/2402.15285
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909118291247104
author Sommer, David
Gruhlke, Robert
Kirstein, Max
Eigel, Martin
Schillings, Claudia
author_facet Sommer, David
Gruhlke, Robert
Kirstein, Max
Eigel, Martin
Schillings, Claudia
contents Sampling from probability densities is a common challenge in fields such as Uncertainty Quantification (UQ) and Generative Modelling (GM). In GM in particular, the use of reverse-time diffusion processes depending on the log-densities of Ornstein-Uhlenbeck forward processes are a popular sampling tool. In Berner et al. [2022] the authors point out that these log-densities can be obtained by solution of a \textit{Hamilton-Jacobi-Bellman} (HJB) equation known from stochastic optimal control. While this HJB equation is usually treated with indirect methods such as policy iteration and unsupervised training of black-box architectures like Neural Networks, we propose instead to solve the HJB equation by direct time integration, using compressed polynomials represented in the Tensor Train (TT) format for spatial discretization. Crucially, this method is sample-free, agnostic to normalization constants and can avoid the curse of dimensionality due to the TT compression. We provide a complete derivation of the HJB equation's action on Tensor Train polynomials and demonstrate the performance of the proposed time-step-, rank- and degree-adaptive integration method on a nonlinear sampling task in 20 dimensions.
format Preprint
id arxiv_https___arxiv_org_abs_2402_15285
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Generative Modelling with Tensor Train approximations of Hamilton--Jacobi--Bellman equations
Sommer, David
Gruhlke, Robert
Kirstein, Max
Eigel, Martin
Schillings, Claudia
Machine Learning
Statistics Theory
35F21, 35Q84, 62F15, 65N75, 65C30
Sampling from probability densities is a common challenge in fields such as Uncertainty Quantification (UQ) and Generative Modelling (GM). In GM in particular, the use of reverse-time diffusion processes depending on the log-densities of Ornstein-Uhlenbeck forward processes are a popular sampling tool. In Berner et al. [2022] the authors point out that these log-densities can be obtained by solution of a \textit{Hamilton-Jacobi-Bellman} (HJB) equation known from stochastic optimal control. While this HJB equation is usually treated with indirect methods such as policy iteration and unsupervised training of black-box architectures like Neural Networks, we propose instead to solve the HJB equation by direct time integration, using compressed polynomials represented in the Tensor Train (TT) format for spatial discretization. Crucially, this method is sample-free, agnostic to normalization constants and can avoid the curse of dimensionality due to the TT compression. We provide a complete derivation of the HJB equation's action on Tensor Train polynomials and demonstrate the performance of the proposed time-step-, rank- and degree-adaptive integration method on a nonlinear sampling task in 20 dimensions.
title Generative Modelling with Tensor Train approximations of Hamilton--Jacobi--Bellman equations
topic Machine Learning
Statistics Theory
35F21, 35Q84, 62F15, 65N75, 65C30
url https://arxiv.org/abs/2402.15285