FFT_SERIAL is a FORTRAN90 program which demonstrates the computation of a Fast Fourier Transform, and is intended as a starting point for developing a parallel version using OpenMP.
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 GNU LGPL license.
FFT_SERIAL is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version.
FFT_OPENMP, a FORTRAN90 program which computes a Fast Fourier Transform using OpenMP.
FFTPACK5, a FORTRAN90 library which implements the Fast Fourier Transform by Paul Swarztrauber and Dick Valent;
FFTPACK5.1, a FORTRAN90 library which implements the Fast Fourier Transform by Paul Swarztrauber and Dick Valent;
FFTW3, FORTRAN90 programs which illustrate the use of the FFTW3 library for Fast Fourier Transforms, by Matteo Frigo and Steven Johnson.
FIRE_SERIAL, a FORTRAN90 program 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 program 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 program 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 program which sets up a matrix multiplication problem A=B*C, intended as a starting point for implementing a parallel version.
OPENMP, FORTRAN90 programs which illustrate the use of the OpenMP application program interface for carrying out parallel computations in a shared memory environment.
POISSON_SERIAL, a FORTRAN90 program 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 program 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 FORTRAN90 program which approximates an integral using a quadrature rule, and is intended as a starting point for parallelization exercises.
SEARCH_SERIAL, a FORTRAN90 program 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.
SFTPACK, a FORTRAN90 library which implements the "slow" Fourier transform, intended as a teaching tool and comparison with the fast Fourier transform.
SINE_TRANSFORM, a FORTRAN90 library which demonstrates some simple properties of the discrete sine transform.
SUBSET_SUM_SERIAL, a FORTRAN90 program 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 program is intended as a starting point for a parallel approach.
You can go up one level to the FORTRAN90 source codes.