Contents
zgssin - Initialize the general sparse solver.
SUBROUTINE ZGSSIN ( 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)
ZGSSIN - Initialize the sparse solver and input the matrix
structure.
MTXTYP - CHARACTER*2
On entry, MTXTYP specifies the coefficient matrix
type. Specifically, the valid options are:
'sp' or 'SP' - symmetric structure, Hermitian positive definite values
'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
the coefficient matrix. NEQNS must be at least one.
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.
-105 : NEQNS < 1