Related papers: Sharp hypocoercive convergence estimates for under…

In molecular dynamics, penalized overdamped Langevin dynamics are used to model the motion of a set of particles that follow constraints up to a parameter $\varepsilon$. The most used schemes for simulating these dynamics are the Euler…

In this article we consider sampling from log concave distributions in Hamiltonian setting, without assuming that the objective gradient is globally Lipschitz. We propose two algorithms based on monotone polygonal (tamed) Euler schemes, to…

In this paper, we prove the existence of solutions of the Poincar\'e-Lelong equation $\sqrt{-1}\partial\bar{\partial}u=f$ on a strictly convex bounded domain $\Omega\subset\mathbb{C}^n$ $(n\geq1)$, where $f$ is a $d$-closed $(1,1)$ form and…

We consider numerical approximations of stochastic Langevin equations by implicit methods. We show a weak backward error analysis result in the sense that the generator associated with the numerical solution coincides with the solution of a…

We study the problem of sampling from a distribution $\mu$ with density $\propto e^{-V}$ for some potential function $V:\mathbb R^d\to \mathbb R$ with query access to $V$ and $\nabla V$. We start with the following standard assumptions: (1)…

This paper is concerned with the multi-dimensional compressible Euler equations with time-dependent over-damping of the form $-\frac{\mu}{(1+t)^\lambda}\rho\boldsymbol u$ in $\mathbb R^n$, where $n\ge2$, $\mu>0$, and $\lambda\in[-1,0)$.…

We propose an approach to obtaining explicit estimates on the resolvent of hypocoercive operators by using Schur complements, rather than from an exponential decay of the evolution semigroup combined with a time integral. We present…

Understanding the dimension dependency of computational complexity in high-dimensional sampling problem is a fundamental problem, both from a practical and theoretical perspective. Compared with samplers with unbiased stationary…

This work introduces an unconventional inexact augmented Lagrangian method where the augmenting term is a Euclidean norm raised to a power between one and two. The proposed algorithm is applicable to a broad class of constrained nonconvex…

This paper considers smooth convex optimization problems with many functional constraints. To solve this general class of problems we propose a new stochastic perturbed augmented Lagrangian method, called SGDPA, where a perturbation is…

We propose first-order pressure-correction scheme for the incompressible Navier-Stokes equations, incorporating the recently developed the Dynamically Regularized Lagrange Multiplier (DRLM) methods. The resulting algorithms are fully…

Approximating the invariant measure and the expectation of the functionals for parabolic stochastic partial differential equations (SPDEs) with non-globally Lipschitz coefficients is an active research area and is far from being well…

This article deals with the realisation of constraints in underdamped Langevin dynamics via soft-constrained dynamics. Specifically, we study systems with a large (or small) parameter that controls the constraint mechanisms, e.g. the…

We study high-dimensional Ornstein--Uhlenbeck processes driven by L\'evy noise and consider drift matrices that decompose into a low-rank plus sparse component, capturing a few latent factors together with a sparse network of direct…

Potential functions in highly pertinent applications, such as deep learning in over-parameterized regime, are empirically observed to admit non-isolated minima. To understand the convergence behavior of stochastic dynamics in such…

In this paper, we quantitative convergence in $W_2$ for a family of Langevin-like stochastic processes that includes stochastic gradient descent and related gradient-based algorithms. Under certain regularity assumptions, we show that the…

Laplace-type results characterize the limit of sequence of measures $(\pi_\varepsilon)_{\varepsilon >0}$ with density w.r.t the Lebesgue measure $(\mathrm{d} \pi_\varepsilon / \mathrm{d} \mathrm{Leb})(x) \propto \exp[-U(x)/\varepsilon]$…

We introduce new second-order adaptive low-dissipation central-upwind (LDCU) schemes for the one- and two-dimensional hyperbolic systems of conservation laws. The new adaptive LDCU schemes employ the LDCU numerical fluxes (recently proposed…

We study the implicit Langevin Monte Carlo (iLMC) method, which simulates the overdamped Langevin equation via an implicit iteration rule. In many applications, iLMC is favored over other explicit schemes such as the (explicit) Langevin…

Let $v \ne 0$ be a vector in $\R^n$. Consider the Laplacian on $\R^n$ with drift $\Delta_{v} = \Delta + 2v\cdot \nabla$ and the measure $d\mu(x) = e^{2 \langle v, x \rangle} dx$, with respect to which $\Delta_{v}$ is self-adjoint. This…