Source author record

Marc Arnaudon

Marc Arnaudon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

Construction of set-valued dual processes on manifolds

The purpose of this paper is to construct a Brownian motion $X := (X_t)_{t\geq 0}$ taking values in a Riemannian manifold $M$, together with a compact valued process $D:= (D_t)_{t\geq 0}$ such that, at least for small enough ${\mathscr F}^D$-stopping time $τ> 0$ and conditioned by ${\mathscr F}_τ^D$, the law of $X_τ$ is the normalized Lebesgue measure on $D_τ$. This intertwining result is a generalization of Pitman theorem. We first construct regular intertwined processes related to Stokes' theorem. Then using several limiting procedures we construct synchronous intertwined, free intertwined, mirror intertwined processes. The local times of the Brownian motion on the (morphological) skeleton or the boundary of $D$ plays an important role. Several examples with moving intervals, discs, annulus, symmetric convex sets are investigated. KEYWORDS: Brownian motions on Riemannian manifolds, intertwining relations, set-valued dual processes, couplings of primal and dual processes, stochastic mean curvature evolutions, boundary and skeleton local times, generalized Pitman theorem.

preprint2022arXiv

On the separation cut-off phenomenon for Brownian motions on high dimensional spheres

This note proves that the separation convergence towards the uniform distribution abruptly occurs at times around ln(n)/n for the (time-accelerated by 2) Brownian motion on the sphere with a high dimension n. The arguments are based on a new and elementary perturbative approach for estimating hitting times in a small noise context. The quantitative estimates thus obtained are applied to the strong stationary times constructed in a privious article by the authors to deduce the wanted cut-off phenomenon.

preprint2022arXiv

Riemannian data-dependent randomized smoothing for neural networks certification

Certification of neural networks is an important and challenging problem that has been attracting the attention of the machine learning community since few years. In this paper, we focus on randomized smoothing (RS) which is considered as the state-of-the-art method to obtain certifiably robust neural networks. In particular, a new data-dependent RS technique called ANCER introduced recently can be used to certify ellipses with orthogonal axis near each input data of the neural network. In this work, we remark that ANCER is not invariant under rotation of input data and propose a new rotationally-invariant formulation of it which can certify ellipses without constraints on their axis. Our approach called Riemannian Data Dependant Randomized Smoothing (RDDRS) relies on information geometry techniques on the manifold of covariance matrices and can certify bigger regions than ANCER based on our experiments on the MNIST dataset.

preprint2020arXiv

A duality formula and a particle Gibbs sampler for continuous time Feynman-Kac measures on path spaces

Continuous time Feynman-Kac measures on path spaces are central in applied probability, partial differential equation theory, as well as in quantum physics. This article presents a new duality formula between normalized Feynman-Kac distribution and their mean field particle interpretations. Among others, this formula allows us to design a reversible particle Gibbs-Glauber sampler for continuous time Feynman-Kac integration on path spaces. This result extends the particle Gibbs samplers introduced by Andrieu-Doucet-Holenstein [2] in the context of discrete generation models to continuous time Feynman-Kac models and their interacting jump particle interpretations. We also provide new propagation of chaos estimates for continuous time genealogical tree based particle models with respect to the time horizon and the size of the systems. These results allow to obtain sharp quantitative estimates of the convergence rate to equilibrium of particle Gibbs-Glauber samplers. To the best of our knowledge these results are the first of this kind for continuous time Feynman-Kac measures.

preprint2017arXiv

An entropic interpolation problem for incompressible viscid fluids

In view of studying incompressible inviscid fluids, Brenier introduced in the late 80's a relaxation of a geodesic problem addressed by Arnold in 1966. Instead of inviscid fluids, the present paper is devoted to incompressible viscid fluids. A natural analogue of Brenier's problem is introduced, where generalized flows are no more supported by absolutely continuous paths, but by Brownian sample paths. It turns out that this new variational problem is an entropy minimization problem with marginal constraints entering the class of convex minimization problems. This paper explores the connection between this variational problem and Brenier's original problem. Its dual problem is derived and the general shape of its solution is described. Under the restrictive assumption that the pressure is a nice function, the kinematics of its solution is made explicit and its connection with the Navier-Stokes equation is established.

preprint2016arXiv

A stochastic algorithm finding $p$-means on the circle

A stochastic algorithm is proposed, finding some elements from the set of intrinsic $p$-mean(s) associated to a probability measure $ν$ on a compact Riemannian manifold and to $p\in[1,\infty)$. It is fed sequentially with independent random variables $(Y_n)_{n\in \mathbb{N}}$ distributed according to $ν$, which is often the only available knowledge of $ν$. Furthermore, the algorithm is easy to implement, because it evolves like a Brownian motion between the random times when it jumps in direction of one of the $Y_n$, $n\in\mathbb{N}$. Its principle is based on simulated annealing and homogenization, so that temperature and approximations schemes must be tuned up (plus a regularizing scheme if $ν$ does not admit a Hölderian density). The analysis of the convergence is restricted to the case where the state space is a circle. In its principle, the proof relies on the investigation of the evolution of a time-inhomogeneous $\mathbb{L}^2$ functional and on the corresponding spectral gap estimates due to Holley, Kusuoka and Stroock. But it requires new estimates on the discrepancies between the unknown instantaneous invariant measures and some convenient Gibbs measures.

preprint2016arXiv

Generalized stochastic Lagrangian paths for the Navier-Stokes equation

In the note added in proof of the seminal paper [Groups of diffeomorphisms andthe motion of an incompressible fluid, Ann. of Math. 92 (1970), 102-163], Ebinand Marsden introduced the so-called correct Laplacian for the Navier-Stokes equationon a compact Riemannian manifold. In the spirit of Brenier's generalized flows forthe Euler equation, we introduce a class of semimartingales on a compact Riemannianmanifold. We prove that these semimartingales are critical points to the correspondingkinetic energy if and only if its drift term solves weakly the Navier-Stokes equationdefined with Ebin-Marsden's Laplacian. We also show that for the torus case,classical solutions of the Navier-Stokes equation realize the minimum of the kineticenergy in a suitable class.

preprint2015arXiv

Reparameterization invariant metric on the space of curves

This paper focuses on the study of open curves in a manifold M, and proposes a reparameterization invariant metric on the space of such paths. We use the square root velocity function (SRVF) introduced by Srivastava et al. in [11] to define a reparameterization invariant metric on the space of immersions M' = Imm([0,1], M) by pullback of a metric on the tangent bundle TM' derived from the Sasaki metric. We observe that such a natural choice of Riemannian metric on TM' induces a first-order Sobolev metric on M' with an extra term involving the origins, and leads to a distance which takes into account the distance between the origins and the distance between the SRV representations of the curves. The geodesic equations for this metric are given, as well as an idea of how to compute the exponential map for observed trajectories in applications. This provides a generalized theoretical SRV framework for curves lying in a general manifold M .

preprint2013arXiv

A stochastic algorithm finding generalized means on compact manifolds

A stochastic algorithm is proposed, finding the set of generalized means associated to a probability measure on a compact Riemannian manifold M and a continuous cost function on the product of M by itself. Generalized means include p-means for p>0, computed with any continuous distance function, not necessarily the Riemannian distance. They also include means for lengths computed from Finsler metrics, or for divergences. The algorithm is fed sequentially with independent random variables Y_n distributed according to the probability measure on the manifold and this is the only knowledge of this measure required. It evolves like a Brownian motion between the times it jumps in direction of the Y_n. Its principle is based on simulated annealing and homogenization, so that temperature and approximations schemes must be tuned up. The proof relies on the investigation of the evolution of a time-inhomogeneous L^2 functional and on the corresponding spectral gap estimates due to Holley, Kusuoka and Stroock.

preprint2013arXiv

Generalized Navier-Stokes flows and applications to incompressible viscous fluids

We introduce a notion of generalized stochastic flows on mani- folds, that extends to the viscous case the one defined by Brenier for perfect fluids. Their kinetic energy extends the classical kinetic energy to Brownian flows, defined as the L2 norm of their drift. We prove that there exists a generalized flow which realizes the infimum of the kinetic energy among all generalized flows with prescribed initial and final configuration. We also con- struct generalized flows with prescribed drift and kinetic energy smaller than the L2 norm of the drift. The results are actually presented for general Lq norms, thus including not only the Navier-Stokes equations but also other equations such as the porous media.

preprint2012arXiv

Equivalent Harnack and Gradient Inequalities for Pointwise Curvature Lower Bound

By using a coupling method, an explicit log-Harnack inequality with local geometry quantities is established for (sub-Markovian) diffusion semigroups on a Riemannian manifold (possibly with boundary). This inequality as well as the consequent $L^2$-gradient inequality, are proved to be equivalent to the pointwise curvature lower bound condition together with the convexity or absence of the boundary. Some applications of the log-Harnack inequality are also introduced.

preprint2012arXiv

Means in complete manifolds: uniqueness and approximation

Let $M$ be a complete Riemannian manifold, $N\in \NN$ and $p\ge 1$. We prove that almost everywhere on $x=(x_1,...,x_N)\in M^N$ for Lebesgue measure in $M^N$, the measure $\di μ(x)=\f1N\sum_{k=1}^N\d_{x_k}$ has a unique $p$-mean $e_p(x)$. As a consequence, if $X=(X_1,...,X_N)$ is a $M^N$-valued random variable with absolutely continuous law, then almost surely $μ(X(\om))$ has a unique $p$-mean. In particular if $(X_n)_{n\ge 1}$ is an independent sample of an absolutely continuous law in $M$, then the process $e_{p,n}(\om)=e_p(X_1(\om),..., X_n(\om))$ is well-defined. Assume $M$ is compact and consider a probability measure $ν$ in $M$. Using partial simulated annealing, we define a continuous semimartingale which converges to the set of minimizers of the integral of distance at power $p$ with respect to $ν$. When the set is a singleton, it converges to the $p$-mean.

preprint2011arXiv

Medians and means in Riemannian geometry: existence, uniqueness and computation

This paper is a short summary of our recent work on the medians and means of probability measures in Riemannian manifolds. Firstly, the existence and uniqueness results of local medians are given. In order to compute medians in practical cases, we propose a subgradient algorithm and prove its convergence. After that, Fréchet medians are considered. We prove their statistical consistency and give some quantitative estimations of their robustness with the aid of upper curvature bounds. We also show that, in compact Riemannian manifolds, the Fréchet medians of generic data points are always unique. Stochastic and deterministic algorithms are proposed for computing Riemannian p-means. The rate of convergence and error estimates of these algorithms are also obtained. Finally, we apply the medians and the Riemannian geometry of Toeplitz covariance matrices to radar target detection.

preprint2011arXiv

Stochastic algorithms for computing means of probability measures

Consider a probability measure supported by a regular geodesic ball in a manifold. For any p larger than or equal to 1 we define a stochastic algorithm which converges almost surely to the p-mean of the measure. Assuming furthermore that the functional to minimize is regular around the p-mean, we prove that a natural renormalization of the inhomogeneous Markov chain converges in law into an inhomogeneous diffusion process. We give an explicit expression of this process, as well as its local characteristic.

preprint2010arXiv

The differentiation of hypoelliptic diffusion semigroups

Basic derivative formulas are presented for hypoelliptic heat semigroups and harmonic functions extending earlier work in the elliptic case. Emphasis is placed on developing integration by parts formulas at the level of local martingales. Combined with the optional sampling theorem, this turns out to be an efficient way of dealing with boundary conditions, as well as with finite lifetime of the underlying diffusion. Our formulas require hypoellipticity of the diffusion in the sense of Malliavin calculus (integrability of the inverse Malliavin covariance) and are formulated in terms of the derivative flow, the Malliavin covariance and its inverse. Finally some extensions to the nonlinear setting of harmonic mappings are discussed.