md, a C++ code which carries out a molecular dynamics simulation.
The computation involves following the paths of particles which exert a distance-dependent force on each other. The particles are not constrained by any walls; if particles meet, they simply pass through each other.
The problem is treated as a coupled set of differential equations. The system of differential equation is discretized by choosing a discrete time step. Given the position and velocity of each particle at one time step, the algorithm estimates these values at the next time step.To compute the next position of each particle requires the evaluation of the right hand side of its corresponding differential equation. Since each of these calculations is independent, there is a potential speedup if the program can take advantage of parallel computing.
This version of the program is intended as a starting point; students are asked to examine the program and add the appropriate OpenMP directives to run it in parallel.
Note that this program exhibits a remarkable speedup if we simply compile it with the "-O2" optimization option.
md nd np step_num dtwhere
The computer code and data files described and made available on this web page are distributed under the MIT license
md is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version and a Python version.
FFT_SERIAL, a C++ code which demonstrates the computation of a Fast Fourier Transform, and is intended as a starting point for developing a parallel version using OpenMP.
HEATED_PLATE, a C++ code which solves the steady state heat equation in a 2D rectangular region, and is intended as a starting point for implementing an OpenMP parallel version.
MD_OPENMP, a C++ code which carries out a molecular dynamics simulation, using OpenMP for parallel execution.
MXM_SERIAL, a C++ code which sets up a matrix multiplication problem A=B*C, intended as a starting point for implementing a parallel version.
POISSON_SERIAL, a C++ code which computes an approximate solution to the Poisson equation in a rectangle, and is intended as the starting point for the creation of a parallel version.
PRIME_SERIAL, a C++ code which counts the number of primes between 1 and N, intended as a starting point for the creation of a parallel version.
QUAD_SERIAL, a C++ code which approximates an integral using a quadrature rule, and is intended as a starting point for parallelization exercises.
SEARCH_SERIAL, a C++ code which searches the integers from A to B for a value J such that F(J) = C. this version of the program is intended as a starting point for a parallel approach.