We discuss how to simulate bridge processes by conditioning a stochastic process on a manifold whose generator is a hypo-elliptic operator. This operator is, up to a drift-term, the sub-Laplacian of a bracketgenerating sub-Riemannian structure, meaning in particular that it has positive smooth density everywhere. The logarithmic gradient of this density is called the score, and we show that it is needed to describe the generator of the bridge process. We therefore discuss several methods for how we can estimate the score using a neural network, with examples. The results are from a joint work with Stefan Sommer (Copenhagen) and Karen Habermann (Warwick).[-]

Simulation of conditioned diffusion processes is an essential tool in inference for stochastic processes, data imputation, generative modelling, and geometric statistics. Whilst simulating diffusion bridge processes is already difficult on Euclidean spaces, when considering diffusion processes on Riemannian manifolds the geometry brings in further complications. In even higher generality, advancing from Riemannian to sub-Riemannian geometries introduces hypoellipticity, and the possibility of finding appropriate explicit approximations for the score, the logarithmic gradient of the density, of the diffusion process is removed. We handle these challenges and construct a method for bridge simulation on sub-Riemannian manifolds by demonstrating how recent progress in machine learning can be modified to allow for training of score approximators on sub-Riemannian manifolds. Since gradients dependent on the horizontal distribution, we generalise the usual notion of denoising loss to work with non-holonomic frames using a stochastic Taylor expansion, and we demonstrate the resulting scheme both explicitly on the Heisenberg group and more generally using adapted coordinates. Joint work with Erlend Grong (Bergen) and Stefan Sommer (Copenhagen).[-]

In computational anatomy and, more generally, shape analysis, the Large Deformation Diffeomorphic Metric Mapping framework models shape variations as diffeomorphic deformations. An important shape space within this framework is the space consisting of shapes characterised by $n \geq 2$ distinct landmark points in $\mathbb{R}^d$. In diffeomorphic landmark matching, two landmark configurations are compared by solving an optimization problem which minimizes a suitable energy functional associated with flows of compactly supported diffeomorphisms transforming one landmark configuration into the other one. The landmark manifold $Q$ of $n$ distinct landmark points in $\mathbb{R}^d$ can be endowed with a Riemannian metric $g$ such that the above optimization problem is equivalent to the geodesic boundary value problem for $g$ on $Q$. Despite its importance for modeling stochastic shape evolutions, no general result concerning long-time existence of Brownian motion on the Riemannian manifold $(Q, g)$ is known. I will present joint work with Philipp Harms and Stefan Sommer on first progress in this direction which provides a full characterization of long-time existence of Brownian motion for configurations of exactly two landmarks, governed by a radial kernel.[-]