Trailing-Edge
-
PDP-10 Archives
-
DEC-20-OT20A-D-MC9
-
2-exercise-libary/aritop.for
There are no other files named aritop.for in the archive.
SUBROUTINE ARITOP
COMMON B1,B2,B3,B4,B5,B6,B7,B8,B9,B10,B11,B12,B13
COMMON D1,D2,D3,D4,D5,D6,D7,D8,D9,D10
COMMON E1,E2,E3,E4,E5,E6,E7,E8,E9,E10
COMMON IFA(10),IFB(50),IFC(100),IFD(10,10),IFE(50,50),
1 IFF(10,10,10)
COMMON IP0,IP1,IP2,IP3,IP5,IP6,IP7,IP8,IP9,IP10,IP11,IP12,IP13,
1 IP14,IP15,IP16
COMMON I1,I2,I3,I4,I5,I6,I7,I8,I9,I10,I11,I12,I13,I14,I15,I16,
1 I17,I18,I19,I20,I21,I22,I23,I24,I25,I26,I27,I28,I29,I30,
2 I31,I32,I33,I34,I35,I36,I37,I38,I39,I40
COMMON L1,L2,L3,L4,L5,L6,L7,L8,L9,L10
COMMON R1,R2,R3,R4,R5,R6,R7,R8,R9,R10
COMMON S1,S2,S3,S4,S5,S6,S7,S8,S9,S10
COMMON CRA,PRA,FL1,FL2,PR
DOUBLE PRECISION S1,S2,S3,S4,S5,S6,S7,S8,S9,S10
DOUBLE PRECISION D1,D2,D3,D4,D5,D6,D7,D8,D9,D10
LOGICAL B1,B2,B3,B4,B5,B6,B7,B8,B9,B10,B11,B12,B13
200 IC=0
201 IF(IC-IP8) 202,203,203
202 R1=E1*E1+E2*E2+E3*E3+E4*E4-E6/E2+E5
R2=E1-E8-E2/E4-E3/E1-E4*E5+E6*E6+E7
R3=E1*E2+E3*E4+E5*E6+E7*E8-E9/E1+R1
R4=E1-E2-E3/E2-E4/E1-E5*E5+E6*E7+E8
R5=E1*E3+E2*E4+E5*E7+E8*E8-E9/E1+E6
R6=E1-E3-E4/E2-E2/E1-E5*E8+E6*E7+E9
R7=E9*E1+E8*E2+E7*E3+E6*E4-E5/E1+R1
R8=E9-E8-E7/E1-E6/E3-E5*E3+E4*E1+E2
R9=E8*E1+E9*E2+E7*E4+E3*E6-E5/E2+R1
R10=E6-E2-E3/E4-E1/E2-E5*E8+E4*E9+E7
IC=IC+1
GOTO 201
203 IC=0
204 IF(IC-IP9) 205,206,206
205 S1=D1*D1+D2*D2+D3*D3+D4*D4-D6/D2+D5
S2=D1-D8-D2/D4-D3/D1-D4*D5+D6*D6+D7
S3=D1*D2+D3*D4+D5*D6+D7*D8-D9/D1+S1
S4=D1-D2-D3/D2-D4/D1-D5*D5+D6*D7+D8
S5=D1*D3+D2*D4+D5*D7+D8*D8-D9/D1+D6
S6=D1-D3-D4/D2-D2/D1-D5*D8+D6*D7+D9
S7=D9*D1+D8*D2+D7*D3+D6*D4-D5/D1+S1
S8=D9-D8-D7/D1-D6/D3-D5*D3+D4*D1+D2
S9=D8*D1+D9*D2+D7*D4+D3*D6-D5/D2+S1
S10=D6-D2-D3/D4-D1/D2-D5*D8+D4*D9+D7
IC=IC+1
GOTO 204
206 IC=0
207 IF(IC-IP10) 208,209,209
208 J1=E1
J2=E2
J3=E3
J4=E4
J5=E5
J6=E6
J7=E7
T1=J1
T2=J2
T3=J3
T4=J4
T5=J5
T6=J6
T7=J7
J1=E1
J2=E2
J3=E3
J4=E4
J5=E5
J6=E6
J7=E7
T1=J1
T2=J2
T3=J3
T4=J4
T5=J5
T6=J6
T7=J7
IC=IC+1
GOTO 207
209 IC=0
210 IF(IC-IP11) 211,212,212
211 L1=I1+I2-I3+I4*I5-I6*I7+I8*I9+I10+I11
L2=I1-I2+I3+I4*I7-I5*I6-I8+I9+I11+I12
L3=I1+I4-I5+I2*I3-I8*I11+I10*I9+I6+I7
L4=I1-I4+I5+I7*I9-I11*I10-I2+I8+I3+I6
L5=I11+I12-I13+I14*I15-I16*I17+I18*I19+I20+I21
L6=I11-I12+I13+I14*I17-I15*I16-I18+I19+I21+I22
L7=I11+I14-I15+I12*I13-I18*I21+I20*I19+I16+I17
L8=I10-I14+I15+I17*I19-I21*I20-I12+I18+I13+I16
L9=I21+I22-I23+I24*I25-I26*I27+I28*I29+I30+I31
L10=I21-I22+I23+I24*I27-I25*I26-I28+I29+I31+I32
IC=IC+1
GOTO 210
212 IC=0
213 IF(IC-IP12) 214,241,241
214 ID=1
215 ID1=11-ID
IFA(ID1)=IFA(ID)
ID=ID+1
IF(ID-11) 215,216,216
216 ID=1
217 ID1=51-ID
IFB(ID1)=IFB(ID)
ID=ID+1
IF(ID-51) 217,218,218
218 JD=1
219 ID=1
220 ID1=101-ID
IFC(ID1)=IFC(ID)
ID=ID+1
IF(ID-101) 220,221,221
221 JD=JD+1
IF(JD-3) 219,222,222
222 JD=1
223 ID=1
224 ID1=11-ID
KD=1
225 KD1=11-KD
IFD(ID1,KD1)=IFD(ID,KD)
KD=KD+1
IF(KD-11) 225,226,226
226 ID=ID+1
IF(ID-11) 224,227,227
227 JD=JD+1
IF(JD-3) 223,228,228
228 JD=1
229 ID=1
230 ID1=51-ID
KD=1
231 KD1=51-KD
IFE(ID1,KD1)=IFE(ID,KD)
KD=KD+1
IF(KD-51) 231,232,232
232 ID=ID+1
IF(ID-51) 230,233,233
233 JD=JD+1
IF(JD-4) 229,234,234
234 ID=1
235 ID1=11-ID
JD=1
236 JD1=11-JD
KD=1
237 KD1=11-KD
IFF(ID1,JD1,KD1)=IFF(ID,JD,KD)
KD=KD+1
IF(KD-11) 237,238,238
238 JD=JD+1
IF(JD-11) 236,239,239
239 ID=ID+1
IF(ID-11) 235,240,240
240 IC=IC+1
GOTO 213
241 IC=0
242 IF(IC-IP13) 243,244,244
243 B7=(B1.AND.(.NOT.B2)).OR.((B1.OR.B2.OR.B3).OR.((.NOT.B1).AND.
1(.NOT.B2).AND.(.NOT.B3)).AND.((.NOT.B4).OR.(B5.AND.B6)))
B8=(.NOT.(B3.AND.(.NOT.B4))).AND.((B2.OR.B3.OR.B4).OR.(B1.OR.(
1.NOT.B3).AND.(.NOT.B4)).AND.((.NOT.B5).AND.(B6.OR.B1)))
B9=(B3.AND.(.NOT.B2)).OR.((B3.OR.B1.OR.B4).OR.((.NOT.B3).AND.
1(.NOT.B1).AND.(.NOT.B4)).AND.((.NOT.B6).OR.(B2.AND.B4)))
B10=(.NOT.(B5.AND.(.NOT.B2))).AND.((B1.OR.B2.OR.B5).OR.(B1.OR.
1(.NOT.B5).AND.(.NOT.B3)).AND.((.NOT.B6).AND.(B3.OR.B2)))
B11=(B5.AND.(.NOT.B4)).OR.((B6.OR.B2.OR.B4).OR.((.NOT.B2).AND.