summaryrefslogtreecommitdiff
path: root/math/lapack/files/patch-ag
blob: 835a5a64eaa6a03eb3a491e82adde312cadcd787 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
*** 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