Data Local Iterative Methods For The Efficient Solution of Partial Differential Equations

A cooperation
between

and
.

Funded by
.

2D Red-Black Gauss Seidel Smoothers

Description:

The tar file contains nine FORTRAN77 variants of a two-dimensional red-black Gauss Seidel relaxation algorithm for a structured grid with a 5-point as well as 9-point discretization of the Laplacian operator. The variants only differ in their data access order and produce bitwise identical results. The file also contains UNIX-shell script files to create executables of all implementation variants and perform benchmarks with them.

Program Variants:

For a detailed description of the data locality optimization techniques used for the data access reordering you might want to look into the article Memory Characteristics of Iterative Methods.

• rb1.F: a straightforward implementation of red-black Gauss Seidel
• rb2.F and rb3.F: implement the fusion data locality optimization technique
• rb4.F - rb8.F: implement the blocking data locality optimization technique
• rb9.F: implements a 2D-blocking technique
  

Comments and Bug Reports:

If you have any comments or bug reports feel free to send an email to jan.treibig@informatik.uni-erlangen.de.

Download:

Download tar file (6K)