diff options
Diffstat (limited to 'math/lapack/files/patch-ag')
| -rw-r--r-- | math/lapack/files/patch-ag | 132 |
1 files changed, 0 insertions, 132 deletions
diff --git a/math/lapack/files/patch-ag b/math/lapack/files/patch-ag deleted file mode 100644 index 835a5a64eaa6..000000000000 --- a/math/lapack/files/patch-ag +++ /dev/null @@ -1,132 +0,0 @@ -*** 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 |