2D Red-Black Gauss Seidel
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.
For a detailed description of the data locality
optimization techniques used for the data access
reordering you might want to look into the article
of Iterative Methods.
- rb1.F: a straightforward
implementation of red-black Gauss Seidel
- rb2.F and rb3.F: implement the
fusion data locality optimization
- rb4.F - rb8.F: implement the
blocking data locality optimization
- rb9.F: implements a
Comments and Bug Reports:
If you have any comments or bug reports feel free
to send an email to firstname.lastname@example.org.