Contents
     dgsssl - Solve routine for the general sparse solver.
      SUBROUTINE DGSSSL ( NRHS, RHS, LDRHS, HANDLE, IER )
      INTEGER           NRHS, LDRHS, IER
      DOUBLE PRECISION  RHS(LDRHS,NRHS)
      DOUBLE PRECISION  HANDLE(150)
     DGSSSL - Triangular solve of a factored sparse matrix.
     NRHS      - INTEGER
             On entry, NRHS specifies the number of right hand
             sides to solve for.  Unchanged on exit.
     RHS(LDRHS,*) - DOUBLE PRECISION array
             On entry, RHS(LDRHS,NRHS) contains the NRHS right
             hand sides.  On exit, it contains the solutions.
     LDRHS     - INTEGER
             On entry, LDRHS specifies the leading dimension of
             the RHS array. 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.
     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:
              -400 : Invalid calling sequence - need to call DGSSFA first.
              -401 : Failure to dynamically allocate memory.
              -402 : NRHS < 1
              -403 : NEQNS > LDRHS