dgssin


NAME

dgssin - Initialize the sparse solver.


SYNOPSIS

 SUBROUTINE DGSSIN ( MTXTYP, PIVOT, NEQNS, COLSTR, ROWIND, OUTUNT,
                     MSGLVL, HANDLE, IER )
 CHARACTER*2       MTXTYP
 CHARACTER*1       PIVOT
 INTEGER           NEQNS, COLSTR(*), ROWIND(*), OUTUNT, MSGLVL, IER
 DOUBLE PRECISION  HANDLE(150)


PURPOSE

DGSSIN - Initialize the sparse solver and input the matrix structure.


PARAMETERS

MTXTYP - CHARACTER*2
On entry, MTXTYP specifies the coefficient matrix type. Specifically, the valid options are:

 'ss' or 'SS' - symmetric structure, symmetric values
 'su' or 'SU' - symmetric structure, unsymmetric values
 'uu' or 'UU' - unsymmetric structure, unsymmetric values

Unchanged on exit.

PIVOT - CHARACTER*1
On entry, PIVOT specifies whether or not pivoting is used in the course of the numeric factorization. The valid options are:

 'n' or 'N' - no pivoting is used
 (Pivoting is not supported for this release).

Unchanged on exit.

NEQNS - INTEGER
On entry, NEQNS specifies the number of equations in coefficient matrix. Unchanged on exit.

COLSTR(*) - INTEGER array
On entry, COLSTR(*) is an array of size (NEQNS+1), containing the pointers of the matrix structure. Unchanged on exit.

ROWIND(*) - INTEGER array
On entry, ROWIND(*) is an array of size COLSTR(NEQNS+1)-1, containing the indices of the matrix structure. Unchanged on exit.

HANDLE(150) - DOUBLE PRECISION array
On entry, HANDLE(*) is an array containing information needed by the solver, and must be passed unchanged to each sparse solver subroutine. Modified on exit.

OUTUNT - INTEGER
Output unit. Unchanged on exit.

MSGLVL - INTEGER
Message level.

 0 - no output from solver.
 (No messages supported for this release.)

Unchanged on exit.

IER - INTEGER
Error number. If no error encountered, unchanged on exit. If error encountered, it is set to a non-zero integer. Error numbers set by this subroutine:

 -101 : Failure to dynamically allocate memory.
 -102 : Invalid matrix type.
 -103 : Invalid pivot option.
 -104 : Number of nonzeros less than NEQNS.

.