22.12.2013 Views

The HMC Algorithm with Overrelaxation and Adaptive--Step ...

The HMC Algorithm with Overrelaxation and Adaptive--Step ...

The HMC Algorithm with Overrelaxation and Adaptive--Step ...

SHOW MORE
SHOW LESS

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>

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!