mc.cpp
void mc::propose_hamilton()
{
    /* Draw random prior momenta. */
    for (int i=0; i<Nq; i++) p[i]=randn(0.0,sqrt(m[i]));
    
    /* Integrate Hamilton's equations. */
    leap_frog(verbose);
}
/* Proposal based on the solution of Hamilton's equation. -------------------------*/
Proposal based on the solution of Hamilton's equation.
Proposal based on the solution of Hamilton ' s equation .
void mc::propose_hamilton()
{
    for (int i=0; i<Nq; i++) p[i]=randn(0.0,sqrt(m[i]));
    leap_frog(verbose);
}
void mc :: propose_hamilton ( ) { for ( int i = 0 ; i < Nq ; i ++ ) p [ i ] = randn ( 0.0 , sqrt ( m [ i ] ) ) ; leap_frog ( verbose ) ; }
Proposal based on the solution of Hamilton's equation.
Proposal based on the solution of Hamilton ' s equation .
/* Draw random prior momenta. */ /* Integrate Hamilton's equations. */