+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+
+ MGD2: PARALLEL 2-D MULTIGRID PACKAGE
+
+ Author: Bernard Bunner (bunner@engin.umich.edu)
+ January 1998
+
+ 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

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

d(cof(x,y)*d(phi)/dx)/dx+d(cof(x,y)*d(phi)/dy=rhs(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

  [ cof(i+1/2,j)*(phi(i+1,j)-phi(i,j))
   -cof(i-1/2,j)*(phi(i,j)-phi(i-1,j)) ] / (dx*dx)
+ [ cof(i,j+1/2)*(phi(i,j+1)-phi(i,j))
   -cof(i,j-1/2)*(phi(i,j)-phi(i,j-1)) ] / (dy*dy)
=  rhs(i,j)

Periodic, Neumann and Dirichlet are possible. The code is
written in Fortran 77 with the MPI library. A Makefile is
provided for the IBM-SP2.