The jacobi.f program solves the Helmholtz equation on a regular mesh, using an iterative Jacobi method with over-relaxation. Parallelism is exploited in both the solver and the numerical error checking
The program md.f implements a simple molecular dynamics simulation in continuous real space. The velocity Verlet algorithm is used to implement the time stepping. The force and energy computations are performed in parallel, as is the time integration. (The program uses some Fortran90 features, so an F90 compiler may be needed.)
Hybrid Monte Carlo sampling.SAMPLES = HMC(F, X, OPTIONS, GRADF) uses a hybrid Monte Carlo
algorithm to sample from the distribution P ~ EXP(-F), where F is the
first argument to HMC. The Markov chain starts at the point X, and
the function GRADF is the gradient of the `energy function F.
