MAKMAT (L)

NAME

MAKMAT, CHMAKMAT

SYNOPSIS

CALL MAKMAT  (lda, n, type, fn, u, ier)

CALL CHMAKMAT (lda, n, type, fn, u, ier)

IMPLEMENTATION

Cray MPP systems

DESCRIPTION

MAKMAT/CHMAKMAT allocates a global array of real/complex type dimensioned (lda,n) on each calling PE. The array is initialized with respect to the indicated types noted below and written to disk as an unformatted, direct-access file with record length 8*lda*n (16*lda*n). Generally, a single PE would call this routine.

Files produced by MAKMAT/CHMAKMAT may be used as input by the Parallel Jacobian eigensolver library routines RSPJAC/CHPJAC.

RESTRICTIONS

For the current T3D memory configuration (8 MWords per PE) it is possible to allocate an array dimensioned 1024x1024 on a single PE. Larger dimensions may generate out of resource errors.

lda (IN)

INTEGER. Leading dimension of the global matrix array.

n (IN)

INTEGER. Number of rows/columns in global matrix array.

type (IN)

INTEGER. Type of test problem as follows. All test problems produce real/complex, symmetric/Hermitian systems.

            type=1 : pentadiagonal system, bandwidth=n/2
                 2 : tridiagonal system, bandwidth=3
                 3 : dense randomized system
                 4 : n/2 x n/2 dense randomized system with 
                     diagonal
                 5 : n-3 x n-3 dense randomized system with 
                     diagonal

The following additional types are valid for the complex case (CHMAKMAT) only:

            type=6 : tridiagonal system, bandwidth=3
                 7 : 3x3 verification problem
                 8 : dense randomized system

fn (IN)

CHARACTER*(*) array. Name of global output file.

u (IN)

INTEGER. Unit number associated with file fn.

ier (OUT)

INTEGER. Return code from an OPEN statement on direct unformatted file fn with status='unknown' and record length 8*lda*n (16*lda*n).

SEE ALSO

RSPJAC(L), SHOMAT(L), FGET(L), FPUT(L)

REFERENCES

O'Neal, D. C., and Reddy, N. R., "Solving Symmetric Eigenvalue Problems on Distributed Memory Machines," 1994 Fall Proceedings, Cray User Group, Tours, France, p. 76-96.

MORE INFORMATION

http://www.psc.edu/general/software/packages/pjac/pjac.html