GAP: ASSEMBLE REVISION 072 OF AGC PROGRAM ARTEMIS BY NASA 2021114-011 11:40 FEB. 26, 1971Note that the date is the date of the printout, not the date of the program revision.
Source Code
|
|
These source-code files were obtained by digitally photographing
an Artemis 72 (Apollo 15-17 Command Module) program listing from a private
collection and then
modifying pre-existing Comanche 55 (Apollo 11) source files to
incorporate changes. Photography was by Ron Burkey, and transcription
was performed by a team of volunteers. Note that the page images
presented online are of reduced quality, and that additional images
may be available. Report any conversion errors or legibility
problems in page images to info@sandroid.org.
Notations on the program listing read, in part:GAP: ASSEMBLE REVISION 072 OF AGC PROGRAM ARTEMIS BY NASA 2021114-011 11:40 FEB. 26, 1971Note that the date is the date of the printout, not the date of the program revision. |
048385,000002: ## Copyright: Public domain.
048386,000003: ## Filename: INTERPRETER.agc
048387,000004: ## Purpose: Part of the source code for Artemis (i.e., Colossus 3),
048388,000005: ## build 072. This is for the Command Module's (CM)
048389,000006: ## Apollo Guidance Computer (AGC), for
048390,000007: ## Apollo 15-17.
048391,000008: ## Assembler: yaYUL
048392,000009: ## Contact: Jim Lawton <jim DOT lawton AT gmail DOT com>
048393,000010: ## Website: www.ibiblio.org/apollo/index.html
048394,000011: ## Page Scans: www.ibiblio.org/apollo/ScansForConversion/Artemis072/
048395,000012: ## Mod history: 2009-09-02 JL Adapted from corresponding Comanche 055 file.
048396,000013: ## 2009-09-16 JL Fixed line comment alignment in top section of file.
048397,000014: ## 2010-02-20 RSB Un-##'d this header.
048398,000015: ## 2011-05-07 JL Removed workarounds.
048399,000016: ## 2017-02-05 RSB Proofed comment text by diff'ing vs Comanche 55
048400,000017: ## and corrected errors found.
048401,000018: ## 2017-03-15 RSB Comment-text fixes identified in 5-way
048402,000019: ## side-by-side diff of Luminary 69/99/116/131/210.
048403,000020:
|
|
Page 1108 |
048405,000022:
048406,000023: # SECTION 1 DISPATCHER
048407,000024:
048408,000025: # ENTRY TO THE INTERPRETER. INTPRET SETS LOC TO THE FIRST INSTRUCTION, BANKSET TO THE BBANK OF THE
048409,000026: # OBJECT INTERPRETIVE PROGRAM, AND INTBIT15 TO THE BIT15 CONTENTS OF FBANK. INTERPRETIVE PROGRAMS MAY BE IN
048410,000027: # VIRTUALLY ALL BANKS PRESENT UNDER ANY SUPER-BANK SETTING, WITH THE RESTRICTION THAT PROGRAMS IN HIGH BANKS
048411,000028: # (BIT15 OF FBANK = 1) DO NOT REFER TO LOWBANKS, AND VICE-VERSA. THE INTERPRETER DOES NOT SWITCH SUPERBANKS.
048412,000029: # E-BANK SWITCHING OCCURS WHENEVER GENERAL ERASABLE (100-3777) IS ADDRESSED.
048413,000030:
048414,000031: 6000 SETLOC INTERF3
048415,000032: 6000 BANK
048416,000033:
048417,000034: 6006 COUNT* $$/INTER
048418,000035: 6006 00003 INTPRET RELINT
048419,000036: 6007 00006 EXTEND # SET LOC TO THE WORD FOLLOWING THE TC.
048420,000037: 6010 22164 +2 QXCH LOC
048421,000038:
048422,000039: 6011 30006 +3 CA BBANK # INTERPRETIVE BRANCHES FINISH HERE.
048423,000040: 6012 54165 TS BANKSET
048424,000041: 6013 74752 MASK BIT15 # GET 15TH BIT FOR INDEXABLE ADDRESSES.
048425,000042: 6014 54115 TS INTBIT15
048426,000043:
048427,000044: 6015 54023 TS EDOP # MAKE SURE NO INSTRUCTIONS LEFT OVER
048428,000045:
048429,000046: 6016 16037 TCF NEWOPS # PICK UP OP CODE PAIR AND BEGIN.
048430,000047:
048431,000048: 6017 22006 INTRSM LXCH BBANK # RESUME SUSPENDED INTERPRETIVE JOB
048432,000049: 6020 16011 TCF INTPRET +3
048433,000050:
048434,000051: # DLOAD LOADS MPAC, MPAC +1, LEAVING ZERO IN MPAC +2.
048435,000052:
048436,000053: 6021 00006 DLOAD EXTEND
048437,000054: 6022 50116 INDEX ADDRWD
048438,000055: 6023 30001 DCA 0 # LOAD DP C(C(ADDRWD)) INTO MPAC,MPAX +1
048439,000056: 6024 52155 SLOAD2 DXCH MPAC
048440,000057: 6025 34772 CAF ZERO # ZERO MPAC +2
048441,000058:
|
|
Page 1109 |
048443,000060: # AT THE END OF MOST INSTRUCTIONS, CONTROL IS GIVEN TO DANZIG TO DISPATCH THE NEXT OPERATION.
048444,000061:
048445,000062: 6026 54156 +2 TS MPAC +2 # AND DECLARE DP MODE
048446,000063: 6027 54163 NEWMODE TS MODE # PROLOGUE FOR MODE-CHANGING INSTRUCTIONS.
048447,000064:
048448,000065: 6030 30165 DANZIG CA BANKSET # SET BBANK BEFORE TESTING NEWJOB SO THAT
048449,000066: 6031 54006 TS BBANK # IT MAY BE SAVED DIRECTLY BY CHANJOB.
048450,000067:
048451,000068: 6032 10023 NOIBNKSW CCS EDOP # SEE IF AN ORDER CODE IS LEFT OVER FROM
048452,000069: 6033 16046 TCF OPJUMP # THE LAST PAIR RETRIEVED. IF SO, EXECUTE.
048453,000070: # EDOP IS SET TO ZERO ON ITS RE-EDITING.
048454,000071:
048455,000072: 6034 10067 CCS NEWJOB # SEE IF A JOB OF HIGHER PRIORITY IS
048456,000073: 6035 15156 TCF CHANG2 # PRESENT, AND IF SO, CHANGE JOBS.
048457,000074:
048458,000075: 6036 24164 INCR LOC # ADVANCE THE LOCATION COUNTER.
048459,000076: # ITRACE (1) REFERS TO "NEWOPS".
048460,000077: 6037 50164 NEWOPS INDEX LOC # ENTRY TO BEGIN BY PICKING OP CODE PAIR.
048461,000078: 6040 30000 CA 0 # MAY BE AN OPCODE PAIR OR A STORE CODE.
048462,000079: 6041 10000 CCS A # TEST SIGN AND GET DABS(A).
048463,000080: 6042 16331 TCF DOSTORE # PROCESS STORE CODE.
048464,000081:
048465,000082: 6043 00177 LOW7 OCT 177
048466,000083:
048467,000084: 6044 54023 TS EDOP # OP CODE PAIR. LEAVE THE OTHER IN EDOP
048468,000085: 6045 76043 MASK LOW7 # WHERE CCS EDOP WILL HONOR IT NEXT.
048469,000086:
048470,000087: 6046 54020 OPJUMP TS CYR # LOWWD ENTERS HERE IF A RIGHT-HAND OP
048471,000088: 6047 10020 CCS CYR # CODE IS TO BE PROCESSED. TEST PREFICES.
048472,000089: 6050 16216 TCF OPJUMP2 # TEST SECOND PREFIX BIT.
048473,000090:
048474,000091: 6051 16712 TCF EXIT # +0 OP CODE IS EXIT.
|
|
Page 1110 |
048476,000093: # PROCESS ADDRESSES WHICH MAY BE DIRECT, INDEXED, OR REFERENCE THE PUSHDOWN LIST.
048477,000094:
048478,000095: 6052 74770 ADDRESS MASK BIT1 # SEE IF ADDRESS IS INDEXED. CYR CONTAINED
048479,000096: 6053 10000 CCS A # 400XX, SO BIT 1 IS NOW AS IT WAS IN CYR.
048480,000097: 6054 16115 TCF INDEX # FORM INDEXED ADDRESS.
048481,000098:
048482,000099: 6055 50164 DIRADRES INDEX LOC # LOOK AHEAD TO NEXT WORD TO SEE IF
048483,000100: 6056 40001 OCT40001 CS 1 # ADDRESS IS GIVEN.
048484,000101: 6057 10000 CCS A
048485,000102: 6060 16164 TCF PUSHUP # IF NOT.
048486,000103:
048487,000104: 6061 77773 NEG4 DEC -4 B-14
048488,000105:
048489,000106: 6062 24164 INCR LOC # IF SO, TO SHOW WE PICKED UP A WORD.
048490,000107: 6063 54116 TS ADDRWD
048491,000108:
|
|
Page 1111 |
048493,000110: # FINAL DIGESTION OF DIRECT ADDRESSES OF OP CODES WITH 01 PREFIX IS DONE HERE. IN EACH CASE, THE
048494,000111: # REQUIRED 12-BIT SUB-ADDRESS IS LEFT IN ADDRWD, WITH ANY REQUIRED E OR F BANK SWITCHING DONE. ADDRESSES LESS
048495,000112: # THAN 45D ARE TAKEN TO BE RELATIVE TO THE WORK AREA. THE OP CODE IS NOW IN BITS 1-5 OF CYR WITH BIT 14 = 1.
048496,000113:
048497,000114: 6064 66220 AD -ENDVAC # SEE IF ADDRESS RELATIVE TO WORK AREA.
048498,000115: 6065 10000 CCS A
048499,000116: 6066 67676 AD -ENDERAS # IF NOT, SEE IF IN GENERAL ERASABLE.
048500,000117: 6067 16074 TCF IERASTST
048501,000118:
048502,000119: 6070 30120 NETZERO CA FIXLOC # IF SO, LEAVE THE MODIFIED ADDRESS IN
048503,000120: 6071 26116 ADS ADDRWD # ADDRWD AND DISPATCH.
048504,000121: 6072 50020 ITR15 INDEX CYR # THIS INDEX MAKES THE NEXT INSTRUCTION
048505,000122: 6073 76242 7 INDJUMP -1 # TCF INDJUMP + OP, EDITING CYR.
048506,000123:
048507,000124: 6074 00006 IERASTST EXTEND
048508,000125: 6075 66105 BZMF GEADDR # GO PROCESS GENERAL-ERASABLE ADDRESS.
048509,000126:
048510,000127: 6076 75041 MASK LOW10 # FIXED BANK ADDRESS. RESTORE AND ADD B15.
048511,000128: 6077 65041 AD LOW10 # SWITCH BANKS AND LEAVE SUBADDRESS IN
048512,000129: 6100 56116 XCH ADDRWD # ADDRWD FOR OPERAND RETRIEVAL. (THIS
048513,000130: 6101 60115 AD INTBIT15 # METHOD PRECLUDES USE OF THE LAST
048514,000131: 6102 54004 TS FBANK # LOCATION IN EACH FBANK.)
048515,000132: 6103 50020 ITR12 INDEX CYR
048516,000133: 6104 76242 7 INDJUMP -1
048517,000134:
048518,000135: 6105 74350 GEADDR MASK LOW8
048519,000136: 6106 65024 AD OCT1400
048520,000137: 6107 56116 XCH ADDRWD
048521,000138: 6110 54003 TS EBANK
048522,000139: 6111 50020 ITR10 INDEX CYR
048523,000140: 6112 76242 7 INDJUMP -1
048524,000141:
|
|
Page 1112 |
048526,000143: # THE FOLLOWING ROUTINE PROCESSES INTERPRETIVE INDEXED ADDRESSES. AN INTERPRETER INDEX REGISTER MAY
048527,000144: # CONTAIN THE ADDRESS OF ANY ERASABLE REGISTER (0-42 BEING RELATIVE TO THE VAC AREA) OR ANY INTERPRETIVE PROGRAM
048528,000145: # BANK, OR ANY INTEGER IN THAT RANGE.
048529,000146:
048530,000147: 6113 37666 DODLOAD* CAF DLOAD* # STODL* COMES HERE TO PROCESS LOAD ADR.
048531,000148: 6114 54020 +1 TS CYR # (STOVL* ENTERS HERE).
048532,000149: 6115 30120 INDEX CA FIXLOC # SET UP INDEX LOCATION.
048533,000150: 6116 54130 TS INDEXLOC
048534,000151: 6117 24164 INCR LOC # (ADDRESS ALWAYS GIVEN).
048535,000152: 6120 50164 INDEX LOC
048536,000153: 6121 40000 CS 0
048537,000154: 6122 10000 CCS A # INDEX 2 IF ADDRESS STORED COMPLEMENTED.
048538,000155: 6123 24130 INCR INDEXLOC
048539,000156: 6124 16125 NOOP
048540,000157:
048541,000158: 6125 54116 TS ADDRWD # 14 BIT ADDRESS TO ADDRWD.
048542,000159: 6126 77675 MASK HIGH4 # IF ADDRESS GREATER THAN 2K, ADD INTBIT15
048543,000160: 6127 00006 EXTEND
048544,000161: 6130 16133 BZF INDEX2
048545,000162: 6131 30115 CA INTBIT15
048546,000163: 6132 26116 ADS ADDRWD
048547,000164:
048548,000165: 6133 50130 INDEX2 INDEX INDEXLOC
048549,000166: 6134 40046 CS X1
048550,000167: 6135 26116 ADS ADDRWD # DO AUGMENT, IGNORING AND CORRECTING OVF.
048551,000168:
048552,000169: 6136 77677 MASK HIGH9 # SEE IF ADDRESS IS IN WORK AREA.
048553,000170: 6137 00006 EXTEND
048554,000171: 6140 16153 BZF INDWORK
048555,000172: 6141 77675 MASK HIGH4 # SEE IF IN FIXED BANK.
048556,000173: 6142 00006 EXTEND
048557,000174: 6143 16155 BZF INDERASE
048558,000175:
048559,000176: 6144 30116 CA ADDRWD # IN FIXED - SWITCH BANKS AND CREATE
048560,000177: 6145 54004 TS FBANK # SUB-ADDRESS.
048561,000178: 6146 75041 MASK LOW10
048562,000179: 6147 64756 AD 2K
048563,000180: 6150 54116 TS ADDRWD
048564,000181: 6151 50020 ITR11 INDEX CYR
048565,000182: 6152 36242 3 INDJUMP -1
048566,000183:
048567,000184: 6153 30120 INDWORK CA FIXLOC # MAKE ADDRWD RELATIVE TO WORK AREA.
048568,000185: 6154 16161 TCF ITR13 -1
048569,000186:
048570,000187: 6155 35024 INDERASE CA OCT1400
048571,000188: 6156 56116 XCH ADDRWD
048572,000189: 6157 54003 TS EBANK
048573,000190: 6160 74350 MASK LOW8
048574,000191: 6161 26116 -1 ADS ADDRWD
048575,000192: 6162 50020 ITR13 INDEX CYR
|
|
Page 1113 |
048577,000194: 6163 36242 3 INDJUMP -1
048578,000195:
|
|
Page 1114 |
048580,000197: # PUSH-UP ROUTINES. WHEN NO OPERAND ADDRESS IS GIVEN, THE APPROPRIATE OPERAND IS TAKEN FROM THE PUSH-DOWN
048581,000198: # LIST. IN MOST CASES THE MODE OF THE RESULT (VECTOR OR SCALAR) OF THE LAST ARITHMETIC OPERATION PERFORMED
048582,000199: # IS THE SAME AS THE TYPE OF OPERAND DESIRED (ALL ADD/SUBTRACT ETC.). EXCEPTIONS TO THIS GENERAL RULE ARE LISTED
048583,000200: # BELOW (NOTE THAT IN EVERY CASE THE MODE REGISTER IS LEFT INTACT):
048584,000201:
048585,000202: # 1. VXSC AND V/SC WANT THE OPPOSITE TYPE OF OPERAND, E.G., IF THE LAST OPERATION YIELDED A VECTOR
048586,000203: # RESULT, VXSC WANTS A SCALAR.
048587,000204:
048588,000205: # 2. THE LOAD CODES SHOULD LOAD THE ACCUMULATOR INDEPENDENT OF THE RESULT OF THE LAST OPERATION. THIS
048589,000206: # INCLUDES VLOAD, DLOAD, TLOAD, PDDL, AND PDVL (NO PUSHUP WITH SLOAD).
048590,000207:
048591,000208: # 3. SOME ARITHMETIC OPERATIONS REQUIRE A STANDARD TYPE OF OPERAND REGARDLESS OF THE PREVIOUS OPERATION.
048592,000209: # THIS INCLUDES SIGN WANTING DP AND TAD REQUIRING TP.
048593,000210:
048594,000211: 6164 34351 PUSHUP CAF OCT23 # IF THE LOW 5 BITS OF CYR ARE LESS THAN
048595,000212: 6165 70020 MASK CYR # 20, THIS OP REQUIRES SPECIAL ATTENTION.
048596,000213: 6166 66171 AD -OCT10 # (NO -0).
048597,000214: 6167 10000 CCS A
048598,000215: 6170 16202 TCF REGUP # FOR ALL CODES GREATER THAN OCT 7.
048599,000216:
048600,000217: 6171 77767 -OCT10 OCT -10
048601,000218:
048602,000219: 6172 66061 AD NEG4 # WE NOW HAVE 7 - OP CODE(MOD4). SEE IF
048603,000220: 6173 10000 CCS A # THE OP CODE (MOD4) IS THREE (REVERSE).
048604,000221: 6174 50000 INDEX A # NO - THE MODE IS DEFINITE. PICK UP THE
048605,000222: 6175 46213 CS NO.WDS
048606,000223: 6176 16204 TCF REGUP +2
048607,000224:
048608,000225: 6177 50163 INDEX MODE # FOR VXSC AND V/SC WE WANT THE REQUIRED
048609,000226: 6200 46211 CS REVCNT # PUSHLOC DECREMENT WITHOUT CHANGING THE
048610,000227: 6201 16204 TCF REGUP +2 # MODE AT THIS TIME.
048611,000228:
048612,000229: 6202 50163 REGUP INDEX MODE # MOST ALL OP CODES PUSHUP HERE.
048613,000230: 6203 46213 CS NO.WDS
048614,000231: 6204 26166 +2 ADS PUSHLOC
048615,000232: 6205 54116 TS ADDRWD
048616,000233: 6206 50020 ITR14 INDEX CYR
048617,000234: 6207 76242 7 INDJUMP -1 # (THE INDEX MAKES THIS A TCF.)
048618,000235:
048619,000236: 6210 00002 OCT 2 # REVERSE PUSHUP DECREMENT. VECTOR TAKES 2
048620,000237: 6211 00006 REVCNT OCT 6 # WORDS, SCALAR TAKES 6.
048621,000238: 6212 00006 OCT 6
048622,000239: 6213 00002 NO.WDS OCT 2 # CONVENTIONAL DECREMENT IS 6 WORDS VECTOR
048623,000240: 6214 00003 OCTAL3 OCT 3 # 2 IN DP, AND 3 IN TP.
048624,000241: 6215 00006 OCT 6
048625,000242:
|
|
Page 1115 |
048627,000244: # TEST THE SECOND PREFIX BIT TO SEE IF THIS IS A MISCELLANEOUS OR A UNARY/SHORT SHIFT OPERATION.
048628,000245:
048629,000246: 6216 10020 OPJUMP2 CCS CYR # TEST SECOND PREFIX BIT.
048630,000247: 6217 16232 TCF OPJUMP3 # TEST THIRD BIT TO SEE IF UNARY OR SHIFT.
048631,000248:
048632,000249: 6220 77722 -ENDVAC DEC -45 B-14
048633,000250:
048634,000251: # THE FOLLOWING ROUTINE PROCESSES ADDRESSES OF SUFFIX CLASS 10. THEY ARE BASICALLY WORK AREA ADDRESSES
048635,000252: # IN THE RANGE 0 - 52, ERASABLE ECADR CONSTANTS FROM 100 - 3777, AND FCADRS ABOVE THAT. ALL 15 BITS ARE AVAILABLE
048636,000253: # IN CONTRAST TO SUFFIX 1, IN WHICH ONLY THE LOW ORDER 14 ARE AVAILABLE.
048637,000254:
048638,000255: 6221 24164 15BITADR INCR LOC # (ENTRY HERE FROM STCALL).
048639,000256: 6222 50164 INDEX LOC # PICK UP ADDRESS WORD.
048640,000257: 6223 30000 CA 0
048641,000258: 6224 54117 TS POLISH # WE MAY NEED A SUBADDRESS LATER.
048642,000259:
048643,000260: 6225 35042 CAF LOW7+2K # THESE INSTRUCTIONS ARE IN BANK 1.
048644,000261: 6226 54004 TS FBANK
048645,000262: 6227 70020 MASK CYR
048646,000263: 6230 50000 ITR7 INDEX A
048647,000264: 6231 16303 TCF MISCJUMP
048648,000265:
|
|
Page 1116 |
048650,000267: # COMPLETE THE DISPATCHING OF UNARY AND SHORT SHIFT OPERATIONS.
048651,000268:
048652,000269: 6232 54004 OPJUMP3 TS FBANK # CALL IN BANK 0 (BITS 11-15 OF A ARE 0.)
048653,000270:
048654,000271: 6233 10020 CCS CYR # TEST THIRD PREFIX BIT.
048655,000272: 6234 50000 INDEX A # THE DECREMENTED UNARY CODE IS IN BITS
048656,000273: 6235 12000 TCF UNAJUMP # 1-4 OF A (ZERO, EXIT, HAS BEEN DETECTED)
048657,000274:
048658,000275: 6236 10163 CCS MODE # ITS A SHORT SHIFT CODE. SEE IF PRESENT
048659,000276: 6237 12017 TCF SHORTT # SCALAR OR VECTOR.
048660,000277: 6240 12017 TCF SHORTT
048661,000278: 6241 12121 TCF SHORTV # CALLS THE APPROPRIATE ROUTINE.
048662,000279:
048663,000280: 6242 4344 FBANKMSK EQUALS BANKMASK
048664,000281: 6242 00122 LVBUF ADRES VBUF
048665,000282:
|
|
Page 1117 |
048667,000284: # THE FOLLOWING IS THE JUMP TABLE FOR OP CODES WHICH MAY HAVE INDEXABLE ADDRESSES OR MAY PUSH UP.
048668,000285:
048669,000286: 6243 16454 INDJUMP TCF VLOAD # 00 - LOAD MPAC WITH A VECTOR.
048670,000287: 6244 17040 TCF TAD # 01 - TRIPLE PRECISION ADD TO MPAC.
048671,000288: 6245 17611 TCF SIGN # 02 - COMPLEMENT MPAC (V OR SC) IF X NEG.
048672,000289: 6246 17335 TCF VXSC # 03 - VECTOR TIMES SCALAR.
048673,000290: 6247 16652 TCF CGOTO # 04 - COMPUTED GO TO.
048674,000291: 6250 16437 TCF TLOAD # 05 - LOAD MPAC WITH TRIPLE PRECISION.
048675,000292: 6251 16021 TCF DLOAD # 06 - LOAD MPAC WITH A DP SCALAR.
048676,000293: 6252 17560 TCF V/SC # 07 - VECTOR DIVIDED BY SCALAR.
048677,000294:
048678,000295: 6253 16450 TCF SLOAD # 10 - LOAD MPAC IN SINGLE PRECISION.
048679,000296: 6254 16567 TCF SSP # 11 - SET SINGLE PRECISION INTO X.
048680,000297: 6255 16472 TCF PDDL # 12 - PUSH DOWN MPAC AND RE-LOAD IN DP.
048681,000298: 6256 17270 TCF MXV # 13 - MATRIX POST-MULTIPLIED BY VECTOR.
048682,000299: 6257 16526 TCF PDVL # 14 - PUSH DOWN AND VECTOR LOAD.
048683,000300: 6260 16575 TCF CCALL # 15 - COMPUTED CALL.
048684,000301: 6261 17273 TCF VXM # 16 - MATRIX PRE-MULTIPLIED BY VECTOR.
048685,000302: 6262 17552 TCF TSLC # 17 - NORMALIZE MPAC (SCALAR ONLY).
048686,000303:
048687,000304: 6263 17530 TCF DMPR # 20 - DP MULTIPLY AND ROUND.
048688,000305: 6264 17533 TCF DDV # 21 - DP DIVIDE BY.
048689,000306: 6265 17537 TCF BDDV # 22 - DP DIVIDE INTO.
048690,000307: 6266 17555 TCF GSHIFT # 23 - GENERAL SHIFT INSTRUCTION.
048691,000308: 6267 16720 TCF VAD # 24 - VECTOR ADD.
048692,000309: 6270 16716 TCF VSU # 25 - VECTOR SUBTRACT.
048693,000310: 6271 17005 TCF BVSU # 26 - VECTOR SUBTRACT FROM.
048694,000311: 6272 17265 TCF DOT # 27 - VECTOR DOT PRODUCT.
048695,000312:
048696,000313: 6273 17414 TCF VXV # 30 - VECTOR CROSS PRODUCT.
048697,000314: 6274 17361 TCF VPROJ # 31 - VECTOR PROJECTION.
048698,000315: 6275 16754 TCF DSU # 32 - DP SUBTRACT.
048699,000316: 6276 17031 TCF BDSU # 33 - DP SUBTRACT FROM.
048700,000317: 6277 16744 TCF DAD # 34 - DP ADD.
048701,000318: 6300 16300 TCF +0 # 35 - AVAILABLE
048702,000319: 6301 17526 TCF DMP1 # 36 - DP MULTIPLY.
048703,000320: 6302 17547 TCF SETPD # 37 - SET PUSH DOWN POINTER (DIRECT ONLY)
048704,000321:
048705,000322: # CODES 10 AND 14 MUST NOT PUSH UP. CODE 04 MAY BE USED FOR VECTOR DECLARE BEFORE PUSHUP IF DESIRED.
048706,000323:
|
|
Page 1118 |
048708,000325: # THE FOLLOWING JUMP TABLE APPLIES TO INDEX, BRANCH, AND MISCELLANEOUS INSTRUCTIONS.
048709,000326:
048710,000327: 6303 12341 MISCJUMP TCF AXT # 00 - ADDRESS TO INDEX TRUE.
048711,000328: 6304 12346 TCF AXC # 01 - ADDRESS TO INDEX COMPLEMENTED.
048712,000329: 6305 12351 TCF LXA # 02 - LOAD INDEX FROM ERASABLE.
048713,000330: 6306 12355 TCF LXC # 03 - LOAD INDEX FROM COMPLEMENT OF ERAS.
048714,000331: 6307 12361 TCF SXA # 04 - STORE INDEX IN ERASABLE.
048715,000332: 6310 12367 TCF XCHX # 05 - EXCHANGE INDEX WITH ERASABLE.
048716,000333: 6311 12403 TCF INCR # 06 - INCREMENT INDEX REGISTER.
048717,000334: 6312 12412 TCF TIX # 07 - TRANSFER ON INDEX.
048718,000335:
048719,000336: 6313 12375 TCF XAD # 10 - INDEX REGISTER ADD FROM ERASABLE.
048720,000337: 6314 12406 TCF XSU # 11 - INDEX SUBTRACT FROM ERASABLE.
048721,000338: 6315 12464 TCF BZE/GOTO # 12 - BRANCH ZERO AND GOTO.
048722,000339: 6316 12471 TCF BPL/BMN # 13 - BRANCH PLUS AND BRANCH MINUS.
048723,000340: 6317 12444 TCF RTB/BHIZ # 14 - RETURN TO BASIC AND BRANCH HI ZERO.
048724,000341: 6320 12504 TCF CALL/ITA # 15 - CALL AND STORE QPRET.
048725,000342: 6321 12513 TCF SW/ # 16 - SWITCH INSTRUCTIONS AND AVAILABLE.
048726,000343: 6322 12454 TCF BOV(B) # 17 - BRANCH ON OVERFLOW TO BASIC OR INT.
048727,000344:
|
|
Page 1119 |
048729,000346: # THE FOLLOWING JUMP TABLE APPIES TO UNARY INSTRUCTIONS.
048730,000347:
048731,000348: 00,2000 SETLOC INTERB0
048732,000349: 00,2000 BANK
048733,000350:
048734,000351: 00,2000 COUNT* $$/INTER
048735,000352: 00,2000 13206 UNAJUMP TCF SQRT # 01 - SQUARE ROOT.
048736,000353: 00,2001 13525 TCF SINE # 02 - SIN.
048737,000354: 00,2002 13514 TCF COSINE # 03 - COS.
048738,000355: 00,2003 13605 TCF ARCSIN # 04 - ARC SIN.
048739,000356: 00,2004 13607 TCF ARCCOS # 05 - ARC COS.
048740,000357: 00,2005 13173 TCF DSQ # 06 - DP SQUARE.
048741,000358: 00,2006 12116 TCF ROUND # 07 - ROUND TO DP.
048742,000359:
048743,000360: 00,2007 17624 TCF COMP # 10 - COMPLEMENT VECTOR OR SCALAR.
048744,000361: 00,2010 13231 TCF VDEF # 11 - VECTOR DEFINE.
048745,000362: 00,2011 13022 TCF UNIT # 12 - UNIT VECTOR.
048746,000363: 00,2012 13175 TCF ABVALABS # 13 - LENGTH OF VECTOR OR MAG OF SCALAR.
048747,000364: 00,2013 13244 TCF VSQ # 14 - SQUARE OF LENGTH OF VECTOR.
048748,000365: 00,2014 16323 TCF STADR # 15 - PUSH UP ON STORE CODE.
048749,000366: 00,2015 13273 TCF RVQ # 16 - RETURN VIA QPRET.
048750,000367: 00,2016 13246 TCF PUSH # 17 - PUSH MPAC DOWN.
048751,000368:
|
|
Page 1120 |
048753,000370: # SECTION 2 LOAD AND STORE PACKAGE.
048754,000371:
048755,000372: # A SET OF EIGHT STORE CODES IS PROVIDED AS THE PRIMARY METHOD OF STORING THE MULTI-PURPOSE
048756,000373: # ACCUMULATOR (MPAC). IF IN THE DANZIG SECTION LOC REFERS TO AN ALGEBRAICALLY POSITIVE WORD, IT IS TAKEN AS A
048757,000374: # STORE CODE WITH A CORRESPONDING ERASABLE ADDRESS. MOST OF THESE CODES ARE TWO ADDRESS, SPECIFYING THAT THE WORD
048758,000375: # FOLLOWING THE STORE CODE IS TO BE USED AS AN ADDRESS FROM WHICH TO RE-LOAD MPAC. FOUR OPTIONS ARE AVAILABLE:
048759,000376:
048760,000377: # 1. STORE STORE MPAC. THE E ADDRESS MAY BE INDEXED.
048761,000378: # 2. STODL STORE MPAC AND RE-LOAD IT IN DP WITH THE NEXT ADDRESS (THE LOAD MAY BE INDEXED).
048762,000379: # 3. STOVL STORE MPAC AND RE-LOAD A VECTOR (AS ABOVE).
048763,000380: # 4. STCALL STORE AND DO A CALL (BOTH ADDRESSES MUST BE DIRECT HERE).
048764,000381:
048765,000382: # STODL AND STOVL WILL TAKE FROM THE PUSH-DOWN LIST IF NO LOAD ADDRESS IS GIVEN.
048766,000383:
048767,000384: 6000 SETLOC INTERF3
048768,000385: 6000 BANK
048769,000386:
048770,000387: 6323 COUNT* $$/INTER
048771,000388: 6323 30165 STADR CA BANKSET # THE STADR CODE (PUSHUP UP ON STORE
048772,000389: 6324 54004 TS FBANK # ADDRESS) ENTERS HERE.
048773,000390: 6325 24164 INCR LOC
048774,000391: 6326 50164 ITR1 INDEX LOC # THE STORECODE WAS STORED COMPLEMENTED TO
048775,000392: 6327 40000 CS 0 # MAKE IT LOOK LIKE AN OPCODE PAIR.
048776,000393: 6330 67702 AD NEGONE # (YUL CANT REMOVE 1 BECAUSE OF EARLY CCS)
048777,000394:
048778,000395: 6331 54116 DOSTORE TS ADDRWD
048779,000396: 6332 74347 MASK LOW11 # ENTRY FROM DISPATCHER. SAVE THE ERASABLE
048780,000397: 6333 56116 XCH ADDRWD # ADDRESS AND JUMP ON THE STORE CODE NO.
048781,000398: 6334 77656 MASK B12T14
048782,000399: 6335 00006 EXTEND
048783,000400: 6336 74764 MP BIT5 # EACH TRANSFER VECTOR ENTRY IS TWO WORDS.
048784,000401: 6337 50000 ITR0 INDEX A
048785,000402: 6340 16341 TCF STORJUMP
048786,000403:
|
|
Page 1121 |
048788,000405: # STORE CODE JUMP TABLE. CALLS THE APPROPRIATE STORING ROUTINE AND EXITS TO DANZIG OR TO ADDRESS WITH
048789,000406: # A SUPPLIED OPERATION CODE.
048790,000407:
048791,000408: # STORE STORE,1 AND STORE,2 RETURN TO DANZIG, THUS RESETTING THE EBANK TO ITS STATE AT INTPRET.
048792,000409:
048793,000410: 6341 06371 STORJUMP TC STORE # STORE.
048794,000411: 6342 16030 TCF DANZIG # PICK UP NEW OP CODE(S).
048795,000412: 6343 06363 TC STORE,1
048796,000413: 6344 16030 TCF DANZIG
048797,000414: 6345 06366 TC STORE,2
048798,000415: 6346 16030 TCF DANZIG
048799,000416:
048800,000417: 6347 06371 TC STORE # STODL.
048801,000418: 6350 16427 TCF DODLOAD
048802,000419:
048803,000420: 6351 06371 TC STORE # STODL WITH INDEXED LOAD ADDRESS.
048804,000421: 6352 16113 TCF DODLOAD*
048805,000422:
048806,000423: 6353 06371 TC STORE # STOVL.
048807,000424: 6354 16432 TCF DOVLOAD
048808,000425:
048809,000426: 6355 06371 TC STORE # STOVL WITH INDEXED LOAD ADDRESS.
048810,000427: 6356 16435 TCF DOVLOAD*
048811,000428: 6357 06371 TC STORE # STOTC.
048812,000429: 6360 35002 CAF CALLCODE
048813,000430: 6361 54020 TS CYR
048814,000431: 6362 16221 TCF 15BITADR # GET A 15 BIT ADDRESS.
048815,000432:
|
|
Page 1122 |
048817,000434: # STORE CODE ADDRESS PROCESSOR.
048818,000435:
048819,000436: 6363 50120 STORE,1 INDEX FIXLOC
048820,000437: 6364 40046 CS X1
048821,000438: 6365 16370 TCF PRESTORE
048822,000439:
048823,000440: 6366 50120 STORE,2 INDEX FIXLOC
048824,000441: 6367 40047 CS X2
048825,000442: 6370 26116 PRESTORE ADS ADDRWD # RESULTANT ADDRESS IS IN ERASABLE.
048826,000443:
048827,000444: 6371 40116 STORE CS ADDRWD
048828,000445: 6372 65006 AD DEC45
048829,000446: 6373 10000 CCS A # DOES THE ADDRESS POINT TO THE WORK AREA?
048830,000447: 6374 30120 CA FIXLOC # YES.
048831,000448: 6375 16402 TCF AHEAD5
048832,000449: 6376 35024 CA OCT1400 # NO. SET EBANK & MAKE UP SUBADDRESS.
048833,000450: 6377 56116 XCH ADDRWD
048834,000451: 6400 54003 TS EBANK
048835,000452: 6401 74350 MASK LOW8
048836,000453: 6402 26116 AHEAD5 ADS ADDRWD
048837,000454:
|
|
Page 1123 |
048839,000456: # STORING ROUTINES. STORE DP, TP, OR VECTOR AS INDICATED BY MODE.
048840,000457:
048841,000458: 6403 00006 STARTSTO EXTEND # MPAC,+1 MUST BE STORED IN ANY EVENT.
048842,000459: # ITRACE (5) REFERS TO "STARTSTO".
048843,000460: 6404 30155 DCA MPAC
048844,000461: 6405 50116 INDEX ADDRWD
048845,000462: 6406 52001 DXCH 0
048846,000463:
048847,000464: 6407 10163 CCS MODE
048848,000465: 6410 16423 TCF TSTORE
048849,000466: 6411 00002 TC Q
048850,000467:
048851,000468: 6412 00006 VSTORE EXTEND
048852,000469: 6413 30160 DCA MPAC +3
048853,000470: 6414 50116 INDEX ADDRWD
048854,000471: 6415 52003 DXCH 2
048855,000472:
048856,000473: 6416 00006 EXTEND
048857,000474: 6417 30162 DCA MPAC +5
048858,000475: 6420 50116 INDEX ADDRWD
048859,000476: 6421 52005 DXCH 4
048860,000477: 6422 00002 TC Q
048861,000478:
048862,000479: 6423 30156 TSTORE CA MPAC +2
048863,000480: 6424 50116 INDEX ADDRWD
048864,000481: 6425 54002 TS 2
048865,000482: 6426 00002 TC Q
048866,000483:
|
|
Page 1124 |
048868,000485: # ROUTINES TO BEGIN PROCESSING OF THE SECOND ADDRESS ASSOCIATED WITH ALL STORE-TYPE CODES EXCEPT STORE
048869,000486: # ITSELF.
048870,000487:
048871,000488: 6427 37665 DODLOAD CAF DLOADCOD
048872,000489: 6430 54020 TS CYR
048873,000490: 6431 16055 TCF DIRADRES # GO GET A DIRECT ADDRESS.
048874,000491:
048875,000492: 6432 34752 DOVLOAD CAF VLOADCOD
048876,000493: 6433 54020 TS CYR
048877,000494: 6434 16055 TCF DIRADRES
048878,000495:
048879,000496: 6435 36056 DOVLOAD* CAF VLOAD*
048880,000497: 6436 16114 TCF DODLOAD* +1 # PROLOGUE TO INDEX ROUTINE.
048881,000498:
|
|
Page 1125 |
048883,000500: # THE FOLLOWING LOAD INSTRUCTIONS ARE PROVIDED FOR LOADING THE MULTI-PURPOSE ACCUMULATOR MPAC.
048884,000501:
048885,000502: 6437 50116 TLOAD INDEX ADDRWD
048886,000503: 6440 30002 CA 2 # LOAD A TRIPLE PRECISION ARGUMENT INTO
048887,000504: 6441 54156 TS MPAC +2 # THE FIRST THREE MPAC REGISTERS, WITH THE
048888,000505: 6442 00006 EXTEND # CONTENTS OF THE OTHER FOUR IRRELEVANT.
048889,000506: 6443 50116 INDEX ADDRWD
048890,000507: 6444 30001 DCA 0
048891,000508: 6445 52155 DXCH MPAC
048892,000509: 6446 34770 TMODE CAF ONE
048893,000510: 6447 16027 TCF NEWMODE # DECLARE TRIPLE PRECISION MODE.
048894,000511:
048895,000512: 6450 22007 SLOAD ZL # LOAD A SINGLE PRECISION NUMBER INTO
048896,000513: 6451 50116 INDEX ADDRWD # MPAC, SETTING MPAC+1,2 TO ZERO. THE
048897,000514: 6452 30000 CA 0 # CONTENTS OF THE REMAINING MPAC REGISTERS
048898,000515: 6453 16024 TCF SLOAD2 # ARE IRRELEVANT.
048899,000516:
048900,000517: 6454 00006 VLOAD EXTEND # LOAD A DOUBLE PRECISION VECTOR INTO
048901,000518: 6455 50116 INDEX ADDRWD # MPAC,+1, MPAC+3,4, AND MPAC+5,6. THE
048902,000519: 6456 30001 DCA 0 # CONTENTS OF MPAC +2 ARE IRRELEVANT.
048903,000520: 6457 52155 DXCH MPAC
048904,000521:
048905,000522: 6460 00006 ENDVLOAD EXTEND # PDVL COMES HERE TO FINISH UP FOR DP, TP.
048906,000523: 6461 50116 INDEX ADDRWD
048907,000524: 6462 30003 DCA 2
048908,000525: 6463 52160 DXCH MPAC +3
048909,000526:
048910,000527: 6464 00006 +4 EXTEND # TPDVL FINISHES HERE.
048911,000528: 6465 50116 INDEX ADDRWD
048912,000529: 6466 30005 DCA 4
048913,000530: 6467 52162 DXCH MPAC +5
048914,000531:
048915,000532: 6470 44770 VMODE CS ONE # DECLARE VECTOR MODE.
048916,000533: 6471 16027 TCF NEWMODE
048917,000534:
|
|
Page 1126 |
048919,000536: # THE FOLLOWING INSTRUCTIONS ARE PROVIDED FOR STORING OPERANDS IN THE PUSHDOWN LIST:
048920,000537: # 1. PUSH PUSHDOWN AND NO LOAD.
048921,000538: # 2. PDDL PUSHDOWN AND DOUBLE PRECISION LOAD.
048922,000539: # 3. PDVL PUSHDOWN AND VECTOR LOAD.
048923,000540:
048924,000541: 6472 00006 PDDL EXTEND
048925,000542: 6473 50116 INDEX ADDRWD # LOAD MPAC,+1, PUSHING THE FORMER
048926,000543: 6474 30001 DCA 0 # CONTENTS DOWN.
048927,000544: 6475 52155 DXCH MPAC
048928,000545: 6476 50166 INDEX PUSHLOC
048929,000546: 6477 52001 DXCH 0
048930,000547:
048931,000548: 6500 50163 INDEX MODE # ADVANCE THE PUSHDOWN POINTER APPRO-
048932,000549: 6501 36213 CAF NO.WDS # PRIATELY.
048933,000550: 6502 26166 ADS PUSHLOC
048934,000551:
048935,000552: 6503 10163 CCS MODE
048936,000553: 6504 16521 TCF ENDTPUSH
048937,000554: 6505 16517 TCF ENDDPUSH
048938,000555:
048939,000556: 6506 54163 TS MODE # NOW DP.
048940,000557: 6507 54156 ENDVPUSH TS MPAC +2
048941,000558: 6510 52160 DXCH MPAC +3 # PUSH DOWN THE REST OF THE VECTOR HERE.
048942,000559: 6511 50166 INDEX PUSHLOC
048943,000560: 6512 51775 DXCH 0 -4
048944,000561:
048945,000562: 6513 52162 DXCH MPAC +5
048946,000563: 6514 50166 INDEX PUSHLOC
048947,000564: 6515 51777 DXCH 0 -2
048948,000565:
048949,000566: 6516 16030 TCF DANZIG
048950,000567:
048951,000568: 6517 54156 ENDDPUSH TS MPAC +2 # SET MPAC +2 TO ZERO AND EXIT ON DP.
048952,000569: 6520 16030 TCF DANZIG
048953,000570:
048954,000571: 6521 54163 ENDTPUSH TS MODE
048955,000572: 6522 56156 XCH MPAC +2 # ON TRIPLE, SET MPAC +2 TO ZERO, PUSHING
048956,000573: 6523 50166 +2 INDEX PUSHLOC # DOWN THE OLD CONTENTS
048957,000574: 6524 53777 TS 0 -1
048958,000575: 6525 16030 TCF DANZIG
048959,000576:
|
|
Page 1127 |
048961,000578: # PDVL - PUSHDOWN AND VECTOR LOAD.
048962,000579:
048963,000580: 6526 00006 PDVL EXTEND # RELOAD MPAC AND PUSH DOWN ITS CONTENTS.
048964,000581: 6527 50116 INDEX ADDRWD
048965,000582: 6530 30001 DCA 0
048966,000583: 6531 52155 DXCH MPAC
048967,000584: 6532 50166 INDEX PUSHLOC
048968,000585: 6533 52001 DXCH 0
048969,000586:
048970,000587: 6534 50163 INDEX MODE # ADVANCE THE PUSHDOWN POINTER.
048971,000588: 6535 36213 CAF NO.WDS
048972,000589: 6536 26166 ADS PUSHLOC
048973,000590:
048974,000591: 6537 10163 CCS MODE # TEST PAST MODE.
048975,000592: 6540 16557 TCF TPDVL
048976,000593: 6541 16460 TCF ENDVLOAD # JUST LOAD LAST FOUR REGISTERS ON DP.
048977,000594:
048978,000595: 6542 00006 VPDVL EXTEND # PUSHDOWN AND RE-LOAD LAST TWO COMPONENTS
048979,000596: 6543 50116 INDEX ADDRWD
048980,000597: 6544 30003 DCA 2
048981,000598: 6545 52160 DXCH MPAC +3
048982,000599: 6546 50166 INDEX PUSHLOC
048983,000600: 6547 51775 DXCH 0 -4
048984,000601:
048985,000602: 6550 00006 EXTEND
048986,000603: 6551 50116 INDEX ADDRWD
048987,000604: 6552 30005 DCA 4
048988,000605: 6553 52162 DXCH MPAC +5
048989,000606: 6554 50166 INDEX PUSHLOC
048990,000607: 6555 51777 DXCH 0 -2
048991,000608:
048992,000609: 6556 16030 TCF DANZIG
048993,000610:
048994,000611: 6557 00006 TPDVL EXTEND # ON TP, WE MUST LOAD THE Y COMPONENT
048995,000612: 6560 50116 INDEX ADDRWD # BEFORE STORING MPAC +2 IN CASE THIS IS A
048996,000613: 6561 30003 DCA 2 # PUSHUP.
048997,000614: 6562 52160 DXCH MPAC +3
048998,000615:
048999,000616: 6563 30156 CA MPAC +2
049000,000617: 6564 50166 INDEX PUSHLOC # IN DP.
049001,000618: 6565 53777 TS 0 -1
049002,000619: 6566 16464 TCF ENDVLOAD +4
049003,000620:
049004,000621: # SSP (STORE SINGLE PRECISION) IS EXECUTED HERE.
049005,000622:
049006,000623: 6567 24164 SSP INCR LOC # PICK UP THE WORD FOLLOWING THE GIVEN
049007,000624: 6570 50164 INDEX LOC # ADDRESS AND STORE IT AT X.
049008,000625: 6571 30000 CA 0
049009,000626: 6572 50116 STORE1 INDEX ADDRWD # SOME INDEX AND MISCELLANEOUS OPS END
049010,000627: 6573 54000 TS 0 # HERE.
|
|
Page 1128 |
049012,000629: 6574 16030 TCF DANZIG
049013,000630:
|
|
Page 1129 |
049015,000632: # SEQUENCE CHANGING AND SUBROUTINE CALLING OPTIONS.
049016,000633:
049017,000634: # THE FOLLOWING OPERATIONS ARE AVAILABLE FOR SEQUENCING CHANGING, BRANCHING, AND CALLING SUBROUTINES:
049018,000635:
049019,000636: # 1. GOTO GO TO.
049020,000637: # 2. CALL CALL SUBROUTINE SETTING QPRET.
049021,000638: # 3. CGOTO COMPUTED GO TO.
049022,000639: # 4. CCALL COMPUTED CALL.
049023,000640: # 7. BPL BRANCH IF MPAC POSITIVE OR ZERO.
049024,000641: # 8. BZE BRANCH IF MPAC ZERO.
049025,000642: # 9. BMN BRANCH IF MPAC NEGATIVE NON-ZERO.
049026,000643:
049027,000644: 6575 24164 CCALL INCR LOC # MAINTAIN LOC FOR QPRET COMPUTATION.
049028,000645: 6576 50164 INDEX LOC
049029,000646: 6577 30000 CAF 0 # GET BASE ADDRESS OF CADR LIST.
049030,000647: 6600 50116 INDEX ADDRWD
049031,000648: 6601 60000 AD 0 # ADD INCREMENT.
049032,000649: 6602 54004 TS FBANK # SELECT DESIRED CADR.
049033,000650: 6603 75041 MASK LOW10
049034,000651: 6604 50000 INDEX A
049035,000652: 6605 32000 CAF 10000
049036,000653: 6606 54117 TS POLISH
049037,000654:
049038,000655: 6607 30165 CALL CA BANKSET # FOR ANY OF THE CALL OPTIONS, MAKE UP THE
049039,000656: 6610 74344 MASK BANKMASK # ADDRESS OF THE NEXT OP-CODE PAIR/STORE
049040,000657: 6611 64344 AD BANKMASK # CODE AND LEAVE IT IN QPRET. NOTE THAT
049041,000658: 6612 60164 AD LOC # BANKMASK = -(2000 - 1).
049042,000659: 6613 50120 INDEX FIXLOC
049043,000660: 6614 54052 TS QPRET
049044,000661:
049045,000662: 6615 30117 GOTO CA POLISH # BASIC BRANCHING SEQUENCE.
049046,000663: 6616 77675 +1 MASK HIGH4
049047,000664: 6617 00006 EXTEND
049048,000665: 6620 16631 BZF GOTOERS # SEE IF ADDRESS POINTS TO FIXED OR ERAS.
049049,000666: 6621 30165 +4 CA BANKSET # SET EBANK PART OF BBANK. NEXT, SET UP
049050,000667: 6622 54006 TS BBANK # FBANK. THE COMBINATION IS PICKED UP &
049051,000668: 6623 30117 CA POLISH # PUT INTO BANKSET AT INTPRET +2.
049052,000669: 6624 54004 TS FBANK
049053,000670: 6625 75041 MASK LOW10
049054,000671: 6626 64756 AD 2K
049055,000672: 6627 54164 TS LOC
049056,000673: 6630 16011 TCF INTPRET +3
049057,000674:
049058,000675: 6631 E3,1400 EBANK= 1400 # SO YUL DOESN'T CUSS THE "CA 1400" BELOW.
049059,000676:
049060,000677: 6631 30117 GOTOERS CA POLISH # THE GIVEN ADDRESS IS IN ERASABLE - SEE
049061,000678: 6632 66220 AD -ENDVAC # IF RELATIVE TO THE WORK AREA.
049062,000679: 6633 10000 CCS A
049063,000680: 6634 30117 CA POLISH # GENERAL ERASABLE.
049064,000681: 6635 16644 TCF GOTOGE
|
|
Page 1130 |
049066,000683: 6636 30120 CA FIXLOC # WORK AREA.
049067,000684: 6637 60117 AD POLISH
049068,000685: 6640 50000 INDEX A # USE THE GIVEN ADDRESS AS THE ADDRESS OF
049069,000686: 6641 30000 CA 0 # THE BRANCH ADDRESS.
049070,000687: 6642 54117 TS POLISH
049071,000688: 6643 16616 TCF GOTO +1 # ALLOWS ARBITRARY INDIRECTNESS LEVELS.
049072,000689:
049073,000690: 6644 54003 GOTOGE TS EBANK
049074,000691: 6645 74350 MASK LOW8
049075,000692: 6646 50000 INDEX A # USE THE GIVEN ADDRESS AS THE ADDRESS OF
049076,000693: 6647 31400 CA 1400 # THE BRANCH ADDRESS.
049077,000694: 6650 54117 TS POLISH
049078,000695: 6651 16616 TCF GOTO +1
049079,000696:
049080,000697: 6652 50164 CGOTO INDEX LOC # COMPUTED GO TO. PICK UP ADDRESS OF CADR
049081,000698: 6653 30001 CA 1 # LIST.
049082,000699: 6654 50116 INDEX ADDRWD # ADD MODIFIER.
049083,000700: 6655 60000 AD 0
049084,000701: 6656 54004 TS FBANK # SELECT GOTO ADDRESS.
049085,000702: 6657 75041 MASK LOW10
049086,000703: 6660 50000 INDEX A
049087,000704: 6661 32000 CA 10000
049088,000705: 6662 54117 TS POLISH
049089,000706: 6663 16616 TCF GOTO +1 # WITH ADDRESS IN A.
049090,000707:
049091,000708: 6664 30165 SWBRANCH CA BANKSET # SWITCH INSTRUCTIONS WHICH ELECT TO
049092,000709: 6665 54004 TS FBANK # BRANCH COME HERE TO DO SO.
049093,000710: 6666 50164 INDEX LOC
049094,000711: 6667 30001 CA 1
049095,000712: 6670 54117 TS POLISH
049096,000713: 6671 16616 TCF GOTO +1
049097,000714:
|
|
Page 1131 |
049099,000716: # TRIPLE PRECISION BRANCHING ROUTINE. IF CALLING TC IS AT L, RETURN IS AS FOLLOWS:
049100,000717:
049101,000718: # L+1 IF MPAC IS GREATER THAN ZERO.
049102,000719: # L+2 IF MPAC IS EQUAL TO +0 OR -0.
049103,000720: # L+3 IF MPAC IS LESS THAN ZERO.
049104,000721:
049105,000722: 6672 10154 BRANCH CCS MPAC
049106,000723: 6673 00002 TC Q
049107,000724: 6674 16676 TCF +2 # ON ZERO.
049108,000725: 6675 16710 TCF NEG
049109,000726:
049110,000727: 6676 10155 CCS MPAC +1
049111,000728: 6677 00002 TC Q
049112,000729: 6700 16702 TCF +2
049113,000730: 6701 16710 TCF NEG
049114,000731:
049115,000732: 6702 10156 CCS MPAC +2
049116,000733: 6703 00002 TC Q
049117,000734: 6704 16706 TCF +2
049118,000735: 6705 16710 TCF NEG
049119,000736:
049120,000737: 6706 50002 Q+1 INDEX Q
049121,000738: 6707 00001 TC 1
049122,000739:
049123,000740: 6710 50002 NEG INDEX Q # IF FIRST NON-ZERO REGISTER WAS NEGATIVE.
049124,000741: 6711 00002 TCQ TC 2
049125,000742:
049126,000743: 6712 6710 Q+2 = NEG
049127,000744:
049128,000745: # ITRACE (3) REFERS TO "EXIT".
049129,000746:
049130,000747: 6712 30165 EXIT CA BANKSET # RESTORE USER'S BANK SETTING, AND LEAVE
049131,000748: 6713 54006 TS BBANK # INTERPRETIVE MODE.
049132,000749: 6714 50164 INDEX LOC
049133,000750: 6715 00001 TC 1
049134,000751:
|
|
Page 1132 |
049136,000753: # SECTION 3 - ADD/SUBTRACT PACKAGE.
049137,000754:
049138,000755: # THE FOLLOWING OPERATIONS ARE PROVIDED FOR ADDING TO AND SUBTRACTING FROM THE MULTI-PURPOSE ACCUMULATOR
049139,000756: # MPAC:
049140,000757:
049141,000758: # 1. DAD DOUBLE PRECISION ADD.
049142,000759: # 2. DSU DOUBLE PRECISION SUBTRACT.
049143,000760: # 3. BDSU DOUBLE PRECISION SUBTRACT FROM.
049144,000761:
049145,000762: # 4. TAD TRIPLE PRECISION ADD.
049146,000763:
049147,000764: # 5. VAD VECTOR ADD.
049148,000765: # 6. VSU VECTOR SUBTRACT.
049149,000766: # 7. BVSU VECTOR SUBTRACT FROM.
049150,000767:
049151,000768: # THE INTERPRETIVE OVERFLOW INDICATOR OVFIND IS SET NON-ZERO IF OVERFLOW OCCURS IN ANY OF THE ABOVE.
049152,000769:
049153,000770: 6716 34752 VSU CAF BIT15 # CHANGES 0 TO DCS.
049154,000771: 6717 16721 TCF +2
049155,000772:
049156,000773: 6720 34346 VAD CAF PRIO30 # CHANGES 0 TO DCA.
049157,000774: 6721 26116 ADS ADDRWD
049158,000775: 6722 00006 EXTEND
049159,000776: 6723 50116 INDEX ADDRWD
049160,000777: 6724 00003 READ HISCALAR # DCA 2 OR DCS 2
049161,000778: 6725 20160 DAS MPAC +3
049162,000779: 6726 00006 EXTEND # CHECK OVERFLOW.
049163,000780: 6727 16731 BZF +2
049164,000781: 6730 06763 TC OVERFLWY
049165,000782:
049166,000783: 6731 00006 EXTEND
049167,000784: 6732 50116 INDEX ADDRWD
049168,000785: 6733 00005 READ CHAN5 # DCA 4 OR DCS 4
049169,000786: 6734 20162 DAS MPAC +5
049170,000787: 6735 00006 EXTEND
049171,000788: 6736 16740 BZF +2
049172,000789: 6737 06760 TC OVERFLWZ
049173,000790:
049174,000791: 6740 00006 EXTEND
049175,000792: 6741 50116 INDEX ADDRWD
049176,000793: 6742 00001 READ LCHAN # DCA 0 OR DCS 0
049177,000794: 6743 16747 TCF ENDVXV
049178,000795:
049179,000796: 6744 00006 DAD EXTEND
049180,000797: 6745 50116 INDEX ADDRWD
049181,000798: 6746 30001 DCA 0
049182,000799: 6747 20155 ENDVXV DAS MPAC # VXV FINISHES HERE.
049183,000800: 6750 00006 EXTEND
049184,000801: 6751 16030 BZF DANZIG
049185,000802:
|
|
Page 1133 |
049187,000804: 6752 06766 SETOVF TC OVERFLOW
049188,000805: 6753 16030 TCF DANZIG
049189,000806:
|
|
Page 1134 |
049191,000808: 6754 00006 DSU EXTEND
049192,000809: 6755 50116 INDEX ADDRWD
049193,000810: 6756 40001 DCS 0
049194,000811: 6757 16747 TCF ENDVXV
049195,000812:
049196,000813: 6760 54001 OVERFLWZ TS L # ENTRY FOR THIRD COMPONENT.
049197,000814: 6761 34773 CAF FIVE
049198,000815: 6762 16765 TCF +3
049199,000816:
049200,000817: 6763 54001 OVERFLWY TS L # ENTRY FOR SECOND COMPONENT.
049201,000818: 6764 36214 CAF THREE
049202,000819: 6765 56001 XCH L
049203,000820:
049204,000821: 6766 50000 OVERFLOW INDEX A # ENTRY FOR 1ST COMP OR DP (L=0).
049205,000822: 6767 44751 CS LIMITS # PICK UP POSMAX OR NEGMAX.
049206,000823: 6770 54130 TS BUF
049207,000824: 6771 00006 EXTEND
049208,000825: 6772 24000 AUG A
049209,000826: 6773 50001 INDEX L
049210,000827: 6774 26155 ADS MPAC +1
049211,000828: 6775 54007 TS 7
049212,000829: 6776 34772 CAF ZERO
049213,000830: 6777 60130 AD BUF
049214,000831: 7000 50001 INDEX L
049215,000832: 7001 26154 ADS MPAC
049216,000833: 7002 54007 TS 7
049217,000834: 7003 00002 TC Q # NO OVERFLOW EXIT.
049218,000835: 7004 17121 TCF SETOVF2 # SET OVFIND AND EXIT.
049219,000836:
049220,000837: 7005 00006 BVSU EXTEND
049221,000838: 7006 50116 INDEX ADDRWD
049222,000839: 7007 30003 DCA 2
049223,000840: 7010 52160 DXCH MPAC +3
049224,000841: 7011 00006 EXTEND
049225,000842: 7012 40001 DCOM
049226,000843: 7013 20160 DAS MPAC +3
049227,000844: 7014 00006 EXTEND
049228,000845: 7015 17017 BZF +2
049229,000846: 7016 06763 TC OVERFLWY
049230,000847:
049231,000848: 7017 00006 EXTEND
049232,000849: 7020 50116 INDEX ADDRWD
049233,000850: 7021 30005 DCA 4
049234,000851: 7022 52162 DXCH MPAC +5
049235,000852: 7023 00006 EXTEND
049236,000853: 7024 40001 DCOM
049237,000854: 7025 20162 DAS MPAC +5
049238,000855: 7026 00006 EXTEND
049239,000856: 7027 17031 BZF +2
049240,000857: 7030 06760 TC OVERFLWZ
049241,000858:
|
|
Page 1135 |
049243,000860: 7031 00006 BDSU EXTEND
049244,000861: 7032 50116 INDEX ADDRWD
049245,000862: 7033 30001 DCA 0
049246,000863: 7034 52155 DXCH MPAC
049247,000864: 7035 00006 EXTEND
049248,000865: 7036 40001 DCOM
049249,000866: 7037 16747 TCF ENDVXV
049250,000867:
|
|
Page 1136 |
049252,000869: # TRIPLE PRECISION ADD ROUTINE.
049253,000870:
049254,000871: 7040 00006 TAD EXTEND
049255,000872: 7041 50116 INDEX ADDRWD
049256,000873: 7042 30002 DCA 1 # ADD MINOR PARTS FIRST.
049257,000874: 7043 20156 DAS MPAC +1
049258,000875: 7044 50116 INDEX ADDRWD
049259,000876: 7045 60000 AD 0
049260,000877: 7046 60154 AD MPAC
049261,000878: 7047 54154 TS MPAC
049262,000879: 7050 16030 TCF DANZIG
049263,000880:
049264,000881: 7051 16752 TCF SETOVF # SET OVFIND IF SUCH OCCURS.
049265,000882:
|
|
Page 1137 |
049267,000884: # ARITHMETIC SUBROUTINES REQUIRED IN FIXED-FIXED.
049268,000885:
049269,000886: # 1. DMPSUB DOUBLE PRECISION MULTIPLY. MULTIPLY THE CONTENTS OF MPAC,+1 BY THE DP WORD WHOSE ADDRESS
049270,000887: # IS IN ADDRWD AND LEAVE A TRIPLE PRECISION RESULT IN MPAC.
049271,000888: # 2. ROUNDSUB ROUND THE TRIPLE PRECISON CONTENTS OF MPAC TO DOUBLE PRECISION.
049272,000889: # 3. DOTSUB TAKE THE DOT PRODUCT OF THE VECTOR IN MPAC AND THE VECTOR WHOSE ADDRESS IS IN ADDRWD
049273,000890: # AND LEAVE THE TRIPLE PRECISION RESULT IN MPAC.
049274,000891: # 4. POLY USING THE CONTENTS OF MPAC AS A DP ARGUMENT, EVALUATE THE POLYNOMIAL WHOSE DEGREE AND
049275,000892: # COEFFICIENTS IMMEDIATELY FOLLOW THE TC POLY INSTRUCTION (SEE ROUTINE FOR DETAILS.)
049276,000893:
049277,000894: 7052 50002 DMP INDEX Q # BASIC SUBROUTINE FOR USE BY PINBALL, ETC
049278,000895: 7053 30000 CAF 0 # ADRES OF ARGUMENT FOLLOWS TC DMP .
049279,000896: 7054 24002 INCR Q
049280,000897: 7055 54116 -1 TS ADDRWD # (PROLOGUE FOR SETTING ADDRWD.)
049281,000898:
049282,000899: 7056 50116 DMPSUB INDEX ADDRWD # GET MINOR PART OF OPERAND AT C(ADDRWD).
049283,000900: 7057 30001 CA 1
049284,000901: 7060 54156 TS MPAC +2 # THIS WORKS FOR SQUARING MPAC AS WELL.
049285,000902: 7061 34772 CAF ZERO # SET MPAC +1 TO ZERO SO WE CAN ACCUMULATE
049286,000903: 7062 56155 XCH MPAC +1 # THE PARTIAL PRODUCTS WITH DAS
049287,000904: 7063 54135 TS MPTEMP # INSTRUCTIONS.
049288,000905: 7064 00006 EXTEND
049289,000906: 7065 70156 MP MPAC +2 # MINOR OF MPAC X MINOR OF C(ADDRWD).
049290,000907:
049291,000908: 7066 56156 XCH MPAC +2 # DISCARD MINOR PART OF ABOVE RESULT AND
049292,000909: 7067 00006 EXTEND # FORM MAJOR OF MPAC X MINOR OF C(ADDRWD).
049293,000910: 7070 70154 MP MPAC
049294,000911: 7071 20156 DAS MPAC +1 # GUARANTEED NO OVERFLOW.
049295,000912:
049296,000913: 7072 50116 INDEX ADDRWD # GET MAJOR PART OF ARGUMENT AT C(ADDRWD).
049297,000914: 7073 30000 CA 0
049298,000915: 7074 56135 XCH MPTEMP # SAVE AND BRING OUT MINOR OF MPAC.
049299,000916: 7075 00006 DMPSUB2 EXTEND
049300,000917: 7076 70135 MP MPTEMP # MAJOR OF C(ADDRWD) X MINOR OF MPAC.
049301,000918: 7077 20156 DAS MPAC +1 # ACCUMULATE, SETTING A TO NET OVERFLOW.
049302,000919:
049303,000920: 7100 56154 XCH MPAC # SETTING MPAC TO 0 OR +-1.
049304,000921: 7101 00006 EXTEND
049305,000922: 7102 70135 MP MPTEMP # MAJOR OF MPAC X MAJOR OF C(ADDRWD).
049306,000923: 7103 20155 DAS MPAC # GUARANTEED NO OVERFLOW.
049307,000924: 7104 00002 TC Q # 49 MCT = .573 MS. INCLUDING RETURN.
049308,000925:
|
|
Page 1138 |
049310,000927: # ROUND MPAC TO DOUBLE PRECISION, SETTING OVFIND ON THE RARE EVENT OF OVERFLOW.
049311,000928:
049312,000929: 7105 34772 ROUNDSUB CAF ZERO # SET MPAC +2 = 0 FOR SCALARS AND CHANGE
049313,000930: 7106 54163 +1 TS MODE # MODE TO DP.
049314,000931:
049315,000932: 7107 56156 VROUND XCH MPAC +2 # BUT WE NEEDNT TAKE THE TIME FOR VECTORS.
049316,000933: 7110 60000 DOUBLE
049317,000934: 7111 54001 TS L
049318,000935: 7112 00002 TC Q
049319,000936:
049320,000937: 7113 60155 AD MPAC +1 # ADD ROUNDING BIT IF MPAC +2 WAS GREATER
049321,000938: 7114 54155 TS MPAC +1 # THAN .5 IN MAGNITUDE.
049322,000939: 7115 00002 TC Q
049323,000940:
049324,000941: 7116 60154 AD MPAC # PROPAGATE INTERFLOW.
049325,000942: 7117 54154 TS MPAC
049326,000943: 7120 00002 TC Q
049327,000944:
049328,000945: 7121 54121 SETOVF2 TS OVFIND # (RARE).
049329,000946: 7122 00002 TC Q
049330,000947:
|
|
Page 1139 |
049332,000949: # THE DOT PRODUCT SUBROUTINE USUALLY FORMS THE DOT PRODUCT OF THE VECTOR IN MPAC WITH A STANDARD SIX
049333,000950: # REGISTER VECTOR WHOSE ADDRESS IS IN ADDRWD. IN THIS CASE C(DOTINC) ARE SET TO 2. VXM, HOWEVER, SETS C(DOTINC) TO
049334,000951: # 6 SO THAT DOTSUB DOTS MPAC WITH A COLUMN VECTOR OF THE MATRIX IN QUESTION IN THIS CASE.
049335,000952:
049336,000953: 7123 34767 PREDOT CAF TWO # PROLOGUE TO SET DOTINC TO 2.
049337,000954: 7124 54136 TS DOTINC
049338,000955:
049339,000956: 7125 00006 DOTSUB EXTEND
049340,000957: 7126 22137 QXCH DOTRET # SAVE RETURN.
049341,000958: 7127 07056 TC DMPSUB # DOT X COMPONENTS.
049342,000959: 7130 52160 DXCH MPAC +3 # POSITION Y COMPONENT OF MPAC FOR
049343,000960: 7131 52155 DXCH MPAC # MULTIPLICATION WHILE SAVING RESULT IN
049344,000961: 7132 52131 DXCH BUF # THREE WORD BUFFER, BUF.
049345,000962: 7133 30156 CA MPAC +2
049346,000963: 7134 54132 TS BUF +2
049347,000964:
049348,000965: 7135 30136 CA DOTINC # ADVANCE ADDRWD TO Y COMPONENT OF
049349,000966: 7136 26116 ADS ADDRWD # OTHER ARGUMENT.
049350,000967: 7137 07056 TC DMPSUB
049351,000968: 7140 52156 DXCH MPAC +1 # ACCUMULATE PARTIAL PRODUCTS.
049352,000969: 7141 20132 DAS BUF +1
049353,000970: 7142 60154 AD MPAC
049354,000971: 7143 60130 AD BUF
049355,000972: 7144 54130 TS BUF
049356,000973: 7145 17147 TCF +2
049357,000974: 7146 54121 TS OVFIND # IF OVERFLOW OCCURS.
049358,000975:
049359,000976: 7147 52162 DXCH MPAC +5 # MULTIPLY Z COMPONENTS.
049360,000977: 7150 52155 DXCH MPAC
049361,000978: 7151 30136 CA DOTINC
049362,000979: 7152 26116 ADS ADDRWD
049363,000980: 7153 07056 TC DMPSUB
049364,000981: 7154 52132 ENDDOT DXCH BUF +1 # LEAVE FINAL ACCUMULATION IN MPAC.
049365,000982: 7155 20156 DAS MPAC +1
049366,000983: 7156 60154 AD MPAC
049367,000984: 7157 60130 AD BUF
049368,000985: 7160 54154 TS MPAC
049369,000986: 7161 00137 TC DOTRET
049370,000987:
049371,000988: 7162 06766 TC OVERFLOW # ON OVERFLOW HERE.
049372,000989: 7163 00137 TC DOTRET
049373,000990:
|
|
Page 1140 |
049375,000992: # DOUBLE PRECISION POLYNOMIAL EVALUATOR
049376,000993:
049377,000994: # N N-1
049378,000995: # THIS ROUTINE EVALUATES A X + A X + ... + A X + A LEAVING THE DP RESULT IN MPAC ON EXIT.
049379,000996: # N N-1 1 0
049380,000997:
049381,000998: # THE ROUTINE HAS TWO ENTRIES
049382,000999:
049383,001000: # 1. ENTRY THRU POWRSERS. THE COEFFICIENTS MAY BE EITHER IN FIXED OR ERASABLE, THE CALL IS BY
049384,001001: # TC POWRSERS, AND THE RETURN IS TO LOC(TC POWRSERS)+1. THE ENTERING DATA MUST BE AS FOLLOWS
049385,001002:
049386,001003: # A SP LOC-3 ADDRESS FOR REFERENCING COEF TABLE
049387,001004: # L SP N-1 N IS THE DEGREE OF THE POWER SERIES
049388,001005: # MPAC DP X ARGUMENT
049389,001006:
049390,001007: # LOC-2N DP A(0)
049391,001008: # ...
049392,001009: # LOC DP A(N)
049393,001010:
049394,001011: # 2. ENTRY THRU POLY. THE CALL TO POLY AND THE ENTERING DATA MUST BE AS FOLLOWS
049395,001012:
049396,001013: # MPAC DP X ARGUMENT
049397,001014:
049398,001015: # LOC TC POLY
049399,001016: # LOC+1 SP N-1
049400,001017: # LOC+2 DP A(0)
049401,001018: # ...
049402,001019: # LOC+2N+2 DP A(N) RETURN IS TO LOC+2N+4
049403,001020:
049404,001021: 7164 00006 POWRSERS EXTEND
049405,001022: 7165 22141 QXCH POLYRET # RETURN ADDRESS
049406,001023: 7166 54117 TS POLISH # POWER SERIES ADDRESS
049407,001024: 7167 22140 LXCH POLYCNT # N-1 TO COUNTER
049408,001025: 7170 17201 TCF POLYCOM # SKIP SET UP BY POLY
049409,001026:
049410,001027: 7171 50002 POLY INDEX Q
049411,001028: 7172 30000 CAF 0
049412,001029: 7173 54140 TS POLYCNT # N-1 TO COUNTER
049413,001030: 7174 60000 DOUBLE
049414,001031: 7175 60002 AD Q
049415,001032: 7176 54117 TS POLISH # L(A(N))-3 TO POLISH
049416,001033: 7177 64773 AD FIVE
049417,001034: 7200 54141 TS POLYRET # STORE RETURN ADDRESS
049418,001035:
049419,001036: 7201 36242 POLYCOM CAF LVBUF # INCOMING X WILL BE MOVED TO VBUF, SO
049420,001037: 7202 54116 TS ADDRWD # SET ADDRWD SO DMPSUB WILL MPY BY VBUF.
049421,001038:
049422,001039: 7203 00006 EXTEND
049423,001040: 7204 50117 INDEX POLISH
049424,001041: 7205 30004 DCA 3
|
|
Page 1141 |
049426,001043: 7206 52155 DXCH MPAC # LOAD A(N) INTO MPAC,
049427,001044: 7207 52123 DXCH VBUF # SAVING X IN VBUF
049428,001045: 7210 17214 TCF POLY2
049429,001046:
049430,001047: 7211 54140 POLYLOOP TS POLYCNT # SAVE DECREMENTED LOOP COUNTER
049431,001048: 7212 44767 CS TWO
049432,001049: 7213 26117 ADS POLISH # REGRESS COEFFICIENT POINTER
049433,001050:
049434,001051: 7214 07056 POLY2 TC DMPSUB # MULTIPLY BY X
049435,001052: 7215 00006 EXTEND
049436,001053: 7216 50117 INDEX POLISH
049437,001054: 7217 30002 DCA 1 # ADD IN NEXT COEFFICIENT
049438,001055: 7220 20155 DAS MPAC # USERS RESPONSIBILITY TO ASSURE NO OVFLOW
049439,001056:
049440,001057: 7221 10140 CCS POLYCNT
049441,001058: 7222 17211 TCF POLYLOOP
049442,001059: 7223 00141 TC POLYRET # RETURN CALLER
049443,001060:
|
|
Page 1142 |
049445,001062: # MISCELLANEOUS MULTI-PRECISION ROUTINES REQUIRED IN FIXED-FIXED BUT NOT USED BY THE INTERPRETER.
049446,001063:
049447,001064: 7224 34772 DPAGREE CAF ZERO # DOUBLE PRECISION ENTRY -
049448,001065: 7225 54156 BUFAGREE TS MPAC +2 # ZERO LOW-ORDER WORD
049449,001066: 7226 22002 TPAGREE LXCH Q # FORCE SIGN AGREEMENT AMONG THE TRIPLE
049450,001067: 7227 06672 TC BRANCH # PRECISION CONTENTS OF MPAC. RETURNING
049451,001068: 7230 17234 TCF ARG+ # WITH SIGNUM OF THE INPUT IN A.
049452,001069: 7231 17254 TCF ARGZERO
049453,001070:
049454,001071: 7232 44750 CS POSMAX # IF NEGATIVE.
049455,001072: 7233 17235 TCF +2
049456,001073:
049457,001074: 7234 34750 ARG+ CAF POSMAX
049458,001075: 7235 54002 TS Q
049459,001076: 7236 00006 EXTEND
049460,001077: 7237 24000 AUG A # FORMS +-1.0.
049461,001078: 7240 60156 AD MPAC +2
049462,001079: 7241 54156 TS MPAC +2
049463,001080: 7242 34772 CAF ZERO
049464,001081: 7243 60002 AD Q
049465,001082: 7244 60155 AD MPAC +1
049466,001083: 7245 54155 TS MPAC +1
049467,001084: 7246 34772 CAF ZERO
049468,001085: 7247 60002 AD Q # Q STILL HAS POSMAX OR NEGMAX IN IT.
049469,001086: 7250 60154 AD MPAC
049470,001087: 7251 54154 ARGZERO2 TS MPAC # ALWAYS SKIPPING UNLESS ARGZERO.
049471,001088: 7252 54155 TS MPAC +1
049472,001089: 7253 00001 TC L # RETURN VIA L.
049473,001090:
049474,001091: 7254 54156 ARGZERO TS MPAC +2 # SET ALL THREE MPAC REGISTERS TO ZERO.
049475,001092: 7255 17251 TCF ARGZERO2
049476,001093:
049477,001094: # SHORTMP MULTIPLIES THE TP CONTENTS OF MPAC BY THE SINGLE PRECISION NUMBER ARRIVING IN A.
049478,001095:
049479,001096: 7256 54135 SHORTMP TS MPTEMP
049480,001097: 7257 00006 +1 EXTEND
049481,001098: 7260 70156 MP MPAC +2
049482,001099: 7261 54156 TS MPAC +2
049483,001100: 7262 34772 SHORTMP2 CAF ZERO # SO SUBSEQUENT DAS WILL WORK.
049484,001101: 7263 56155 XCH MPAC +1
049485,001102: 7264 17075 TCF DMPSUB2
049486,001103:
|
|
Page 1143 |
049488,001105: # MISCELLANEOUS VECTOR OPERATIONS. INCLUDED HERE ARE THE FOLLOWING:
049489,001106:
049490,001107: # 1. DOT DP VECTOR DOT PRODUCT.
049491,001108: # 2. VXV DP VECTOR CROSS PRODUCT.
049492,001109: # 3. VXSC DP VECTOR TIMES SCALAR.
049493,001110: # 4. V/SC DP VECTOR DIVIDED BY SCALAR.
049494,001111: # 5. VPROJ DP VECTOR PROJECTION. ( (MPAC.X)MPAC ).
049495,001112: # 6. VXM DP VECTOR POST-MULTIPLIED BY MATRIX.
049496,001113: # 7. MXV DP VECTOR PRE-MULTIPLIED BY MATRIX.
049497,001114:
049498,001115: 7265 07123 DOT TC PREDOT # DO THE DOT PRODUCT AND EXIT, CHANGING
049499,001116: 7266 34772 DMODE CAF ZERO # THE MODE TO DP SCALAR.
049500,001117: 7267 16027 TCF NEWMODE
049501,001118:
049502,001119: 7270 34767 MXV CAF TWO # SET UP MATINC AND DOTINC FOR ROW
049503,001120: 7271 54140 TS MATINC # VECTORS.
049504,001121: 7272 17276 TCF VXM/MXV # GO TO COMMON PORTION.
049505,001122:
049506,001123: 7273 44355 VXM CS TEN # SET MATINC AND DOTINC TO REFER TO MATRIX
049507,001124: 7274 54140 TS MATINC # AS THREE COLUMN VECTORS.
049508,001125: 7275 36211 CAF SIX
049509,001126:
|
|
Page 1144 |
049511,001128: # COMMON PORTION OF MXV AND VXM.
049512,001129:
049513,001130: 7276 54136 VXM/MXV TS DOTINC
049514,001131:
049515,001132: # ITRACE (2) REFERS TO "VXM/MXV".
049516,001133:
049517,001134: 7277 07466 TC MPACVBUF # SAVE VECTOR IN MPAC FOR FURTHER USE.
049518,001135:
049519,001136: 7300 07125 TC DOTSUB # GO DOT TO GET X COMPONENT OF ANSWER.
049520,001137: 7301 00006 EXTEND
049521,001138: 7302 30123 DCA VBUF # MOVE MPAC VECTOR BACK INTO MPAC, SAVING
049522,001139: 7303 52155 DXCH MPAC # NEW X COMPONENT IN BUF2.
049523,001140: 7304 52134 DXCH BUF2
049524,001141: 7305 00006 EXTEND
049525,001142: 7306 30125 DCA VBUF +2
049526,001143: 7307 52160 DXCH MPAC +3
049527,001144: 7310 00006 EXTEND
049528,001145: 7311 30127 DCA VBUF +4
049529,001146: 7312 52162 DXCH MPAC +5
049530,001147: 7313 30140 CA MATINC # INITIALIZE ADDRWD FOR NEXT DOT PRODUCT.
049531,001148: 7314 26116 ADS ADDRWD # FORMS BASE ADDRESS OF NEXT COLUMN(ROW).
049532,001149:
049533,001150: 7315 07125 TC DOTSUB
049534,001151: 7316 52123 DXCH VBUF # MOVE GIVEN VECTOR BACK TO MPAC, SAVING Y
049535,001152: 7317 52155 DXCH MPAC # COMPONENT OF ANSWER IN VBUF +2.
049536,001153: 7320 52125 DXCH VBUF +2
049537,001154: 7321 52160 DXCH MPAC +3
049538,001155: 7322 52127 DXCH VBUF +4
049539,001156: 7323 52162 DXCH MPAC +5
049540,001157: 7324 30140 CA MATINC # FORM ADDRESS OF LAST COLUMN OR ROW.
049541,001158: 7325 26116 ADS ADDRWD
049542,001159:
049543,001160: 7326 07125 TC DOTSUB
049544,001161: 7327 52134 DXCH BUF2 # ANSWER NOW COMPLETE. PUT COMPONENTS INTO
049545,001162: 7330 52155 DXCH MPAC # PROPER MPAC REGISTERS.
049546,001163: 7331 52162 DXCH MPAC +5
049547,001164: 7332 52125 DXCH VBUF +2
049548,001165: 7333 52160 DXCH MPAC +3
049549,001166: 7334 16030 TCF DANZIG # EXIT.
049550,001167:
|
|
Page 1145 |
049552,001169: # VXSC - VECTOR TIMES SCALAR.
049553,001170:
049554,001171: 7335 10163 VXSC CCS MODE # TEST PRESENT MODE.
049555,001172: 7336 17364 TCF DVXSC # SEPARATE ROUTINE WHEN SCALAR IS IN MPAC.
049556,001173: 7337 17364 TCF DVXSC
049557,001174:
049558,001175: 7340 07056 VVXSC TC DMPSUB # COMPUTE X COMPONENT
049559,001176: 7341 07107 TC VROUND # AND ROUND IT.
049560,001177: 7342 52160 DXCH MPAC +3 # PUT Y COMPONENT INTO MPAC SAVING MPAC IN
049561,001178: 7343 52155 DXCH MPAC # MPAC +3.
049562,001179: 7344 52160 DXCH MPAC +3
049563,001180:
049564,001181: 7345 07056 TC DMPSUB # DO SAME FOR Y AND Z COMPONENTS.
049565,001182: 7346 07107 TC VROUND
049566,001183: 7347 52162 DXCH MPAC +5
049567,001184: 7350 52155 DXCH MPAC
049568,001185: 7351 52162 DXCH MPAC +5
049569,001186:
049570,001187: 7352 07056 TC DMPSUB
049571,001188: 7353 07107 TC VROUND
049572,001189: 7354 52155 VROTATEX DXCH MPAC # EXIT USED TO RESTORE MPAC AFTER THIS
049573,001190: 7355 52162 DXCH MPAC +5 # TYPE OF ROTATION. CALLED BY VECTOR SHIFT
049574,001191: 7356 52160 DXCH MPAC +3 # RIGHT, V/SC, ETC.
049575,001192: 7357 52155 DXCH MPAC
049576,001193: 7360 16030 TCF DANZIG
049577,001194:
|
|
Page 1146 |
049579,001196: # DP VECTOR PROJECTION ROUTINE.
049580,001197:
049581,001198: 7361 07123 VPROJ TC PREDOT # (MPAC.X)MPAC IS COMPUTED AND LEFT IN
049582,001199: 7362 44766 CS FOUR # MPAC. DO DOT AND FALL INTO DVXSC.
049583,001200: 7363 26116 ADS ADDRWD
049584,001201:
049585,001202: # VXSC WHEN SCALAR ARRIVES IN MPAC AND VECTOR IS AT X.
049586,001203:
049587,001204: 7364 00006 DVXSC EXTEND # SAVE SCALAR IN MPAC +3 AND GET X
049588,001205: 7365 30155 DCA MPAC # COMPONENT OF ANSWER.
049589,001206: 7366 52160 DXCH MPAC +3
049590,001207: 7367 07056 TC DMPSUB
049591,001208: 7370 07107 TC VROUND
049592,001209:
049593,001210: 7371 34767 CAF TWO # ADVANCE ADDRWD TO Y COMPONENT OF X.
049594,001211: 7372 26116 ADS ADDRWD
049595,001212: 7373 00006 EXTEND
049596,001213: 7374 30160 DCA MPAC +3 # PUT SCALAR BACK INTO MPAC AND SAVE
049597,001214: 7375 52155 DXCH MPAC # X RESULT IN MPAC +5.
049598,001215: 7376 52162 DXCH MPAC +5
049599,001216: 7377 07056 TC DMPSUB
049600,001217: 7400 07107 TC VROUND
049601,001218:
049602,001219: 7401 34767 CAF TWO
049603,001220: 7402 26116 ADS ADDRWD # TO Z COMPONENT.
049604,001221: 7403 52160 DXCH MPAC +3 # BRING SCALAR BACK, PUTTING Y RESULT IN
049605,001222: 7404 52155 DXCH MPAC # THE PROPER PLACE.
049606,001223: 7405 52160 DXCH MPAC +3
049607,001224: 7406 07056 TC DMPSUB
049608,001225: 7407 07107 TC VROUND
049609,001226:
049610,001227: 7410 52155 DXCH MPAC # PUT Z COMPONENT IN PROPER PLACE, ALSO
049611,001228: 7411 52162 DXCH MPAC +5 # POSITIONING X.
049612,001229: 7412 52155 DXCH MPAC
049613,001230:
049614,001231: 7413 16470 TCF VMODE # MODE HAS CHANGED TO VECTOR.
049615,001232:
|
|
Page 1147 |
049617,001234: # THE VECTOR CROSS PRODUCT ROUTINE CALCULATES (X M -X M ,X M -X M ,X M -X M ) WHERE M IS THE VECTOR IN
049618,001235: # 3 2 2 3 1 3 3 1 2 1 1 2
049619,001236: # MPAC AND X THE VECTOR AT THE GIVEN ADDRESS.
049620,001237:
049621,001238: 7414 00006 VXV EXTEND
049622,001239: 7415 30162 DCA MPAC +5 # FORM UP M3X1, LEAVING M1 IN VBUF.
049623,001240: 7416 52155 DXCH MPAC
049624,001241: 7417 52123 DXCH VBUF
049625,001242: 7420 07056 TC DMPSUB # BY X1.
049626,001243:
049627,001244: 7421 00006 EXTEND
049628,001245: 7422 40160 DCS MPAC +3 # CALCULATE -X1M2, SAVING X1M3 IN VBUF +2.
049629,001246: 7423 52155 DXCH MPAC
049630,001247: 7424 52125 DXCH VBUF +2
049631,001248: 7425 07056 TC DMPSUB
049632,001249:
049633,001250: 7426 34767 CAF TWO # ADVANCE ADDRWD TO X2.
049634,001251: 7427 26116 ADS ADDRWD
049635,001252: 7430 00006 EXTEND
049636,001253: 7431 40162 DCS MPAC +5 # PREPARE TO GET -X2M3, SAVING -X1M2 IN
049637,001254: 7432 52155 DXCH MPAC # MPAC +5.
049638,001255: 7433 52162 DXCH MPAC +5
049639,001256: 7434 07056 TC DMPSUB
049640,001257:
049641,001258: 7435 00006 EXTEND
049642,001259: 7436 30123 DCA VBUF # GET X2M1, SAVING -X2M3 IN VBUF +4.
049643,001260: 7437 52155 DXCH MPAC
049644,001261: 7440 52127 DXCH VBUF +4
049645,001262: 7441 07056 TC DMPSUB
049646,001263:
049647,001264: 7442 34767 CAF TWO # ADVANCE ADDRWD TO X3.
049648,001265: 7443 26116 ADS ADDRWD
049649,001266: 7444 00006 EXTEND
049650,001267: 7445 40123 DCS VBUF # GET -X3M1, ADDING X2M1 TO MPAC +5 TO
049651,001268: 7446 52155 DXCH MPAC # COMPLETE THE Z COMPONENT OF THE ANSWER.
049652,001269: 7447 20162 DAS MPAC +5
049653,001270:
049654,001271: 7450 00006 EXTEND
049655,001272: 7451 17453 BZF +2
049656,001273: 7452 06760 TC OVERFLWZ
049657,001274:
049658,001275: 7453 07056 TC DMPSUB
049659,001276: 7454 52125 DXCH VBUF +2 # MOVE X1M3 TO MPAC +3 SETTING UP FOR X3M2
049660,001277: 7455 52160 DXCH MPAC +3 # AND ADD -X3M1 TO MPAC +3 TO COMPLETE THE
049661,001278: 7456 52155 DXCH MPAC # Y COMPONENT OF THE RESULT.
049662,001279: 7457 20160 DAS MPAC +3
049663,001280:
049664,001281: 7460 00006 EXTEND
049665,001282: 7461 17463 BZF +2
|
|
Page 1148 |
049667,001284: 7462 06763 TC OVERFLWY
049668,001285:
049669,001286: 7463 07056 TC DMPSUB
049670,001287: 7464 52127 DXCH VBUF +4 # GO ADD -X2M3 TO X3M2 TO COMPLETE THE X
049671,001288: 7465 16747 TCF ENDVXV # COMPONENT (TAIL END OF DAD).
049672,001289:
049673,001290: # THE MPACVBUF SUBROUTINE SAVES THE VECTOR IN MPAC IN VBUF WITHOUT CLOBBERING MPAC.
049674,001291:
049675,001292: 7466 00006 MPACVBUF EXTEND # CALLED BY MXV, VXM, AND UNIT.
049676,001293: 7467 30155 DCA MPAC
049677,001294: 7470 52123 DXCH VBUF
049678,001295: 7471 00006 EXTEND
049679,001296: 7472 30160 DCA MPAC +3
049680,001297: 7473 52125 DXCH VBUF +2
049681,001298: 7474 00006 EXTEND
049682,001299: 7475 30162 DCA MPAC +5
049683,001300: 7476 52127 DXCH VBUF +4
049684,001301: 7477 00002 TC Q # RETURN TO CALLER.
049685,001302:
049686,001303: # DOUBLE PRECISION SIGN AGREE ROUTINE. ARRIVE WITH INPUT IN A+L. OUTPUT IS IN A + L.
049687,001304:
049688,001305: 7500 10000 ALSIGNAG CCS A # TEST UPPER PART.
049689,001306: 7501 17505 TCF UPPOS # IT IS POSITIVE
049690,001307: 7502 00002 TC Q # ZERO
049691,001308: 7503 17515 TCF UPNEG # NEGATIVE
049692,001309: 7504 00002 TC Q # ZERO
049693,001310:
049694,001311: 7505 56001 UPPOS XCH L # SAVE DECREMENTED UPPER PART.
049695,001312: 7506 64753 AD HALF
049696,001313: 7507 64753 AD HALF
049697,001314: 7510 54000 TS A # SKIPS ON OVERFLOW
049698,001315: 7511 17513 TCF +2
049699,001316: 7512 24001 INCR L # RESTORE UPPER TO ROIGINAL VALUE
049700,001317: 7513 56001 XCH L # SWAP A + L BACK.
049701,001318: 7514 00002 TC Q
049702,001319:
049703,001320: 7515 56001 UPNEG XCH L # SAVE COMPLEMENTED + DECREMENTED UPPER PT
049704,001321: 7516 64752 AD NEGMAX
049705,001322: 7517 67702 AD NEGONE
049706,001323: 7520 54000 TS A
049707,001324: 7521 17523 TCF +2 # DONT INCREMENT IF NO OVERFLOW.
049708,001325: 7522 24001 INCR L
049709,001326: 7523 56001 XCH L
049710,001327: 7524 40000 COM # MAKE NEGATIVE AGAIN.
049711,001328: 7525 00002 TC Q
049712,001329:
|
|
Page 1149 |
049714,001331: # INTERPRETIVE INSTRUCTIONS WHOSE EXECUTION CONSISTS OF PRINCIPALLY CALLING SUBROUTINES.
049715,001332:
049716,001333: 7526 07056 DMP1 TC DMPSUB # DMP INSTRUCTION.
049717,001334: 7527 16030 TCF DANZIG
049718,001335:
049719,001336: 7530 07056 DMPR TC DMPSUB
049720,001337: 7531 07106 TC ROUNDSUB +1 # (C(A) = +0).
049721,001338: 7532 16030 TCF DANZIG
049722,001339:
049723,001340: 7533 00006 DDV EXTEND
049724,001341: 7534 50116 INDEX ADDRWD # MOVE DIVIDEND INTO BUF.
049725,001342: 7535 30001 DCA 0
049726,001343: 7536 17543 TCF BDDV +4
049727,001344:
049728,001345: 7537 00006 BDDV EXTEND # MOVE DIVISOR INTO MPAC SAVING MPAC, THE
049729,001346: 7540 50116 INDEX ADDRWD # DIVIDEND, IN BUF.
049730,001347: 7541 30001 DCA 0
049731,001348: 7542 52155 DXCH MPAC
049732,001349: 7543 52131 +4 DXCH BUF
049733,001350: 7544 34772 CAF ZERO # DIVIDE ROUTINES IN BANK 0.
049734,001351: 7545 54004 TS FBANK
049735,001352: 7546 12353 TCF DDV/BDDV
049736,001353:
049737,001354: 7547 30116 SETPD CA ADDRWD # MUST SET TO WORK AREA, OR EBANK TROUBLE.
049738,001355: 7550 54166 TS PUSHLOC
049739,001356: 7551 16032 TCF NOIBNKSW # NO FBANK SWITCH REQUIRED.
049740,001357:
049741,001358: 7552 34772 TSLC CAF ZERO # SHIFTING ROUTINES LOCATED IN BANK 00.
049742,001359: 7553 54004 TS FBANK
049743,001360: 7554 12172 TCF TSLC2
049744,001361:
049745,001362: 7555 36043 GSHIFT CAF LOW7 # USED AS MASK AT GENSHIFT. THIS PROCESSES
049746,001363: 7556 54004 TS FBANK # ANY SHIFT INSTRUCTION (EXCEPT TSLC) WITH
049747,001364: 7557 12214 TCF GENSHIFT # AN ADDRESS (ROUTINES IN BANK 0).
049748,001365:
|
|
Page 1150 |
049750,001367: # THE FOLLOWING IS THE PROLOGUE TO V/SC. IF THE PRESENT MODE IS VECTOR, IT SAVES THE SCALAR AT X IN BUF
049751,001368: # AND CALLS THE V/SC ROUTINE IN BANK 0. IF THE PRESENT MODE IS SCALAR, IT MOVES THE VECTOR AT X INTO MPAC, SAVING
049752,001369: # THE SCALAR IN MPAC IN BUF BEFORE CALLING THE V/SC ROUTINE IN BANK 0.
049753,001370:
049754,001371: 7560 10163 V/SC CCS MODE
049755,001372: 7561 17572 TCF DV/SC # MOVE VECTOR INTO MPAC.
049756,001373: 7562 17572 TCF DV/SC
049757,001374:
049758,001375: 7563 00006 VV/SC EXTEND
049759,001376: 7564 50116 INDEX ADDRWD
049760,001377: 7565 30001 DCA 0
049761,001378: 7566 52131 V/SC1 DXCH BUF # IN BOTH CASES, VECTOR IS NOW IN MPAC AND
049762,001379: 7567 34772 CAF ZERO # SCALAR IN BUF.
049763,001380: 7570 54004 TS FBANK
049764,001381: 7571 12653 TCF V/SC2
049765,001382:
049766,001383: 7572 00006 DV/SC EXTEND
049767,001384: 7573 50116 INDEX ADDRWD
049768,001385: 7574 30003 DCA 2
049769,001386: 7575 52160 DXCH MPAC +3
049770,001387: 7576 00006 EXTEND
049771,001388: 7577 50116 INDEX ADDRWD
049772,001389: 7600 30005 DCA 4
049773,001390: 7601 52162 DXCH MPAC +5
049774,001391:
049775,001392: 7602 44770 CS ONE # CHANGE MODE TO VECTOR.
049776,001393: 7603 54163 TS MODE
049777,001394:
049778,001395: 7604 00006 EXTEND
049779,001396: 7605 50116 INDEX ADDRWD
049780,001397: 7606 30001 DCA 0
049781,001398: 7607 52155 DXCH MPAC
049782,001399: 7610 17566 TCF V/SC1 # FINISH PROLOGUE AT COMMON SECTION.
049783,001400:
|
|
Page 1151 |
049785,001402: # SIGN AND COMPLEMENT INSTRUCTIONS.
049786,001403:
049787,001404: 7611 50116 SIGN INDEX ADDRWD # CALL COMP INSTRUCTION IF WORD AT X IS
049788,001405: 7612 10000 CCS 0 # NEGATIVE NON-ZERO.
049789,001406: 7613 16030 TCF DANZIG
049790,001407: 7614 17616 TCF +2
049791,001408: 7615 17624 TCF COMP # DO THE COMPLEMENT.
049792,001409:
049793,001410: 7616 50116 INDEX ADDRWD
049794,001411: 7617 10001 CCSL CCS 1
049795,001412: 7620 16030 TCF DANZIG
049796,001413: 7621 16030 TCF DANZIG
049797,001414: 7622 17624 TCF COMP
049798,001415: 7623 16030 TCF DANZIG
049799,001416: 7624 00006 COMP EXTEND # COMPLEMENT DP MPAC IN EVERY CASE.
049800,001417: 7625 40155 DCS MPAC
049801,001418: 7626 52155 DXCH MPAC
049802,001419:
049803,001420: 7627 10163 CCS MODE # EITHER COMPLEMENT MPAC +3 OR THE REST OF
049804,001421: 7630 17641 TCF DCOMP # THE VECTOR ACCUMULATOR.
049805,001422: 7631 17641 TCF DCOMP
049806,001423:
049807,001424: 7632 00006 EXTEND # VECTOR COMPLEMENT.
049808,001425: 7633 40160 DCS MPAC +3
049809,001426: 7634 52160 DXCH MPAC +3
049810,001427: 7635 00006 EXTEND
049811,001428: 7636 40162 DCS MPAC +5
049812,001429: 7637 52162 DXCH MPAC +5
049813,001430: 7640 16030 TCF DANZIG
049814,001431:
049815,001432: 7641 40156 DCOMP CS MPAC +2
049816,001433: 7642 54156 TS MPAC +2
049817,001434: 7643 16030 TCF DANZIG
049818,001435:
|
|
Page 1152 |
049820,001437: # THE FOLLOWING SHORT SHIFT CODES REQUIRE NO ADDRESS WORD:
049821,001438:
049822,001439: # 1. SR1 TO SR4 SCALAR SHIFT RIGHT.
049823,001440: # 2. SR1R TO SR4R SCALAR SHIFT RIGHT AND ROUND.
049824,001441: # 3. SL1 TO SL4 SCALAR SHIFT LEFT.
049825,001442: # 4. SL1R TO SL4R SCALAR SHIFT LEFT AND ROUND.
049826,001443:
049827,001444: # 5. VSR1 TO VSR8 VECTOR SHIFT RIGHT (ALWAYS ROUNDS).
049828,001445: # 6. VSL1 TO VSL8 VECTOR SHIFT LEFT (NEVER ROUNDS).
049829,001446:
049830,001447: # THE FOLLOWING CODES REQUIRE AN ADDRESS WHICH MAY BE INDEXED:*
049831,001448:
049832,001449: # 1. SR SCALAR SHIFT RIGHT.
049833,001450: # 2. SRR SCALAR SHIFT RIGHT AND ROUND.
049834,001451: # 3. SL SCALAR SHIFT LEFT.
049835,001452: # 4. SLR SCALAR SHIFT LEFT AND ROUND.
049836,001453:
049837,001454: # 5. VSR VECTOR SHIFT RIGHT.
049838,001455: # 6. VSL VECTOR SHIFT LEFT.
049839,001456:
049840,001457: # * IF THE ADDRESS IS INDEXED, AND THE INDEX MODIFICATION RESULTS IN A NEGATIVE SHIFT COUNT, A SHIFT OF THE
049841,001458: # ABSOLUTE VALUE OF THE COUNT IS DONE IN THE OPPOSITE DIRECTION.
049842,001459:
049843,001460: 00,2000 SETLOC INTERB0
049844,001461: 00,2000 BANK
049845,001462:
049846,001463: 00,2017 COUNT* $$/INTER
049847,001464: 00,2017 36211 SHORTT CAF SIX # SCALAR SHORT SHIFTS COME HERE. THE SHIFT
049848,001465: 00,2020 70020 MASK CYR # COUNT-1 IS NOW IN BITS 2-3 OF CYR. THE
049849,001466: 00,2021 54021 TS SR # ROUNDING BIT IS IN BIT1 AT THIS POINT.
049850,001467:
049851,001468: 00,2022 10020 CCS CYR # SEE IF RIGHT OR LEFT SHIFT DESIRED.
049852,001469: 00,2023 12101 TCF TSSL # SHIFT LEFT.
049853,001470:
049854,001471: 00,2024 00024 SRDDV DEC 20 B-14 # MPTEMP SETTING FOR SR BEFORE DDV.
049855,001472:
049856,001473: 00,2025 50021 TSSR INDEX SR # GET SHIFTING BIT.
049857,001474: 00,2026 34753 CAF BIT14
049858,001475: 00,2027 54135 TS MPTEMP
049859,001476:
049860,001477: 00,2030 10020 CCS CYR # SEE IF A ROUND IS DESIRED.
049861,001478: 00,2031 02050 RIGHTR TC MPACSRND # YES - SHIFT RIGHT AND ROUND.
049862,001479: 00,2032 16027 TCF NEWMODE # SET MODE TO DP (C(A) = 0).
049863,001480: 00,2033 30135 MPACSHR CA MPTEMP # DO A TRIPLE PRECISION SHIFT RIGHT.
049864,001481: 00,2034 00006 EXTEND
049865,001482: 00,2035 70156 MP MPAC +2
049866,001483: 00,2036 54156 +3 TS MPAC +2 # (EXIT FROM SQRT AND ABVAL).
049867,001484: 00,2037 30135 CA MPTEMP
049868,001485: 00,2040 00006 EXTEND
|
|
Page 1153 |
049870,001487: 00,2041 70154 MP MPAC # SHIFT MAJOR PART INTO A,L AND PLACE IN
049871,001488: 00,2042 52155 DXCH MPAC # MPAC,+1.
049872,001489: 00,2043 30135 CA MPTEMP
049873,001490: 00,2044 00006 EXTEND
049874,001491: 00,2045 70001 MP L # ORIGINAL C(MPAC +1).
049875,001492: 00,2046 20156 DAS MPAC +1 # GUARANTEED NO OVERFLOW.
049876,001493: 00,2047 16030 TCF DANZIG
049877,001494:
049878,001495: # MPAC SHIFT RIGHT AND ROUND SUBROUTINES.
049879,001496:
049880,001497: 00,2050 30156 MPACSRND CA MPAC +2 # WE HAVE TO DO ALL THREE MULTIPLIES SINCE
049881,001498: 00,2051 00006 EXTEND # MPAC +1 AND MPAC +2 MIGHT HAVE SIGN
049882,001499: 00,2052 70135 MP MPTEMP # DISAGREEMENT WITH A SHIFT RIGHT OF 1.
049883,001500: 00,2053 56155 XCH MPAC +1
049884,001501: 00,2054 00006 EXTEND
049885,001502: 00,2055 70135 MP MPTEMP
049886,001503: 00,2056 56155 XCH MPAC +1 # TRIAL MINOR PART.
049887,001504: 00,2057 60001 AD L
049888,001505:
049889,001506: 00,2060 60000 VSHR2 DOUBLE # (FINISH VECTOR COMPONENT SHIFT RIGHT
049890,001507: 00,2061 54156 TS MPAC +2 # AND ROUND.
049891,001508: 00,2062 12064 TCF +2
049892,001509: 00,2063 26155 ADS MPAC +1 # GUARANTEED NO OVERFLOW.
049893,001510:
049894,001511: 00,2064 34772 CAF ZERO
049895,001512: 00,2065 54156 TS MPAC +2
049896,001513: 00,2066 56154 XCH MPAC # SETTING TO ZERO SO FOLLOWING DAS WORKS.
049897,001514: 00,2067 00006 EXTEND
049898,001515: 00,2070 70135 MP MPTEMP
049899,001516: 00,2071 20155 DAS MPAC # AGAIN NO OVERFLOW.
049900,001517: 00,2072 00002 TC Q
049901,001518:
049902,001519: 00,2073 30135 VSHRRND CA MPTEMP # ENTRY TO SHIFT RIGHT AND ROUND MPAC WHEN
049903,001520: 00,2074 00006 EXTEND # MPAC CONTAINS A VECTOR COMPONENT.
049904,001521: 00,2075 70155 MP MPAC +1
049905,001522: 00,2076 54155 TS MPAC +1
049906,001523: 00,2077 56001 XCH L
049907,001524: 00,2100 12060 TCF VSHR2 # GO ADD ONE IF NECESSARY AND FINISH.
049908,001525:
|
|
Page 1154 |
049910,001527: # ROUTINE FOR SHORT SCALAR SHIFT LEFT (AND MAYBE ROUND).
049911,001528:
049912,001529: 00,2101 30021 TSSL CA SR # GET SHIFT COUNT FOR SR.
049913,001530: 00,2102 54135 +1 TS MPTEMP
049914,001531:
049915,001532: 00,2103 00006 +2 EXTEND # ENTRY HERE FROM SL FOR SCALARS.
049916,001533: 00,2104 30156 DCA MPAC +1 # SHIFTING LEFT ONE PLACE AT A TIME IS
049917,001534: 00,2105 20156 DAS MPAC +1 # FASTER THAN DOING THE WHOLE SHIFT WITH
049918,001535: 00,2106 60154 AD MPAC # MULTIPLIES ASSUMING THAT FREQUENCY OF
049919,001536: 00,2107 60154 AD MPAC # SHIFT COUNTS GOES DOWN RAPIDLY AS A
049920,001537: 00,2110 54154 TS MPAC # FUNCTION OF THEIR MAGNITUDE.
049921,001538: 00,2111 12113 TCF +2
049922,001539: 00,2112 54121 TS OVFIND # OVERFLOW. (LEAVES OVERFLOW-CORRECTED
049923,001540: # RESULT ANYWAY).
049924,001541: 00,2113 10135 CCS MPTEMP # LOOP ON DECREMENTED SHIFT COUNT.
049925,001542: 00,2114 12102 TCF TSSL +1
049926,001543:
049927,001544: 00,2115 10020 CCS CYR # SEE IF ROUND WANTED.
049928,001545: 00,2116 07105 ROUND TC ROUNDSUB # YES - ROUND AND EXIT.
049929,001546: 00,2117 16030 TCF DANZIG # SL LEAVES A ZERO IN CYR FOR NO ROUND.
049930,001547: 00,2120 16030 TCF DANZIG # NO - EXIT IMMEDIATL
049931,001548:
|
|
Page 1155 |
049933,001550: # VECTOR SHIFTING ROUTINES.
049934,001551:
049935,001552: 00,2121 34774 SHORTV CAF LOW3 # SAVE 3 BIT SHIFT COUNT - 1 WITHOUT
049936,001553: 00,2122 70020 MASK CYR # EDITING CYR.
049937,001554: 00,2123 54135 TS MPTEMP
049938,001555: 00,2124 10020 CCS CYR # SEE IF LEFT OR RIGHT SHIFT.
049939,001556: 00,2125 12145 TCF VSSL # VECTOR SHIFT LEFT.
049940,001557: 00,2126 00176 OCT176 OCT 176 # USED IN PROCESSED SHIFTS WITH - COUNT.
049941,001558:
049942,001559: 00,2127 50135 VSSR INDEX MPTEMP # (ENTRY FROM SR). PICK UP SHIFTING BIT.
049943,001560: 00,2130 34753 CAF BIT14 # MPTEMP CONTAINS THE SHIFT COUNT - 1.
049944,001561: 00,2131 54135 TS MPTEMP
049945,001562: 00,2132 02073 TC VSHRRND # SHIFT X COMPONENT.
049946,001563:
049947,001564: 00,2133 52155 DXCH MPAC # SWAP X AND Y COMPONENTS.
049948,001565: 00,2134 52160 DXCH MPAC +3
049949,001566: 00,2135 52155 DXCH MPAC
049950,001567: 00,2136 02073 TC VSHRRND # SHIFT Y COMPONENT.
049951,001568:
049952,001569: 00,2137 52155 DXCH MPAC # SWAP Y AND Z COMPONENTS.
049953,001570: 00,2140 52162 DXCH MPAC +5
049954,001571: 00,2141 52155 DXCH MPAC
049955,001572: 00,2142 02073 TC VSHRRND # SHIFT Z COMPONENT.
049956,001573:
049957,001574: 00,2143 17354 TCF VROTATEX # RESTORE COMPONENTS TO PROPER PLACES.
049958,001575:
|
|
Page 1156 |
049960,001577: # VECTOR SHIFT LEFT - DONE ONE PLACE AT A TIME.
049961,001578:
049962,001579: 00,2144 54135 -1 TS MPTEMP # SHIFTING LOOP.
049963,001580:
049964,001581: 00,2145 00006 VSSL EXTEND
049965,001582: 00,2146 30155 DCA MPAC
049966,001583: 00,2147 20155 DAS MPAC
049967,001584: 00,2150 00006 EXTEND
049968,001585: 00,2151 12153 BZF +2
049969,001586: 00,2152 06766 TC OVERFLOW
049970,001587:
049971,001588: 00,2153 00006 EXTEND
049972,001589: 00,2154 30160 DCA MPAC +3
049973,001590: 00,2155 20160 DAS MPAC +3
049974,001591: 00,2156 00006 EXTEND
049975,001592: 00,2157 12161 BZF +2
049976,001593: 00,2160 06763 TC OVERFLWY
049977,001594:
049978,001595: 00,2161 00006 EXTEND
049979,001596: 00,2162 30162 DCA MPAC +5
049980,001597: 00,2163 20162 DAS MPAC +5
049981,001598: 00,2164 00006 EXTEND
049982,001599: 00,2165 12167 BZF +2
049983,001600: 00,2166 06760 TC OVERFLWZ
049984,001601:
049985,001602: 00,2167 10135 CCS MPTEMP # LOOP ON DECREMENTED SHIFT COUNTER.
049986,001603: 00,2170 12144 TCF VSSL -1
049987,001604: 00,2171 16030 TCF DANZIG # EXIT.
049988,001605:
|
|
Page 1157 |
049990,001607: # TSLC - TRIPLE SHIFT LEFT AND COUNT. SHIFTS MPAC LEFT UNTIL GREATER THAN .5 IN MAGNITUDE, LEAVING
049991,001608: # THE COMPLEMENT OF THE NUMBER OF SHIFTS REQUIRED IN X.
049992,001609:
049993,001610: 00,2172 54135 TSLC2 TS MPTEMP # START BY ZEROING SHIFT COUNT (IN A NOW).
049994,001611: 00,2173 06672 TC BRANCH # EXIT WITH NO SHIFTING IF ARGUMENT ZERO.
049995,001612: 00,2174 12176 TCF +2
049996,001613: 00,2175 12212 TCF ENDTSLC # STORES ZERO SHIFT COUNT IN THIS CASE.
049997,001614:
049998,001615: 00,2176 07226 TC TPAGREE # MAY CAUSE UPSHIFT OF ONE EXTRA PLACE.
049999,001616:
050000,001617: 00,2177 30154 CA MPAC # BEGIN NORMALIZATION LOOP.
050001,001618: 00,2200 12207 TCF TSLCTEST
050002,001619:
050003,001620: 00,2201 24135 TSLCLOOP INCR MPTEMP # INCREMENT SHIFT COUNTER.
050004,001621: 00,2202 00006 EXTEND
050005,001622: 00,2203 30156 DCA MPAC +1
050006,001623: 00,2204 20156 DAS MPAC +1
050007,001624: 00,2205 60154 AD MPAC
050008,001625: 00,2206 26154 ADS MPAC
050009,001626: 00,2207 60000 TSLCTEST DOUBLE # SEE IF (ANOTHER) SHIFT IS REQUIRED.
050010,001627: 00,2210 54000 OVSK
050011,001628: 00,2211 12201 TCF TSLCLOOP # YES - INCREMENT COUNT AND SHIFT AGAIN.
050012,001629:
050013,001630: 00,2212 40135 ENDTSLC CS MPTEMP
050014,001631: 00,2213 16572 TCF STORE1 # STORE SHIFT COUNT AND RETURN TO DANZIG.
050015,001632:
|
|
Page 1158 |
050017,001634: # THE FOLLOWING ROUTINES PROCESSES THE GENERAL SHIFT INSTRUCTIONS SR, SRR, SL, AND SLR.
050018,001635: # THE GIVEN ADDRESS IS DECODED AS FOLLOWS:
050019,001636:
050020,001637: # BITS 1-7 SHIFT COUNT (SUBADDRESS) LESS THAN 125 DECIMAL.
050021,001638: # BIT 8 PSEUDO SIGN BIT (DETECTS CHANGE IN SIGN IN INDEXED SHIFTS).
050022,001639: # BIT 9 0 FOR LEFT SHIFT, AND 1 FOR RIGHT SHIFT.
050023,001640: # BIT 10 1 FOR TERMINAL ROUND ON SCALAR SHIFTS, 0 OTHERWISE.
050024,001641: # BITS 11-13 0.
050025,001642: # BIT 14 1.
050026,001643: # BIT 15 0.
050027,001644:
050028,001645: # THE ABOVE ENCODING IS DONE BY THE YUL SYSTEM.
050029,001646:
050030,001647: 00,2214 70116 GENSHIFT MASK ADDRWD # GET SHIFT COUNT, TESTING FOR ZERO.
050031,001648: 00,2215 10000 CCS A # (ARRIVES WITH C(A) = LOW7).
050032,001649: 00,2216 12224 TCF GENSHFT2 # IF NON-ZERO, PROCEED WITH DECREMENTED CT
050033,001650:
050034,001651: 00,2217 34757 CAF BIT10 # ZERO SHIFT COUNT. NO SHIFTS NEEDED BUT
050035,001652: 00,2220 70116 MASK ADDRWD # WE MIGHT HAVE TO ROUND MPAC ON SLR AND
050036,001653: 00,2221 10000 CCS A # SRR (SCALAR ONLY).
050037,001654: 00,2222 07105 TC ROUNDSUB
050038,001655: 00,2223 16030 TCF DANZIG
050039,001656:
050040,001657: 00,2224 54135 GENSHFT2 TS MPTEMP # DECREMENTED SHIFT COUNT TO MPTEMP.
050041,001658: 00,2225 34761 CAF BIT8 # TEST MEANING OF LOW SEVEN BIT COUNT IN
050042,001659: 00,2226 00006 EXTEND # MPTEMP NOW.
050043,001660: 00,2227 70116 MP ADDRWD
050044,001661: 00,2230 76214 MASK LOW2 # JUMPS ON SHIFT DIRECTION (BIT8) AND
050045,001662: 00,2231 50000 INDEX A
050046,001663: 00,2232 12233 TCF +1 # ORIGINAL SHIFT DIRECTION (BIT 9).
050047,001664: 00,2233 12332 TCF RIGHT- # NEGATIVE SHIFT COUNT FOR SL OR SLR.
050048,001665: 00,2234 12342 TCF LEFT # SL OR SLR.
050049,001666: 00,2235 12336 TCF LEFT- # NEGATIVE SHIFT COUNT WITH SR OR SRR.
050050,001667:
|
|
Page 1159 |
050052,001669: # GENERAL SHIFT RIGHT.
050053,001670:
050054,001671: 00,2236 10163 RIGHT CCS MODE # SEE IF VECTOR OR SCALAR.
050055,001672: 00,2237 12277 TCF GENSCR
050056,001673: 00,2240 12277 TCF GENSCR
050057,001674:
050058,001675: 00,2241 30135 CA MPTEMP # SEE IF SHIFT COUNT LESS THAN 14D.
050059,001676: 00,2242 63726 VRIGHT2 AD NEG12
050060,001677: 00,2243 00006 EXTEND
050061,001678: 00,2244 62127 BZMF VSSR # IF SO, BRANCH AND SHIFT IMMEDIATELY.
050062,001679:
050063,001680: 00,2245 67702 AD NEGONE # IF NOT, REDUCE MPTEMP BY A TOTAL OF 14,
050064,001681: 00,2246 54135 TS MPTEMP # AND DO A SHIFT RIGHT AND ROUND BY 14.
050065,001682: 00,2247 34772 CAF ZERO # THE ROUND AT THIS STAGE MAY INTRODUCE A
050066,001683: 00,2250 54001 TS L # ONE BIT ERROR IN A SHIFT RIGHT 15D.
050067,001684: 00,2251 56154 XCH MPAC
050068,001685: 00,2252 56155 XCH MPAC +1
050069,001686: 00,2253 02272 TC SETROUND # X COMPONENT NOW SHIFTED, SO MAKE UP THE
050070,001687: 00,2254 20155 DAS MPAC # ROUNDING QUANTITY (0 IN A AND 0 OR +-1
050071,001688: # IN L).
050072,001689: 00,2255 56157 XCH MPAC +3 # REPEAT THE ABOVE PROCESS FOR Y AND Z.
050073,001690: 00,2256 56160 XCH MPAC +4
050074,001691: 00,2257 02272 TC SETROUND
050075,001692: 00,2260 20160 DAS MPAC +3 # NO OVERFLOW ON THESE ADDS.
050076,001693:
050077,001694: 00,2261 56161 XCH MPAC +5
050078,001695: 00,2262 56162 XCH MPAC +6
050079,001696: 00,2263 02272 TC SETROUND
050080,001697: 00,2264 20162 DAS MPAC +5
050081,001698:
050082,001699: 00,2265 10135 CCS MPTEMP # SEE IF DONE, DOING FINAL DECREMENT.
050083,001700: 00,2266 54135 TS MPTEMP
050084,001701: 00,2267 12242 TCF VRIGHT2
050085,001702: 00,2270 04604 BIASLO DEC .2974 B-1 # SQRT CONSTANT
050086,001703:
050087,001704: 00,2271 16030 TCF DANZIG
050088,001705:
050089,001706: 00,2272 60000 SETROUND DOUBLE # MAKES UP ROUNDING QUANTITY FROM ARRIVING
050090,001707: 00,2273 54156 TS MPAC +2 # C(A). L IS ZERO INITIALLY.
050091,001708: 00,2274 34772 CAF ZERO
050092,001709: 00,2275 56001 XCH L
050093,001710: 00,2276 00002 TC Q # RETURN AND DO THE DAS, RESETTING L TO 0.
050094,001711:
|
|
Page 1160 |
050096,001713: # PROCESS SR AND SRR FOR SCALARS.
050097,001714:
050098,001715: 00,2277 30135 GENSCR CA MPTEMP # SEE IF THE ORIGINAL SHIFT COUNT WAS LESS
050099,001716: 00,2300 63726 +1 AD NEG12 # THAN 14D.
050100,001717: 00,2301 00006 EXTEND
050101,001718: 00,2302 62322 BZMF DOSSHFT # DO THE SHIFT IMMEDIATELY IF SO.
050102,001719:
050103,001720: 00,2303 67702 +4 AD NEGONE # IF NOT, DECREMENT SHIFT COUNT BY 14D AND
050104,001721: 00,2304 54135 TS MPTEMP # SHIFT MPAC RIGHT 14 PLACES.
050105,001722: 00,2305 34772 CAF ZERO
050106,001723: 00,2306 56154 XCH MPAC
050107,001724: 00,2307 56155 XCH MPAC +1
050108,001725: 00,2310 54156 TS MPAC +2
050109,001726: 00,2311 10135 CCS MPTEMP # SEE IF FINISHED, DO FINAL DECREMENT.
050110,001727: 00,2312 54135 TS MPTEMP
050111,001728: 00,2313 02300 TC GENSCR +1
050112,001729: 00,2314 22650 SLOPEHI DEC .5884 # SQRT CONSTANT.
050113,001730: 00,2315 34757 CAF BIT10 # FINISHED WITH SHIFT. SEE IF ROUND
050114,001731: 00,2316 70116 MASK ADDRWD # WANTED.
050115,001732: 00,2317 10000 CCS A
050116,001733: 00,2320 07105 TC ROUNDSUB
050117,001734: 00,2321 16030 TCF DANZIG # DO SO AND/OR EXIT.
050118,001735:
050119,001736: 00,2322 50135 DOSSHFT INDEX MPTEMP # PICK UP SHIFTING BIT.
050120,001737: 00,2323 34753 CAF BIT14
050121,001738: 00,2324 54135 TS MPTEMP
050122,001739: 00,2325 34757 CAF BIT10 # SEE IF TERMINAL ROUND DESIRED.
050123,001740: 00,2326 70116 MASK ADDRWD
050124,001741: 00,2327 10000 CCS A
050125,001742: 00,2330 12031 TCF RIGHTR # YES.
050126,001743: 00,2331 12033 TCF MPACSHR # JUST SHIFT RIGHT.
050127,001744:
|
|
Page 1161 |
050129,001746: # PROCESS THE RIGHT- (SL(R) WITH A NEGATIVE COUNT), LEFT-, AND LEFT OPTIONS.
050130,001747:
050131,001748: 00,2332 40135 RIGHT- CS MPTEMP # GET ABSOLUTE VALUE - 1 OF SHIFT COUNT
050132,001749: 00,2333 62126 AD OCT176 # UNDERSTANDING THAT BIT8 (PSEUDO-SIGN)
050133,001750: 00,2334 54135 TS MPTEMP # WAS 1 INITIALLY.
050134,001751: 00,2335 12236 TCF RIGHT # DO NORMAL SHIFT RIGHT.
050135,001752:
050136,001753: 00,2336 42126 LEFT- CS OCT176 # SAME PROLOGUE TO LEFT FOR INDEXED RIGHT
050137,001754: 00,2337 60135 AD MPTEMP # SHIFTS WHOSE NET SHIFT COUNT IS NEGATIVE
050138,001755: 00,2340 40000 COM
050139,001756: 00,2341 54135 TS MPTEMP
050140,001757:
050141,001758: 00,2342 10163 LEFT CCS MODE # SINCE LEFT SHIFTING IS SONE ONE PLACE AT
050142,001759: 00,2343 12346 TCF GENSCL # A TIME, NO COMPARISON WITH 14 NEED BE
050143,001760: 00,2344 12346 TCF GENSCL # DONE. FOR SCALARS, SEE IF TERMINAL ROUND
050144,001761: 00,2345 12145 TCF VSSL # DESIRED. FOR VECTORS, SHIFT IMMEDIATELY.
050145,001762:
050146,001763: 00,2346 40116 GENSCL CS ADDRWD # PUT ROUNDING BIT (BIT 10 OF ADDRWD) INTO
050147,001764: 00,2347 00006 EXTEND # BIT 15 OF CYR WHERE THE ROUNDING BIT OF
050148,001765: 00,2350 74763 MP BIT6 # A SHORT SHIFT LEFT WOULD BE
050149,001766: 00,2351 54020 TS CYR
050150,001767: 00,2352 12103 TCF TSSL +2 # DO THE SHIFT.
050151,001768:
|
|
Page 1162 |
050153,001770: # SCALAR DIVISION INSTRUCTIONS, DDV AND BDDV, ARE EXECUTED HERE. AT THIS POINT, THE DIVIDEND IS IN MPAC
050154,001771: # AND THE DIVISOR IN BUF.
050155,001772:
050156,001773: 00,2353 44770 DDV/BDDV CS ONE # INITIALIZATION.
050157,001774: 00,2354 54136 TS DVSIGN # +-1 FOR POSITIVE QUOTIENT - -0 FOR NEG.
050158,001775: 00,2355 54137 TS DVNORMCT # DIVIDEND NORMALIZATION COUNT.
050159,001776: 00,2356 54140 TS MAXDVSW # NEAR-ONE DIVIDE FLAG.
050160,001777:
050161,001778: 00,2357 10130 CCS BUF # FORCE BUF POSITIVE WITH THE MAJOR PART
050162,001779: 00,2360 12515 TCF BUFPOS # NON-ZERO.
050163,001780: 00,2361 12363 TCF +2
050164,001781: 00,2362 12530 TCF BUFNEG
050165,001782:
050166,001783: 00,2363 07225 BUFZERO TC BUFAGREE
050167,001784: 00,2364 10154 CCS MPAC # TEST TO SEE IF MPAC NON-ZERO. (TOO BIG)
050168,001785: 00,2365 12413 TCF OVF+ # MAJOR PART OF DIVIDEND IS POSITIVE NON-0
050169,001786: 00,2366 12370 TCF +2
050170,001787: 00,2367 12412 TCF OVF+ -1 # MAJOR PART OF DIVIDEND IS NEG. NON-ZERO
050171,001788:
050172,001789: 00,2370 56131 XCH BUF +1 # SHIFT DIVIDEND AND DIVISOR LEFT 14.
050173,001790: 00,2371 56130 XCH BUF
050174,001791: 00,2372 56155 XCH MPAC +1
050175,001792: 00,2373 56154 XCH MPAC
050176,001793: 00,2374 10130 CCS BUF # TRY AGAIN ON FORMER MINOR PART.
050177,001794: 00,2375 12421 TCF BUF+
050178,001795: 00,2376 12400 TCF +2 # OVERFLOW ON ZERO DIVISOR.
050179,001796: 00,2377 12415 TCF BUF-
050180,001797:
050181,001798: 00,2400 40154 CS MPAC # SIGN OF MPAC DETERMINES SIGN OF RESULT.
050182,001799: 00,2401 00006 SGNDVOVF EXTEND
050183,001800: 00,2402 62404 BZMF +2
050184,001801: 00,2403 24136 INCR DVSIGN # NEGMAX IN MPAC PERHAPS.
050185,001802: 00,2404 34750 DVOVF CAF POSMAX # ON DIVISION OVERFLOW OF ANY SORT, SET
050186,001803: 00,2405 54154 TS MPAC # SET DP MPAC TO +-POSMAX.
050187,001804: 00,2406 02627 TC FINALDV +3
050188,001805: 00,2407 34770 CAF ONE # SET OVERFLOW INDICATOR AND EXIT.
050189,001806: 00,2410 54121 TS OVFIND
050190,001807: 00,2411 06030 TC DANZIG
050191,001808:
050192,001809: 00,2412 24136 -1 INCR DVSIGN
050193,001810: 00,2413 40131 OVF+ CS BUF +1 # LOAD LOWER ORDER PART OF DIVISOR.
050194,001811: 00,2414 12401 TCF SGNDVOVF # GET SIGN OF RESULT.
050195,001812:
050196,001813: 00,2415 00006 BUF- EXTEND # IF BUF IS NEGATIVE, COMPLEMENT IT AND
050197,001814: 00,2416 40131 DCS BUF # MAINTAIN DVSIGN FOR FINAL QUOTIENT SIGN.
050198,001815: 00,2417 52131 DXCH BUF
050199,001816: 00,2420 24136 INCR DVSIGN # NOW -0.
050200,001817: 00,2421 10154 BUF+ CCS MPAC # FORCE MPAC POSITIVE, CHECKING FOR ZERO
|
|
Page 1163 |
050202,001819: 00,2422 12436 TCF MPAC+ # DIVIDEND IN THE PROCESS.
050203,001820: 00,2423 12425 TCF +2
050204,001821: 00,2424 12432 TCF MPAC-
050205,001822: 00,2425 10155 CCS MPAC +1
050206,001823: 00,2426 12436 TCF MPAC+
050207,001824: 00,2427 16030 TCF DANZIG # EXIT IMMEDIATELY ON ZERO DIVIDEND.
050208,001825: 00,2430 12432 TCF MPAC-
050209,001826: 00,2431 16030 TCF DANZIG
050210,001827:
050211,001828: 00,2432 00006 MPAC- EXTEND # FORCE MPAC POSITIVE AS BUF IN BUF-.
050212,001829: 00,2433 40155 DCS MPAC
050213,001830: 00,2434 52155 DXCH MPAC
050214,001831: 00,2435 24136 INCR DVSIGN # NOW +1 OR -0.
050215,001832:
|
|
Page 1164 |
050217,001834: 00,2436 40154 MPAC+ CS MPAC # CHECK FOR DIVISION OVERFLOW. IF THE
050218,001835: 00,2437 67702 AD NEGONE # MAJOR PART OF THE DIVIDEND IS LESS THAN
050219,001836: 00,2440 60130 AD BUF # THE MAJOR PART OF THE DIVISOR BY AT
050220,001837: 00,2441 10000 CCS A # LEAST TWO, WE CAN PROCEED IMMEDIATELY
050221,001838: 00,2442 12504 TCF DVNORM # WITHOUT NORMALIZATION PRODUCING A DVMAX.
050222,001839: 00,2443 60001 -1/2+2 OCT 60001 # USED IN SQRTSUB.
050223,001840:
050224,001841: 00,2444 12445 TCF +1 # IF THE ABOVE DOES NOT HOLD, FORCE SIGN
050225,001842: 00,2445 34753 CAF HALF # AGREEMENT IN NUMERATOR AND DENOMINATOR
050226,001843: 00,2446 60000 DOUBLE # TO FACILITATE OVERFLOW AND NEAR-ONE
050227,001844: 00,2447 60155 AD MPAC +1 # CHECKING.
050228,001845: 00,2450 54155 TS MPAC +1
050229,001846: 00,2451 34772 CAF ZERO
050230,001847: 00,2452 64750 AD POSMAX
050231,001848: 00,2453 26154 ADS MPAC
050232,001849:
050233,001850: 00,2454 34753 CAF HALF # SAME FOR BUF.
050234,001851: 00,2455 60000 DOUBLE
050235,001852: 00,2456 60131 AD BUF +1
050236,001853: 00,2457 54131 TS BUF +1
050237,001854: 00,2460 34772 CAF ZERO
050238,001855: 00,2461 64750 AD POSMAX
050239,001856: 00,2462 26130 ADS BUF
050240,001857:
050241,001858: 00,2463 40154 CS MPAC # CHECK MAGNITUDE OF SIGN-CORRECTED
050242,001859: 00,2464 60130 AD BUF # OPERANDS.
050243,001860: 00,2465 10000 CCS A
050244,001861: 00,2466 12504 TCF DVNORM # DIVIDE OK - WILL NOT BECOME MAXDV CASE.
050245,001862: 00,2467 00133 LBUF2 ADRES BUF2
050246,001863: 00,2470 12404 TCF DVOVF # DIVISOR NOT LESS THAN DIVIDEND - OVF.
050247,001864:
050248,001865: 00,2471 54140 TS MAXDVSW # IF THE MAJOR PARTS OF THE DIVIDEND AND
050249,001866: 00,2472 40155 CS MPAC +1 # DIVISOR ARE EQUAL, A SPECIAL APPROXIMA-
050250,001867: 00,2473 60131 AD BUF +1 # TION IS USED (PROVIDED THE DIVISION IS
050251,001868: 00,2474 00006 EXTEND # POSSIBLE, OF COURSE).
050252,001869: 00,2475 62404 BZMF DVOVF
050253,001870: 00,2476 12504 TCF DVNORM # IF NO OVERFLOW.
050254,001871:
|
|
Page 1165 |
050256,001873: 00,2477 00006 BUFNORM EXTEND # ADD -1 TO AUGMENT SHIFT COUNT AND SHIFT
050257,001874: 00,2500 24137 AUG DVNORMCT # LEFT ONE PLACE.
050258,001875: 00,2501 00006 EXTEND
050259,001876: 00,2502 30131 DCA BUF
050260,001877: 00,2503 20131 DAS BUF
050261,001878:
050262,001879: 00,2504 30130 DVNORM CA BUF # SEE IF DIVISOR NORMALIZED YET.
050263,001880: 00,2505 60000 DOUBLE
050264,001881: 00,2506 54000 OVSK
050265,001882: 00,2507 12477 TCF BUFNORM # NO - SHIFT LEFT ONE AND TRY AGAIN.
050266,001883:
050267,001884: 00,2510 52155 DXCH MPAC # CALL DIVIDEND NORMALIZATION SEQUENCE
050268,001885: 00,2511 50137 INDEX DVNORMCT # PRIOR TO DOING THE DIVIDE.
050269,001886: 00,2512 02564 TC MAXTEST
050270,001887:
050271,001888: 00,2513 54156 TS MPAC +2 # RETURNS WITH DIVISION DONE AND C(A) = 0.
050272,001889: 00,2514 16030 TCF DANZIG
050273,001890:
050274,001891: 00,2515 10000 BUFPOS CCS A
050275,001892: 00,2516 12421 TCF BUF+ # TO BUF+ IF BUF IS GREATER THAN +1.
050276,001893:
050277,001894: 00,2517 40131 CS BUF +1 # IF BUF IS +1, FORCING SIGN AGREEMENT
050278,001895: 00,2520 00006 EXTEND # MAY CAUSE BUF TO BECOME ZERO.
050279,001896: 00,2521 62421 BZMF BUF+ # BRANCH IF SIGNS AGREE.
050280,001897:
050281,001898: 00,2522 34753 CA HALF # SIGNS DISAGREE. FORCE AGREEMENT.
050282,001899: 00,2523 60000 +6 DOUBLE
050283,001900: 00,2524 26131 ADS BUF +1
050284,001901: 00,2525 34772 CA ZERO
050285,001902: 00,2526 54130 TS BUF
050286,001903: 00,2527 12363 TCF BUFZERO
050287,001904:
050288,001905: 00,2530 10000 BUFNEG CCS A
050289,001906: 00,2531 12415 TCF BUF- # TO BUF- IF BUF IS LESS THAN -1.
050290,001907:
050291,001908: 00,2532 30131 CA BUF +1 # IF BUF IS -1, FORCING SIGN AGREEMENT
050292,001909: 00,2533 00006 EXTEND # MAY CAUSE BUF TO BECOME ZERO.
050293,001910: 00,2534 62415 BZMF BUF- # BRANCH IF SIGNS AGREE.
050294,001911:
050295,001912: 00,2535 44753 CS HALF # SIGNS DISAGREE. FORCE AGREEMENT.
050296,001913: 00,2536 12523 TCF BUFPOS +6
050297,001914:
|
|
Page 1166 |
050299,001916: # THE FOLLOWING ARE PROLOGUES TO SHIFT THE DIVIDEND ARRIVING IN A AND L BEFORE THE DIVIDE.
050300,001917:
050301,001918: 00,2537 22021 -21D LXCH SR # SPECIAL PROLOGUE FOR UNIT WHEN THE
050302,001919: 00,2540 00006 EXTEND # LENGTH OF THE ARGUMENT WAS NOT LESS THAN
050303,001920: 00,2541 74753 MP HALF # .5. IN THIS CASE, EACH COMPONENT MUST BE
050304,001921: 00,2542 56001 XCH L # SHIFTED RIGHT ONE TO PRODUCE A HALF-UNIT
050305,001922: 00,2543 60021 AD SR # VECTOR.
050306,001923: 00,2544 56001 XCH L
050307,001924: 00,2545 12570 TCF GENDDV +1 # WITH DP DIVIDEND IN A,L.
050308,001925:
050309,001926: 00,2546 20001 DDOUBL # PROLOGUE WHICH NORMALIZES THE DIVIDEND
050310,001927: 00,2547 20001 DDOUBL # WHEN IT IS KNOWN THAT NO DIVISION
050311,001928: 00,2550 20001 DDOUBL # OVERFLOW WILL OCCUR.
050312,001929: 00,2551 20001 DDOUBL
050313,001930: 00,2552 20001 DDOUBL
050314,001931: 00,2553 20001 DDOUBL
050315,001932: 00,2554 20001 DDOUBL
050316,001933: 00,2555 20001 DDOUBL
050317,001934: 00,2556 20001 DDOUBL
050318,001935: 00,2557 20001 DDOUBL
050319,001936: 00,2560 20001 DDOUBL
050320,001937: 00,2561 20001 DDOUBL
050321,001938: 00,2562 20001 DDOUBL
050322,001939: 00,2563 52155 -1 DXCH MPAC
050323,001940: 00,2564 10140 MAXTEST CCS MAXDVSW # 0 IF MAJORS MIGHT BE =, -1 OTHERWISE.
050324,001941: 00,2565 06552 BIASHI DEC .4192 B-1 # SQRT CONSTANTS
050325,001942:
050326,001943: 00,2566 12641 TCF MAXDV # CHECK TO SEE IF THEY ARE NOW EQUAL.
050327,001944:
|
|
Page 1167 |
050329,001946: # THE FOLLOWING IS A GENERAL PURPOSE DOUBLE PRECISION DIVISION ROUTINE. IT DIVIDES MPAC BY BUF AND LEAVES
050330,001947: # THE RESULT IN MPAC. THE FOLLOWING CONDITIONS MUST BE SATISFIED:
050331,001948:
050332,001949: # 1. THE DIVISOR (BUF) MUST BE POSITIVE AND NOT LESS THAN .5.
050333,001950:
050334,001951: # 2. THE DIVIDEND (MPAC) MUST BE POSITIVE WITH THE MAJOR PART OF MPAC STRICTLY LESS THAN THAT OF BUF
050335,001952: # (A SPECIAL APPROXIMATION, MAXDV, IS USED WHEN THE MAJOR PARTS ARE EQUAL).
050336,001953:
050337,001954: # UNDERSTANDING THAT A/B = Q + S(R/B) WHERE S = 2(-14) AND Q AND R ARE QUOTIENT AND REMAINDER, RESPEC-
050338,001955: # TIVELY, THE FOLLOWING APPROXIMATION IS OBTAINED BY MULTIPLYING ABOVE AND BELOW BY C - SD AND NEGLECTING TERMS OF
050339,001956: # ORDER S-SQUARED (POSSIBLY INTRODUCING ERROR INTO THE LOW TWO BITS OF THE RESULT). SIGN AGREEMENT IS UNNECESSARY.
050340,001957:
050341,001958: # A + SB . (R - QD) A + SB
050342,001959: # ------ = Q + S(------) WHERE Q AND R ARE QUOTIENT AND REMAINDER OF ------ RESPECTIVELY.
050343,001960: # C + SD ( C ) C
050344,001961:
050345,001962: 00,2567 52155 GENDDV DXCH MPAC # WE NEED A AND B ONLY FOR FIRST DV.
050346,001963: 00,2570 00006 +1 EXTEND # (SPECIAL UNIT PROLOGUE ENTERS HERE).
050347,001964: 00,2571 10130 DV BUF # A NOW CONTAINS Q AND L, R.
050348,001965: 00,2572 52155 DXCH MPAC
050349,001966:
050350,001967: 00,2573 40154 CS MPAC # FORM DIVIDEND FOR MINOR PART OF RESULT.
050351,001968: 00,2574 00006 EXTEND
050352,001969: 00,2575 70131 MP BUF +1
050353,001970: 00,2576 60155 AD MPAC +1 # OVERFLOW AT THIS POINT IS POSITIVE SINCE
050354,001971: 00,2577 54000 OVSK # R IS POSITIVE IN EVERY CASE.
050355,001972: 00,2600 12605 TCF +5
050356,001973:
050357,001974: 00,2601 00006 EXTEND # OVERFLOW CAN BE REMOVED BY SUBTRACTING C
050358,001975: 00,2602 60130 SU BUF # (BUF) ONCE SINCE R IS ALWAYS LESS THAN C
050359,001976: 00,2603 24154 INCR MPAC # IN THIS CASE. INCR COMPENSATES SUBTRACT.
050360,001977: 00,2604 12607 TCF +DOWN # (SINCE C(A) IS STILL POSITIVE).
050361,001978:
050362,001979: 00,2605 00006 +5 EXTEND # C(A) CAN BE MADE LESS THAN C IN MAGNI-
050363,001980: 00,2606 62617 BZMF -UP # TUDE BY DIMINISHING IT BY C (SINCE C IS
050364,001981: # NOT LESS THAN .5) UNLESS C(A) = 0.
050365,001982:
|
|
Page 1168 |
050367,001984: 00,2607 00006 +DOWN EXTEND
050368,001985: 00,2610 60130 SU BUF # IF POSITIVE, REDUCE ONLY IF NECESSARY
050369,001986: 00,2611 00006 EXTEND # SINCE THE COMPENSATING INCR MIGHT CAUSE
050370,001987: 00,2612 12615 BZF +3 # OVERFLOW.
050371,001988: 00,2613 00006 EXTEND # DONT SUBTRACT UNLESS RESULT IS POSITIVE
050372,001989: 00,2614 62623 BZMF ENDMAXDV # OR ZERO.
050373,001990:
050374,001991: 00,2615 24154 +3 INCR MPAC # KEEP SUBTRACT HERE AND COMPENSATE.
050375,001992: 00,2616 12624 TCF FINALDV
050376,001993:
050377,001994: 00,2617 00006 -UP EXTEND # IF ZERO, SET MINOR PART OF RESULT TO
050378,001995: 00,2620 12627 BZF FINALDV +3 # ZERO.
050379,001996:
050380,001997: 00,2621 00006 EXTEND # IF NEGATIVE, ADD C TO A, SUBTRACTING ONE
050381,001998: 00,2622 26154 DIM MPAC # TO COMPENSATE. DIM IS OK HERE SINCE THE
050382,001999: 00,2623 60130 ENDMAXDV AD BUF # MAJOR PART NEVER GOES NEGATIVE.
050383,002000:
|
|
Page 1169 |
050385,002002: 00,2624 22007 FINALDV ZL # DO DV TO OBTAIN MINOR PART OF RESULT.
050386,002003: 00,2625 00006 EXTEND
050387,002004: 00,2626 10130 DV BUF
050388,002005: 00,2627 54155 +3 TS MPAC +1
050389,002006:
050390,002007: 00,2630 10136 CCS DVSIGN # LEAVE RESULT POSITIVE UNLESS C(DVSIGN)=
050391,002008: 00,2631 00002 TC Q # -0.
050392,002009: 00,2632 00002 TC Q
050393,002010: 00,2633 00002 TC Q
050394,002011:
050395,002012: 00,2634 00006 EXTEND
050396,002013: 00,2635 40155 DCS MPAC
050397,002014: 00,2636 52155 DXCH MPAC
050398,002015: 00,2637 34772 CAF ZERO # SO WE ALWAYS RETURN WITH C(A) = 0.
050399,002016: 00,2640 00002 TC Q
050400,002017:
|
|
Page 1170 |
050402,002019: # IF THE MAJOR PARTS OF THE DIVISOR AND DIVIDEND ARE EQUAL, BUT THE MINOR PARTS ARE SUCH THAT THE
050403,002020: # DIVIDEND IS STRICTLY LESS THAN THE DIVISOR IN MAGNITUDE, THE FOLLOWING APPROXIMATION IS USED. THE ASSUMPTIONS
050404,002021: # ARE THE SAME AS THE GENERAL ROUTINE WITH THE ADDITION THAT SIGN AGREEMENT IS NECESSARY (B, C, & D POSITIVE).
050405,002022:
050406,002023: # C + SB . (C + B - D)
050407,002024: # ------ = 37777 + S(---------)
050408,002025: # C + SD ( C )
050409,002026:
050410,002027: # THE DIVISION MAY BE PERFORMED IMMEDIATELY SINCE B IS STRICTLY LESS THAN D AND C IS NOT LESS THAN .5.
050411,002028:
050412,002029: 00,2641 40154 MAXDV CS MPAC # SEE IF MAXDV CASE STILL HOLDS AFTER
050413,002030: 00,2642 60130 AD BUF # NORMALIZATION.
050414,002031: 00,2643 00006 EXTEND
050415,002032: 00,2644 12646 BZF +2
050416,002033: 00,2645 12567 TCF GENDDV # MPAC NOW LESS THAN BUF - DIVIDE AS USUAL
050417,002034:
050418,002035: 00,2646 34750 +2 CAF POSMAX # SET MAJOR PART OF RESULT.
050419,002036: 00,2647 54154 TS MPAC
050420,002037:
050421,002038: 00,2650 40131 CS BUF +1 # FORM DIVIDEND OF MINOR PART OF RESULT.
050422,002039: 00,2651 60155 AD MPAC +1
050423,002040: 00,2652 12623 TCF ENDMAXDV # GO ADD C AND DO DIVIDE, ATTACHING SIGN
050424,002041: # BEFORE EXITING.
050425,002042:
|
|
Page 1171 |
050427,002044: # VECTOR DIVIDED BY SCALAR, V/SC, IS EXECUTED HERE. THE VECTOR IS NOW IN MPAC WITH SCALAR IN BUF.
050428,002045:
050429,002046: 00,2653 44770 V/SC2 CS ONE # INITIALIZE DIVIDEND NORMALIZATION COUNT
050430,002047: 00,2654 54137 TS DVNORMCT # AND DIVISION SIGN REGISTER.
050431,002048: 00,2655 54127 TS VBUF +5
050432,002049:
050433,002050: 00,2656 03007 TC VECAGREE # FORCE SIGN AGREEMENT IN VECTOR
050434,002051:
050435,002052: 00,2657 52131 DXCH BUF
050436,002053: 00,2660 07500 TC ALSIGNAG # SIGN AGREE BUF
050437,002054: 00,2661 52131 DXCH BUF
050438,002055: 00,2662 10130 CCS BUF # FORCE DIVISOR POSITIVE WITH MAJOR PART
050439,002056: 00,2663 12720 TCF /BUF+ # NON-ZERO (IF POSSIBLE).
050440,002057: 00,2664 12666 TCF +2
050441,002058: 00,2665 12714 TCF /BUF-
050442,002059:
050443,002060: 00,2666 56131 XCH BUF +1 # SHIFT VECTOR AND SCALAR LEFT 14.
050444,002061: 00,2667 56130 XCH BUF
050445,002062: 00,2670 56155 XCH MPAC +1
050446,002063: 00,2671 56154 XCH MPAC
050447,002064: 00,2672 00006 EXTEND # CHECK FOR OVERFLOW IN EACH CASE.
050448,002065: 00,2673 12675 BZF +2
050449,002066: 00,2674 12404 TCF DVOVF
050450,002067:
050451,002068: 00,2675 56160 XCH MPAC +4
050452,002069: 00,2676 56157 XCH MPAC +3
050453,002070: 00,2677 00006 EXTEND
050454,002071: 00,2700 12702 BZF +2
050455,002072: 00,2701 12404 TCF DVOVF
050456,002073:
050457,002074: 00,2702 56162 XCH MPAC +6
050458,002075: 00,2703 56161 XCH MPAC +5
050459,002076: 00,2704 00006 EXTEND
050460,002077: 00,2705 12707 BZF +2
050461,002078: 00,2706 12404 TCF DVOVF
050462,002079:
050463,002080: 00,2707 10130 CCS BUF
050464,002081: 00,2710 12720 TCF /BUF+
050465,002082: 00,2711 12404 TCF DVOVF # ZERO DIVISOR - OVERFLOW.
050466,002083: 00,2712 12714 TCF /BUF-
050467,002084: 00,2713 12404 TCF DVOVF
050468,002085:
050469,002086: 00,2714 00006 /BUF- EXTEND # ON NEGATIVE, COMPLEMENT BUF AND MAINTAIN
050470,002087: 00,2715 40131 DCS BUF # DVSIGN IN VBUF +5.
050471,002088: 00,2716 52131 DXCH BUF
050472,002089: 00,2717 24127 INCR VBUF +5
050473,002090:
|
|
Page 1172 |
050475,002092: 00,2720 00006 /BUF+ EXTEND
050476,002093: 00,2721 30131 DCA BUF # LEAVE ABS(ORIG DIVISOR) IN BUF2
050477,002094: 00,2722 52134 DXCH BUF2 # FOR OVERFLOW TESTING
050478,002095: 00,2723 12731 TCF /NORM # NORMALIZE DIVISOR IN BUF.
050479,002096:
050480,002097: 00,2724 00006 /NORM2 EXTEND # IF LESS THAN .5, AUGMENT DVNORMCT AND
050481,002098: 00,2725 24137 AUG DVNORMCT # DOUBLE DIVISOR.
050482,002099: 00,2726 00006 EXTEND
050483,002100: 00,2727 30131 DCA BUF
050484,002101: 00,2730 20131 DAS BUF
050485,002102:
050486,002103: 00,2731 30130 /NORM CA BUF # SEE IF DIVISOR NORMALIZED.
050487,002104: 00,2732 60000 DOUBLE
050488,002105: 00,2733 54000 OVSK
050489,002106: 00,2734 12724 TCF /NORM2 # DOUBLE AND TRY AGAIN IF NOT.
050490,002107:
050491,002108: 00,2735 02747 TC V/SCDV # DO X COMPONENT DIVIDE.
050492,002109: 00,2736 52160 DXCH MPAC +3 # SUPPLY ARGUMENTS IN USUAL SEQUENCE.
050493,002110: 00,2737 52155 DXCH MPAC
050494,002111: 00,2740 52160 DXCH MPAC +3
050495,002112:
050496,002113: 00,2741 02747 TC V/SCDV # Y COMPONENT.
050497,002114: 00,2742 52162 DXCH MPAC +5
050498,002115: 00,2743 52155 DXCH MPAC
050499,002116: 00,2744 52162 DXCH MPAC +5
050500,002117:
050501,002118: 00,2745 02747 TC V/SCDV # Z COMPONENT.
050502,002119: 00,2746 17354 TCF VROTATEX # GO RE-ARRANGE COMPONENTS BEFORE EXIT.
050503,002120:
|
|
Page 1173 |
050505,002122: # SUBROUTINE USED BY V/SC TO DIVIDE VECTOR COMPONENT IN MPAC,+1 BY THE SCALAR GIVEN IN BUF.
050506,002123:
050507,002124: 00,2747 30127 V/SCDV CA VBUF +5 # REFLECTS SIGN OF SCALAR.
050508,002125: 00,2750 54136 TS DVSIGN
050509,002126:
050510,002127: 00,2751 10154 CCS MPAC # FORCE MPAC POSITIVE, EXITING ON ZERO.
050511,002128: 00,2752 12766 TCF /MPAC+
050512,002129: 00,2753 12755 TCF +2
050513,002130: 00,2754 12762 TCF /MPAC-
050514,002131:
050515,002132: 00,2755 10155 CCS MPAC +1
050516,002133: 00,2756 12766 TCF /MPAC+
050517,002134: 00,2757 00002 TC Q
050518,002135: 00,2760 12762 TCF /MPAC-
050519,002136: 00,2761 00002 TC Q
050520,002137:
050521,002138: 00,2762 00006 /MPAC- EXTEND # USUAL COMPLEMENTING AND SETTING OF SIGN.
050522,002139: 00,2763 40155 DCS MPAC
050523,002140: 00,2764 52155 DXCH MPAC
050524,002141: 00,2765 24136 INCR DVSIGN
050525,002142:
050526,002143: 00,2766 44770 /MPAC+ CS ONE # INITIALIZE NEAR-ONE SWITCH.
050527,002144: 00,2767 54140 TS MAXDVSW
050528,002145:
050529,002146: 00,2770 40154 CS MPAC # CHECK POSSIBLE OVERFLOW.
050530,002147: 00,2771 60133 AD BUF2 # UNNORMALIZED INPUT DIVISOR.
050531,002148: 00,2772 10000 CCS A
050532,002149: 00,2773 13003 TCF DDVCALL # NOT NEAR-ONE
050533,002150: 00,2774 12776 TCF +2 # +0 IS JUST POSSIBLE
050534,002151: 00,2775 12404 TCF DVOVF # NO HOPE
050535,002152: 00,2776 54140 TS MAXDVSW # SIGNAL POSSIBLE NEAR-ONE CASE
050536,002153: 00,2777 40155 CS MPAC +1 # SEE IF DIVISION CAN BE DONE
050537,002154: 00,3000 60134 AD BUF2 +1
050538,002155: 00,3001 00006 EXTEND
050539,002156: 00,3002 62404 BZMF DVOVF
050540,002157:
050541,002158: 00,3003 52155 DDVCALL DXCH MPAC # CALL PRE-DIVIDE NORMALIZATION.
050542,002159: 00,3004 50137 INDEX DVNORMCT
050543,002160: 00,3005 12564 TCF MAXTEST
050544,002161:
|
|
Page 1174 |
050546,002163: 00,3006 32506 SLOPELO DEC .8324
050547,002164:
050548,002165: 00,3007 56002 VECAGREE XCH Q # SAVE Q IN A
050549,002166: 00,3010 52155 DXCH MPAC
050550,002167: 00,3011 07500 TC ALSIGNAG # SIGNAGREE MPAC
050551,002168: 00,3012 52155 DXCH MPAC
050552,002169: 00,3013 52160 DXCH MPAC +3
050553,002170: 00,3014 07500 TC ALSIGNAG # SIGN AGREE MPAC +3
050554,002171: 00,3015 52160 DXCH MPAC +3
050555,002172: 00,3016 52162 DXCH MPAC +5
050556,002173: 00,3017 07500 TC ALSIGNAG # SIGNAGREE MPAC +5
050557,002174: 00,3020 52162 DXCH MPAC +5
050558,002175: 00,3021 00000 TC A
050559,002176:
|
|
Page 1175 |
050561,002178: # THE FOLLOWING ROUTINE EXECUTES THE UNIT INSTRUCTION, WHICH TAKES THE UNIT OF THE VECTOR IN MPAC.
050562,002179:
050563,002180: 00,3022 03007 UNIT TC VECAGREE # FORCE SIGN AGREEMENT IN VECTOR
050564,002181: 00,3023 07466 TC MPACVBUF # SAVE ARGUMENT IN VBUF
050565,002182: 00,3024 34772 CAF ZERO # MUST SENSE OVERFLOW IN FOLLOWING DOT.
050566,002183: 00,3025 56121 XCH OVFIND
050567,002184: 00,3026 54141 TS TEM1
050568,002185: 00,3027 03316 TC VSQSUB # DOT MPAC WITH ITSELF.
050569,002186: 00,3030 30141 CA TEM1
050570,002187: 00,3031 56121 XCH OVFIND
050571,002188: 00,3032 00006 EXTEND
050572,002189: 00,3033 13035 BZF +2
050573,002190: 00,3034 12404 TCF DVOVF
050574,002191: 00,3035 00006 EXTEND
050575,002192: 00,3036 30155 DCA MPAC # LEAVE THE SQUARE OF THE LENGTH OF THE
050576,002193: 00,3037 50120 INDEX FIXLOC # ARGUMENT IN LVSQUARE.
050577,002194: 00,3040 52043 DXCH LVSQUARE
050578,002195:
050579,002196: 00,3041 03342 TC SQRTSUB # GO TAKE THE NORMALIZED SQUARE ROOT.
050580,002197:
050581,002198: 00,3042 10154 CCS MPAC # CHECK FOR UNIT OVERFLOW.
050582,002199: 00,3043 13050 TCF +5 # MPAC IS NOT LESS THAN .5 UNLESS
050583,002200: 00,3044 54001 TS L
050584,002201: 00,3045 50120 INDEX FIXLOC
050585,002202: 00,3046 52045 DXCH LV
050586,002203: 00,3047 12404 TCF DVOVF # INPUT TO SQRTSUB WAS 0.
050587,002204:
050588,002205: 00,3050 44321 CS FOURTEEN # SEE IF THE INPUT WAS SO SMALL THE THE
050589,002206: 00,3051 60135 AD MPTEMP # FIRST TWO REGISTERS OF THE SQUARE WERE 0
050590,002207: 00,3052 10000 CCS A
050591,002208: 00,3053 40000 COM # IF SO, SAVE THE NEGATIVE OF THE SHIFT
050592,002209: 00,3054 13132 TCF SMALL # COUNT -15D.
050593,002210:
050594,002211: 00,3055 13064 TCF LARGE # (THIS IS USUALLY THE CASE.)
050595,002212:
050596,002213: 00,3056 44776 CS THIRTEEN # IF THE SHIFT COUNT WAS EXACTLY 14, SET
050597,002214: 00,3057 54135 TS MPTEMP # THE PRE-DIVIDE NORM COUNT TO -13D.
050598,002215:
050599,002216: 00,3060 30154 CA MPAC # SHIFT THE LENGTH RIGHT 14 BEFORE STORING
050600,002217: 00,3061 54001 SMALL2 TS L # (SMALL EXITS TO THIS POINT).
050601,002218: 00,3062 34772 CAF ZERO
050602,002219: 00,3063 13111 TCF LARGE2 # GO TO STORE LENGTH AND PROCEED.
050603,002220:
050604,002221: 00,3064 10135 LARGE CCS MPTEMP # MOST ALL CASES COME HERE.
050605,002222: 00,3065 13073 TCF LARGE3 # SEE IF NO NORMALIZATION WAS REQUIRED BY
050606,002223: 00,3066 42024 CS SRDDV # SQRT, AND IF SO, SET UP FOR A SHIFT
050607,002224: 00,3067 54135 TS MPTEMP # RIGHT 1 BEFORE DIVIDING TO PRODUCE
050608,002225: 00,3070 00006 EXTEND # THE DESIRED HALF UNIT VECTOR.
050609,002226: 00,3071 30155 DCA MPAC
|
|
Page 1176 |
050611,002228: 00,3072 13111 TCF LARGE2
050612,002229:
|
|
Page 1177 |
050614,002231: 00,3073 40000 LARGE3 COM # LEAVE NEGATIVE OF SHIFT COUNT-1 FOR
050615,002232: 00,3074 54135 TS MPTEMP # PREDIVIDE LEFT SHIFT.
050616,002233:
050617,002234: 00,3075 40000 COM # PICK UP REQUIRED SHIFTING BIT TO UNNORM-
050618,002235: 00,3076 50000 INDEX A # ALIZE THE SQRT RESULT.
050619,002236: 00,3077 34753 CAF BIT14
050620,002237: 00,3100 54130 TS BUF
050621,002238: 00,3101 00006 EXTEND
050622,002239: 00,3102 70155 MP MPAC +1
050623,002240: 00,3103 56130 XCH BUF
050624,002241: 00,3104 00006 EXTEND # (UNNORMALIZE THE SQRT FOR LV).
050625,002242: 00,3105 70154 MP MPAC
050626,002243: 00,3106 56001 XCH L
050627,002244: 00,3107 60130 AD BUF
050628,002245: 00,3110 56001 XCH L
050629,002246:
050630,002247: 00,3111 50120 LARGE2 INDEX FIXLOC
050631,002248: 00,3112 52045 DXCH LV # LENGTH NOW STORED IN WORK AREA.
050632,002249:
050633,002250: 00,3113 44770 CS ONE
050634,002251: 00,3114 54140 TS MAXDVSW # NO MAXDV CASES IN UNIT.
050635,002252:
050636,002253: 00,3115 52123 DXCH VBUF # PREPARE X COMPONENT FOR DIVIDE, SETTING
050637,002254: 00,3116 52155 DXCH MPAC # LENGTH OF VECTOR AS DIVISOR IN BUF.
050638,002255: 00,3117 52131 DXCH BUF
050639,002256: 00,3120 03150 TC UNITDV
050640,002257:
050641,002258: 00,3121 52125 DXCH VBUF +2 # DO Y AND Z IN USUAL FASHION SO WE CAN
050642,002259: 00,3122 52155 DXCH MPAC # EXIT THROUGH VROTATEX.
050643,002260: 00,3123 52160 DXCH MPAC +3
050644,002261: 00,3124 03150 TC UNITDV
050645,002262:
050646,002263: 00,3125 52127 DXCH VBUF +4
050647,002264: 00,3126 52155 DXCH MPAC
050648,002265: 00,3127 52162 DXCH MPAC +5
050649,002266: 00,3130 03150 TC UNITDV
050650,002267: 00,3131 17354 TCF VROTATEX # AND EXIT.
050651,002268:
|
|
Page 1178 |
050653,002270: # IF THE LENGTH OF THE ARGUMENT VECTOR WAS LESS THAN 2(-28), EACH COMPONENT MUST BE SHIFTED LEFT AT LEAST
050654,002271: # 14 PLACES BEFORE THE DIVIDE. NOTE THAT IN THIS CASE, THE MAJOR PART OF EACH COMPONENT IS ZERO.
050655,002272:
050656,002273: 00,3132 54135 SMALL TS MPTEMP # NEGATIVE OF PRE-DIVIDE SHIFT COUNT.
050657,002274:
050658,002275: 00,3133 34772 CAF ZERO # SHIFT EACH COMPONENT LEFT 14.
050659,002276: 00,3134 56123 XCH VBUF +1
050660,002277: 00,3135 56122 XCH VBUF
050661,002278: 00,3136 56125 XCH VBUF +3
050662,002279: 00,3137 56124 XCH VBUF +2
050663,002280: 00,3140 56127 XCH VBUF +5
050664,002281: 00,3141 56126 XCH VBUF +4
050665,002282:
050666,002283: 00,3142 40135 CS MPTEMP
050667,002284: 00,3143 50000 INDEX A
050668,002285: 00,3144 34753 CAF BIT14
050669,002286: 00,3145 00006 EXTEND
050670,002287: 00,3146 70154 MP MPAC
050671,002288: 00,3147 13061 TCF SMALL2
050672,002289:
050673,002290: 00,3150 4776 THIRTEEN = OCT15
050674,002291: 00,3150 4321 FOURTEEN = OCT16
050675,002292: 00,3150 4321 OCT16 = R1D1
050676,002293:
|
|
Page 1179 |
050678,002295: # THE FOLLOWING ROUTINE SETS UP THE CALL TO THE DIVIDE ROUTINES.
050679,002296:
050680,002297: 00,3150 10154 UNITDV CCS MPAC # FORCE MPAC POSITIVE IF POSSIBLE, SETTING
050681,002298: 00,3151 13167 TCF UMPAC+ # DVSIGN ACCORDING TO THE SIGN OF MPAC
050682,002299: 00,3152 13154 TCF +2 # SINCE THE DIVISOR IS ALWAYS POSITIVE
050683,002300: 00,3153 13161 TCF UMPAC- # HERE.
050684,002301:
050685,002302: 00,3154 10155 CCS MPAC +1
050686,002303: 00,3155 13167 TCF UMPAC+
050687,002304: 00,3156 00002 TC Q # EXIT IMMEDIATELY ON ZERO.
050688,002305: 00,3157 13161 TCF UMPAC-
050689,002306: 00,3160 00002 TC Q
050690,002307:
050691,002308: 00,3161 44772 UMPAC- CS ZERO # IF NEGATIVE, SET -0 IN DVSIGN FOR FINAL
050692,002309: 00,3162 54136 TS DVSIGN # COMPLEMENT.
050693,002310: 00,3163 00006 EXTEND
050694,002311: 00,3164 40155 DCS MPAC # PICK UP ABSOLUTE VALUE OF ARG AND JUMP.
050695,002312: 00,3165 50135 INDEX MPTEMP
050696,002313: 00,3166 12563 TCF MAXTEST -1
050697,002314:
050698,002315: 00,3167 54136 UMPAC+ TS DVSIGN # SET DVSIGN FOR POSITIVE QUOTIENT.
050699,002316: 00,3170 52155 DXCH MPAC
050700,002317: 00,3171 50135 INDEX MPTEMP
050701,002318: 00,3172 12563 TCF MAXTEST -1
050702,002319:
|
|
Page 1180 |
050704,002321: # MISCELLANEOUS UNARY OPERATIONS.
050705,002322:
050706,002323: 00,3173 03277 DSQ TC DSQSUB # SQUARE THE DP CONTENTS OF MPAC.
050707,002324: 00,3174 16030 TCF DANZIG
050708,002325:
050709,002326: 00,3175 10163 ABVALABS CCS MODE # ABVAL OR ABS INSTRUCTION.
050710,002327: 00,3176 13225 TCF ABS # DO ABS ON SCALAR.
050711,002328: 00,3177 13225 TCF ABS
050712,002329:
050713,002330: 00,3200 03316 ABVAL TC VSQSUB # DOT MPAC WITH ITSELF.
050714,002331: 00,3201 22163 LXCH MODE # MODE IS NOW DP (L ZERO AFTER DAS).
050715,002332:
050716,002333: 00,3202 00006 EXTEND # STORE SQUARE OF LENGTH IN WORK AREA.
050717,002334: 00,3203 30155 DCA MPAC
050718,002335: 00,3204 50120 INDEX FIXLOC
050719,002336: 00,3205 52043 DXCH LVSQUARE
050720,002337:
|
|
Page 1181 |
050722,002339: # PROGRAM DESCRIPTION - SUBROUTINE SQRT
050723,002340: #
050724,002341: # FUNCTIONAL DESCRIPTION - DOUBLE PRECISION SQUARE ROOT ROUTINE
050725,002342: #
050726,002343: # THIS PROGRAM TAKES THE SQUARE ROOT OF THE 27 OR 28 MOST SIGNIFICANT BITS IN THE TRIPLE PRECISION SET OF
050727,002344: # NUMBERS - MPAC, MPAC+1, AND MPAC+2. THE ROOT IS RETURNED DOUBLE PRECISION IN MPAC AND MPAC+1.
050728,002345: #
050729,002346: # WARNING - THIS SUBROUTINE USES A TRIPLE PRECISION INPUT. THE PROGRAMMER MUST ASSURE THE CONTENTS OF MPAC+2
050730,002347: # ESPECIALLY IF THE CONTENTS OF MPAC IS SMALL OR ZERO. FOR DETAILS SEE STG MEMO NO.949.
050731,002348: #
050732,002349: # CALLING SEQUENCE - IN INTERPRETIVE MODE I.E., FOLLOWING TC INTPRET, SQRT NO ADDRESS IS ALLOWED
050733,002350:
050734,002351: # INPUT SCALING THE BINARY POINT IS ASSUMED TO THE RIGHT OF BIT 15. THE ANSWER IS RETURNED WITH THE SAME SCALING
050735,002352:
050736,002353: # SUBROUTINES - GENSCR, MPACSHR, SQRTSUB, ABORT
050737,002354:
050738,002355: # ABORT EXIT MODE - ABORTS ON NEGATIVE INPUT -1.2X10E-4 (77775 OCTAL) OR LESS.
050739,002356: # DISPLAYS ERROR CODE 1302
050740,002357: # TC ABORT
050741,002358: # OCT 1302
050742,002359:
050743,002360: # DEBRIS - LOCATIONS BUF, MPTEMP, ADDRWD ARE USED
050744,002361:
050745,002362: 00,3206 03342 SQRT TC SQRTSUB # TAKE THE SQUARE ROOT OF MPAC.
050746,002363: 00,3207 10135 CCS MPTEMP # RETURNED NORMALIZED SQUARE ROOT. SEE IF
050747,002364: 00,3210 13212 TCF +2 # ANY UN-NORMALIZATION REQUIRED AND EXIT
050748,002365: 00,3211 16030 TCF DANZIG # IF NOT.
050749,002366:
050750,002367: 00,3212 63726 AD NEG12 # A RIGHT SHIFT OF MORE THAN 13 COULD BE
050751,002368: 00,3213 00006 EXTEND # REQUIRED IF INPUT WAS ZERO IN MPAC,+1.
050752,002369: 00,3214 63220 BZMF SQRTSHFT # GOES HERE IN MOST CASES.
050753,002370: 00,3215 22007 ZL # IF A LONG SHIFT IS REQUIRED, GO TO
050754,002371: 00,3216 22116 LXCH ADDRWD # GENERAL RIGHT SHIFT ROUTINES.
050755,002372: 00,3217 12303 TCF GENSCR +4 # ADDRWD WAS ZERO TO PREVENT ROUND.
050756,002373:
050757,002374: 00,3220 50135 SQRTSHFT INDEX MPTEMP # SELECT SHIFTING BIT AND EXIT THROUGH
050758,002375: 00,3221 34752 CAF BIT15 # SHIFT ROUTINES.
050759,002376: 00,3222 54135 TS MPTEMP
050760,002377: 00,3223 34772 CAF ZERO # TO ZERO MPAC +2 IN THE PROCESS.
050761,002378: 00,3224 12036 TCF MPACSHR +3
050762,002379:
050763,002380: 00,3225 06672 ABS TC BRANCH # TEST SIGN OF MPAC AND COMPLEMENT IF
050764,002381: 00,3226 16030 TCF DANZIG
050765,002382: 00,3227 16030 TCF DANZIG
050766,002383: 00,3230 17624 TCF COMP
050767,002384:
|
|
Page 1182 |
050769,002386: 00,3231 44766 VDEF CS FOUR # VECTOR DEFINE - ESSENTIALLY TREATS
050770,002387: 00,3232 26166 ADS PUSHLOC # SCALAR IN MPAC AS X COMPONENT, PUSHES UP
050771,002388: 00,3233 00006 EXTEND # FOR Y AND THEN AGAIN FOR Z.
050772,002389: 00,3234 50000 INDEX A
050773,002390: 00,3235 30003 DCA 2
050774,002391: 00,3236 52160 DXCH MPAC +3
050775,002392: 00,3237 00006 EXTEND
050776,002393: 00,3240 50166 INDEX PUSHLOC
050777,002394: 00,3241 30001 DCA 0
050778,002395: 00,3242 52162 DXCH MPAC +5
050779,002396: 00,3243 16470 TCF VMODE # MODE IS NON VECTOR.
050780,002397:
050781,002398: 00,3244 03316 VSQ TC VSQSUB # DOT MPAC WITH ITSELF.
050782,002399: 00,3245 17266 TCF DMODE # MODE IS NOW DP.
050783,002400:
050784,002401: 00,3246 00006 PUSH EXTEND # PUSH DOWN MPAC LEAVING IT LOADED.
050785,002402: 00,3247 30155 DCA MPAC
050786,002403: 00,3250 50166 INDEX PUSHLOC # PUSH DOWN FIRST TWO REGISTERS IN EACH
050787,002404: 00,3251 52001 DXCH 0
050788,002405:
050789,002406: 00,3252 50163 INDEX MODE # INCREMENT PUSHDOWN POINTER.
050790,002407: 00,3253 36213 CAF NO.WDS
050791,002408: 00,3254 26166 ADS PUSHLOC
050792,002409:
050793,002410: 00,3255 10163 CCS MODE
050794,002411: 00,3256 13271 TCF TPUSH # PUSH DOWN MPAC +2.
050795,002412: 00,3257 16030 TCF DANZIG # DONE FOR DP.
050796,002413:
050797,002414: 00,3260 00006 EXTEND # ON VECTOR, PUSH DOWN Y AND Z COMPONENTS.
050798,002415: 00,3261 30160 DCA MPAC +3
050799,002416: 00,3262 50166 INDEX PUSHLOC
050800,002417: 00,3263 51775 DXCH 0 -4
050801,002418: 00,3264 00006 EXTEND
050802,002419: 00,3265 30162 DCA MPAC +5
050803,002420: 00,3266 50166 INDEX PUSHLOC
050804,002421: 00,3267 51777 DXCH 0 -2
050805,002422: 00,3270 16030 TCF DANZIG
050806,002423:
050807,002424: 00,3271 30156 TPUSH CA MPAC +2
050808,002425: 00,3272 16523 TCF ENDTPUSH +2
050809,002426:
050810,002427: 00,3273 50120 RVQ INDEX FIXLOC # RVQ - RETURN IVA QPRET.
050811,002428: 00,3274 30052 CA QPRET
050812,002429: 00,3275 54117 TS POLISH
050813,002430: 00,3276 16621 TCF GOTO +4 # (ASSUME QPRET POINTS TO FIXED ONLY.)
050814,002431:
|
|
Page 1183 |
050816,002433: # THE FOLLOWING SUBROUTINES ARE USED IN SQUARING MPAC, IN BOTH THE SCALAR AND VECTOR SENSE. THEY ARE
050817,002434: # SPECIAL CASES OF DMPSUB AND DOTSUB, PUT IN TO SAVE SOME TIME.
050818,002435:
050819,002436: 00,3277 30155 DSQSUB CA MPAC +1 # SQUARES THE SCALAR CONTENTS OF MPAC.
050820,002437: 00,3300 00006 EXTEND
050821,002438: 00,3301 70000 SQUARE
050822,002439: 00,3302 54156 TS MPAC +2
050823,002440: 00,3303 34772 CAF ZERO # FORM 2(CROSS TERM).
050824,002441: 00,3304 56155 XCH MPAC +1
050825,002442: 00,3305 00006 EXTEND
050826,002443: 00,3306 70154 MP MPAC
050827,002444: 00,3307 20001 DDOUBL # AND MAYBE OVERFLOW.
050828,002445: 00,3310 20156 DAS MPAC +1 # AND SET A TO NET OVERFLOW.
050829,002446: 00,3311 56154 XCH MPAC
050830,002447: 00,3312 00006 EXTEND
050831,002448: 00,3313 70000 SQUARE
050832,002449: 00,3314 20155 DAS MPAC
050833,002450: 00,3315 00002 TC Q
050834,002451:
050835,002452: 00,3316 00006 VSQSUB EXTEND # DOTS THE VECTOR IN MPAC WITH ITSELF.
050836,002453: 00,3317 22137 QXCH DOTRET
050837,002454: 00,3320 03277 TC DSQSUB # SQUARE THE X COMPONENT.
050838,002455: 00,3321 52160 DXCH MPAC +3
050839,002456: 00,3322 52155 DXCH MPAC
050840,002457: 00,3323 52131 DXCH BUF # SO WE CAN END IN DOTSUB.
050841,002458: 00,3324 30156 CA MPAC +2
050842,002459: 00,3325 54132 TS BUF +2
050843,002460:
050844,002461: 00,3326 03277 TC DSQSUB # SQUARE Y COMPONENT.
050845,002462: 00,3327 52156 DXCH MPAC +1
050846,002463: 00,3330 20132 DAS BUF +1
050847,002464: 00,3331 60154 AD MPAC
050848,002465: 00,3332 60130 AD BUF
050849,002466: 00,3333 54130 TS BUF
050850,002467: 00,3334 13336 TCF +2
050851,002468: 00,3335 54121 TS OVFIND # IF OVERFLOW.
050852,002469:
050853,002470: 00,3336 52162 DXCH MPAC +5
050854,002471: 00,3337 52155 DXCH MPAC
050855,002472: 00,3340 03277 TC DSQSUB # SQUARE Z COMPONENT.
050856,002473: 00,3341 17154 TCF ENDDOT # END AS IN DOTSUB.
050857,002474:
|
|
Page 1184 |
050859,002476: # DOUBLE PRECISION SQUARE ROOT ROUTINE. TAKE THE SQUARE ROOT OF THE TRIPLE PRECISION (MPAC +2 USED ONLY
050860,002477: # IN NORMALIZATION) CONTENTS OF MPAC AND LEAVE THE NORMALIZED RESULT IN MPAC (C(MPAC) GREATER THAN OR EQUAL TO
050861,002478: # .5). THE RIGHT SHIFT COUNT (TO UNNORMALIZE) IS LEFT IN MPTEMP.
050862,002479:
050863,002480: 00,3342 34772 SQRTSUB CAF ZERO # START BY ZEROING RIGHT SHIFT COUNT.
050864,002481: 00,3343 54135 TS MPTEMP
050865,002482:
050866,002483: 00,3344 10154 CCS MPAC # CHECK FOR POSITIVE ARGUMENT, SHIFTING
050867,002484: 00,3345 13403 TCF SMPAC+ # FIRST SIGNIFICANT MPAC REGISTER INTO
050868,002485: 00,3346 13350 TCF +2 # MPAC ITSELF.
050869,002486: 00,3347 13372 TCF SQRTNEG # SEE IF MAG OF ARGUMENT LESS THAN 10(-4).
050870,002487:
050871,002488: 00,3350 56156 XCH MPAC +2 # MPAC IS ZERO - SHIFT LEFT 14.
050872,002489: 00,3351 56155 XCH MPAC +1
050873,002490: 00,3352 54154 TS MPAC
050874,002491: 00,3353 34774 CAF SEVEN # AUGMENT RIGHT SHIFT COUNTER.
050875,002492: 00,3354 54135 TS MPTEMP
050876,002493:
050877,002494: 00,3355 10154 CCS MPAC # SEE IF MPAC NOW PNZ.
050878,002495: 00,3356 13403 TCF SMPAC+
050879,002496: 00,3357 13361 TCF +2
050880,002497: 00,3360 13375 TCF ZEROANS # NEGATIVE BUT LESS THAN 10(-4) IN MAG.
050881,002498:
050882,002499: 00,3361 56155 XCH MPAC +1 # ZERO - SHIFT LEFT 14 AGAIN.
050883,002500: 00,3362 54154 TS MPAC
050884,002501: 00,3363 34774 CAF SEVEN # AUGMENT RIGHT SHIFT COUNTER.
050885,002502: 00,3364 26135 ADS MPTEMP
050886,002503:
050887,002504: 00,3365 10154 CCS MPAC
050888,002505: 00,3366 13403 TCF SMPAC+
050889,002506: 00,3367 00002 TC Q # SQRT(0) = 0.
050890,002507: 00,3370 13375 TCF ZEROANS
050891,002508: 00,3371 13450 TCF FIXROOT # DO NOT LEAVE SQRTSUB WITH -0 IN MPAC.
050892,002509:
050893,002510: 00,3372 10000 SQRTNEG CCS A # ARGUMENT IS NEGATIVE, BUT SEE IF SIGN-
050894,002511: 00,3373 13401 TCF SQRTABRT # CORRECTED ARGUMENT IS LESS THAN 10(-4)
050895,002512:
050896,002513: 00,3374 10155 CCS MPAC +1 # IN MAGNITUDE. IF SO, CALL ANSWER ZERO.
050897,002514: 00,3375 34772 ZEROANS CAF ZERO # FORCE ANSWER TO ZERO HERE.
050898,002515: 00,3376 13450 TCF FIXROOT
050899,002516: 00,3377 13401 TCF SQRTABRT
050900,002517: 00,3400 13450 TCF FIXROOT
050901,002518:
050902,002519: 00,3401 05726 SQRTABRT TC POODOO
050903,002520: 00,3402 21302 OCT 21302 # SQRT CALLED WITH NEGATIVE ARGUMENT
050904,002521:
050905,002522: 00,3403 62443 SMPAC+ AD -1/2+2 # SEE IF ARGUMENT GREATER THAN OR EQUAL TO
050906,002523: 00,3404 00006 EXTEND # .5.
|
|
Page 1185 |
050908,002525: 00,3405 63453 BZMF SRTEST # IF SO, SEE IF LESS THAN .25.
050909,002526:
050910,002527: 00,3406 52155 DXCH MPAC # WE WILL TAKE THE SQUARE ROOT OF MPAC/2.
050911,002528: 00,3407 22021 LXCH SR # SHIFT RIGHT 1 AND GO TO THE SQRT ROUTINE
050912,002529: 00,3410 00006 EXTEND
050913,002530: 00,3411 74753 MP HALF
050914,002531: 00,3412 52155 DXCH MPAC
050915,002532: 00,3413 56021 XCH SR
050916,002533: 00,3414 26155 ADS MPAC +1 # GUARANTEED NO OVERFLOW.
050917,002534:
050918,002535: 00,3415 32314 ARGHI CAF SLOPEHI # ARGUMENT BETWEEN .25 AND .5. GET A
050919,002536: 00,3416 00006 EXTEND # LINEAR APPROXIMATION FOR THIS RANGE.
050920,002537: 00,3417 70154 MP MPAC
050921,002538: 00,3420 62565 AD BIASHI # X0/2 = (MPAC/2)(SLOPEHI) + BIASHI/2.
050922,002539:
050923,002540: 00,3421 54130 +4 TS BUF # X0/2 (ARGLO ENTERS HERE).
050924,002541: 00,3422 30154 CA MPAC # SINGLE-PRECISION THROUGHOUT.
050925,002542: 00,3423 00006 EXTEND
050926,002543: 00,3424 74753 MP HALF # .5(MPAC/2)
050927,002544: 00,3425 00006 EXTEND
050928,002545: 00,3426 10130 DV BUF # .5(MPAC/2)/(X0/2)
050929,002546: 00,3427 26130 ADS BUF # X1 = X0/2 + .5(MPAC/2)/(X0/2).
050930,002547:
050931,002548: 00,3430 00006 EXTEND
050932,002549: 00,3431 74753 MP HALF # FORM UP X1/2.
050933,002550: 00,3432 52155 DXCH MPAC # SAVE AND BRING OUT ARGUMENT.
050934,002551: 00,3433 00006 EXTEND # TAKE DP QUOTIENT WITH X1.
050935,002552: 00,3434 10130 DV BUF
050936,002553: 00,3435 54131 TS BUF +1 # SAVE MAJOR PART OF QUOTIENT.
050937,002554: 00,3436 34772 CAF ZERO # FORM MINOR PART OF QUOTIENT USING
050938,002555: 00,3437 56001 XCH L # (REMAINDER,0).
050939,002556: 00,3440 00006 EXTEND
050940,002557: 00,3441 10130 DV BUF
050941,002558: 00,3442 54001 TS L # IN PREPARATION FOR DAS.
050942,002559: 00,3443 30131 CA BUF +1
050943,002560: 00,3444 20155 DAS MPAC # X2 = X1/2 + (MPAC/2)X1
050944,002561:
050945,002562: 00,3445 00006 EXTEND # OVERFLOWS IF ARG. NEAR POSMAX.
050946,002563: 00,3446 13452 BZF TCQBNK00
050947,002564: 00,3447 34750 CAF POSMAX
050948,002565: 00,3450 54154 FIXROOT TS MPAC
050949,002566: 00,3451 54155 TS MPAC +1
050950,002567: 00,3452 00002 TCQBNK00 TC Q # RETURN TO CALLER TO UNNORMALIZE, ETC.
050951,002568:
|
|
Page 1186 |
050953,002570: 00,3453 64754 SRTEST AD QUARTER # ARGUMENT WAS LESS THAN .5, SEE IF LESS
050954,002571: 00,3454 00006 EXTEND # THAN .25.
050955,002572: 00,3455 63477 BZMF SQRTNORM # IF SO, BEGIN NORMALIZATION.
050956,002573:
050957,002574: 00,3456 52155 DXCH MPAC # IF BETWEEN .5 AND .25, SHIFT RIGHT 1 AND
050958,002575: 00,3457 22021 LXCH SR # START AT ARGLO.
050959,002576: 00,3460 00006 EXTEND
050960,002577: 00,3461 74753 MP HALF
050961,002578: 00,3462 52155 DXCH MPAC
050962,002579: 00,3463 56021 XCH SR
050963,002580: 00,3464 26155 ADS MPAC +1 # NO OVERFLOW.
050964,002581:
050965,002582: 00,3465 33006 ARGLO CAF SLOPELO # (NORMALIZED) ARGUMENT BETWEEN .125 AND
050966,002583: 00,3466 00006 EXTEND # .25
050967,002584: 00,3467 70154 MP MPAC
050968,002585: 00,3470 62270 AD BIASLO
050969,002586: 00,3471 13421 TCF ARGHI +4 # BEGIN SQUARE ROOT.
050970,002587:
050971,002588: 00,3472 00006 SQRTNM2 EXTEND # SHIFT LEFT 2 AND INCREMENT RIGHT SHIFT
050972,002589: 00,3473 30156 DCA MPAC +1 # COUNT (FOR TERMINAL UNNORMALIZATION).
050973,002590: 00,3474 20156 DAS MPAC +1
050974,002591: 00,3475 60154 AD MPAC
050975,002592: 00,3476 26154 ADS MPAC # (NO OVERFLOW).
050976,002593:
050977,002594: 00,3477 24135 SQRTNORM INCR MPTEMP # FIRST TIME THROUGH, JUST SHIFT LEFT 1
050978,002595: 00,3500 00006 EXTEND # (PUTS IN EFFECTIVE RIGHT SHIFT SINCE
050979,002596: 00,3501 30156 DCA MPAC +1 # WE WANT MPAC/2).
050980,002597: 00,3502 20156 DAS MPAC +1
050981,002598: 00,3503 60154 AD MPAC
050982,002599: 00,3504 26154 ADS MPAC # (AGAIN NO OVERFLOW).
050983,002600: 00,3505 60000 DOUBLE
050984,002601: 00,3506 54022 TS CYL
050985,002602:
050986,002603: 00,3507 10022 NORMTEST CCS CYL # SEE IF ARGUMENT NOW NORMALIZED AT
050987,002604: 00,3510 10022 CCS CYL # GREATER THAN .125.
050988,002605: 00,3511 13472 TCF SQRTNM2 # NO - SHIFT LEFT 2 MORE AND TRY AGAIN.
050989,002606: 00,3512 13415 TCF ARGHI # YES - NOW BETWEEN .5 AND .25.
050990,002607: 00,3513 13465 TCF ARGLO # ARGUMENT NOW BETWEEN .25 AND .125.
050991,002608:
|
|
Page 1187 |
050993,002610: # TRIGONOMETRIC FUNCTION PACKAGE.
050994,002611:
050995,002612: # THE FOLLOWING TRIGONOMETRIC FUNCTIONS ARE AVAILABLE AS INTERPRETIVE OPERATIONS:
050996,002613:
050997,002614: # 1. SIN COMPUTES (1/2)SINE(2 PI MPAC).
050998,002615: # 2. COS COMPUTES (1/2)COSINE(2 PI MPAC).
050999,002616:
051000,002617: # 3. ASIN COMPUTES (1/2PI)ARCSINE(2 MPAC).
051001,002618: # 4. ACOS COMPUTES (1/2PI)ARCCOSINE(2 MPAC).
051002,002619:
051003,002620: # SIN-ASIN AND COS-ACOS ARE MUTUALLY INVERSE, IE SIN(ASIN(X)) = X.
051004,002621:
051005,002622: 00,3514 06672 COSINE TC BRANCH # FINDS COSINE USING THE IDENTITY
051006,002623: 00,3515 13520 TCF +3 # COS(X) = SIN(PI/2 - ABS(X)).
051007,002624: 00,3516 13523 TCF PRESINE
051008,002625: 00,3517 13523 TCF PRESINE
051009,002626:
051010,002627: 00,3520 00006 +3 EXTEND
051011,002628: 00,3521 40155 DCS MPAC
051012,002629: 00,3522 52155 DXCH MPAC
051013,002630:
051014,002631: 00,3523 34754 PRESINE CAF QUARTER # PI/2 SCALED.
051015,002632: 00,3524 26154 ADS MPAC
051016,002633:
051017,002634: 00,3525 52155 SINE DXCH MPAC # DOUBLE ARGUMENT.
051018,002635: 00,3526 20001 +1 DDOUBL
051019,002636: 00,3527 54000 OVSK # SEE IF OVERFLOW PRESENT.
051020,002637: 00,3530 13533 TCF +3 # IF NOT, ARGUMENT OK AS IS.
051021,002638:
051022,002639: 00,3531 00006 EXTEND # IF SO, WE LOST (OR GAINED) PI, SO
051023,002640: 00,3532 40001 DCOM # COMPLEMENT MPAC USING THE IDENTITY
051024,002641: # SIN(X-(+)PI) = SIN(-X).
051025,002642: 00,3533 52155 +3 DXCH MPAC
051026,002643: 00,3534 30154 CA MPAC # SEE IF ARGUMENT GREATER THAN .5 IN
051027,002644: 00,3535 60000 DOUBLE # MAGNITUDE. IF SO, REDUCE IT TO LESS THAN
051028,002645: 00,3536 54001 TS L # .5 (+-PI/2 SCALED) AS FOLLOWS:
051029,002646: 00,3537 13550 TCF SN1
051030,002647:
051031,002648: 00,3540 50000 INDEX A # IF POSITIVE, FORM PI - X, IF NEGATIVE
051032,002649: 00,3541 34752 CAF NEG1/2 +1 # USE -PI -X.
051033,002650: 00,3542 60000 DOUBLE
051034,002651: 00,3543 00006 EXTEND
051035,002652: 00,3544 60154 SU MPAC # GUARANTEED NO OVERFLOW.
051036,002653: 00,3545 54154 TS MPAC
051037,002654: 00,3546 40155 CS MPAC +1
051038,002655: 00,3547 54155 TS MPAC +1
051039,002656:
|
|
Page 1188 |
051041,002658: 00,3550 00006 SN1 EXTEND # SET UP TO EVALUATE HASTINGS POLYNOMIAL
051042,002659: 00,3551 30155 DCA MPAC
051043,002660: 00,3552 52134 DXCH BUF2
051044,002661: 00,3553 03277 TC DSQSUB # SQUARE MPAC.
051045,002662:
051046,002663: 00,3554 07171 TC POLY # EVALUATE FOURTH ORDER POLYNOMIAL.
051047,002664: 00,3555 00003 DEC 3 B-14
051048,002665: 00,3556 14441 37325 2DEC +.3926990796
051049,002666: 00,3560 53250 60764 2DEC -.6459637111
051050,002667: 00,3562 12146 21276 2DEC +.318758717
051051,002668: 00,3564 75466 71471 2DEC -.074780249
051052,002669: 00,3566 00236 32757 2DEC +.009694988
051053,002670:
051054,002671: 00,3570 32467 CAF LBUF2 # MULTIPLY BY ARGUMENT AND SHIFT LEFT 2.
051055,002672: 00,3571 07055 TC DMPSUB -1
051056,002673:
051057,002674: 00,3572 00006 EXTEND
051058,002675: 00,3573 30156 DCA MPAC +1
051059,002676: 00,3574 20156 DAS MPAC +1
051060,002677: 00,3575 60154 AD MPAC
051061,002678: 00,3576 26154 ADS MPAC # NEITHER SHIFT OVERFLOWS.
051062,002679: 00,3577 00006 EXTEND
051063,002680: 00,3600 30156 DCA MPAC +1
051064,002681: 00,3601 20156 DAS MPAC +1
051065,002682: 00,3602 60154 AD MPAC
051066,002683: 00,3603 26154 ADS MPAC
051067,002684: 00,3604 16030 TCF DANZIG
051068,002685:
|
|
Page 1189 |
051070,002687: # ARCSIN/ARCCOS ROUTINE.
051071,002688:
051072,002689: 00,3605 33626 ARCSIN CAF LASINEX # COMPUTE ARCSIN BY USING THE IDENTITY
051073,002690: 00,3606 13610 TCF +2 # ARCSIN(X) = PI/2 - ARCCOS(X).
051074,002691:
051075,002692: 00,3607 33710 ARCCOS CAF LDANZIG # (EXITS IMMEDIATELY).
051076,002693: 00,3610 54136 TS ESCAPE
051077,002694: 00,3611 06672 TC BRANCH # TEST SIGN OF INPUT.
051078,002695: 00,3612 13622 TCF ACOSST # START IMMEDIATELY IF POSITIVE.
051079,002696: 00,3613 13724 TCF ACOSZERO # ARCCOS(0) = PI/2 = .25.
051080,002697: 00,3614 00006 EXTEND # IF NEGATIVE, USE THE IDENTITY
051081,002698: 00,3615 40155 DCS MPAC # ARCCOS(X) = PI - ARCCOS(-X), FORCING
051082,002699: 00,3616 52155 DXCH MPAC # ARGUMENT POSITIVE.
051083,002700: 00,3617 33727 CAF TCSUBTR # SET EXIT TO DO ABOVE BEFORE
051084,002701: 00,3620 56136 XCH ESCAPE # ARCSIN/ARCCOS CONSIDERATIONS.
051085,002702: 00,3621 54137 TS ESCAPE2
051086,002703:
051087,002704: 00,3622 44753 ACOSST CS HALF # TEST MAGNITUDE OF INPUT.
051088,002705: 00,3623 60154 AD MPAC
051089,002706: 00,3624 10000 CCS A
051090,002707: 00,3625 13716 TCF ACOSOVF # THIS IS PROBABLY AN OVERFLOW CASE.
051091,002708:
051092,002709: 00,3626 13704 LASINEX TCF ASINEX
051093,002710:
051094,002711: 00,3627 13637 TCF ACOSST2 # NO OVERFLOW - PROCEED.
051095,002712:
051096,002713: 00,3630 10155 CCS MPAC +1 # IF MAJOR PART IS .5, CALL ANSWER 0
051097,002714: 00,3631 34772 CAF ZERO # UNLESS MINOR PART NEGATIVE.
051098,002715: 00,3632 13634 TCF ACOS=0
051099,002716:
051100,002717: 00,3633 13637 TCF ACOSST2
051101,002718:
051102,002719: 00,3634 54155 ACOS=0 TS MPAC +1
051103,002720: 00,3635 54154 +1 TS MPAC
051104,002721: 00,3636 00136 TC ESCAPE
051105,002722:
051106,002723: 00,3637 00006 ACOSST2 EXTEND # NOW THAT ARGUMENT IS IN PROPER RANGE,
051107,002724: 00,3640 40155 DCS MPAC # BEGIN COMPUTATION. USE HASTINGS
051108,002725: 00,3641 64753 AD HALF # APPROXIMATION ARCCOS(X) = SQRT(1-X)P(X)
051109,002726: 00,3642 52155 DXCH MPAC # IN A SCALED VERSION WHERE P(X) IS A
051110,002727: 00,3643 52134 DXCH BUF2 # SEVENTH ORDER POLYNOMIAL.
051111,002728:
051112,002729: 00,3644 03342 TC SQRTSUB # RETURNS WITH NORMALIZED SQUARE ROOT.
051113,002730:
051114,002731: 00,3645 10135 CCS MPTEMP # SEE IF UN-NORMALIZATION REQUIRED.
051115,002732: 00,3646 13711 TCF ACOSSHR # IF SO.
051116,002733:
|
|
Page 1190 |
051118,002735: 00,3647 52155 ACOS3 DXCH MPAC # SET UP FOR POLYNOMIAL EVALUATION.
051119,002736: 00,3650 52134 DXCH BUF2
051120,002737: 00,3651 52155 DXCH MPAC
051121,002738:
051122,002739: 00,3652 07171 TC POLY
051123,002740: 00,3653 00006 DEC 6 B-14
051124,002741: 00,3654 13240 23630 2DEC +.353553385 # COEFFICIENTS ARE C 2(+I)/PISQRT(2) WHERE
051125,002742: 00,3656 74721 47775 2DEC* -.0483017006 B+1* # I
051126,002743: 00,3660 02440 20237 2DEC* +.0200273085 B+2* # WHERE C STANDS FOR ORIGINAL COEFFS.
051127,002744: 00,3662 75067 70742 2DEC* -.0112931863 B+3*
051128,002745: 00,3664 03436 26756 2DEC* +.00695311612 B+4*
051129,002746: 00,3666 74037 57640 2DEC* -.00384617957 B+5*
051130,002747: 00,3670 03046 07143 2DEC* +.001501297736 B+6*
051131,002748: 00,3672 76654 42244 2DEC* -.000284160334 B+7*
051132,002749:
051133,002750: 00,3674 32467 CAF LBUF2 # DO FINAL MULTIPLY AND GO TO ANY
051134,002751: 00,3675 07055 TC DMPSUB -1 # EPILOGUE SEQUENCES.
051135,002752: 00,3676 00136 TC ESCAPE
051136,002753:
051137,002754: 00,3677 00006 SUBTR EXTEND # EPILOGUE FOR NEGATIVE INPUTS TO ARCCOS.
051138,002755: 00,3700 40155 DCS MPAC
051139,002756: 00,3701 64753 AD HALF # FORMS PI - ARCCOS(-X) = ARCCOS(X).
051140,002757: 00,3702 52155 DXCH MPAC
051141,002758: 00,3703 00137 TC ESCAPE2 # GO TO POSSIBLE ARCSIN EPILOGUE.
051142,002759:
051143,002760: 00,3704 00006 ASINEX EXTEND
051144,002761: 00,3705 40155 DCS MPAC # ARCSIN EPILOGUE - GET ARCSIN(X)
051145,002762: 00,3706 64754 AD QUARTER # = PI/2 - ARCCOS(X).
051146,002763: 00,3707 52155 DXCH MPAC
051147,002764: 00,3710 16030 LDANZIG TCF DANZIG
051148,002765:
|
|
Page 1191 |
051150,002767: 00,3711 50000 ACOSSHR INDEX A # THE SHIFT RIGHT IS LESS THAN 14 SINCE
051151,002768: 00,3712 34753 CAF BIT14 # THE INPUT WAS NON-ZERO DP.
051152,002769: 00,3713 54135 TS MPTEMP
051153,002770: 00,3714 02073 TC VSHRRND # DP SHIFT RIGHT AND ROUND.
051154,002771: 00,3715 13647 TCF ACOS3 # PROCEED.
051155,002772:
051156,002773: 00,3716 00006 ACOSOVF EXTEND # IF MAJOR PART WAS ONLY 1 MORE THAN .5,
051157,002774: 00,3717 13634 BZF ACOS=0 # CALL ANSWER ZERO.
051158,002775:
051159,002776: 00,3720 05671 ACOSABRT TC ALARM # IF OVERFLOW, CALL ANSWER ZERO BUT
051160,002777: 00,3721 01301 OCT 1301 # SOUND AN ALARM.
051161,002778:
051162,002779: 00,3722 34772 CAF ZERO
051163,002780: 00,3723 13634 TCF ACOS=0
051164,002781:
051165,002782: 00,3724 34754 ACOSZERO CAF QUARTER # ACOS(0) = PI/2.
051166,002783: 00,3725 13635 TCF ACOS=0 +1 # SET MPAC AND EXIT VIA ESCAPE.
051167,002784:
051168,002785: 00,3726 77763 NEG12 DEC -12 B-14
051169,002786: 00,3727 13677 TCSUBTR TCF SUBTR
051170,002787:
|
|
Page 1192 |
051172,002789: # THE FOLLOWING INSTRUCTIONS ARE AVAILABLE FOR SETTING, MODIFYING, AND BRANCHING ON INDEX REGISTERS:
051173,002790:
051174,002791: # 1. AXT ADDRESS TO INDEX TRUE.
051175,002792: # 1. AXC ADDRESS TO INDEX COMPLEMENTED.
051176,002793: # 3. LXA LOAD INDEX FROM ERASABLE.
051177,002794: # 4. LXC LOAD INDEX COMPLEMENTED FROM ERASABLE.
051178,002795: # 5. SXA STORE INDEX IN ERASABLE.
051179,002796: # 6. XCHX EXCHANGE INDEX REIGSTER WITH ERASABLE.
051180,002797:
051181,002798: # 7. INCR INCREMENT INDEX REGISTER.
051182,002799: # 8. XAD ERASABLE ADD TO INDEX REGISTER.
051183,002800: # 9. XSU ERASABLE SUBTRACT FROM INDEX REGISTER.
051184,002801:
051185,002802: # 10. TIX BRANCH ON INDEX REGISTER AND DECREMENT.
051186,002803:
051187,002804: 01,2000 SETLOC INTERB1
051188,002805: 01,2000 BANK
051189,002806:
051190,002807: 01,2341 COUNT* $$/INTER
051191,002808: 01,2341 02436 AXT TC TAGSUB # SELECT APPROPRIATE INDEX REGISTER.
051192,002809: 01,2342 30117 CA POLISH
051193,002810: 01,2343 50130 XSTORE INDEX INDEXLOC # CONTAINS C(FIXLOC) OR C(FIXLOC)+1.
051194,002811: 01,2344 54046 TS X1
051195,002812: 01,2345 16030 TCF DANZIG
051196,002813:
051197,002814: 01,2346 02436 AXC TC TAGSUB
051198,002815: 01,2347 40117 CS POLISH
051199,002816: 01,2350 02343 TC XSTORE
051200,002817:
051201,002818: 01,2351 02424 LXA TC 15ADRERS # LOAD INDEX REGISTER FROM ERASABLE.
051202,002819: 01,2352 50117 INDEX POLISH
051203,002820: 01,2353 30000 CA 0
051204,002821: 01,2354 12343 TCF XSTORE
051205,002822:
051206,002823: 01,2355 02424 LXC TC 15ADRERS # LOAD NDX REG FROM ERASABLE COMPLEMENTED.
051207,002824: 01,2356 50117 INDEX POLISH
051208,002825: 01,2357 40000 CS 0
051209,002826: 01,2360 12343 TCF XSTORE
051210,002827:
051211,002828: 01,2361 02424 SXA TC 15ADRERS # STORE INDEX REGISTER IN ERASABLE.
051212,002829: 01,2362 50130 INDEX INDEXLOC
051213,002830: 01,2363 30046 CA X1
051214,002831: 01,2364 50117 MSTORE1 INDEX POLISH
051215,002832: 01,2365 54000 TS 0
051216,002833: 01,2366 16030 TCF DANZIG
051217,002834:
|
|
Page 1193 |
051219,002836: 01,2367 02424 XCHX TC 15ADRERS # EXCHANGE INDEX REGISTER WITH ERASABLE.
051220,002837: 01,2370 50117 INDEX POLISH
051221,002838: 01,2371 30000 CA 0
051222,002839: 01,2372 50130 INDEX INDEXLOC
051223,002840: 01,2373 56046 XCH X1
051224,002841: 01,2374 12364 TCF MSTORE1
051225,002842:
051226,002843: 01,2375 02424 XAD TC 15ADRERS # ADD ERASABLE TO INDEX REGISTER.
051227,002844: 01,2376 50117 INDEX POLISH
051228,002845: 01,2377 30000 CA 0
051229,002846: 01,2400 50130 XAD2 INDEX INDEXLOC
051230,002847: 01,2401 26046 ADS X1 # IGNORING OVERFLOWS.
051231,002848: 01,2402 16030 TCF DANZIG
051232,002849:
051233,002850: 01,2403 02436 INCR TC TAGSUB # INCREMENT INDEX REGISTER.
051234,002851: 01,2404 30117 CA POLISH
051235,002852: 01,2405 12400 TCF XAD2
051236,002853:
051237,002854: 01,2406 02424 XSU TC 15ADRERS # SUBTRACT ERASABLE FROM INDEX REGISTER.
051238,002855: 01,2407 50117 INDEX POLISH
051239,002856: 01,2410 40000 CS 0
051240,002857: 01,2411 12400 TCF XAD2
051241,002858:
051242,002859: 01,2412 02436 TIX TC TAGSUB # BRANCH AND DECREMENT ON INDEX.
051243,002860: 01,2413 50130 INDEX INDEXLOC
051244,002861: 01,2414 40050 CS S1
051245,002862: 01,2415 50130 INDEX INDEXLOC
051246,002863: 01,2416 60046 AD X1
051247,002864: 01,2417 00006 EXTEND # NO OPERATION IF DECREMENTED INDEX IS
051248,002865: 01,2420 66030 BZMF DANZIG # NEGATIVE OR ZERO.
051249,002866:
051250,002867: 01,2421 50130 DOTIXBR INDEX INDEXLOC
051251,002868: 01,2422 56046 XCH X1 # IGNORING OVERFLOWS.
051252,002869:
051253,002870: 01,2423 16615 TCF GOTO # DO THE BRANCH USING THE CADR IN POLISH.
051254,002871:
|
|
Page 1194 |
051256,002873: # SUBROUTINE TO CONVERT AN ERASABLE ADDRESS (11 BITS) TO AN EBANK SETTING AND SUBADDRESS.
051257,002874:
051258,002875: 01,2424 40117 15ADRERS CS POLISH
051259,002876: 01,2425 65006 AD DEC45
051260,002877: 01,2426 10000 CCS A # DOES THE ADDRESS POINT TO THE WORK AREA?
051261,002878: 01,2427 30120 CA FIXLOC # YES. ADD FIXLOC. EBANK OK AS IS.
051262,002879: 01,2430 12435 TCF +5
051263,002880:
051264,002881: 01,2431 35024 CA OCT1400 # NO. SET EBANK & MAKE UP SUBADDRESS.
051265,002882: 01,2432 56117 XCH POLISH
051266,002883: 01,2433 54003 TS EBANK
051267,002884: 01,2434 74350 MASK LOW8
051268,002885: 01,2435 26117 +5 ADS POLISH # FALL INTO TAGSUB, AND RETURN VIA Q.
051269,002886:
051270,002887: # SUBROUTINE WHICH SETS THE ADDRESS OF THE SPECIFIED INDEX IN INDEXLOC. (ACTUALLY, THE ADDRESS -38D.)
051271,002888:
051272,002889: 01,2436 30120 TAGSUB CA FIXLOC
051273,002890: 01,2437 54130 TS INDEXLOC
051274,002891:
051275,002892: 01,2440 10020 CCS CYR # BIT 15 SPECIFIES INDEX.
051276,002893: 01,2441 24130 INCR INDEXLOC # 0 MEANS USE X2.
051277,002894: 01,2442 00002 TC Q
051278,002895: 01,2443 00002 TC Q # 1 FOR X1.
051279,002896:
|
|
Page 1195 |
051281,002898: # MISCELLANEOUS OPERATION CODES WITH DIRECT ADDRESSES. INCLUDED HERE ARE:
051282,002899:
051283,002900: # 1. ITA STORE QPRET (RETURN ADDRESS) IN ERASABLE.
051284,002901: # 2. CALL CALL A SUBROUTINE, LEAVING RETURN IN QPRET.
051285,002902: # 3. RTB RETURN TO BASIC LANGUAGE AT THE GIVEN ADDRESS.
051286,002903: # 4. BHIZ BRANCH IF THE HIGHORDER OF MPAC IS ZERO (SINGLE PRECISION).
051287,002904: # 5. BOV BRANCH ON OVERFLOW.
051288,002905: # 6. GOTO SIMPLE SEQUENCE CHANGE.
051289,002906:
051290,002907: 01,2444 10020 RTB/BHIZ CCS CYR
051291,002908: 01,2445 30117 RTB CA POLISH
051292,002909: 01,2446 04641 TC SWCALL -1 # SO A "TC Q" FROM ROUTINE LEADS TO DANZIG
051293,002910:
051294,002911: 01,2447 10154 BHIZ CCS MPAC
051295,002912: 01,2450 16030 TCF DANZIG
051296,002913: 01,2451 16615 TCF GOTO
051297,002914: 01,2452 16030 TCF DANZIG
051298,002915: 01,2453 16615 TCF GOTO
051299,002916:
051300,002917: 01,2454 10121 BOV(B) CCS OVFIND # BRANCH ON OVERFLOW TO BASIC OR INTERP.
051301,002918: 01,2455 12457 TCF +2
051302,002919: 01,2456 16030 TCF DANZIG
051303,002920: 01,2457 54121 TS OVFIND
051304,002921: 01,2460 10020 CCS CYR
051305,002922: 01,2461 12445 TCF RTB # IF BASIC.
051306,002923: 01,2462 00360 B5TOB8 OCT 360
051307,002924: 01,2463 16615 TCF GOTO
051308,002925:
|
|
Page 1196 |
051310,002927: 01,2464 10020 BZE/GOTO CCS CYR # SEE WHICH OP-CODE IS DESIRED.
051311,002928: 01,2465 06672 TC BRANCH # DO BZE.
051312,002929: 01,2466 16030 TCF DANZIG
051313,002930: 01,2467 16615 TCF GOTO # DO GOTO.
051314,002931: 01,2470 16030 TCF DANZIG
051315,002932:
051316,002933: 01,2471 10020 BPL/BMN CCS CYR
051317,002934: 01,2472 12500 TCF BPL
051318,002935: 01,2473 12000 5B10 DEC 5 B+10 # SHIFTS OP CODE IN SWITCH INSTRUCTION ADR
051319,002936:
051320,002937: 01,2474 06672 TC BRANCH # DO BMN.
051321,002938: 01,2475 16030 TCF DANZIG
051322,002939: 01,2476 16030 TCF DANZIG
051323,002940: 01,2477 16615 TCF GOTO # ONLY IF NNZ.
051324,002941:
051325,002942: 01,2500 06672 BPL TC BRANCH
051326,002943: 01,2501 16615 TCF GOTO # IF POSITIVE OR ZERO.
051327,002944: 01,2502 16615 TCF GOTO
051328,002945: 01,2503 16030 TCF DANZIG
051329,002946:
051330,002947: 01,2504 10020 CALL/ITA CCS CYR
051331,002948: 01,2505 16607 TCF CALL
051332,002949:
051333,002950: 01,2506 12507 NOOP # CANNOT GET HERE
051334,002951: 01,2507 02424 TC 15ADRERS # STORE QPRET. (TAGSUB AFTER 15ADRERS IS
051335,002952: 01,2510 50120 INDEX FIXLOC # SLOW IN THIS CASE, BUT SAVES STORAGE.)
051336,002953: 01,2511 30052 CA QPRET
051337,002954: 01,2512 12364 TCF MSTORE1
051338,002955:
|
|
Page 1197 |
051340,002957: # THE FOLLOWING OPERATIONS ARE AVAILABLE FOR ALTERING AND TESTING INTERPRETIVE SWITCHES:
051341,002958: #
051342,002959: # 00 BONSET SET A SWITCH AND DO A GOTO IF IT WAS ON.
051343,002960: # 01 SETGO SET A SWITCH AND DO A GOTO.
051344,002961: # 02 BOFSET SET A SWITCH AND DO A GOTO IF IT WAS OFF
051345,002962: # 03 SET SET A SWITCH.
051346,002963:
051347,002964: # 04 BONINV INVERT A SWITCH AND BRANCH IF IT WAS ON.
051348,002965: # 05 INVGO INVERT A SWITCH AND DO A GOTO.
051349,002966: # 06 BOFINV INVERT A SWITCH AND BRANCH IF IT WAS OFF
051350,002967: # 07 INVERT INVERT A SWITCH.
051351,002968:
051352,002969: # 10 BONCLR CLEAR A SWITCH AND BRANCH IF IT WAS ON.
051353,002970: # 11 CLRGO CLEAR A SWITCH AND DO A GOTO.
051354,002971: # 12 BOFCLR CLEAR A SWITCH AND BRANCH IF IT WAS OFF.
051355,002972: # 13 CLEAR CLEAR A SWITCH.
051356,002973:
051357,002974: # 14 BON BRANCH IF A SWITCH WAS ON.
051358,002975: # 16 BOFF BRANCH IF A SWITCH WAS OFF.
051359,002976:
051360,002977: # THE ADDRESS SUPPLIED WITH THE SWITCH INSTRUCTION IS INTERPRETED AS FOLLOWS:
051361,002978:
051362,002979: # BITS 1-4 SWITCH BIT NUMBER (1-15).
051363,002980: # BITS 5-8 SWITCH OPERATION NUMBER.
051364,002981: # BITS 9- SWITCH WORD NUMBER (UP TO 64 SWITCH WORDS).
051365,002982:
051366,002983: # THE ADDRESS ITSELF IS MADE UP BY THE YUL SYSTEM ASSEMBLER. THE BRANCH INSTRUCTIONS REQUIRE TWO
051367,002984: # ADDRESSES, THE SECOND TAKEN AS THE DIRECT (OR INDIRECT IF IN ERASABLE) ADDRESS OF THE BRANCH.
051368,002985:
051369,002986: 01,2513 34777 SWITCHES CAF LOW4 # LEAVE THE SWITCH BIT IN SWBIT .
051370,002987: 01,2514 70117 MASK POLISH
051371,002988: 01,2515 50000 INDEX A
051372,002989: 01,2516 34752 CAF BIT15 # (NUMBER FROM LEFT TO RIGHT.)
051373,002990: 01,2517 54131 TS SWBIT
051374,002991:
051375,002992: 01,2520 34762 CAF BIT7 # LEAVE THE SWITCH NUMBER IN SWWORD.
051376,002993: 01,2521 00006 EXTEND
051377,002994: 01,2522 70117 MP POLISH
051378,002995: 01,2523 54130 TS SWWORD
051379,002996:
051380,002997: 01,2524 00004 INHINT # DURING SWITCH CHANGE SO RUPT CAN USE TOO
051381,002998: 01,2525 50000 INDEX A # LEAVE THE SWITCH WORD ITSELF IN L.
051382,002999: 01,2526 30074 CA STATE
051383,003000: 01,2527 54002 TS Q # Q WILL BE USED AS A CHANNEL.
|
|
Page 1198 |
051385,003002: 01,2530 34756 CAF BIT11
051386,003003: 01,2531 00006 EXTEND # DISPATCH SWITCH BIT OPERATION AS IN BITS
051387,003004: 01,2532 70117 MP POLISH # 7-8 OF POLISH.
051388,003005: 01,2533 73024 MASK B3TOB4 # GETS 4X2-BIT CODE.
051389,003006: 01,2534 50000 INDEX A
051390,003007: 01,2535 12536 TCF +1
051391,003008:
051392,003009: 01,2536 30131 +1 CA SWBIT # 00 - SET SWITCH IN QUESTION.
051393,003010: 01,2537 00006 EXTEND
051394,003011: 01,2540 04002 ROR QCHAN
051395,003012: 01,2541 12550 TCF SWSTORE
051396,003013:
051397,003014: 01,2542 30131 +5 CA SWBIT # 01 - INVERT SWITCH.
051398,003015: 01,2543 00006 EXTEND
051399,003016: 01,2544 06002 RXOR QCHAN
051400,003017: 01,2545 12550 TCF SWSTORE
051401,003018:
051402,003019: 01,2546 40131 +9D CS SWBIT # 10 - CLEAR.
051403,003020: 01,2547 70002 MASK Q
051404,003021: 01,2550 50130 SWSTORE INDEX SWWORD
051405,003022: 01,2551 54074 TS STATE # NEW SWITCH WORD.
|
|
Page 1199 |
051407,003024: 01,2552 00003 +13D RELINT # 11 - NOOP.
051408,003025: 01,2553 34754 CAF BIT13
051409,003026: 01,2554 00006 EXTEND # DISPATCH SEQUCE CHANGING OR BRANCHING
051410,003027: 01,2555 70117 MP POLISH # CODE.
051411,003028: 01,2556 73024 MASK B3TOB4
051412,003029: 01,2557 50000 INDEX A
051413,003030: 01,2560 12561 TCF +1 # ORIGINALLY STORED IN BITS 5-6.
051414,003031:
051415,003032: 01,2561 40002 +1 CS Q # 00 - BRANCH IF ON.
051416,003033: 01,2562 70131 TEST MASK SWBIT
051417,003034: 01,2563 10000 CCS A
051418,003035: 01,2564 12574 TCF SWSKIP
051419,003036:
051420,003037: 01,2565 16664 +5 TCF SWBRANCH # 01 - GO TO.
051421,003038:
051422,003039: 01,2566 12574 TCF SWSKIP # HERE ONLY ON BIT 15.
051423,003040:
051424,003041: 01,2567 12570 NOOP # CANNOT GET HERE
051425,003042: 01,2570 12571 NOOP # CANNOT GET HERE
051426,003043: 01,2571 30002 +9D CA Q # 10 - BRANCH IF OFF.
051427,003044: 01,2572 12562 TCF TEST
051428,003045:
051429,003046: 01,2573 12574 +11D NOOP # CANNOT GET HERE
051430,003047: 01,2574 24164 SWSKIP INCR LOC
051431,003048:
051432,003049: 01,2575 16030 +13D TCF DANZIG # 11 - NOOP.
051433,003050:
051434,003051: 01,2576 01,3024 B3TOB4 EQUALS COREINC
051435,003052: 01,2576 01,2513 SW/ EQUALS SWITCHES
End of include-file INTERPRETER.agc. Parent file is MAIN.agc