Generate Samples from Multivariate Truncated Normal Distributions
Create a Zigzag-HMC engine object
Create a Zigzag-NUTS engine object
Draw a random Laplace momentum
One-step Harmonic HMC Sampler (Whitened Coordinates)
Efficient Cholesky decomposition
Set the mean for the target MTN
Set the precision matrix for the target MTN
Get an eligible initial value for a MTN with given mean and truncation...
Draw one Markovian zigzag sample
Draw one MTN sample with Zigzag-HMC or Zigzag-NUTS
Sample from a truncated Gaussian distribution with the harmonic HMC
Markovian Zigzag Sampler
Sample from a truncated Gaussian distribution
Efficient sampling from high-dimensional truncated Gaussian distributions, or multivariate truncated normal (MTN). Techniques include zigzag Hamiltonian Monte Carlo as in Akihiko Nishimura, Zhenyu Zhang and Marc A. Suchard (2024) <doi:10.1080/01621459.2024.2395587>, and harmonic Monte Carlo in Ari Pakman and Liam Paninski (2014) <doi:10.1080/10618600.2013.788448>.