dgssin(3P) Sun Performance Library dgssin(3P)NAME
dgssin - Initialize the general 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)PURPOSEDGSSIN - Initialize the sparse solver and input the matrix structure.
PARAMETERS
MTXTYP - CHARACTER*2
On entry, MTXTYP specifies the coefficient matrix type. Valid
options are:
'sp' or 'SP' - symmetric structure, positive-definite values
'ss' or 'SS' - symmetric structure, symmetric values
'su' or 'SU' - symmetric structure, unsymmetric values
'uu' or 'UU' - unsymmetric structure, unsymmetric values
'S0' or 's0' - unsymmetric,non-definite structure, use SuperLU
If MTXTYP is 'sp', 'SP', 'ss', or 'SS', a packed column repre‐
sentation of the lower triangular part of the symmetric matrix
is expected in COLSTR and ROWIND. Otherwise, a packed column
representation of the full matrix is expected in COLSTR and
ROWIND. Unchanged on exit.
PIVOT - CHARACTER*1
On entry, PIVOT specifies whether or not pivoting is used in
the course of the numeric factorization. This argument is
ignored if MTXTYP = 'S0' or 's0'. Currently, the only valid
option is:
'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 coef‐
ficient matrix. NEQNS must be at least one. Unchanged on
exit.
COLSTR(*) - INTEGER array
On entry, COLSTR(*) is an array of size (NEQNS+1). If MTXTYP
is 'sp', 'SP', 'ss', or 'SS', it contains the column pointers
of the lower triangular part of the symmetric matrix. Other‐
wise, it contains the column pointers of the full matrix.
Unchanged on exit.
ROWIND(*) - INTEGER array
On entry, ROWIND(*) is an array of size COLSTR(NEQNS+1)-1. If
MTXTYP is 'sp', 'SP', 'ss', or 'SS', it contains the row
indices of the nonzero values in the lower triangular part of
the symmetric matrix. Otherwise, it contains the row indices
of the nonzero values in the full matrix. 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 num‐
bers 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
-106 : structural zero on diagonal.
-107 : Zero-based indexing required for SuperLU.
-108 : Row index out of range.
-109 : Row index unsorted or duplicate entry present.
-110 : Column pointers not in ascending order.
-111 : Fortran (one-based) indexing required.
-666 : Abnormal termination - internal error.
3rd Berkeley Distribution 24 Mar 2009 dgssin(3P)