fft_serial, a Fortran90 code which computes a Fast Fourier Transform, and is intended as a starting point for developing a parallel version.
On an Apple PowerPC G5 with two processors, the following results were observed:
| N | Time | MegaFLOPS |
|---|---|---|
| 2 | 0.1092E-05 | 9.161 |
| 4 | 0.2270E-05 | 17.62 |
| 8 | 0.3381E-05 | 35.49 |
| 16 | 0.4817E-05 | 66.43 |
| 32 | 0.8672E-05 | 92.25 |
| 64 | 0.1849E-04 | 103.8 |
| 128 | 0.3628E-04 | 123.5 |
| 256 | 0.8159E-04 | 125.5 |
| 512 | 0.1651E-03 | 139.6 |
| 1024 | 0.3702E-03 | 138.3 |
| 2048 | 0.7618E-03 | 147.9 |
| 4096 | 0.1704E-02 | 144.3 |
| 8192 | 0.3487E-02 | 152.7 |
| 16384 | 0.8183E-02 | 140.2 |
| 32768 | 0.1712E-01 | 143.6 |
| 65536 | 0.4002E-01 | 131.0 |
| 131072 | 0.8538E-01 | 130.5 |
| 262144 | 0.1901 | 124.1 |
| 524288 | 0.391 | 127.3 |
| 1048576 | 0.8493 | 123.5 |
The computer code and data files described and made available on this web page are distributed under the MIT license.
fft_serial is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version.
fft_openmp, a Fortran90 code which computes a Fast Fourier Transform using OpenMP.
fftpack5, a Fortran90 code which implements the Fast Fourier Transform by Paul Swarztrauber and Dick Valent;
fftpack51, a Fortran90 code which implements the Fast Fourier Transform by Paul Swarztrauber and Dick Valent;
fftw_test, Fortran90 codes which illustrate the use of the FFTW3 code for Fast Fourier Transforms, by Matteo Frigo and Steven Johnson.
fire_serial, a Fortran90 code which simulates a forest fire over a rectangular array of trees, starting at a single random location. It is intended as a starting point for the development of a parallel version.
heated_plate, a Fortran90 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, a Fortran90 code which carries out a molecular dynamics simulation, and is intended as a starting point for developing a parallel version using OpenMP.
mxm_serial, a Fortran90 code which sets up a matrix multiplication problem A=B*C, intended as a starting point for implementing a parallel version.
openmp_test, a Fortran90 code which illustrates the use of the OpenMP application code interface for carrying out parallel computations in a shared memory environment.
poisson_serial, a Fortran90 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 Fortran90 code which counts the primes between 1 and N, intended as a starting point for the creation of a parallel version.
quad_serial, a Fortran90 code which approximates an integral using a quadrature rule, and is intended as a starting point for parallelization exercises.
search_serial, a Fortran90 code which searches the integers from A to B for a value J such that F(J) = C. this version of the code is intended as a starting point for a parallel approach.
sftpack, a Fortran90 code which implements the "slow" Fourier transform, intended as a teaching tool and comparison with the fast Fourier transform.
sine_transform, a Fortran90 code which demonstrates some simple properties of the discrete sine transform.
subset_sum_serial, a Fortran90 code which seeks solutions of the subset sum problem, in which it is desired to find a subset of a set of integers which has a given sum; this version of the code is intended as a starting point for a parallel approach.