Previous: figi Up: ../eispad.html Next: hqr


FIGI2(NM,N,T,D,E,Z,IERR)

       SUBROUTINE FIGI2(NM,N,T,D,E,Z,IERR)
 C
       INTEGER I,J,N,NM,IERR
       DOUBLE PRECISION T(NM,3),D(N),E(N),Z(NM,N)
       DOUBLE PRECISION H
 C
 C     GIVEN A NONSYMMETRIC TRIDIAGONAL MATRIX SUCH THAT THE PRODUCTS
 C     OF CORRESPONDING PAIRS OF OFF-DIAGONAL ELEMENTS ARE ALL
 C     NON-NEGATIVE, AND ZERO ONLY WHEN BOTH FACTORS ARE ZERO, THIS
 C     SUBROUTINE REDUCES IT TO A SYMMETRIC TRIDIAGONAL MATRIX
 C     USING AND ACCUMULATING DIAGONAL SIMILARITY TRANSFORMATIONS.
 C
 C     ON INPUT
 C
 C        NM MUST BE SET TO THE ROW DIMENSION OF TWO-DIMENSIONAL
 C          ARRAY PARAMETERS AS DECLARED IN THE CALLING PROGRAM
 C          DIMENSION STATEMENT.
 C
 C        N IS THE ORDER OF THE MATRIX.
 C
 C        T CONTAINS THE INPUT MATRIX.  ITS SUBDIAGONAL IS
 C          STORED IN THE LAST N-1 POSITIONS OF THE FIRST COLUMN,
 C          ITS DIAGONAL IN THE N POSITIONS OF THE SECOND COLUMN,
 C          AND ITS SUPERDIAGONAL IN THE FIRST N-1 POSITIONS OF
 C          THE THIRD COLUMN.  T(1,1) AND T(N,3) ARE ARBITRARY.
 C
 C     ON OUTPUT
 C
 C        T IS UNALTERED.
 C
 C        D CONTAINS THE DIAGONAL ELEMENTS OF THE SYMMETRIC MATRIX.
 C
 C        E CONTAINS THE SUBDIAGONAL ELEMENTS OF THE SYMMETRIC
 C          MATRIX IN ITS LAST N-1 POSITIONS.  E(1) IS NOT SET.
 C
 C        Z CONTAINS THE TRANSFORMATION MATRIX PRODUCED IN
 C          THE REDUCTION.
 C
 C        IERR IS SET TO
 C          ZERO       FOR NORMAL RETURN,
 C          N+I        IF T(I,1)*T(I-1,3) IS NEGATIVE,
 C          2*N+I      IF T(I,1)*T(I-1,3) IS ZERO WITH
 C                     ONE FACTOR NON-ZERO.
 C
 C     QUESTIONS AND COMMENTS SHOULD BE DIRECTED TO BURTON S. GARBOW,
 C     MATHEMATICS AND COMPUTER SCIENCE DIV, ARGONNE NATIONAL LABORATORY
 C
 C     THIS VERSION DATED AUGUST 1983.
 C
 C     ------------------------------------------------------------------
 C