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
|