BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andreas Eberle

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

ResearcherUnclaimed source record
math.PRComputationMachine LearningData Structures and AlgorithmsDistributed, Parallel, and Cluster Computingphysics.chem-phphysics.comp-ph

Catalog footprint

What is connected

8works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

8 published item(s)

preprint2020arXiv

Couplings for Andersen Dynamics

Andersen dynamics is a standard method for molecular simulations, and a precursor of the Hamiltonian Monte Carlo algorithm used in MCMC inference. The stochastic process corresponding to Andersen dynamics is a PDMP (piecewise deterministic Markov process) that iterates between Hamiltonian flows and velocity randomizations of randomly selected particles. Both from the viewpoint of molecular dynamics and MCMC inference, a basic question is to understand the convergence to equilibrium of this PDMP particularly in high dimension. Here we present couplings to obtain sharp convergence bounds in the Wasserstein sense that do not require global convexity of the underlying potential energy.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Coupling and Convergence for Hamiltonian Monte Carlo

Based on a new coupling approach, we prove that the transition step of the Hamiltonian Monte Carlo algorithm is contractive w.r.t. a carefully designed Kantorovich (L1 Wasserstein) distance. The lower bound for the contraction rate is explicit. Global convexity of the potential is not required, and thus multimodal target distributions are included. Explicit quantitative bounds for the number of steps required to approximate the stationary distribution up to a given error are a direct consequence of contractivity. These bounds show that HMC can overcome diffusive behaviour if the duration of the Hamiltonian dynamics is adjusted appropriately.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Two-scale coupling for preconditioned Hamiltonian Monte Carlo in infinite dimensions

We derive non-asymptotic quantitative bounds for convergence to equilibrium of the exact preconditioned Hamiltonian Monte Carlo algorithm (pHMC) on a Hilbert space. As a consequence, explicit and dimension-free bounds for pHMC applied to high-dimensional distributions arising in transition path sampling and path integral molecular dynamics are given. Global convexity of the underlying potential energies is not required. Our results are based on a two-scale coupling which is contractive in a carefully designed distance.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sticky couplings of multidimensional diffusions with different drifts

We present a novel approach of coupling two multidimensional and non-degenerate Itô processes $(X_t)$ and $(Y_t)$ which follow dynamics with different drifts. Our coupling is sticky in the sense that there is a stochastic process $(r_t)$, which solves a one-dimensional stochastic differential equation with a sticky boundary behavior at zero, such that almost surely $|X_t-Y_t|\leq r_t$ for all $t\geq 0$. The coupling is constructed as a weak limit of Markovian couplings. We provide explicit, non-asymptotic and long-time stable bounds for the probability of the event $\{X_t=Y_t\}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reflection couplings and contraction rates for diffusions

We consider contractivity for diffusion semigroups w.r.t. Kantorovich ($L^1$ Wasserstein) distances based on appropriately chosen concave functions. These distances are inbetween total variation and usual Wasserstein distances. It is shown that by appropriate explicit choices of the underlying distance, contractivity with rates of close to optimal order can be obtained in several fundamental classes of examples where contractivity w.r.t. standard Wasserstein distances fails. Applications include overdamped Langevin diffusions with locally non-convex potentials, products of these processes, and systems of weakly interacting diffusions, both of mean-field and nearest neighbour type.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Engineering Parallel String Sorting

We discuss how string sorting algorithms can be parallelized on modern multi-core shared memory machines. As a synthesis of the best sequential string sorting algorithms and successful parallel sorting algorithms for atomic objects, we first propose string sample sort. The algorithm makes effective use of the memory hierarchy, uses additional word level parallelism, and largely avoids branch mispredictions. Then we focus on NUMA architectures, and develop parallel multiway LCP-merge and -mergesort to reduce the number of random memory accesses to remote nodes. Additionally, we parallelize variants of multikey quicksort and radix sort that are also useful in certain situations. Comprehensive experiments on five current multi-core platforms are then reported and discussed. The experiments show that our implementations scale very well on real-world inputs and modern machines.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Error bounds for Metropolis-Hastings algorithms applied to perturbations of Gaussian measures in high dimensions

The Metropolis-adjusted Langevin algorithm (MALA) is a Metropolis-Hastings method for approximate sampling from continuous distributions. We derive upper bounds for the contraction rate in Kantorovich-Rubinstein-Wasserstein distance of the MALA chain with semi-implicit Euler proposals applied to log-concave probability measures that have a density w.r.t. a Gaussian reference measure. For sufficiently "regular" densities, the estimates are dimension-independent, and they hold for sufficiently small step sizes $h$ that do not depend on the dimension either. In the limit $h\downarrow0$, the bounds approach the known optimal contraction rates for overdamped Langevin diffusions in a convex potential. A similar approach also applies to Metropolis-Hastings chains with Ornstein-Uhlenbeck proposals. In this case, the resulting estimates are still independent of the dimension but less optimal, reflecting the fact that MALA is a higher order approximation of the diffusion limit than Metropolis-Hastings with Ornstein-Uhlenbeck proposals.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantitative approximations of evolving probability measures and sequential Markov Chain Monte Carlo methods

We study approximations of evolving probability measures by an interacting particle system. The particle system dynamics is a combination of independent Markov chain moves and importance sampling/resampling steps. Under global regularity conditions, we derive non-asymptotic error bounds for the particle system approximation. In a few simple examples, including high dimensional product measures, bounds with explicit constants of feasible size are obtained. Our main motivation are applications to sequential MCMC methods for Monte Carlo integral estimation.

0 citations0 reviewsNo code linkedOpen access