Bunner

    A public domain, parallel processor replacement Fortran-77 code for Mudpack written by Bernard Bunner. It requires both MPI and the C preprocessor. This code features staggered grids, rectangular domains, constant mesh spacing in 2 or 3 dimensions, V or W cycling, full weighting for the restriction and bilinear interpolation for the correction, either vertex-centered or cell-centered operations, and periodic, Neumann, or Dirichlet boundary conditions.

    There are four files comprising two packages

File	Bytes	Date	Comment
README_mgd2	910	Feb 4 1998	Two dimensional code
mgd2.tgz	23725	Feb 4 1998
README_mgd3	1071	Feb 4 1998	Three dimensional code
mgd3.tgz	29873	Feb 4 1998

Should you have trouble unpacking one of the .tgz files, the contents should be consulted.  General instructions on how to unpack the files in the .tgz files are in this hyperlink.

    mgd2 is a parallel 2D multigrid program which solves the non-separable Poisson equation:

(A(x,y)u_x)_x + (A(x,y)u_y)_y = f(x,y)

on a staggered grid. The rectangular domain has a constant grid step in both directions and is decomposed into rectangular subdomains. In discretized form, this equation can be written as

[A(i+.5,j)(u(i+1,j)-u(i,j)) - A(i-.5,j)(u(i,j)-u(i-1,j))](dx)^-2 +

[A(i,j+.5)(u(i,j+1)-u(i,j)) - A(i,j-.5)(u(i,j)-u(i,j-1))](dy)^-2 =

f(i,j)

    mgd3 is a parallel 3D multigrid program which solves the non-separable Poisson equation:

(A(x,y,z)u_x)_x + (A(x,y,z)u_y)_y + (A(x,y,z)u_z)_z = f(x,y,z)

on a staggered grid. The rectangular domain has a constant grid step in both directions and is decomposed into rectangular subdomains. In discretized form, this equation can be written as

[A(i+.5,j,k)(u(i+1,j,k)-u(i,j,k)) - A(i-.5,j,k)(u(i,j,k)-u(i-1,j,k))](dx)^-2 +

[A(i,j+.5,k)(u(i,j+1,k)-u(i,j,k)) - A(i,j-.5,k)(u(i,j,k)-u(i,j-1,k))](dy)^-2 +

[A(i,j,k+.5)(u(i,j,k+1)-u(i,j,k)) - A(i,j,k-.5)(u(i,j,k)-u(i,j,k-1))](dz)^-2 =

f(i,j,k)

    For both packages, periodic, Neumann, and Dirichlet boundary conditions are possible.  The codes are written in Fortran-77 and use the MPI library.  A Makefile is provided for the IBM-SP2.

    These packages can be unpacked using either the GNU tar utility and the command

tar zxvf mgd2.tgz

or using the GNU gunzip utility and the commands

gunzip mgd2.tgz
tar xvf mgd2.tar

and similarly for mgd3.tgz.

Cheers,
Craig C. Douglas

Last modified: 

MGNet Sites	Newsletters	Bibliography	Conferences
Free Software	Free Subscription	Preprints	Proceedings
Adding to MGNet	Retrieving Files	Books	Sponsors
Index	Links	Tutorials	Mystery Link