The HMC Algorithm with Overrelaxation and Adaptive--Step ...
The HMC Algorithm with Overrelaxation and Adaptive--Step ...
The HMC Algorithm with Overrelaxation and Adaptive--Step ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
−3<br />
10 10−2 10−1 100 101<br />
Background<br />
Aim of Talk<br />
<strong>The</strong> Hamiltonian Monte Carlo (<strong>HMC</strong>) <strong>Algorithm</strong><br />
Improving Performance of <strong>HMC</strong> <strong>Algorithm</strong><br />
Numerical Experiments & Results<br />
Example results from the improved <strong>HMC</strong> algorithm<br />
Evaluation criteria – Geometric Illustration<br />
Degree of correlation<br />
Spectral analysis<br />
Autocorrelation function<br />
0.9<br />
0.7<br />
0.5<br />
10 0<br />
ρ(τ)<br />
0.3<br />
P(κ)<br />
0.1<br />
10 −1<br />
−0.1<br />
−0.3<br />
−0.5<br />
0 1 2 3 4 5 6 7 8<br />
τ<br />
10 1 κ<br />
10−2<br />
−3<br />
10 10−2 10−1 100 101<br />
1.3<br />
Integrated autocorrelation function<br />
<strong>The</strong> curve turn over<br />
1.1<br />
ρ int<br />
(τ)<br />
0.9<br />
10 0<br />
10 1 κ<br />
P(0)<br />
P(κ)<br />
P(κ)≈ (κ) −α<br />
0.7<br />
0.5<br />
0.3<br />
0 1 2 3 4 5 6 7 8<br />
τ<br />
κ*<br />
M. Alfaki, S. Subbey, <strong>and</strong> D. Haugl<strong>and</strong> <strong>The</strong> Hamiltonian Monte Carlo (<strong>HMC</strong>) <strong>Algorithm</strong>