  We propose a generic approach for numerically efficient simulation from
analytically intractable distributions with constrained support. Our approach
relies upon Generalized Randomized Hamiltonian Monte Carlo (GRHMC) processes
and combines these with a randomized transition kernel that appropriately
adjusts the Hamiltonian flow at the boundary of the constrained domain,
ensuring that it remains within the domain. The numerical implementation of
this constrained GRHMC process exploits the sparsity of the randomized
transition kernel and the specific structure of the constraints so that the
proposed approach is numerically accurate, computationally fast and operational
even in high-dimensional applications. We illustrate this approach with
posterior distributions of several Bayesian models with challenging parameter
domain constraints in applications to real-word data sets. Building on the
capability of GRHMC processes to efficiently explore otherwise challenging and
high-dimensional posteriors, the proposed method expands the set of Bayesian
models that can be analyzed by using the standard Markov-Chain Monte-Carlo
(MCMC) methodology, As such, it can advance the development and use of Bayesian
models with useful constrained priors, which are difficult to handle with
existing methods. The article is accompanied by an R-package
(\url{https://github.com/torekleppe/pdmphmc}), which allows for automatically
implementing GRHMC processes for arbitrary target distributions and domain
constraints.
