version 1.17, 2017/06/17 11:06:23
|
version 1.19, 2020/05/21 21:45:59
|
Line 161
|
Line 161
|
IF( C.EQ.ZERO ) THEN |
IF( C.EQ.ZERO ) THEN |
CS = ONE |
CS = ONE |
SN = ZERO |
SN = ZERO |
GO TO 10 |
|
* |
* |
ELSE IF( B.EQ.ZERO ) THEN |
ELSE IF( B.EQ.ZERO ) THEN |
* |
* |
Line 174
|
Line 173
|
A = TEMP |
A = TEMP |
B = -C |
B = -C |
C = ZERO |
C = ZERO |
GO TO 10 |
* |
ELSE IF( ( A-D ).EQ.ZERO .AND. SIGN( ONE, B ).NE.SIGN( ONE, C ) ) |
ELSE IF( ( A-D ).EQ.ZERO .AND. SIGN( ONE, B ).NE.SIGN( ONE, C ) ) |
$ THEN |
$ THEN |
CS = ONE |
CS = ONE |
SN = ZERO |
SN = ZERO |
GO TO 10 |
* |
ELSE |
ELSE |
* |
* |
TEMP = A - D |
TEMP = A - D |
Line 207
|
Line 206
|
SN = C / TAU |
SN = C / TAU |
B = B - C |
B = B - C |
C = ZERO |
C = ZERO |
|
* |
ELSE |
ELSE |
* |
* |
* Complex eigenvalues, or real (almost) equal eigenvalues. |
* Complex eigenvalues, or real (almost) equal eigenvalues. |
Line 269
|
Line 269
|
* |
* |
END IF |
END IF |
* |
* |
10 CONTINUE |
|
* |
|
* Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I). |
* Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I). |
* |
* |
RT1R = A |
RT1R = A |