Inc. The following example programs demonstrate the use of some OpenMP directives. Helmholtz Equation using Jacobi 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 file permanent.c contains a straightforward implementation of a recursive algorithm to calculate the permanent of a matrix. The algorithm illustrated is not particularly efficient, but the code is extremely effective in demonstrating the use of KAI's WorkQueuing constructs.