version 1.3, 2010/08/06 15:32:19
|
version 1.4, 2010/08/07 13:18:05
|
Line 33
|
Line 33
|
* |
* |
* DD1 (input/output) DOUBLE PRECISION |
* DD1 (input/output) DOUBLE PRECISION |
* |
* |
* DD2 (input/output) DOUBLE PRECISION |
* DD2 (input/output) DOUBLE PRECISION |
* |
* |
* DX1 (input/output) DOUBLE PRECISION |
* DX1 (input/output) DOUBLE PRECISION |
* |
* |
* DY1 (input) DOUBLE PRECISION |
* DY1 (input) DOUBLE PRECISION |
* |
* |
Line 71
|
Line 71
|
IF (.NOT.DP2.EQ.ZERO) GO TO 20 |
IF (.NOT.DP2.EQ.ZERO) GO TO 20 |
DFLAG = -TWO |
DFLAG = -TWO |
GO TO 260 |
GO TO 260 |
* REGULAR-CASE.. |
|
20 CONTINUE |
20 CONTINUE |
|
* REGULAR-CASE.. |
DP1 = DD1*DX1 |
DP1 = DD1*DX1 |
DQ2 = DP2*DY1 |
DQ2 = DP2*DY1 |
DQ1 = DP1*DX1 |
DQ1 = DP1*DX1 |
Line 108
|
Line 108
|
DX1 = DY1*DU |
DX1 = DY1*DU |
* GO SCALE-CHECK |
* GO SCALE-CHECK |
GO TO 100 |
GO TO 100 |
* PROCEDURE..ZERO-H-D-AND-DX1.. |
|
60 CONTINUE |
60 CONTINUE |
|
* PROCEDURE..ZERO-H-D-AND-DX1.. |
DFLAG = -ONE |
DFLAG = -ONE |
DH11 = ZERO |
DH11 = ZERO |
DH12 = ZERO |
DH12 = ZERO |
Line 121
|
Line 121
|
DX1 = ZERO |
DX1 = ZERO |
* RETURN.. |
* RETURN.. |
GO TO 220 |
GO TO 220 |
* PROCEDURE..FIX-H.. |
|
70 CONTINUE |
70 CONTINUE |
|
* PROCEDURE..FIX-H.. |
IF (.NOT.DFLAG.GE.ZERO) GO TO 90 |
IF (.NOT.DFLAG.GE.ZERO) GO TO 90 |
* |
* |
IF (.NOT.DFLAG.EQ.ZERO) GO TO 80 |
IF (.NOT.DFLAG.EQ.ZERO) GO TO 80 |
Line 135
|
Line 135
|
DH12 = ONE |
DH12 = ONE |
DFLAG = -ONE |
DFLAG = -ONE |
90 CONTINUE |
90 CONTINUE |
GO TO IGO(120,150,180,210) |
GO TO (150,180,210) IGO |
* PROCEDURE..SCALE-CHECK |
GO TO 120 |
100 CONTINUE |
100 CONTINUE |
|
* PROCEDURE..SCALE-CHECK |
110 CONTINUE |
110 CONTINUE |
IF (.NOT.DD1.LE.RGAMSQ) GO TO 130 |
IF (.NOT.DD1.LE.RGAMSQ) GO TO 130 |
IF (DD1.EQ.ZERO) GO TO 160 |
IF (DD1.EQ.ZERO) GO TO 160 |
ASSIGN 120 TO IGO |
IGO = 0 |
* FIX-H.. |
* FIX-H.. |
GO TO 70 |
GO TO 70 |
120 CONTINUE |
120 CONTINUE |
Line 153
|
Line 154
|
130 CONTINUE |
130 CONTINUE |
140 CONTINUE |
140 CONTINUE |
IF (.NOT.DD1.GE.GAMSQ) GO TO 160 |
IF (.NOT.DD1.GE.GAMSQ) GO TO 160 |
ASSIGN 150 TO IGO |
IGO = 1 |
* FIX-H.. |
* FIX-H.. |
GO TO 70 |
GO TO 70 |
150 CONTINUE |
150 CONTINUE |
Line 166
|
Line 167
|
170 CONTINUE |
170 CONTINUE |
IF (.NOT.DABS(DD2).LE.RGAMSQ) GO TO 190 |
IF (.NOT.DABS(DD2).LE.RGAMSQ) GO TO 190 |
IF (DD2.EQ.ZERO) GO TO 220 |
IF (DD2.EQ.ZERO) GO TO 220 |
ASSIGN 180 TO IGO |
IGO = 2 |
* FIX-H.. |
* FIX-H.. |
GO TO 70 |
GO TO 70 |
180 CONTINUE |
180 CONTINUE |
Line 177
|
Line 178
|
190 CONTINUE |
190 CONTINUE |
200 CONTINUE |
200 CONTINUE |
IF (.NOT.DABS(DD2).GE.GAMSQ) GO TO 220 |
IF (.NOT.DABS(DD2).GE.GAMSQ) GO TO 220 |
ASSIGN 210 TO IGO |
IGO = 3 |
* FIX-H.. |
* FIX-H.. |
GO TO 70 |
GO TO 70 |
210 CONTINUE |
210 CONTINUE |
Line 186
|
Line 187
|
DH22 = DH22*GAM |
DH22 = DH22*GAM |
GO TO 200 |
GO TO 200 |
220 CONTINUE |
220 CONTINUE |
IF (DFLAG) 250,230,240 |
IF (DFLAG.LT.ZERO) THEN |
|
GO TO 250 |
|
ELSE IF (DFLAG.EQ.ZERO) THEN |
|
GO TO 230 |
|
ELSE |
|
GO TO 240 |
|
END IF |
230 CONTINUE |
230 CONTINUE |
DPARAM(3) = DH21 |
DPARAM(3) = DH21 |
DPARAM(4) = DH12 |
DPARAM(4) = DH12 |