*** TIMING/EIG/zeispack.f~	Fri Sep 23 18:36:18 1994
--- TIMING/EIG/zeispack.f	Wed Jul  5 01:01:51 1995
***************
*** 3062,3068 ****
        DREAL(ZDUMR) = ZDUMR
        DIMAG(ZDUMI) = (0.0D0,-1.0D0)*ZDUMI
        CABS1(ZDUM) = DABS(DREAL(ZDUM)) + DABS(DIMAG(ZDUM))
!       CSIGN(ZDUM1,ZDUM2) = CDABS(ZDUM1)*(ZDUM2/CDABS(ZDUM2))
  *
  *     GET EPS FROM DLAMCH FOR NEW STOPPING CRITERION
        IF (N.LE.0 .OR. P.LE.0) RETURN
--- 3062,3068 ----
        DREAL(ZDUMR) = ZDUMR
        DIMAG(ZDUMI) = (0.0D0,-1.0D0)*ZDUMI
        CABS1(ZDUM) = DABS(DREAL(ZDUM)) + DABS(DIMAG(ZDUM))
!       CSIGN(ZDUM1,ZDUM2) = ABS(ZDUM1)*(ZDUM2/ABS(ZDUM2))
  *
  *     GET EPS FROM DLAMCH FOR NEW STOPPING CRITERION
        IF (N.LE.0 .OR. P.LE.0) RETURN
***************
*** 3281,3287 ****
  *           INCREMENT OP COUNT
              IOPS = IOPS + 23
              IF (WANTU) IOPS = IOPS + 6*N
!             T = DCMPLX(CDABS(S(I)),0.0D0)
              R = S(I)/T
              S(I) = T
              IF (I .LT. M) E(I) = E(I)/R
--- 3281,3287 ----
  *           INCREMENT OP COUNT
              IOPS = IOPS + 23
              IF (WANTU) IOPS = IOPS + 6*N
!             T = DCMPLX(ABS(S(I)),0.0D0)
              R = S(I)/T
              S(I) = T
              IF (I .LT. M) E(I) = E(I)/R
***************
*** 3294,3300 ****
  *           INCREMENT OP COUNT
              IOPS = IOPS + 20
              IF (WANTV) IOPS = IOPS + 6*P
!             T = DCMPLX(CDABS(E(I)),0.0D0)
              R = T/E(I)
              E(I) = T
              S(I+1) = S(I+1)*R
--- 3294,3300 ----
  *           INCREMENT OP COUNT
              IOPS = IOPS + 20
              IF (WANTV) IOPS = IOPS + 6*P
!             T = DCMPLX(ABS(E(I)),0.0D0)
              R = T/E(I)
              E(I) = T
              S(I+1) = S(I+1)*R
***************
*** 3346,3358 ****
  *
  *           INCREMENT OP COUNT
              IOPST = IOPST + 17
!             TEST = CDABS(S(L)) + CDABS(S(L+1))
  *
  *           REPLACE STOPPING CRITERION WITH NEW ONE
  *
! *           ZTEST = TEST + CDABS(E(L))
  *           IF (ZTEST .NE. TEST) GO TO 420
!             IF (CDABS(E(L)) .GT. EPS * TEST) GOTO 420
  *
                 E(L) = (0.0D0,0.0D0)
  C        ......EXIT
--- 3346,3358 ----
  *
  *           INCREMENT OP COUNT
              IOPST = IOPST + 17
!             TEST = ABS(S(L)) + ABS(S(L+1))
  *
  *           REPLACE STOPPING CRITERION WITH NEW ONE
  *
! *           ZTEST = TEST + ABS(E(L))
  *           IF (ZTEST .NE. TEST) GO TO 420
!             IF (ABS(E(L)) .GT. EPS * TEST) GOTO 420
  *
                 E(L) = (0.0D0,0.0D0)
  C        ......EXIT
***************
*** 3374,3387 ****
  *
  *              INCREMENT OP COUNT
                 IOPST = IOPST + 18
!                IF (LS .NE. M) TEST = TEST + CDABS(E(LS))
!                IF (LS .NE. L + 1) TEST = TEST + CDABS(E(LS-1))
  *
  *              REPLACE STOPPING CRITERION WITH NEW ONE AS IN LAPACK
  *
! *              ZTEST = TEST + CDABS(S(LS))
  *              IF (ZTEST .NE. TEST) GO TO 460
!                IF (CDABS(S(LS))  .GT. EPS * TEST) GOTO 460
  *
                    S(LS) = (0.0D0,0.0D0)
  C           ......EXIT
--- 3374,3387 ----
  *
  *              INCREMENT OP COUNT
                 IOPST = IOPST + 18
!                IF (LS .NE. M) TEST = TEST + ABS(E(LS))
!                IF (LS .NE. L + 1) TEST = TEST + ABS(E(LS-1))
  *
  *              REPLACE STOPPING CRITERION WITH NEW ONE AS IN LAPACK
  *
! *              ZTEST = TEST + ABS(S(LS))
  *              IF (ZTEST .NE. TEST) GO TO 460
!                IF (ABS(S(LS))  .GT. EPS * TEST) GOTO 460
  *
                    S(LS) = (0.0D0,0.0D0)
  C           ......EXIT
***************
*** 3458,3465 ****
  *
  *           INCREMENT OP COUNT
              IOPST = IOPST + 48
!             SCALE = DMAX1(CDABS(S(M)),CDABS(S(M-1)),CDABS(E(M-1)),
!      *                    CDABS(S(L)),CDABS(E(L)))
              SM = DREAL(S(M))/SCALE
              SMM1 = DREAL(S(M-1))/SCALE
              EMM1 = DREAL(E(M-1))/SCALE
--- 3458,3465 ----
  *
  *           INCREMENT OP COUNT
              IOPST = IOPST + 48
!             SCALE = DMAX1(ABS(S(M)),ABS(S(M-1)),ABS(E(M-1)),
!      *                    ABS(S(L)),ABS(E(L)))
              SM = DREAL(S(M))/SCALE
              SMM1 = DREAL(S(M-1))/SCALE
              EMM1 = DREAL(E(M-1))/SCALE