Source Code
![]() |
This is a reconstruction of the AGC program Luminary 99 Rev 0. It was the
second release of the Lunar Module flight software targeted for use in Apollo 11, after Luminary 97.
A bug (which had been around since at least Apollo 10, Luminary 69) was
discovered in Rev 0 shortly before the Apollo 11 flight, resulting in a last minute
revision into Rev 1,
which is what actually flew rather than
the Rev 0 presented here. A hardcopy of Rev 0 is known to
exist — it belonged to AGC developer Allan Klumpp for many years — but
unfortunately the Virtual AGC Project has not had access to that hardcopy. Thus
the code you see here had been reconstructed (we believe accurately) rather than
transcribed. The first step of the reconstruction of the Rev 0 source code was
the transcription of the Rev 1 source code from a hardcopy in the MIT Museum collection.
The process of reverting the active portion of the source code (i.e., other than
program comments) from Rev 1 to Rev 0 was very minor, consisting
only of moving the position of the STARTSB1 label in
FRESH START AND RESTART. Allan had previously given
us the checksums of the memory banks of Rev 0, and we have verified the program presented here has checksums
identical to all banks of Allan's listing of Rev 0. The notations on Allan's Rev 0 program listing read,
in part:GAP: ASSEMBLE REVISION 099 OF AGC PROGRAM LUMINARY BY NASA 2021112-051 A single program comment is known to differ between Rev 0 and Rev 1, but these are harder to reconstruct and verify than changes to the active portion of the code. Thus it is possible that there are additional differences between the program comments in Rev 0 and Rev 1 that are unknown to us, and therefore are not reflected in the code presented here. |
044519,000002: ## Copyright: Public domain.
044520,000003: ## Filename: INTERPRETER.agc
044521,000004: ## Purpose: Part of the reconstructed source code for LMY99 Rev 0,
044522,000005: ## otherwise known as Luminary Rev 99, the second release
044523,000006: ## of the Apollo Guidance Computer (AGC) software for Apollo 11.
044524,000007: ## It differs from LMY99 Rev 1 (the flown version) only in the
044525,000008: ## placement of a single label. The corrections shown here have
044526,000009: ## been verified to have the same bank checksums as AGC developer
044527,000010: ## Allan Klumpp's copy of Luminary Rev 99, and so are believed
044528,000011: ## to be accurate. This file is intended to be a faithful
044529,000012: ## recreation, except that the code format has been changed to
044530,000013: ## conform to the requirements of the yaYUL assembler rather than
044531,000014: ## the original YUL assembler.
044532,000015: ##
044533,000016: ## Assembler: yaYUL
044534,000017: ## Contact: Ron Burkey <info@sandroid.org>.
044535,000018: ## Website: www.ibiblio.org/apollo.
044536,000019: ## Pages: 1002-1094
044537,000020: ## Mod history: 2009-05-25 RSB Adapted from the corresponding
044538,000021: ## Luminary131 file, using page
044539,000022: ## images from Luminary 1A.
044540,000023: ## 2011-01-06 JL Fixed pseudo-label indentation.
044541,000024: ## 2011-05-08 JL Removed workarounds.
044542,000025: ## 2016-12-17 RSB Proofed text comments with octopus/ProoferComments
044543,000026: ## and corrected the errors found.
044544,000027: ## 2017-01-28 RSB Back-ported some comment-text fixes made in Luminary 69.
044545,000028: ## 2017-02-05 RSB Back-ported comment corrections
044546,000029: ## identified while proofing Artemis 072.
044547,000030: ## 2017-03-13 RSB Comment-text fixes noted in proofing Luminary 116.
044548,000031: ## 2017-03-15 RSB Comment-text fixes identified in 5-way
044549,000032: ## side-by-side diff of Luminary 69/99/116/131/210.
044550,000033: ## 2017-08-01 MAS Created from LMY99 Rev 1.
044551,000034:
044552,000035: ## This source code has been transcribed or otherwise adapted from
044553,000036: ## digitized images of a hardcopy from the MIT Museum. The digitization
044554,000037: ## was performed by Paul Fjeld, and arranged for by Deborah Douglas of
044555,000038: ## the Museum. Many thanks to both. The images (with suitable reduction
044556,000039: ## in storage size and consequent reduction in image quality as well) are
044557,000040: ## available online at www.ibiblio.org/apollo. If for some reason you
044558,000041: ## find that the images are illegible, contact me at info@sandroid.org
044559,000042: ## about getting access to the (much) higher-quality images which Paul
044560,000043: ## actually created.
044561,000044: ##
044562,000045: ## The code has been modified to match LMY99 Revision 0, otherwise
044563,000046: ## known as Luminary Revision 99, the Apollo 11 software release preceeding
044564,000047: ## the listing from which it was transcribed. It has been verified to
044565,000048: ## contain the same bank checksums as AGC developer Allan Klumpp's listing
044566,000049: ## of Luminary Revision 99 (for which we do not have scans).
044567,000050: ##
044568,000051: ## Notations on Allan Klumpp's listing read, in part:
044569,000052: ##
044570,000053: ## ASSEMBLE REVISION 099 OF AGC PROGRAM LUMINARY BY NASA 2021112-51
044571,000054:
![]() |
Page 1002 |
044573,000056: # SECTION 1 DISPATCHER
044574,000057:
044575,000058: # ENTRY TO THE INTERPRETER. INTPRET SETS LOC TO THE FIRST INSTRUCTION, BANKSET TO THE BBANK OF THE
044576,000059: # OBJECT INTERPRETIVE PROGRAM, AND INTBIT15 TO THE BIT15 CONTENTS OF FBANK. INTERPRETIVE PROGRAMS MAY BE IN
044577,000060: # VIRTUALLY ALL BANKS PRESENT UNDER ANY SUPER-BANK SETTING, WITH THE RESTRICTION THAT PROGRAMS IN HIGH BANKS
044578,000061: # (BIT15 OF FBANK = 1) DO NOT REFER TO LOWBANKS, AND VICE-VERSA. THE INTERPRETER DOES NOT SWITCH SUPERBANKS.
044579,000062: # E-BANK SWITCHING OCCURS WHENEVER GENERAL ERASABLE (100-3777) IS ADDRESSED.
044580,000063:
044581,000064: 6037 BLOCK 03
044582,000065:
044583,000066: 6037 COUNT* $$/INTER
044584,000067: 6037 00003 INTPRET RELINT
044585,000068: 6040 00006 EXTEND # SET LOC TO THE WORD FOLLOWING THE TC.
044586,000069: 6041 22164 QXCH LOC
044587,000070:
044588,000071: 6042 30006 +2 CA BBANK # INTERPRETIVE BRANCHES FINISH HERE.
044589,000072: 6043 54165 TS BANKSET
044590,000073: 6044 74735 MASK BIT15 # GET 15TH BIT FOR INDEXABLE ADDRESSES.
044591,000074: 6045 54115 TS INTBIT15
044592,000075:
044593,000076: 6046 54023 TS EDOP # MAKE SURE NO INSTRUCTIONS LEFT OVER
044594,000077:
044595,000078: 6047 16070 TCF NEWOPS # PICK UP OP CODE PAIR AND BEGIN.
044596,000079:
044597,000080: 6050 22006 INTRSM LXCH BBANK # RESUME SUSPENDED INTERPRETIVE JOB
044598,000081: 6051 16042 TCF INTPRET +3
044599,000082:
044600,000083: # DLOAD LOADS MPAC, MPAC +1, LEAVING ZERO IN MPAC +2.
044601,000084:
044602,000085: 6052 00006 DLOAD EXTEND
044603,000086: 6053 50116 INDEX ADDRWD
044604,000087: 6054 30001 DCA 0 # LOAD DP C(C(ADDRWD)) INTO MPAC,MPAX +1
044605,000088: 6055 52155 SLOAD2 DXCH MPAC
044606,000089: 6056 34755 CAF ZERO # ZERO MPAC +2
044607,000090:
![]() |
Page 1003 |
044609,000092: # AT THE END OF MOST INSTRUCTIONS, CONTROL IS GIVEN TO DANZIG TO DISPATCH THE NEXT OPERATION.
044610,000093:
044611,000094: 6057 54156 TS MPAC +2 # AND DECLARE DP MODE
044612,000095:
044613,000096: 6060 54163 NEWMODE TS MODE # PROLOGUE FOR MODE-CHANGING INSTRUCTIONS.
044614,000097:
044615,000098: 6061 30165 DANZIG CA BANKSET # SET BBANK BEFORE TESTING NEWJOB SO THAT
044616,000099: 6062 54006 TS BBANK # IT MAY BE SAVED DIRECTLY BY CHANJOB.
044617,000100:
044618,000101: 6063 10023 NOIBNKSW CCS EDOP # SEE IF AN ORDER CODE IS LEFT OVER FROM
044619,000102: 6064 16077 TCF OPJUMP # THE LAST PAIR RETRIEVED. IF SO, EXECUTE.
044620,000103: # EDOP IS SET TO ZERO ON ITS RE-EDITING.
044621,000104:
044622,000105: 6065 10067 CCS NEWJOB # SEE IF A JOB OF HIGHER PRIORITY IS
044623,000106: 6066 15126 TCF CHANG2 # PRESENT, AND IF SO, CHANGE JOBS.
044624,000107:
044625,000108: 6067 24164 INCR LOC # ADVANCE THE LOCATION COUNTER.
044626,000109:
044627,000110: # ITRACE (1) REFERS TO "NEWOPS".
044628,000111: 6070 50164 NEWOPS INDEX LOC # ENTRY TO BEGIN BY PICKING OP CODE PAIR.
044629,000112: 6071 30000 CA 0 # MAY BE AN OPCODE PAIR OR A STORE CODE.
044630,000113: 6072 10000 CCS A # TEST SIGN AND GET DABS(A).
044631,000114: 6073 16362 TCF DOSTORE # PROCESS STORE CODE.
044632,000115:
044633,000116: 6074 00177 LOW7 OCT 177
044634,000117:
044635,000118: 6075 54023 TS EDOP # OP CODE PAIR. LEAVE THE OTHER IN EDOP
044636,000119: 6076 76074 MASK LOW7 # WHERE CCS EDOP WILL HONOR IT NEXT.
044637,000120:
044638,000121: 6077 54020 OPJUMP TS CYR # LOWWD ENTERS HERE IF A RIGHT-HAND OP
044639,000122: 6100 10020 CCS CYR # CODE IS TO BE PROCESSED. TEST PREFICES.
044640,000123: 6101 16247 TCF OPJUMP2 # TEST SECOND PREFIX BIT.
044641,000124:
044642,000125: 6102 16743 TCF EXIT # +0 OP CODE IS EXIT.
044643,000126:
![]() |
Page 1004 |
044645,000128: # PROCESS ADDRESSES WHICH MAY BE DIRECT, INDEXED, OR REFERENCE THE PUSHDOWN LIST.
044646,000129:
044647,000130: 6103 74753 ADDRESS MASK BIT1 # SEE IF ADDRESS IS INDEXED. CYR CONTAINED
044648,000131: 6104 10000 CCS A # 400XX, SO BIT 1 IS NOW AS IT WAS IN CYR.
044649,000132: 6105 16146 TCF INDEX # FORM INDEXED ADDRESS.
044650,000133:
044651,000134: 6106 50164 DIRADRES INDEX LOC # LOOK AHEAD TO NEXT WORD TO SEE IF
044652,000135: 6107 40001 OCT40001 CS 1 # ADDRESS IS GIVEN.
044653,000136: 6110 10000 CCS A
044654,000137: 6111 16215 TCF PUSHUP # IF NOT.
044655,000138:
044656,000139: 6112 77773 NEG4 DEC -4 B-14
044657,000140:
044658,000141: 6113 24164 INCR LOC # IF SO, TO SHOW WE PICKED UP A WORD.
044659,000142: 6114 54116 TS ADDRWD
044660,000143:
![]() |
Page 1005 |
044662,000145: # FINAL DIGESTION OF DIRECT ADDRESSES OF OP CODES WITH 01 PREFIX IS DONE HERE. IN EACH CASE, THE
044663,000146: # REQUIRED 12-BIT SUB-ADDRESS IS LEFT IN ADDRWD, WITH ANY REQUIRED E OR F BANK SWITCHING DONE. ADDRESSES LESS
044664,000147: # 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.
044665,000148:
044666,000149: 6115 66251 AD -ENDVAC # SEE IF ADDRESS RELATIVE TO WORK AREA.
044667,000150: 6116 10000 CCS A
044668,000151: 6117 67743 AD -ENDERAS # IF NOT, SEE IF IN GENERAL ERASABLE.
044669,000152: 6120 16125 TCF IERASTST
044670,000153:
044671,000154: 6121 30120 NETZERO CA FIXLOC # IF SO, LEAVE THE MODIFIED ADDRESS IN
044672,000155: 6122 26116 ADS ADDRWD # ADDRWD AND DISPATCH.
044673,000156: 6123 50020 ITR15 INDEX CYR # THIS INDEX MAKES THE NEXT INSTRUCTION
044674,000157: 6124 76273 7 INDJUMP -1 # TCF INDJUMP + OP, EDITING CYR.
044675,000158:
044676,000159: 6125 00006 IERASTST EXTEND
044677,000160: 6126 66136 BZMF GEADDR # GO PROCESS GENERAL-ERASABLE ADDRESS.
044678,000161:
044679,000162: 6127 75012 MASK LOW10 # FIXED BANK ADDRESS. RESTORE AND ADD B15.
044680,000163: 6130 65012 AD LOW10 # SWITCH BANKS AND LEAVE SUBADDRESS IN
044681,000164: 6131 56116 XCH ADDRWD # ADDRWD FOR OPERAND RETRIEVAL. (THIS
044682,000165: 6132 60115 AD INTBIT15 # METHOD PRECLUDES USE OF THE LAST
044683,000166: 6133 54004 TS FBANK # LOCATION IN EACH FBANK.)
044684,000167: 6134 50020 ITR12 INDEX CYR
044685,000168: 6135 76273 7 INDJUMP -1
044686,000169:
044687,000170: 6136 74357 GEADDR MASK LOW8
044688,000171: 6137 65007 AD OCT1400
044689,000172: 6140 56116 XCH ADDRWD
044690,000173: 6141 54003 TS EBANK
044691,000174: 6142 50020 ITR10 INDEX CYR
044692,000175: 6143 76273 7 INDJUMP -1
044693,000176:
![]() |
Page 1006 |
044695,000178: # THE FOLLOWING ROUTINE PROCESSES INTERPRETIVE INDEXED ADDRESSES. AN INTERPRETER INDEX REGISTER MAY
044696,000179: # CONTAIN THE ADDRESS OF ANY ERASABLE REGISTER (0-42 BEING RELATIVE TO THE VAC AREA) OR ANY INTERPRETIVE PROGRAM
044697,000180: # BANK, OR ANY INTEGER IN THAT RANGE.
044698,000181:
044699,000182: 6144 37733 DODLOAD* CAF DLOAD* # STODL* COMES HERE TO PROCESS LOAD ADR.
044700,000183: 6145 54020 TS CYR # (STOVL* ENTERS HERE).
044701,000184:
044702,000185: 6146 30120 INDEX CA FIXLOC # SET UP INDEX LOCATION.
044703,000186: 6147 54130 TS INDEXLOC
044704,000187: 6150 24164 INCR LOC # (ADDRESS ALWAYS GIVEN).
044705,000188: 6151 50164 INDEX LOC
044706,000189: 6152 40000 CS 0
044707,000190: 6153 10000 CCS A # INDEX 2 IF ADDRESS STORED COMPLEMENTED.
044708,000191: 6154 24130 INCR INDEXLOC
044709,000192: 6155 16156 NOOP
044710,000193:
044711,000194: 6156 54116 TS ADDRWD # 14 BIT ADDRESS TO ADDRWD.
044712,000195: 6157 77742 MASK HIGH4 # IF ADDRESS GREATER THAN 2K, ADD INTBIT15
044713,000196: 6160 00006 EXTEND
044714,000197: 6161 16164 BZF INDEX2
044715,000198: 6162 30115 CA INTBIT15
044716,000199: 6163 26116 ADS ADDRWD
044717,000200:
044718,000201: 6164 50130 INDEX2 INDEX INDEXLOC
044719,000202: 6165 40046 CS X1
044720,000203: 6166 26116 ADS ADDRWD # DO AUGMENT, IGNORING AND CORRECTING OVF.
044721,000204:
044722,000205: 6167 77744 MASK HIGH9 # SEE IF ADDRESS IS IN WORK AREA.
044723,000206: 6170 00006 EXTEND
044724,000207: 6171 16204 BZF INDWORK
044725,000208: 6172 77742 MASK HIGH4 # SEE IF IN FIXED BANK.
044726,000209: 6173 00006 EXTEND
044727,000210: 6174 16206 BZF INDERASE
044728,000211:
044729,000212: 6175 30116 CA ADDRWD # IN FIXED - SWITCH BANKS AND CREATE
044730,000213: 6176 54004 TS FBANK # SUB-ADDRESS.
044731,000214: 6177 75012 MASK LOW10
044732,000215: 6200 64741 AD 2K
044733,000216: 6201 54116 TS ADDRWD
044734,000217: 6202 50020 ITR11 INDEX CYR
044735,000218: 6203 36273 3 INDJUMP -1
044736,000219:
044737,000220: 6204 30120 INDWORK CA FIXLOC # MAKE ADDRWD RELATIVE TO WORK AREA.
044738,000221: 6205 16212 TCF ITR13 -1
044739,000222:
044740,000223: 6206 35007 INDERASE CA OCT1400
044741,000224: 6207 56116 XCH ADDRWD
044742,000225: 6210 54003 TS EBANK
044743,000226: 6211 74357 MASK LOW8
044744,000227: 6212 26116 -1 ADS ADDRWD
![]() |
Page 1007 |
044746,000229: 6213 50020 ITR13 INDEX CYR
044747,000230: 6214 36273 3 INDJUMP -1
044748,000231:
![]() |
Page 1008 |
044750,000233: # PUSH-UP ROUTINES. WHEN NO OPERAND ADDRESS IS GIVEN, THE APPROPRIATE OPERAND IS TAKEN FROM THE PUSH-DOWN
044751,000234: # LIST. IN MOST CASES THE MODE OF THE RESULT (VECTOR OR SCALAR) OF THE LAST ARITHMETIC OPERATION PERFORMED
044752,000235: # IS THE SAME AS THE TYPE OF OPERAND DESIRED (ALL ADD/SUBTRACT ETC.). EXCEPTIONS TO THIS GENERAL RULE ARE LISTED
044753,000236: # BELOW (NOTE THAT IN EVERY CASE THE MODE REGISTER IS LEFT INTACT):
044754,000237:
044755,000238: # 1. VXSC AND V/SC WANT THE OPPOSITE TYPE OF OPERAND, E.G., IF THE LAST OPERATION YIELDED A VECTOR
044756,000239: # RESULT, VXSC WANTS A SCALAR.
044757,000240:
044758,000241: # 2. THE LOAD CODES SHOULD LOAD THE ACCUMULATOR INDEPENDENT OF THE RESULT OF THE LAST OPERATION. THIS
044759,000242: # INCLUDES VLOAD, DLOAD, TLOAD, PDDL, AND PDVL (NO PUSHUP WITH SLOAD).
044760,000243:
044761,000244: # 3. SOME ARITHMETIC OPERATIONS REQUIRE A STANDARD TYPE OF OPERAND REGARDLESS OF THE PREVIOUS OPERATION.
044762,000245: # THIS INCLUDES SIGN WANTING DP AND TAD REQUIRING TP.
044763,000246:
044764,000247: 6215 34360 PUSHUP CAF OCT23 # IF THE LOW 5 BITS OF CYR ARE LESS THAN
044765,000248: 6216 70020 MASK CYR # 20, THIS OP REQUIRES SPECIAL ATTENTION.
044766,000249: 6217 66222 AD -OCT10 # (NO -0).
044767,000250: 6220 10000 CCS A
044768,000251: 6221 16233 TCF REGUP # FOR ALL CODES GREATER THAN OCT 7.
044769,000252:
044770,000253: 6222 77767 -OCT10 OCT -10
044771,000254:
044772,000255: 6223 66112 AD NEG4 # WE NOW HAVE 7 - OP CODE (MOD4). SEE IF
044773,000256: 6224 10000 CCS A # THE OP CODE (MOD4) IS THREE (REVERSE).
044774,000257: 6225 50000 INDEX A # NO - THE MODE IS DEFINITE. PICK UP THE
044775,000258: 6226 46244 CS NO.WDS
044776,000259: 6227 16235 TCF REGUP +2
044777,000260:
044778,000261: 6230 50163 INDEX MODE # FOR VXSC AND V/SC WE WANT THE REQUIRED
044779,000262: 6231 46242 CS REVCNT # PUSHLOC DECREMENT WITHOUT CHANGING THE
044780,000263: 6232 16235 TCF REGUP +2 # MODE AT THIS TIME.
044781,000264:
044782,000265: 6233 50163 REGUP INDEX MODE # MOST ALL OP CODES PUSHUP HERE.
044783,000266: 6234 46244 CS NO.WDS
044784,000267: 6235 26166 +2 ADS PUSHLOC
044785,000268: 6236 54116 TS ADDRWD
044786,000269: 6237 50020 ITR14 INDEX CYR
044787,000270: 6240 76273 7 INDJUMP -1 # (THE INDEX MAKES THIS A TCF.)
044788,000271:
044789,000272: 6241 00002 OCT 2 # REVERSE PUSHUP DECREMENT. VECTOR TAKES 2
044790,000273: 6242 00006 REVCNT OCT 6 # WORDS, SCALAR TAKES 6.
044791,000274: 6243 00006 OCT 6
044792,000275: 6244 00002 NO.WDS OCT 2 # CONVENTIONAL DECREMENT IS 6 WORDS VECTOR
044793,000276: 6245 00003 OCTAL3 OCT 3 # 2 IN DP, AND 3 IN TP.
044794,000277: 6246 00006 OCT 6
044795,000278:
![]() |
Page 1009 |
044797,000280: # TEST THE SECOND PREFIX BIT TO SEE IF THIS IS A MISCELLANEOUS OR A UNARY/SHORT SHIFT OPERATION.
044798,000281:
044799,000282: 6247 10020 OPJUMP2 CCS CYR # TEST SECOND PREFIX BIT.
044800,000283: 6250 16263 TCF OPJUMP3 # TEST THIRD BIT TO SEE IF UNARY OR SHIFT.
044801,000284:
044802,000285: 6251 77722 -ENDVAC DEC -45 B-14
044803,000286:
044804,000287: # THE FOLLOWING ROUTINE PROCESSES ADDRESSES OF SUFFIX CLASS 10. THEY ARE BASICALLY WORK AREA ADDRESSES
044805,000288: # IN THE RANGE 0-52, ERASABLE ECADR CONSTANTS FROM 100-3777, AND FCADRS ABOVE THAT. ALL 15 BITS ARE AVAILABLE
044806,000289: # IN CONTRAST TO SUFFIX 1, IN WHICH ONLY THE LOW ORDER 14 ARE AVAILABLE.
044807,000290:
044808,000291: 6252 24164 15BITADR INCR LOC # (ENTRY HERE FROM STCALL).
044809,000292: 6253 50164 INDEX LOC # PICK UP ADDRESS WORD.
044810,000293: 6254 30000 CA 0
044811,000294: 6255 54117 TS POLISH # WE MAY NEED A SUBADDRESS LATER.
044812,000295:
044813,000296: 6256 35013 CAF LOW7+2K # THESE INSTRUCTIONS ARE IN BANK 1.
044814,000297: 6257 54004 TS FBANK
044815,000298: 6260 70020 MASK CYR
044816,000299: 6261 50000 ITR7 INDEX A
044817,000300: 6262 16334 TCF MISCJUMP
044818,000301:
![]() |
Page 1010 |
044820,000303: # COMPLETE THE DISPATCHING OF UNARY AND SHORT SHIFT OPERATIONS.
044821,000304:
044822,000305: 6263 54004 OPJUMP3 TS FBANK # CALL IN BANK 0 (BITS 11-15 OF A ARE 0.)
044823,000306: # ITRACE (6) REFERS TO "OPJUMP3".
044824,000307: 6264 10020 CCS CYR # TEST THIRD PREFIX BIT.
044825,000308: 6265 50000 INDEX A # THE DECREMENTED UNARY CODE IS IN BITS
044826,000309: 6266 12000 TCF UNAJUMP # 1-4 OF A (ZERO, EXIT, HAS BEEN DETECTED)
044827,000310:
044828,000311: 6267 10163 CCS MODE # ITS A SHORT SHIFT CODE. SEE IF PRESENT
044829,000312: 6270 12017 TCF SHORTT # SCALAR OR VECTOR.
044830,000313: 6271 12017 TCF SHORTT
044831,000314: 6272 12121 TCF SHORTV # CALLS THE APPROPRIATE ROUTINE.
044832,000315:
044833,000316: 6273 4350 FBANKMSK EQUALS BANKMASK
044834,000317: 6273 00122 LVBUF ADRES VBUF
044835,000318:
![]() |
Page 1011 |
044837,000320: # THE FOLLOWING IS THE JUMP TABLE FOR OP CODES WHICH MAY HAVE INDEXABLE ADDRESSES OR MAY PUSH UP.
044838,000321:
044839,000322: 6274 16505 INDJUMP TCF VLOAD # 00 - LOAD MPAC WITH A VECTOR.
044840,000323: 6275 17071 TCF TAD # 01 - TRIPLE PRECISION ADD TO MPAC.
044841,000324: 6276 17655 TCF SIGN # 02 - COMPLEMENT MPAC (V OR SC) IF X NEG.
044842,000325: 6277 17401 TCF VXSC # 03 - VECTOR TIMES SCALAR.
044843,000326: 6300 16703 TCF CGOTO # 04 - COMPUTED GO TO.
044844,000327: 6301 16470 TCF TLOAD # 05 - LOAD MPAC WITH TRIPLE PRECISION.
044845,000328: 6302 16052 TCF DLOAD # 06 - LOAD MPAC WITH A DP SCALAR.
044846,000329: 6303 17624 TCF V/SC # 07 - VECTOR DIVIDED BY SCALAR.
044847,000330:
044848,000331: 6304 16501 TCF SLOAD # 10 - LOAD MPACIN SINGLE PRECISION.
044849,000332: 6305 16620 TCF SSP # 11 - SET SINGLE PRECISION INTO X.
044850,000333: 6306 16523 TCF PDDL # 12 - PUSH DOWN MPAC AND RE-LOAD IN DP.
044851,000334: 6307 17334 TCF MXV # 13 - MATRIX POST-MULTIPLIED BY VECTOR.
044852,000335: 6310 16557 TCF PDVL # 14 - PUSH DOWN AND VECTOR LOAD.
044853,000336: 6311 16626 TCF CCALL # 15 - COMPUTED CALL.
044854,000337: 6312 17337 TCF VXM # 16 - MATRIX PRE-MULTIPLIED BY VECTOR.
044855,000338: 6313 17616 TCF TSLC # 17 - NORMALIZE MPAC (SCALAR ONLY).
044856,000339:
044857,000340: 6314 17574 TCF DMPR # 20 - DP MULTIPLY AND ROUND.
044858,000341: 6315 17577 TCF DDV # 21 - DP DIVIDE BY.
044859,000342: 6316 17603 TCF BDDV # 22 - DP DIVIDE INTO.
044860,000343: 6317 17621 TCF GSHIFT # 23 - GENERAL SHIFT INSTRUCTION.
044861,000344: 6320 16751 TCF VAD # 24 - VECTOR ADD.
044862,000345: 6321 16747 TCF VSU # 25 - VECTOR SUBTRACT.
044863,000346: 6322 17036 TCF BVSU # 26 - VECTOR SUBTRACT FROM.
044864,000347: 6323 17331 TCF DOT # 27 - VECTOR DOT PRODUCT.
044865,000348:
044866,000349: 6324 17460 TCF VXV # 30 - VECTOR CROSS PRODUCT.
044867,000350: 6325 17425 TCF VPROJ # 31 - VECTOR PROJECTION.
044868,000351: 6326 17005 TCF DSU # 32 - DP SUBTRACT.
044869,000352: 6327 17062 TCF BDSU # 33 - DP SUBTRACT FROM.
044870,000353: 6330 16775 TCF DAD # 34 - DP ADD.
044871,000354: 6331 16331 TCF +0 # 35 - AVAILABLE
044872,000355: 6332 17572 TCF DMP1 # 36 - DP MULTIPLY.
044873,000356: 6333 17613 TCF SETPD # 37 - SET PUSH DOWN POINTER (DIRECT ONLY)
044874,000357:
044875,000358: # CODES 10 AND 14 MUST NOT PUSH UP. CODE 04 MAY BE USED FOR VECTOR DECLARE BEFORE PUSHUP IF DESIRED.
044876,000359:
![]() |
Page 1012 |
044878,000361: # THE FOLLOWING JUMP TABLE APPLIES TO INDEX, BRANCH, AND MISCELLANEOUS INSTRUCTIONS.
044879,000362:
044880,000363: 6334 12344 MISCJUMP TCF AXT # 00 - ADDRESS TO INDEX TRUE.
044881,000364: 6335 12351 TCF AXC # 01 - ADDRESS TO INDEX COMPLEMENTED.
044882,000365: 6336 12354 TCF LXA # 02 - LOAD INDEX FROM ERASABLE.
044883,000366: 6337 12360 TCF LXC # 03 - LOAD INDEX FROM COMPLEMENT OF ERAS.
044884,000367: 6340 12364 TCF SXA # 04 - STORE INDEX IN ERASABLE.
044885,000368: 6341 12372 TCF XCHX # 05 - EXCHANGE INDEX WITH ERASABLE.
044886,000369: 6342 12406 TCF INCR # 06 - INCREMENT INDEX REGISTER.
044887,000370: 6343 12415 TCF TIX # 07 - TRANSFER ON INDEX.
044888,000371:
044889,000372: 6344 12400 TCF XAD # 10 - INDEX REGISTER ADD FROM ERASABLE.
044890,000373: 6345 12411 TCF XSU # 11 - INDEX SUBTRACT FROM ERASABLE.
044891,000374: 6346 12467 TCF BZE/GOTO # 12 - BRANCH ZERO AND GOTO.
044892,000375: 6347 12474 TCF BPL/BMN # 13 - BRANCH PLUS AND BRANCH MINUS.
044893,000376: 6350 12447 TCF RTB/BHIZ # 14 - RETURN TO BASIC AND BRANCH HI ZERO.
044894,000377: 6351 12507 TCF CALL/ITA # 15 - CALL AND STORE QPRET.
044895,000378: 6352 12516 TCF SW/ # 16 - SWITCH INSTRUCTIONS AND AVAILABLE.
044896,000379: 6353 12457 TCF BOV(B) # 17 - BRANCH ON OVERFLOW TO BASIC OR INT.
044897,000380:
![]() |
Page 1013 |
044899,000382: # THE FOLLOWING JUMP TABLE APPIES TO UNARY INSTRUCTIONS.
044900,000383:
044901,000384: 6354 COUNT* $$/INTER
044902,000385: 00,2000 BANK 0 # 00 - EXIT - DETECTED EARLIER.
044903,000386: 00,2000 13207 UNAJUMP TCF SQRT # 01 - SQUARE ROOT.
044904,000387: 00,2001 13530 TCF SINE # 02 - SIN.
044905,000388: 00,2002 13517 TCF COSINE # 03 - COS.
044906,000389: 00,2003 13610 TCF ARCSIN # 04 - ARC SIN.
044907,000390: 00,2004 13612 TCF ARCCOS # 05 - ARC COS.
044908,000391: 00,2005 13174 TCF DSQ # 06 - DP SQUARE.
044909,000392: 00,2006 12116 TCF ROUND # 07 - ROUND TO DP.
044910,000393:
044911,000394: 00,2007 17670 TCF COMP # 10 - COMPLEMENT VECTOR OR SCALAR.
044912,000395: 00,2010 13232 TCF VDEF # 11 - VECTOR DEFINE.
044913,000396: 00,2011 13023 TCF UNIT # 12 - UNIT VECTOR.
044914,000397: 00,2012 13176 TCF ABVALABS # 13 - LENGTH OF VECTOR OR MAG OF SCALAR.
044915,000398: 00,2013 13245 TCF VSQ # 14 - SQUARE OF LENGTH OF VECTOR.
044916,000399: 00,2014 16354 TCF STADR # 15 - PUSH UP ON STORE CODE.
044917,000400: 00,2015 13274 TCF RVQ # 16 - RETURN VIA QPRET.
044918,000401: 00,2016 13247 TCF PUSH # 17 - PUSH MPAC DOWN.
044919,000402:
![]() |
Page 1014 |
044921,000404: # SECTION 2 LOAD AND STORE PACKAGE.
044922,000405:
044923,000406: # A SET OF EIGHT STORE CODES IS PROVIDED AS THE PRIMARY METHOD OF STORING THE MULTI-PURPOSE
044924,000407: # ACCUMULATOR (MPAC). IF IN THE DANZIG SECTION LOC REFERS TO AN ALGEBRAICALLY POSITIVE WORD, IT IS TAKEN AS A
044925,000408: # STORE CODE WITH A CORRESPONDING ERASABLE ADDRESS. MOST OF THESE CODES ARE TWO ADDRESS, SPECIFYING THAT THE WORD
044926,000409: # FOLLOWING THE STORE CODE IS TO BE USED AS AN ADDRESS FROM WHICH TO RE-LOAD MPAC. FOUR OPTIONS ARE AVAILABLE:
044927,000410:
044928,000411: # 1. STORE STORE MPAC. THE E ADDRESS MAY BE INDEXED.
044929,000412: # 2. STODL STORE MPAC AND RE-LOAD IT IN DP WITH THE NEXT ADDRESS (THE LOAD MAY BE INDEXED).
044930,000413: # 3. STOVL STORE MPAC AND RE-LOAD A VECTOR (AS ABOVE).
044931,000414: # 4. STCALL STORE AND DO A CALL (BOTH ADDRESSES MUST BE DIRECT HERE).
044932,000415:
044933,000416: # STODL AND STOVL WILL TAKE FROM THE PUSH-DOWN LIST IF NO LOAD ADDRESS IS GIVEN.
044934,000417:
044935,000418: 6354 BLOCK 3
044936,000419:
044937,000420: 6354 COUNT* $$/INTER
044938,000421: 6354 30165 STADR CA BANKSET # THE STADR CODE (PUSHUP UP ON STORE
044939,000422: 6355 54004 TS FBANK # ADDRESS) ENTERS HERE.
044940,000423: 6356 24164 INCR LOC
044941,000424: 6357 50164 ITR1 INDEX LOC # THE STORECODE WAS STORED COMPLEMENTED TO
044942,000425: 6360 40000 CS 0 # MAKE IT LOOK LIKE AN OPCODE PAIR.
044943,000426: 6361 67747 AD NEGONE # (YUL CANT REMOVE 1 BECAUSE OF EARLY CCS)
044944,000427:
044945,000428: 6362 54116 DOSTORE TS ADDRWD
044946,000429: 6363 74356 MASK LOW11 # ENTRY FROM DISPATCHER. SAVE THE ERASABLE
044947,000430: 6364 56116 XCH ADDRWD # ADDRESS AND JUMP ON THE STORE CODE NO.
044948,000431: 6365 77722 MASK B12T14
044949,000432: 6366 00006 EXTEND
044950,000433: 6367 74747 MP BIT5 # EACH TRANSFER VECTOR ENTRY IS TWO WORDS.
044951,000434: 6370 50000 INDEX A
044952,000435: 6371 16372 TCF STORJUMP
044953,000436:
![]() |
Page 1015 |
044955,000438: # STORE CODE JUMP TABLE. CALLS THE APPROPRIATE STORING ROUTINE AND EXITS TO DANZIG OR TO ADDRESS WITH
044956,000439: # A SUPPLIED OPERATION CODE.
044957,000440:
044958,000441: # STORE STORE,1 AND STORE,2 RETURN TO DANZIG, THUS RESETTING THE EBANK TO ITS STATE AT INTPRET.
044959,000442:
044960,000443: 6372 06422 STORJUMP TC STORE # STORE.
044961,000444: 6373 16061 TCF DANZIG # PICK UP NEW OP CODE(S).
044962,000445: 6374 06414 TC STORE,1
044963,000446: 6375 16061 TCF DANZIG
044964,000447: 6376 06417 TC STORE,2
044965,000448: 6377 16061 TCF DANZIG
044966,000449:
044967,000450: 6400 06422 TC STORE # STODL.
044968,000451: 6401 16460 TCF DODLOAD
044969,000452:
044970,000453: 6402 06422 TC STORE # STODL WITH INDEXED LOAD ADDRESS.
044971,000454: 6403 16144 TCF DODLOAD*
044972,000455:
044973,000456: 6404 06422 TC STORE # STOVL.
044974,000457: 6405 16463 TCF DOVLOAD
044975,000458:
044976,000459: 6406 06422 TC STORE # STOVL WITH INDEXED LOAD ADDRESS.
044977,000460: 6407 16466 TCF DOVLOAD*
044978,000461: 6410 06422 TC STORE # STOTC.
044979,000462: 6411 34766 CAF CALLCODE
044980,000463: 6412 54020 TS CYR
044981,000464: 6413 16252 TCF 15BITADR # GET A 15 BIT ADDRESS.
044982,000465:
![]() |
Page 1016 |
044984,000467: # STORE CODE ADDRESS PROCESSOR.
044985,000468:
044986,000469: 6414 50120 STORE,1 INDEX FIXLOC
044987,000470: 6415 40046 CS X1
044988,000471: 6416 16421 TCF PRESTORE
044989,000472:
044990,000473: 6417 50120 STORE,2 INDEX FIXLOC
044991,000474: 6420 40047 CS X2
044992,000475: 6421 26116 PRESTORE ADS ADDRWD # RESULTANT ADDRESS IS IN ERASABLE.
044993,000476:
044994,000477: 6422 40116 STORE CS ADDRWD
044995,000478: 6423 64772 AD DEC45
044996,000479: 6424 10000 CCS A # DOES THE ADDRESS POINT TO THE WORK AREA?
044997,000480: 6425 30120 CA FIXLOC # YES.
044998,000481: 6426 16433 TCF AHEAD5
044999,000482: 6427 35007 CA OCT1400 # NO. SET EBANK & MAKE UP SUBADDRESS.
045000,000483: 6430 56116 XCH ADDRWD
045001,000484: 6431 54003 TS EBANK
045002,000485: 6432 74357 MASK LOW8
045003,000486: 6433 26116 AHEAD5 ADS ADDRWD
045004,000487:
![]() |
Page 1017 |
045006,000489: # STORING ROUTINES. STORE DP, TP, OR VECTOR AS INDICATED BY MODE.
045007,000490:
045008,000491: 6434 00006 STARTSTO EXTEND # MPAC,+1 MUST BE STORED IN ANY EVENT.
045009,000492: # ITRACE (5) REFERS TO "STARTSTO".
045010,000493: 6435 30155 DCA MPAC
045011,000494: 6436 50116 INDEX ADDRWD
045012,000495: 6437 52001 DXCH 0
045013,000496:
045014,000497: 6440 10163 CCS MODE
045015,000498: 6441 16454 TCF TSTORE
045016,000499: 6442 00002 TC Q
045017,000500:
045018,000501: 6443 00006 VSTORE EXTEND
045019,000502: 6444 30160 DCA MPAC +3
045020,000503: 6445 50116 INDEX ADDRWD
045021,000504: 6446 52003 DXCH 2
045022,000505:
045023,000506: 6447 00006 EXTEND
045024,000507: 6450 30162 DCA MPAC +5
045025,000508: 6451 50116 INDEX ADDRWD
045026,000509: 6452 52005 DXCH 4
045027,000510: 6453 00002 TC Q
045028,000511:
045029,000512: 6454 30156 TSTORE CA MPAC +2
045030,000513: 6455 50116 INDEX ADDRWD
045031,000514: 6456 54002 TS 2
045032,000515: 6457 00002 TC Q
045033,000516:
![]() |
Page 1018 |
045035,000518: # ROUTINES TO BEGIN PROCESSING OF THE SECOND ADDRESS ASSOCIATED WITH ALL STORE-TYPE CODES EXCEPT STORE
045036,000519: # ITSELF.
045037,000520:
045038,000521: 6460 37732 DODLOAD CAF DLOADCOD
045039,000522: 6461 54020 TS CYR
045040,000523: 6462 16106 TCF DIRADRES # GO GET A DIRECT ADDRESS.
045041,000524:
045042,000525: 6463 34735 DOVLOAD CAF VLOADCOD
045043,000526: 6464 54020 TS CYR
045044,000527: 6465 16106 TCF DIRADRES
045045,000528:
045046,000529: 6466 36107 DOVLOAD* CAF VLOAD*
045047,000530: 6467 16145 TCF DODLOAD* +1 # PROLOGUE TO INDEX ROUTINE.
045048,000531:
![]() |
Page 1019 |
045050,000533: # THE FOLLOWING LOAD INSTRUCTIONS ARE PROVIDED FOR LOADING THE MULTI-PURPOSE ACCUMULATOR MPAC.
045051,000534:
045052,000535: 6470 50116 TLOAD INDEX ADDRWD
045053,000536: 6471 30002 CA 2 # LOAD A TRIPLE PRECISION ARGUMENT INTO
045054,000537: 6472 54156 TS MPAC +2 # THE FIRST THREE MPAC REGISTERS, WITH THE
045055,000538: 6473 00006 EXTEND # CONTENTS OF THE OTHER FOUR IRRELEVANT.
045056,000539: 6474 50116 INDEX ADDRWD
045057,000540: 6475 30001 DCA 0
045058,000541: 6476 52155 DXCH MPAC
045059,000542: 6477 34753 TMODE CAF ONE
045060,000543: 6500 16060 TCF NEWMODE # DECLARE TRIPLE PRECISION MODE.
045061,000544:
045062,000545: 6501 22007 SLOAD ZL # LOAD A SINGLE PRECISION NUMBER INTO
045063,000546: 6502 50116 INDEX ADDRWD # MPAC, SETTING MPAC+1,2 TO ZERO. THE
045064,000547: 6503 30000 CA 0 # CONTENTS OF THE REMAINING MPAC REGISTERS
045065,000548: 6504 16055 TCF SLOAD2 # ARE IRRELEVANT.
045066,000549:
045067,000550: 6505 00006 VLOAD EXTEND # LOAD A DOUBLE PRECISION VECTOR INTO
045068,000551: 6506 50116 INDEX ADDRWD # MPAC,+1, MPAC+3,4, AND MPAC+5,6. THE
045069,000552: 6507 30001 DCA 0 # CONTENTS OF MPAC +2 ARE IRRELEVANT.
045070,000553: 6510 52155 DXCH MPAC
045071,000554:
045072,000555: 6511 00006 ENDVLOAD EXTEND # PDVL COMES HERE TO FINISH UP FOR DP, TP.
045073,000556: 6512 50116 INDEX ADDRWD
045074,000557: 6513 30003 DCA 2
045075,000558: 6514 52160 DXCH MPAC +3
045076,000559:
045077,000560: 6515 00006 +4 EXTEND # TPDVL FINISHES HERE.
045078,000561: 6516 50116 INDEX ADDRWD
045079,000562: 6517 30005 DCA 4
045080,000563: 6520 52162 DXCH MPAC +5
045081,000564:
045082,000565: 6521 44753 VMODE CS ONE # DECLARE VECTOR MODE.
045083,000566: 6522 16060 TCF NEWMODE
045084,000567:
![]() |
Page 1020 |
045086,000569: # THE FOLLOWING INSTRUCTIONS ARE PROVIDED FOR STORING OPERANDS IN THE PUSHDOWN LIST:
045087,000570: # 1. PUSH PUSHDOWN AND NO LOAD.
045088,000571: # 2. PDDL PUSHDOWN AND DOUBLE PRECISION LOAD.
045089,000572: # 3. PDVL PUSHDOWN AND VECTOR LOAD.
045090,000573:
045091,000574: 6523 00006 PDDL EXTEND
045092,000575: 6524 50116 INDEX ADDRWD # LOAD MPAC,+1, PUSHING THE FORMER
045093,000576: 6525 30001 DCA 0 # CONTENTS DOWN.
045094,000577: 6526 52155 DXCH MPAC
045095,000578: 6527 50166 INDEX PUSHLOC
045096,000579: 6530 52001 DXCH 0
045097,000580:
045098,000581: 6531 50163 INDEX MODE # ADVANCE THE PUSHDOWN POINTER APPRO-
045099,000582: 6532 36244 CAF NO.WDS # PRIATELY.
045100,000583: 6533 26166 ADS PUSHLOC
045101,000584:
045102,000585: 6534 10163 CCS MODE
045103,000586: 6535 16552 TCF ENDTPUSH
045104,000587: 6536 16550 TCF ENDDPUSH
045105,000588:
045106,000589: 6537 54163 TS MODE # NOW DP.
045107,000590: 6540 54156 ENDVPUSH TS MPAC +2
045108,000591: 6541 52160 DXCH MPAC +3 # PUSH DOWN THE REST OF THE VECTOR HERE.
045109,000592: 6542 50166 INDEX PUSHLOC
045110,000593: 6543 51775 DXCH 0 -4
045111,000594:
045112,000595: 6544 52162 DXCH MPAC +5
045113,000596: 6545 50166 INDEX PUSHLOC
045114,000597: 6546 51777 DXCH 0 -2
045115,000598:
045116,000599: 6547 16061 TCF DANZIG
045117,000600:
045118,000601: 6550 54156 ENDDPUSH TS MPAC +2 # SET MPAC +2 TO ZERO AND EXIT ON DP.
045119,000602: 6551 16061 TCF DANZIG
045120,000603:
045121,000604: 6552 54163 ENDTPUSH TS MODE
045122,000605: 6553 56156 XCH MPAC +2 # ON TRIPLE, SET MPAC +2 TO ZERO, PUSHING
045123,000606: 6554 50166 +2 INDEX PUSHLOC # DOWN THE OLD CONTENTS
045124,000607: 6555 53777 TS 0 -1
045125,000608: 6556 16061 TCF DANZIG
045126,000609:
![]() |
Page 1021 |
045128,000611: # PDVL - PUSHDOWN AND VECTOR LOAD.
045129,000612:
045130,000613: 6557 00006 PDVL EXTEND # RELOAD MPAC AND PUSH DOWN ITS CONTENTS.
045131,000614: 6560 50116 INDEX ADDRWD
045132,000615: 6561 30001 DCA 0
045133,000616: 6562 52155 DXCH MPAC
045134,000617: 6563 50166 INDEX PUSHLOC
045135,000618: 6564 52001 DXCH 0
045136,000619:
045137,000620: 6565 50163 INDEX MODE # ADVANCE THE PUSHDOWN POINTER.
045138,000621: 6566 36244 CAF NO.WDS
045139,000622: 6567 26166 ADS PUSHLOC
045140,000623:
045141,000624: 6570 10163 CCS MODE # TEST PAST MODE.
045142,000625: 6571 16610 TCF TPDVL
045143,000626: 6572 16511 TCF ENDVLOAD # JUST LOAD LAST FOUR REGISTERS ON DP.
045144,000627:
045145,000628: 6573 00006 VPDVL EXTEND # PUSHDOWN AND RE-LOAD LAST TWO COMPONENTS
045146,000629: 6574 50116 INDEX ADDRWD
045147,000630: 6575 30003 DCA 2
045148,000631: 6576 52160 DXCH MPAC +3
045149,000632: 6577 50166 INDEX PUSHLOC
045150,000633: 6600 51775 DXCH 0 -4
045151,000634:
045152,000635: 6601 00006 EXTEND
045153,000636: 6602 50116 INDEX ADDRWD
045154,000637: 6603 30005 DCA 4
045155,000638: 6604 52162 DXCH MPAC +5
045156,000639: 6605 50166 INDEX PUSHLOC
045157,000640: 6606 51777 DXCH 0 -2
045158,000641:
045159,000642: 6607 16061 TCF DANZIG
045160,000643:
045161,000644: 6610 00006 TPDVL EXTEND # ON TP, WE MUST LOAD THE Y COMPONENT
045162,000645: 6611 50116 INDEX ADDRWD # BEFORE STORING MPAC +2 IN CASE THIS IS A
045163,000646: 6612 30003 DCA 2 # PUSHUP.
045164,000647: 6613 52160 DXCH MPAC +3
045165,000648:
045166,000649: 6614 30156 CA MPAC +2
045167,000650: 6615 50166 INDEX PUSHLOC # IN DP.
045168,000651: 6616 53777 TS 0 -1
045169,000652: 6617 16515 TCF ENDVLOAD +4
045170,000653:
045171,000654: # SSP (STORE SINGLE PRECISION) IS EXECUTED HERE.
045172,000655:
045173,000656: 6620 24164 SSP INCR LOC # PICK UP THE WORD FOLLOWING THE GIVEN
045174,000657: 6621 50164 INDEX LOC # ADDRESS AND STORE IT AT X.
045175,000658: 6622 30000 CA 0
045176,000659: 6623 50116 STORE1 INDEX ADDRWD # SOME INDEX AND MISCELLANEOUS OPS END
045177,000660: 6624 54000 TS 0 # HERE.
![]() |
Page 1022 |
045179,000662: 6625 16061 TCF DANZIG
045180,000663:
![]() |
Page 1023 |
045182,000665: # SEQUENCE CHANGING AND SUBROUTINE CALLING OPTIONS.
045183,000666:
045184,000667: # THE FOLLOWING OPERATIONS ARE AVAILABLE FOR SEQUENCING CHANGING, BRANCHING, AND CALLING SUBROUTINES:
045185,000668: # 1. GOTO GO TO.
045186,000669: # 2. CALL CALL SUBROUTINE SETTING QPRET.
045187,000670: # 3. CGOTO COMPUTED GO TO.
045188,000671: # 4. CCALL COMPUTED CALL.
045189,000672: # 7. BPL BRANCH IF MPAC POSITIVE OR ZERO.
045190,000673: # 8. BZE BRANCH IF MPAC ZERO.
045191,000674: # 9. BMN BRANCH IF MPAC NEGATIVE NON-ZERO.
045192,000675:
045193,000676: 6626 24164 CCALL INCR LOC # MAINTAIN LOC FOR QPRET COMPUTATION.
045194,000677: 6627 50164 INDEX LOC
045195,000678: 6630 30000 CAF 0 # GET BASE ADDRESS OF CADR LIST.
045196,000679: 6631 50116 INDEX ADDRWD
045197,000680: 6632 60000 AD 0 # ADD INCREMENT.
045198,000681: 6633 54004 TS FBANK # SELECT DESIRED CADR.
045199,000682: 6634 75012 MASK LOW10
045200,000683: 6635 50000 INDEX A
045201,000684: 6636 32000 CAF 10000
045202,000685: 6637 54117 TS POLISH
045203,000686:
045204,000687: 6640 30165 CALL CA BANKSET # FOR ANY OF THE CALL OPTIONS, MAKE UP THE
045205,000688: 6641 74350 MASK BANKMASK # ADDRESS OF THE NEXT OP-CODE PAIR/STORE
045206,000689: 6642 64350 AD BANKMASK # CODE AND LEAVE IT IN QPRET. NOTE THAT
045207,000690: 6643 60164 AD LOC # BANKMASK = -(2000 - 1).
045208,000691: 6644 50120 INDEX FIXLOC
045209,000692: 6645 54052 TS QPRET
045210,000693:
045211,000694: 6646 30117 GOTO CA POLISH # BASIC BRANCHING SEQUENCE.
045212,000695: 6647 77742 +1 MASK HIGH4
045213,000696: 6650 00006 EXTEND
045214,000697: 6651 16662 BZF GOTOERS # SEE IF ADDRESS POINTS TO FIXED OR ERAS.
045215,000698: 6652 30165 +4 CA BANKSET # SET EBANK PART OF BBANK. NEXT, SET UP
045216,000699: 6653 54006 TS BBANK # FBANK. THE COMBINATION IS PICKED UP &
045217,000700: 6654 30117 CA POLISH # PUT INTO BANKSET AT INTPRET +2.
045218,000701: 6655 54004 TS FBANK
045219,000702: 6656 75012 MASK LOW10
045220,000703: 6657 64741 AD 2K
045221,000704: 6660 54164 TS LOC
045222,000705: 6661 16042 TCF INTPRET +3
045223,000706:
045224,000707: 6662 E3,1400 EBANK= 1400 # SO YUL DOESN'T CUSS THE "CA 1400" BELOW.
045225,000708:
045226,000709: 6662 30117 GOTOERS CA POLISH # THE GIVEN ADDRESS IS IN ERASABLE - SEE
045227,000710: 6663 66251 AD -ENDVAC # IF RELATIVE TO THE WORK AREA.
045228,000711: 6664 10000 CCS A
045229,000712: 6665 30117 CA POLISH # GENERAL ERASABLE.
045230,000713: 6666 16675 TCF GOTOGE
045231,000714:
![]() |
Page 1024 |
045233,000716: 6667 30120 CA FIXLOC # WORK AREA.
045234,000717: 6670 60117 AD POLISH
045235,000718: 6671 50000 INDEX A # USE THE GIVEN ADDRESS AS THE ADDRESS OF
045236,000719: 6672 30000 CA 0 # THE BRANCH ADDRESS.
045237,000720: 6673 54117 TS POLISH
045238,000721: 6674 16647 TCF GOTO +1 # ALLOWS ARBITRARY INDIRECTNESS LEVELS.
045239,000722:
045240,000723: 6675 54003 GOTOGE TS EBANK
045241,000724: 6676 74357 MASK LOW8
045242,000725: 6677 50000 INDEX A # USE THE GIVEN ADDRESS AS THE ADDRESS OF
045243,000726: 6700 31400 CA 1400 # THE BRANCH ADDRESS.
045244,000727: 6701 54117 TS POLISH
045245,000728: 6702 16647 TCF GOTO +1
045246,000729:
045247,000730: 6703 50164 CGOTO INDEX LOC # COMPUTED GO TO. PICK UP ADDRESS OF CADR
045248,000731: 6704 30001 CA 1 # LIST.
045249,000732: 6705 50116 INDEX ADDRWD # ADD MODIFIER.
045250,000733: 6706 60000 AD 0
045251,000734: 6707 54004 TS FBANK # SELECT GOTO ADDRESS.
045252,000735: 6710 75012 MASK LOW10
045253,000736: 6711 50000 INDEX A
045254,000737: 6712 32000 CA 10000
045255,000738: 6713 54117 TS POLISH
045256,000739: 6714 16647 TCF GOTO +1 # WITH ADDRESS IN A.
045257,000740:
045258,000741: 6715 30165 SWBRANCH CA BANKSET # SWITCH INSTRUCTIONS WHICH ELECT TO
045259,000742: 6716 54004 TS FBANK # BRANCH COME HERE TO DO SO.
045260,000743: 6717 50164 INDEX LOC
045261,000744: 6720 30001 CA 1
045262,000745: 6721 54117 TS POLISH
045263,000746: 6722 16647 TCF GOTO +1
045264,000747:
![]() |
Page 1025 |
045266,000749: # TRIPLE PRECISION BRANCHING ROUTINE. IF CALLING TC IS AT L, RETURN IS AS FOLLOWS:
045267,000750: # L+1 IF MPAC IS GREATER THAN ZERO.
045268,000751: # L+2 IF MPAC IS EQUAL TO +0 OR -0.
045269,000752: # L+3 IF MPAC IS LESS THAN ZERO.
045270,000753:
045271,000754: 6723 10154 BRANCH CCS MPAC
045272,000755: 6724 00002 TC Q
045273,000756: 6725 16727 TCF +2 # ON ZERO.
045274,000757: 6726 16741 TCF NEG
045275,000758:
045276,000759: 6727 10155 CCS MPAC +1
045277,000760: 6730 00002 TC Q
045278,000761: 6731 16733 TCF +2
045279,000762: 6732 16741 TCF NEG
045280,000763:
045281,000764: 6733 10156 CCS MPAC +2
045282,000765: 6734 00002 TC Q
045283,000766: 6735 16737 TCF +2
045284,000767: 6736 16741 TCF NEG
045285,000768:
045286,000769: 6737 50002 Q+1 INDEX Q
045287,000770: 6740 00001 TC 1
045288,000771:
045289,000772: 6741 50002 NEG INDEX Q # IF FIRST NON-ZERO REGISTER WAS NEGATIVE.
045290,000773: 6742 00002 TC 2
045291,000774:
045292,000775: 6743 6741 Q+2 = NEG
045293,000776:
045294,000777: # ITRACE (3) REFERS TO "EXIT".
045295,000778:
045296,000779: 6743 30165 EXIT CA BANKSET # RESTORE USER'S BANK SETTING, AND LEAVE
045297,000780: 6744 54006 TS BBANK # INTERPRETIVE MODE.
045298,000781: 6745 50164 INDEX LOC
045299,000782: 6746 00001 TC 1
045300,000783:
![]() |
Page 1026 |
045302,000785: # SECTION 3 - ADD/SUBTRACT PACKAGE.
045303,000786:
045304,000787: # THE FOLLOWING OPERATIONS ARE PROVIDED FOR ADDING TO AND SUBTRACTING FROM THE MULTI-PURPOSE ACCUMULATOR
045305,000788: # MPAC:
045306,000789:
045307,000790: # 1. DAD DOUBLE PRECISION ADD.
045308,000791: # 2. DSU DOUBLE PRECISION SUBTRACT.
045309,000792: # 3. BDSU DOUBLE PRECISION SUBTRACT FROM.
045310,000793: # 4. TAD TRIPLE PRECISION ADD.
045311,000794: # 5. VAD VECTOR ADD.
045312,000795: # 6. VSU VECTOR SUBTRACT.
045313,000796: # 7. BVSU VECTOR SUBTRACT FROM.
045314,000797: # THE INTERPRETIVE OVERFLOW INDICATOR OVFIND IS SET NON-ZERO IF OVERFLOW OCCURS IN ANY OF THE ABOVE.
045315,000798:
045316,000799: 6747 34735 VSU CAF BIT15 # CHANGES 0 TO DCS.
045317,000800: 6750 16752 TCF +2
045318,000801:
045319,000802: 6751 34355 VAD CAF PRIO30 # CHANGES 0 TO DCA.
045320,000803: 6752 26116 ADS ADDRWD
045321,000804: 6753 00006 EXTEND
045322,000805: 6754 50116 INDEX ADDRWD
045323,000806: 6755 00003 READ HISCALAR # DCA 2 OR DCS 2
045324,000807: 6756 20160 DAS MPAC +3
045325,000808: 6757 00006 EXTEND # CHECK OVERFLOW.
045326,000809: 6760 16762 BZF +2
045327,000810: 6761 07014 TC OVERFLWY
045328,000811:
045329,000812: 6762 00006 EXTEND
045330,000813: 6763 50116 INDEX ADDRWD
045331,000814: 6764 00005 READ CHAN5 # DCA 4 OR DCS 4
045332,000815: 6765 20162 DAS MPAC +5
045333,000816: 6766 00006 EXTEND
045334,000817: 6767 16771 BZF +2
045335,000818: 6770 07011 TC OVERFLWZ
045336,000819:
045337,000820: 6771 00006 EXTEND
045338,000821: 6772 50116 INDEX ADDRWD
045339,000822: 6773 00001 READ LCHAN # DCA 0 OR DCS 0
045340,000823: 6774 17000 TCF ENDVXV
045341,000824:
045342,000825: 6775 00006 DAD EXTEND
045343,000826: 6776 50116 INDEX ADDRWD
045344,000827: 6777 30001 DCA 0
045345,000828: 7000 20155 ENDVXV DAS MPAC # VXV FINISHES HERE.
045346,000829: 7001 00006 EXTEND
045347,000830: 7002 16061 BZF DANZIG
045348,000831:
![]() |
Page 1027 |
045350,000833: 7003 07017 SETOVF TC OVERFLOW
045351,000834: 7004 16061 TCF DANZIG
045352,000835:
![]() |
Page 1028 |
045354,000837: 7005 00006 DSU EXTEND
045355,000838: 7006 50116 INDEX ADDRWD
045356,000839: 7007 40001 DCS 0
045357,000840: 7010 17000 TCF ENDVXV
045358,000841:
045359,000842: 7011 54001 OVERFLWZ TS L # ENTRY FOR THIRD COMPONENT.
045360,000843: 7012 34756 CAF FIVE
045361,000844: 7013 17016 TCF +3
045362,000845:
045363,000846: 7014 54001 OVERFLWY TS L # ENTRY FOR SECOND COMPONENT.
045364,000847: 7015 36245 CAF THREE
045365,000848: 7016 56001 XCH L
045366,000849:
045367,000850: 7017 50000 OVERFLOW INDEX A # ENTRY FOR 1ST COMP OR DP (L=0).
045368,000851: 7020 44734 CS LIMITS # PICK UP POSMAX OR NEGMAX.
045369,000852: 7021 54130 TS BUF
045370,000853: 7022 00006 EXTEND
045371,000854: 7023 24000 AUG A # FORCE OVERFLOW.
045372,000855: 7024 50001 INDEX L
045373,000856: 7025 26155 ADS MPAC +1
045374,000857: 7026 54007 TS 7
045375,000858: 7027 34755 CAF ZERO
045376,000859: 7030 60130 AD BUF
045377,000860: 7031 50001 INDEX L
045378,000861: 7032 26154 ADS MPAC
045379,000862: 7033 54007 TS 7
045380,000863: 7034 00002 TC Q # NO OVERFLOW EXIT.
045381,000864: 7035 17152 TCF SETOVF2 # SET OVFIND AND EXIT.
045382,000865:
045383,000866: 7036 00006 BVSU EXTEND
045384,000867: 7037 50116 INDEX ADDRWD
045385,000868: 7040 30003 DCA 2
045386,000869: 7041 52160 DXCH MPAC +3
045387,000870: 7042 00006 EXTEND
045388,000871: 7043 40001 DCOM
045389,000872: 7044 20160 DAS MPAC +3
045390,000873: 7045 00006 EXTEND
045391,000874: 7046 17050 BZF +2
045392,000875: 7047 07014 TC OVERFLWY
045393,000876:
045394,000877: 7050 00006 EXTEND
045395,000878: 7051 50116 INDEX ADDRWD
045396,000879: 7052 30005 DCA 4
045397,000880: 7053 52162 DXCH MPAC +5
045398,000881: 7054 00006 EXTEND
045399,000882: 7055 40001 DCOM
045400,000883: 7056 20162 DAS MPAC +5
045401,000884: 7057 00006 EXTEND
045402,000885: 7060 17062 BZF +2
045403,000886: 7061 07011 TC OVERFLWZ
045404,000887:
![]() |
Page 1029 |
045406,000889: 7062 00006 BDSU EXTEND
045407,000890: 7063 50116 INDEX ADDRWD
045408,000891: 7064 30001 DCA 0
045409,000892: 7065 52155 DXCH MPAC
045410,000893: 7066 00006 EXTEND
045411,000894: 7067 40001 DCOM
045412,000895: 7070 17000 TCF ENDVXV
045413,000896:
![]() |
Page 1030 |
045415,000898: # TRIPLE PRECISION ADD ROUTINE.
045416,000899:
045417,000900: 7071 00006 TAD EXTEND
045418,000901: 7072 50116 INDEX ADDRWD
045419,000902: 7073 30002 DCA 1 # ADD MINOR PARTS FIRST.
045420,000903: 7074 20156 DAS MPAC +1
045421,000904: 7075 50116 INDEX ADDRWD
045422,000905: 7076 60000 AD 0
045423,000906: 7077 60154 AD MPAC
045424,000907: 7100 54154 TS MPAC
045425,000908: 7101 16061 TCF DANZIG
045426,000909:
045427,000910: 7102 17003 TCF SETOVF # SET OVFIND IF SUCH OCCURS.
045428,000911:
![]() |
Page 1031 |
045430,000913: # ARITHMETIC SUBROUTINES REQUIRED IN FIXED-FIXED.
045431,000914: # 1. DMPSUB DOUBLE PRECISION MULTIPLY. MULTIPLY THE CONTENTS OF MPAC,+1 BY THE DP WORD WHOSE ADDRESS
045432,000915: # IS IN ADDRWD AND LEAVE A TRIPLE PRECISION RESULT IN MPAC.
045433,000916: # 2. ROUNDSUB ROUND THE TRIPLE PRECISON CONTENTS OF MPAC TO DOUBLE PRECISION.
045434,000917: # 3. DOTSUB TAKE THE DOT PRODUCT OF THE VECTOR IN MPAC AND THE VECTOR WHOSE ADDRESS IS IN ADDRWD
045435,000918: # AND LEAVE THE TRIPLE PRECISION RESULT IN MPAC.
045436,000919: # 4. POLY USING THE CONTENTS OF MPAC AS A DP ARGUMENT, EVALUATE THE POLYNOMIAL WHOSE DEGREE AND
045437,000920: # COEFFICIENTS IMMEDIATELY FOLLOW THE TC POLY INSTRUCTION (SEE ROUTINE FOR DETAILS.)
045438,000921:
045439,000922: 7103 50002 DMP INDEX Q # BASIC SUBROUTINE FOR USE BY PINBALL, ETC
045440,000923: 7104 30000 CAF 0 # ADRES OF ARGUMENT FOLLOWS TC DMP .
045441,000924: 7105 24002 INCR Q
045442,000925: 7106 54116 -1 TS ADDRWD # (PROLOGUE FOR SETTING ADDRWD.)
045443,000926:
045444,000927: 7107 50116 DMPSUB INDEX ADDRWD # GET MINOR PART OF OPERAND AT C(ADDRWD).
045445,000928: 7110 30001 CA 1
045446,000929: 7111 54156 TS MPAC +2 # THIS WORKS FOR SQUARING MPAC AS WELL.
045447,000930: 7112 34755 CAF ZERO # SET MPAC +1 TO ZERO SO WE CAN ACCUMULATE
045448,000931: 7113 56155 XCH MPAC +1 # THE PARTIAL PRODUCTS WITH DAS
045449,000932: 7114 54135 TS MPTEMP # INSTRUCTIONS.
045450,000933: 7115 00006 EXTEND
045451,000934: 7116 70156 MP MPAC +2 # MINOR OF MPAC X MINOR OF C(ADDRWD).
045452,000935:
045453,000936: 7117 56156 XCH MPAC +2 # DISCARD MINOR PART OF ABOVE RESULT AND
045454,000937: 7120 00006 EXTEND # FORM MAJOR OF MPAC X MINOR OF C(ADDRWD).
045455,000938: 7121 70154 MP MPAC
045456,000939: 7122 20156 DAS MPAC +1 # GUARANTEED NO OVERFLOW.
045457,000940:
045458,000941: 7123 50116 INDEX ADDRWD # GET MAJOR PART OF ARGUMENT AT C(ADDRWD).
045459,000942: 7124 30000 CA 0
045460,000943: 7125 56135 XCH MPTEMP # SAVE AND BRING OUT MINOR OF MPAC.
045461,000944: 7126 00006 DMPSUB2 EXTEND
045462,000945: 7127 70135 MP MPTEMP # MAJOR OF C(ADDRWD) X MINOR OF MPAC.
045463,000946: 7130 20156 DAS MPAC +1 # ACCUMULATE, SETTING A TO NET OVERFLOW.
045464,000947:
045465,000948: 7131 56154 XCH MPAC # SETTING MPAC TO 0 OR +-1.
045466,000949: 7132 00006 EXTEND
045467,000950: 7133 70135 MP MPTEMP # MAJOR OF MPAC X MAJOR OF C(ADDRWD).
045468,000951: 7134 20155 DAS MPAC # GUARANTEED NO OVERFLOW.
045469,000952: 7135 00002 TC Q # 49 MCT = .573 MS. INCLUDING RETURN.
045470,000953:
![]() |
Page 1032 |
045472,000955: # ROUND MPAC TO DOUBLE PRECISION, SETTING OVFIND ON THE RARE EVENT OF OVERFLOW.
045473,000956:
045474,000957: 7136 34755 ROUNDSUB CAF ZERO # SET MPAC +2 = 0 FOR SCALARS AND CHANGE
045475,000958: 7137 54163 +1 TS MODE # MODE TO DP.
045476,000959:
045477,000960: 7140 56156 VROUND XCH MPAC +2 # BUT WE NEEDNT TAKE THE TIME FOR VECTORS.
045478,000961: 7141 60000 DOUBLE
045479,000962: 7142 54001 TS L
045480,000963: 7143 00002 TC Q
045481,000964:
045482,000965: 7144 60155 AD MPAC +1 # ADD ROUNDING BIT IF MPAC +2 WAS GREATER
045483,000966: 7145 54155 TS MPAC +1 # THAN .5 IN MAGNITUDE.
045484,000967: 7146 00002 TC Q
045485,000968:
045486,000969: 7147 60154 AD MPAC # PROPAGATE INTERFLOW.
045487,000970: 7150 54154 TS MPAC
045488,000971: 7151 00002 TC Q
045489,000972:
045490,000973: 7152 54121 SETOVF2 TS OVFIND # (RARE).
045491,000974: 7153 00002 TC Q
045492,000975:
![]() |
Page 1033 |
045494,000977: # THE DOT PRODUCT SUBROUTINE USUALLY FORMS THE DOT PRODUCT OF THE VECTOR IN MPAC WITH A STANDARD SIX
045495,000978: # REGISTER VECTOR WHOSE ADDRESS IS IN ADDRWD. IN THIS CASE C(DOTINC) ARE SET TO 2. VXM, HOWEVER, SETS C(DOTINC) TO
045496,000979: # 6 SO THAT DOTSUB DOTS MPAC WITH A COLUMN VECTOR OF THE MATRIX IN QUESTION IN THIS CASE.
045497,000980:
045498,000981: 7154 34752 PREDOT CAF TWO # PROLOGUE TO SET DOTINC TO 2.
045499,000982: 7155 54136 TS DOTINC
045500,000983:
045501,000984: 7156 00006 DOTSUB EXTEND
045502,000985: 7157 22137 QXCH DOTRET # SAVE RETURN.
045503,000986: 7160 07107 TC DMPSUB # DOT X COMPONENTS.
045504,000987: 7161 52160 DXCH MPAC +3 # POSITION Y COMPONENT OF MPAC FOR
045505,000988: 7162 52155 DXCH MPAC # MULTIPLICATION WHILE SAVING RESULT IN
045506,000989: 7163 52131 DXCH BUF # THREE WORD BUFFER, BUF.
045507,000990: 7164 30156 CA MPAC +2
045508,000991: 7165 54132 TS BUF +2
045509,000992:
045510,000993: 7166 30136 CA DOTINC # ADVANCE ADDRWD TO Y COMPONENT OF
045511,000994: 7167 26116 ADS ADDRWD # OTHER ARGUMENT.
045512,000995: 7170 07107 TC DMPSUB
045513,000996: 7171 52156 DXCH MPAC +1 # ACCUMULATE PARTIAL PRODUCTS.
045514,000997: 7172 20132 DAS BUF +1
045515,000998: 7173 60154 AD MPAC
045516,000999: 7174 60130 AD BUF
045517,001000: 7175 54130 TS BUF
045518,001001: 7176 17200 TCF +2
045519,001002: 7177 54121 TS OVFIND # IF OVERFLOW OCCURS.
045520,001003:
045521,001004: 7200 52162 DXCH MPAC +5 # MULTIPLY Z COMPONENTS.
045522,001005: 7201 52155 DXCH MPAC
045523,001006: 7202 30136 CA DOTINC
045524,001007: 7203 26116 ADS ADDRWD
045525,001008: 7204 07107 TC DMPSUB
045526,001009: 7205 52132 ENDDOT DXCH BUF +1 # LEAVE FINAL ACCUMULATION IN MPAC.
045527,001010: 7206 20156 DAS MPAC +1
045528,001011: 7207 60154 AD MPAC
045529,001012: 7210 60130 AD BUF
045530,001013: 7211 54154 TS MPAC
045531,001014: 7212 00137 TC DOTRET
045532,001015:
045533,001016: 7213 07017 TC OVERFLOW # ON OVERFLOW HERE.
045534,001017: 7214 00137 TC DOTRET
045535,001018:
![]() |
Page 1034 |
045537,001020: # DOUBLE PRECISION POLYNOMIAL EVALUATOR
045538,001021: # N N-1
045539,001022: # THIS ROUTINE EVALUATES A X + A X + ... + A X + A LEAVING THE DP RESULT IN MPAC ON EXIT.
045540,001023: # N N-1 1 0
045541,001024:
045542,001025: # THE ROUTINE HAS TWO ENTRIES
045543,001026:
045544,001027: # 1. ENTRY THRU POWRSERS. THE COEFFICIENTS MAY BE EITHER IN FIXED OR ERASABLE, THE CALL IS BY
045545,001028: # TC POWRSERS, AND THE RETURN IS TO LOC(TC POWRSERS)+1. THE ENTERING DATA MUST BE AS FOLLOWS
045546,001029:
045547,001030: # A SP LOC-3 ADDRESS FOR REFERENCING COEF TABLE
045548,001031: # L SP N-1 N IS THE DEGREE OF THE POWER SERIES
045549,001032: # MPAC DP X ARGUMENT
045550,001033: # LOC-2N DP A(0)
045551,001034: # ...
045552,001035: # LOC DP A(N)
045553,001036:
045554,001037: # 2. ENTRY THRU POLY. THE CALL TO POLY AND THE ENTERING DATA MUST BE AS FOLLOWS
045555,001038:
045556,001039: # MPAC DP X ARGUMENT
045557,001040: # LOC TC POLY
045558,001041: # LOC+1 SP N-1
045559,001042: # LOC+2 DP A(0)
045560,001043: # ...
045561,001044: # LOC+2N+2 DP A(N) RETURN IS TO LOC+2N+4
045562,001045:
045563,001046: 7215 00006 POWRSERS EXTEND
045564,001047: 7216 22141 QXCH POLYRET # RETURN ADDRESS
045565,001048: 7217 54117 TS POLISH # POWER SERIES ADDRESS
045566,001049: 7220 22140 LXCH POLYCNT # N-1 TO COUNTER
045567,001050: 7221 17232 TCF POLYCOM # SKIP SET UP BY POLY
045568,001051:
045569,001052: 7222 50002 POLY INDEX Q
045570,001053: 7223 30000 CAF 0
045571,001054: 7224 54140 TS POLYCNT # N-1 TO COUNTER
045572,001055: 7225 60000 DOUBLE
045573,001056: 7226 60002 AD Q
045574,001057: 7227 54117 TS POLISH # L(A(N))-3 TO POLISH
045575,001058: 7230 64756 AD FIVE
045576,001059: 7231 54141 TS POLYRET # STORE RETURN ADDRESS
045577,001060:
045578,001061: 7232 36273 POLYCOM CAF LVBUF # INCOMING X WILL BE MOVED TO VBUF, SO
045579,001062: 7233 54116 TS ADDRWD # SET ADDRWD SO DMPSUB WILL MPY BY VBUF.
045580,001063:
045581,001064: 7234 00006 EXTEND
045582,001065: 7235 50117 INDEX POLISH
045583,001066: 7236 30004 DCA 3
045584,001067:
![]() |
Page 1035 |
045586,001069: 7237 52155 DXCH MPAC # LOAD A(N) INTO MPAC,
045587,001070: 7240 52123 DXCH VBUF # SAVING X IN VBUF
045588,001071: 7241 17245 TCF POLY2
045589,001072:
045590,001073: 7242 54140 POLYLOOP TS POLYCNT # SAVE DECREMENTED LOOP COUNTER
045591,001074: 7243 44752 CS TWO
045592,001075: 7244 26117 ADS POLISH # REGRESS COEFFICIENT POINTER
045593,001076:
045594,001077: 7245 07107 POLY2 TC DMPSUB # MULTIPLY BY X
045595,001078: 7246 00006 EXTEND
045596,001079: 7247 50117 INDEX POLISH
045597,001080: 7250 30002 DCA 1 # ADD IN NEXT COEFFICIENT
045598,001081: 7251 20155 DAS MPAC # USERS RESPONSIBILITY TO ASSURE NO OVFLOW
045599,001082:
045600,001083: 7252 10140 CCS POLYCNT
045601,001084: 7253 17242 TCF POLYLOOP
045602,001085: 7254 00141 TC POLYRET # RETURN CALLER
045603,001086:
![]() |
Page 1036 |
045605,001088: # MISCELLANEOUS MULTI-PRECISION ROUTINES REQUIRED IN FIXED-FIXED BUT NOT USED BY THE INTERPRETER.
045606,001089:
045607,001090: 7255 34755 DPAGREE CAF ZERO # DOUBLE PRECISION ENTRY -
045608,001091: 7256 54156 TS MPAC +2 # ZERO LOW-ORDER WORD
045609,001092:
045610,001093: 7257 22002 TPAGREE LXCH Q # FORCE SIGN AGREEMENT AMONG THE TRIPLE
045611,001094: 7260 06723 TC BRANCH # PRECISION CONTENTS OF MPAC. RETURNING
045612,001095: 7261 17265 TCF ARG+ # WITH SIGNUM OF THE INPUT IN A.
045613,001096: 7262 17305 TCF ARGZERO
045614,001097:
045615,001098: 7263 44733 CS POSMAX # IF NEGATIVE.
045616,001099: 7264 17266 TCF +2
045617,001100:
045618,001101: 7265 34733 ARG+ CAF POSMAX
045619,001102: 7266 54002 TS Q
045620,001103: 7267 00006 EXTEND
045621,001104: 7270 24000 AUG A # FORMS +-1.0.
045622,001105: 7271 60156 AD MPAC +2
045623,001106: 7272 54156 TS MPAC +2
045624,001107: 7273 34755 CAF ZERO
045625,001108: 7274 60002 AD Q
045626,001109: 7275 60155 AD MPAC +1
045627,001110: 7276 54155 TS MPAC +1
045628,001111: 7277 34755 CAF ZERO
045629,001112: 7300 60002 AD Q # Q STILL HAS POSMAX OR NEGMAX IN IT.
045630,001113: 7301 60154 AD MPAC
045631,001114: 7302 54154 ARGZERO2 TS MPAC # ALWAYS SKIPPING UNLESS ARGZERO.
045632,001115: 7303 54155 TS MPAC +1
045633,001116: 7304 00001 TC L # RETURN VIA L.
045634,001117:
045635,001118: 7305 54156 ARGZERO TS MPAC +2 # SET ALL THREE MPAC REGISTERS TO ZERO.
045636,001119: 7306 17302 TCF ARGZERO2
045637,001120:
045638,001121: # SHORTMP MULTIPLIES THE TP CONTENTS OF MPAC BY THE SINGLE PRECISION NUMBER ARRIVING IN A.
045639,001122:
045640,001123: 7307 54135 SHORTMP TS MPTEMP
045641,001124: 7310 00006 EXTEND
045642,001125: 7311 70156 MP MPAC +2
045643,001126: 7312 54156 TS MPAC +2
045644,001127: 7313 34755 SHORTMP2 CAF ZERO # SO SUBSEQUENT DAS WILL WORK.
045645,001128: 7314 56155 XCH MPAC +1
045646,001129: 7315 17126 TCF DMPSUB2
045647,001130:
![]() |
Page 1037 |
045649,001132: # DMPNSUB MULTIPLIES THE DP FRACTION ARRIVING IN MPAC BY THE SP
045650,001133: # INTEGER ARRIVING IN A. THE DP PRODUCT DEPARTS BOTH IN MPAC AND IN
045651,001134: # A AND L. NOTE THAT DMPNSUB NORMALLY INCREASES THE MAGNITUDE OF THE
045652,001135: # CONTENTS OF MPAC. THE CUSTOMER MUST INSURE THAT B(A) X B(MPAC,MPAC+1)
045653,001136: # AND B(A) X B(MPAC) ARE LESS THAN 1 IN MAGNITUDE, WHERE B, AS IS OBVIOUS,
045654,001137: # INDICATES THE ARRIVING CONTENTS.
045655,001138:
045656,001139: 7316 54135 DMPNSUB TS DMPNTEMP
045657,001140: 7317 00006 EXTEND
045658,001141: 7320 70155 MP MPAC +1
045659,001142: 7321 52155 DXCH MPAC # LOW PRODUCT TO MPAC, HIGH FACTOR TO A
045660,001143: 7322 00006 EXTEND
045661,001144: 7323 70135 MP DMPNTEMP
045662,001145: 7324 30001 CA L
045663,001146: 7325 26154 ADS MPAC # COMPLETING THE PRODUCT IN MPAC
045664,001147: 7326 00006 EXTEND
045665,001148: 7327 30155 DCA MPAC # BRINGING THE PRODUCT INTO A AND L
045666,001149: 7330 00002 TC Q
045667,001150:
![]() |
Page 1038 |
045669,001152: # MISCELLANEOUS VECTOR OPERATIONS. INCLUDED HERE ARE THE FOLLOWING:
045670,001153: # 1. DOT DP VECTOR DOT PRODUCT.
045671,001154: # 2. VXV DP VECTOR CROSS PRODUCT.
045672,001155: # 3. VXSC DP VECTOR TIMES SCALAR.
045673,001156: # 4. V/SC DP VECTOR DIVIDED BY SCALAR.
045674,001157: # 5. VPROJ DP VECTOR PROJECTION. ( (MPAC.X)MPAC ).
045675,001158: # 6. VXM DP VECTOR POST-MULTIPLIED BY MATRIX.
045676,001159: # 7. MXV DP VECTOR PRE-MULTIPLIED BY MATRIX.
045677,001160:
045678,001161: 7331 07154 DOT TC PREDOT # DO THE DOT PRODUCT AND EXIT, CHANGING
045679,001162: 7332 34755 DMODE CAF ZERO # THE MODE TO DP SCALAR.
045680,001163: 7333 16060 TCF NEWMODE
045681,001164:
045682,001165: 7334 34752 MXV CAF TWO # SET UP MATINC AND DOTINC FOR ROW
045683,001166: 7335 54140 TS MATINC # VECTORS.
045684,001167: 7336 17342 TCF VXM/MXV # GO TO COMMON PORTION.
045685,001168:
045686,001169: 7337 44363 VXM CS TEN # SET MATINC AND DOTINC TO REFER TO MATRIX
045687,001170: 7340 54140 TS MATINC # AS THREE COLUMN VECTORS.
045688,001171: 7341 36242 CAF SIX
045689,001172:
![]() |
Page 1039 |
045691,001174: # COMMON PORTION OF MXV AND VXM.
045692,001175:
045693,001176: 7342 54136 VXM/MXV TS DOTINC
045694,001177: # ITRACE (2) REFERS TO "VXM/MXV".
045695,001178: 7343 07532 TC MPACVBUF # SAVE VECTOR IN MPAC FOR FURTHER USE.
045696,001179:
045697,001180: 7344 07156 TC DOTSUB # GO DOT TO GET X COMPONENT OF ANSWER.
045698,001181: 7345 00006 EXTEND
045699,001182: 7346 30123 DCA VBUF # MOVE MPAC VECTOR BACK INTO MPAC, SAVING
045700,001183: 7347 52155 DXCH MPAC # NEW X COMPONENT IN BUF2.
045701,001184: 7350 52134 DXCH BUF2
045702,001185: 7351 00006 EXTEND
045703,001186: 7352 30125 DCA VBUF +2
045704,001187: 7353 52160 DXCH MPAC +3
045705,001188: 7354 00006 EXTEND
045706,001189: 7355 30127 DCA VBUF +4
045707,001190: 7356 52162 DXCH MPAC +5
045708,001191: 7357 30140 CA MATINC # INITIALIZE ADDRWD FOR NEXT DOT PRODUCT.
045709,001192: 7360 26116 ADS ADDRWD # FORMS BASE ADDRESS OF NEXT COLUMN(ROW).
045710,001193:
045711,001194: 7361 07156 TC DOTSUB
045712,001195: 7362 52123 DXCH VBUF # MOVE GIVEN VECTOR BACK TO MPAC, SAVING Y
045713,001196: 7363 52155 DXCH MPAC # COMPONENT OF ANSWER IN VBUF +2.
045714,001197: 7364 52125 DXCH VBUF +2
045715,001198: 7365 52160 DXCH MPAC +3
045716,001199: 7366 52127 DXCH VBUF +4
045717,001200: 7367 52162 DXCH MPAC +5
045718,001201: 7370 30140 CA MATINC # FORM ADDRESS OF LAST COLUMN OR ROW.
045719,001202: 7371 26116 ADS ADDRWD
045720,001203:
045721,001204: 7372 07156 TC DOTSUB
045722,001205: 7373 52134 DXCH BUF2 # ANSWER NOW COMPLETE. PUT COMPONENTS INTO
045723,001206: 7374 52155 DXCH MPAC # PROPER MPAC REGISTERS.
045724,001207: 7375 52162 DXCH MPAC +5
045725,001208: 7376 52125 DXCH VBUF +2
045726,001209: 7377 52160 DXCH MPAC +3
045727,001210: 7400 16061 TCF DANZIG # EXIT.
045728,001211:
![]() |
Page 1040 |
045730,001213: # VXSC - VECTOR TIMES SCALAR.
045731,001214:
045732,001215: 7401 10163 VXSC CCS MODE # TEST PRESENT MODE.
045733,001216: 7402 17430 TCF DVXSC # SEPARATE ROUTINE WHEN SCALAR IS IN MPAC.
045734,001217: 7403 17430 TCF DVXSC
045735,001218:
045736,001219: 7404 07107 VVXSC TC DMPSUB # COMPUTE X COMPONENT
045737,001220: 7405 07140 TC VROUND # AND ROUND IT.
045738,001221: 7406 52160 DXCH MPAC +3 # PUT Y COMPONENT INTO MPAC SAVING MPAC IN
045739,001222: 7407 52155 DXCH MPAC # MPAC +3.
045740,001223: 7410 52160 DXCH MPAC +3
045741,001224:
045742,001225: 7411 07107 TC DMPSUB # DO SAME FOR Y AND Z COMPONENTS.
045743,001226: 7412 07140 TC VROUND
045744,001227: 7413 52162 DXCH MPAC +5
045745,001228: 7414 52155 DXCH MPAC
045746,001229: 7415 52162 DXCH MPAC +5
045747,001230:
045748,001231: 7416 07107 TC DMPSUB
045749,001232: 7417 07140 TC VROUND
045750,001233: 7420 52155 VROTATEX DXCH MPAC # EXIT USED TO RESTORE MPAC AFTER THIS
045751,001234: 7421 52162 DXCH MPAC +5 # TYPE OF ROTATION. CALLED BY VECTOR SHIFT
045752,001235: 7422 52160 DXCH MPAC +3 # RIGHT, V/SC, ETC.
045753,001236: 7423 52155 DXCH MPAC
045754,001237: 7424 16061 TCF DANZIG
045755,001238:
![]() |
Page 1041 |
045757,001240: # DP VECTOR PROJECTION ROUTINE.
045758,001241:
045759,001242: 7425 07154 VPROJ TC PREDOT # (MPAC.X)MPAC IS COMPUTED AND LEFT IN
045760,001243: 7426 44751 CS FOUR # MPAC. DO DOT AND FALL INTO DVXSC.
045761,001244: 7427 26116 ADS ADDRWD
045762,001245:
045763,001246: # VXSC WHEN SCALAR ARRIVES IN MPAC AND VECTOR IS AT X.
045764,001247:
045765,001248: 7430 00006 DVXSC EXTEND # SAVE SCALAR IN MPAC +3 AND GET X
045766,001249: 7431 30155 DCA MPAC # COMPONENT OF ANSWER.
045767,001250: 7432 52160 DXCH MPAC +3
045768,001251: 7433 07107 TC DMPSUB
045769,001252: 7434 07140 TC VROUND
045770,001253:
045771,001254: 7435 34752 CAF TWO # ADVANCE ADDRWD TO Y COMPONENT OF X.
045772,001255: 7436 26116 ADS ADDRWD
045773,001256: 7437 00006 EXTEND
045774,001257: 7440 30160 DCA MPAC +3 # PUT SCALAR BACK INTO MPAC AND SAVE
045775,001258: 7441 52155 DXCH MPAC # X RESULT IN MPAC +5.
045776,001259: 7442 52162 DXCH MPAC +5
045777,001260: 7443 07107 TC DMPSUB
045778,001261: 7444 07140 TC VROUND
045779,001262:
045780,001263: 7445 34752 CAF TWO
045781,001264: 7446 26116 ADS ADDRWD # TO Z COMPONENT.
045782,001265: 7447 52160 DXCH MPAC +3 # BRING SCALAR BACK, PUTTING Y RESULT IN
045783,001266: 7450 52155 DXCH MPAC # THE PROPER PLACE.
045784,001267: 7451 52160 DXCH MPAC +3
045785,001268: 7452 07107 TC DMPSUB
045786,001269: 7453 07140 TC VROUND
045787,001270:
045788,001271: 7454 52155 DXCH MPAC # PUT Z COMPONENT IN PROPER PLACE, ALSO
045789,001272: 7455 52162 DXCH MPAC +5 # POSITIONING X.
045790,001273: 7456 52155 DXCH MPAC
045791,001274:
045792,001275: 7457 16521 TCF VMODE # MODE HAS CHANGED TO VECTOR.
045793,001276:
![]() |
Page 1042 |
045795,001278: # THE VECTOR CROSS PRODUCT ROUTINE CALCULATES (X M -X M ,X M -X M ,X M -X M ) WHERE M IS THE VECTOR IN
045796,001279: # 3 2 2 3 1 3 3 1 2 1 1 2
045797,001280: # MPAC AND X THE VECTOR AT THE GIVEN ADDRESS.
045798,001281:
045799,001282: 7460 00006 VXV EXTEND
045800,001283: 7461 30162 DCA MPAC +5 # FORM UP M3X1, LEAVING M1 IN VBUF.
045801,001284: 7462 52155 DXCH MPAC
045802,001285: 7463 52123 DXCH VBUF
045803,001286: 7464 07107 TC DMPSUB # BY X1.
045804,001287:
045805,001288: 7465 00006 EXTEND
045806,001289: 7466 40160 DCS MPAC +3 # CALCULATE -X1M2, SAVING X1M3 IN VBUF +2.
045807,001290: 7467 52155 DXCH MPAC
045808,001291: 7470 52125 DXCH VBUF +2
045809,001292: 7471 07107 TC DMPSUB
045810,001293:
045811,001294: 7472 34752 CAF TWO # ADVANCE ADDRWD TO X2.
045812,001295: 7473 26116 ADS ADDRWD
045813,001296: 7474 00006 EXTEND
045814,001297: 7475 40162 DCS MPAC +5 # PREPARE TO GET -X2M3, SAVING -X1M2 IN
045815,001298: 7476 52155 DXCH MPAC # MPAC +5.
045816,001299: 7477 52162 DXCH MPAC +5
045817,001300: 7500 07107 TC DMPSUB
045818,001301:
045819,001302: 7501 00006 EXTEND
045820,001303: 7502 30123 DCA VBUF # GET X2M1, SAVING -X2M3 IN VBUF +4.
045821,001304: 7503 52155 DXCH MPAC
045822,001305: 7504 52127 DXCH VBUF +4
045823,001306: 7505 07107 TC DMPSUB
045824,001307:
045825,001308: 7506 34752 CAF TWO # ADVANCE ADDRWD TO X3.
045826,001309: 7507 26116 ADS ADDRWD
045827,001310: 7510 00006 EXTEND
045828,001311: 7511 40123 DCS VBUF # GET -X3M1, ADDING X2M1 TO MPAC +5 TO
045829,001312: 7512 52155 DXCH MPAC # COMPLETE THE Z COMPONENT OF THE ANSWER.
045830,001313: 7513 20162 DAS MPAC +5
045831,001314:
045832,001315: 7514 00006 EXTEND
045833,001316: 7515 17517 BZF +2
045834,001317: 7516 07011 TC OVERFLWZ
045835,001318:
045836,001319: 7517 07107 TC DMPSUB
045837,001320: 7520 52125 DXCH VBUF +2 # MOVE X1M3 TO MPAC +3 SETTING UP FOR X3M2
045838,001321: 7521 52160 DXCH MPAC +3 # AND ADD -X3M1 TO MPAC +3 TO COMPLETE THE
045839,001322: 7522 52155 DXCH MPAC # Y COMPONENT OF THE RESULT.
045840,001323: 7523 20160 DAS MPAC +3
045841,001324:
045842,001325: 7524 00006 EXTEND
045843,001326: 7525 17527 BZF +2
![]() |
Page 1043 |
045845,001328: 7526 07014 TC OVERFLWY
045846,001329:
045847,001330: 7527 07107 TC DMPSUB
045848,001331: 7530 52127 DXCH VBUF +4 # GO ADD -X2M3 TO X3M2 TO COMPLETE THE X
045849,001332: 7531 17000 TCF ENDVXV # COMPONENT (TAIL END OF DAD).
045850,001333:
045851,001334: # THE MPACVBUF SUBROUTINE SAVES THE VECTOR IN MPAC IN VBUF WITHOUT CLOBBERING MPAC.
045852,001335:
045853,001336: 7532 00006 MPACVBUF EXTEND # CALLED BY MXV, VXM, AND UNIT.
045854,001337: 7533 30155 DCA MPAC
045855,001338: 7534 52123 DXCH VBUF
045856,001339: 7535 00006 EXTEND
045857,001340: 7536 30160 DCA MPAC +3
045858,001341: 7537 52125 DXCH VBUF +2
045859,001342: 7540 00006 EXTEND
045860,001343: 7541 30162 DCA MPAC +5
045861,001344: 7542 52127 DXCH VBUF +4
045862,001345: 7543 00002 TC Q # RETURN TO CALLER.
045863,001346:
045864,001347: # DOUBLE PRECISION SIGN AGREE ROUTINE. ARRIVE WITH INPUT IN A+L. OUTPUT IS IN A + L.
045865,001348:
045866,001349: 7544 10000 ALSIGNAG CCS A # TEST UPPER PART.
045867,001350: 7545 17551 TCF UPPOS # IT IS POSITIVE
045868,001351: 7546 00002 TC Q # ZERO
045869,001352: 7547 17561 TCF UPNEG # NEGATIVE
045870,001353: 7550 00002 TC Q # ZERO
045871,001354:
045872,001355: 7551 56001 UPPOS XCH L # SAVE DECREMENTED UPPER PART.
045873,001356: 7552 64736 AD HALF
045874,001357: 7553 64736 AD HALF
045875,001358: 7554 54000 TS A # SKIPS ON OVERFLOW
045876,001359: 7555 17557 TCF +2
045877,001360: 7556 24001 INCR L # RESTORE UPPER TO ROIGINAL VALUE
045878,001361: 7557 56001 XCH L # SWAP A + L BACK.
045879,001362: 7560 00002 TC Q
045880,001363:
045881,001364: 7561 56001 UPNEG XCH L # SAVE COMPLEMENTED + DECREMENTED UPPER PT
045882,001365: 7562 64735 AD NEGMAX
045883,001366: 7563 67747 AD NEGONE
045884,001367: 7564 54000 TS A
045885,001368: 7565 17567 TCF +2 # DONT INCREMENT IF NO OVERFLOW.
045886,001369: 7566 24001 INCR L
045887,001370: 7567 56001 XCH L
045888,001371: 7570 40000 COM # MAKE NEGATIVE AGAIN.
045889,001372: 7571 00002 TC Q
045890,001373:
![]() |
Page 1044 |
045892,001375: # INTERPRETIVE INSTRUCTIONS WHOSE EXECUTION CONSISTS OF PRINCIPALLY CALLING SUBROUTINES.
045893,001376:
045894,001377: 7572 07107 DMP1 TC DMPSUB # DMP INSTRUCTION.
045895,001378: 7573 16061 TCF DANZIG
045896,001379:
045897,001380: 7574 07107 DMPR TC DMPSUB
045898,001381: 7575 07137 TC ROUNDSUB +1 # (C(A) = +0).
045899,001382: 7576 16061 TCF DANZIG
045900,001383:
045901,001384: 7577 00006 DDV EXTEND
045902,001385: 7600 50116 INDEX ADDRWD # MOVE DIVIDEND INTO BUF.
045903,001386: 7601 30001 DCA 0
045904,001387: 7602 17607 TCF BDDV +4
045905,001388:
045906,001389: 7603 00006 BDDV EXTEND # MOVE DIVISOR INTO MPAC SAVING MPAC, THE
045907,001390: 7604 50116 INDEX ADDRWD # DIVIDEND, IN BUF.
045908,001391: 7605 30001 DCA 0
045909,001392: 7606 52155 DXCH MPAC
045910,001393: 7607 52131 +4 DXCH BUF
045911,001394: 7610 34755 CAF ZERO # DIVIDE ROUTINES IN BANK 0.
045912,001395: 7611 54004 TS FBANK
045913,001396: 7612 12353 TCF DDV/BDDV
045914,001397:
045915,001398: 7613 30116 SETPD CA ADDRWD # MUST SET TO WORK AREA, OR EBANK TROUBLE.
045916,001399: 7614 54166 TS PUSHLOC
045917,001400: 7615 16063 TCF NOIBNKSW # NO FBANK SWITCH REQUIRED.
045918,001401:
045919,001402: 7616 34755 TSLC CAF ZERO # SHIFTING ROUTINES LOCATED IN BANK 00.
045920,001403: 7617 54004 TS FBANK
045921,001404: 7620 12172 TCF TSLC2
045922,001405:
045923,001406: 7621 36074 GSHIFT CAF LOW7 # USED AS MASK AT GENSHIFT. THIS PROCESSES
045924,001407: 7622 54004 TS FBANK # ANY SHIFT INSTRUCTION (EXCEPT TSLC) WITH
045925,001408: 7623 12214 TCF GENSHIFT # AN ADDRESS (ROUTINES IN BANK 0).
045926,001409:
![]() |
Page 1045 |
045928,001411: # THE FOLLOWING IS THE PROLOGUE TO V/SC. IF THE PRESENT MODE IS VECTOR, IT SAVES THE SCALAR AT X IN BUF
045929,001412: # AND CALLS THE V/SC ROUTINE IN BANK 0. IF THE PRESENT MODE IS SCALAR, IT MOVES THE VECTOR AT X INTO MPAC, SAVING
045930,001413: # THE SCALAR IN MPAC IN BUF BEFORE CALLING THE V/SC ROUTINE IN BANK 0.
045931,001414:
045932,001415: 7624 10163 V/SC CCS MODE
045933,001416: 7625 17636 TCF DV/SC # MOVE VECTOR INTO MPAC.
045934,001417: 7626 17636 TCF DV/SC
045935,001418:
045936,001419: 7627 00006 VV/SC EXTEND
045937,001420: 7630 50116 INDEX ADDRWD
045938,001421: 7631 30001 DCA 0
045939,001422: 7632 52131 V/SC1 DXCH BUF # IN BOTH CASES, VECTOR IS NOW IN MPAC AND
045940,001423: 7633 34755 CAF ZERO # SCALAR IN BUF.
045941,001424: 7634 54004 TS FBANK
045942,001425: 7635 12654 TCF V/SC2
045943,001426:
045944,001427: 7636 00006 DV/SC EXTEND
045945,001428: 7637 50116 INDEX ADDRWD
045946,001429: 7640 30003 DCA 2
045947,001430: 7641 52160 DXCH MPAC +3
045948,001431: 7642 00006 EXTEND
045949,001432: 7643 50116 INDEX ADDRWD
045950,001433: 7644 30005 DCA 4
045951,001434: 7645 52162 DXCH MPAC +5
045952,001435:
045953,001436: 7646 44753 CS ONE # CHANGE MODE TO VECTOR.
045954,001437: 7647 54163 TS MODE
045955,001438:
045956,001439: 7650 00006 EXTEND
045957,001440: 7651 50116 INDEX ADDRWD
045958,001441: 7652 30001 DCA 0
045959,001442: 7653 52155 DXCH MPAC
045960,001443: 7654 17632 TCF V/SC1 # FINISH PROLOGUE AT COMMON SECTION.
045961,001444:
![]() |
Page 1046 |
045963,001446: # SIGN AND COMPLEMENT INSTRUCTIONS.
045964,001447:
045965,001448: 7655 50116 SIGN INDEX ADDRWD # CALL COMP INSTRUCTION IF WORD AT X IS
045966,001449: 7656 10000 CCS 0 # NEGATIVE NON-ZERO.
045967,001450: 7657 16061 TCF DANZIG
045968,001451: 7660 17662 TCF +2
045969,001452: 7661 17670 TCF COMP # DO THE COMPLEMENT.
045970,001453:
045971,001454: 7662 50116 INDEX ADDRWD
045972,001455: 7663 10001 CCSL CCS 1
045973,001456: 7664 16061 TCF DANZIG
045974,001457: 7665 16061 TCF DANZIG
045975,001458: 7666 17670 TCF COMP
045976,001459: 7667 16061 TCF DANZIG
045977,001460: 7670 00006 COMP EXTEND # COMPLEMENT DP MPAC IN EVERY CASE.
045978,001461: 7671 40155 DCS MPAC
045979,001462: 7672 52155 DXCH MPAC
045980,001463:
045981,001464: 7673 10163 CCS MODE # EITHER COMPLEMENT MPAC +3 OR THE REST OF
045982,001465: 7674 17705 TCF DCOMP # THE VECTOR ACCUMULATOR.
045983,001466: 7675 17705 TCF DCOMP
045984,001467:
045985,001468: 7676 00006 EXTEND # VECTOR COMPLEMENT.
045986,001469: 7677 40160 DCS MPAC +3
045987,001470: 7700 52160 DXCH MPAC +3
045988,001471: 7701 00006 EXTEND
045989,001472: 7702 40162 DCS MPAC +5
045990,001473: 7703 52162 DXCH MPAC +5
045991,001474: 7704 16061 TCF DANZIG
045992,001475:
045993,001476: 7705 40156 DCOMP CS MPAC +2
045994,001477: 7706 54156 TS MPAC +2
045995,001478: 7707 16061 TCF DANZIG
045996,001479:
![]() |
Page 1047 |
045998,001481: # THE FOLLOWING SHORT SHIFT CODES REQUIRE NO ADDRESS WORD:
045999,001482: # 1. SR1 TO SR4 SCALAR SHIFT RIGHT.
046000,001483: # 2. SR1R TO SR4R SCALAR SHIFT RIGHT AND ROUND.
046001,001484: # 3. SL1 TO SL4 SCALAR SHIFT LEFT.
046002,001485: # 4. SL1R TO SL4R SCALAR SHIFT LEFT AND ROUND.
046003,001486: # 5. VSR1 TO VSR8 VECTOR SHIFT RIGHT (ALWAYS ROUNDS).
046004,001487: # 6. VSL1 TO VSL8 VECTOR SHIFT LEFT (NEVER ROUNDS).
046005,001488: # THE FOLLOWING CODES REQUIRE AN ADDRESS WHICH MAY BE INDEXED:*
046006,001489: # 1. SR SCALAR SHIFT RIGHT.
046007,001490: # 2. SRR SCALAR SHIFT RIGHT AND ROUND.
046008,001491: # 3. SL SCALAR SHIFT LEFT.
046009,001492: # 4. SLR SCALAR SHIFT LEFT AND ROUND.
046010,001493: # 5. VSR VECTOR SHIFT RIGHT.
046011,001494: # 6. VSL VECTOR SHIFT LEFT.
046012,001495: # * IF THE ADDRESS IS INDEXED, AND THE INDEX MODIFICATION RESULTS IN A NEGATIVE SHIFT COUNT, A SHIFT OF THE
046013,001496: # ABSOLUTE VALUE OF THE COUNT IS DONE IN THE OPPOSITE DIRECTION.
046014,001497:
046015,001498: 00,2017 BANK 00
046016,001499:
046017,001500: 00,2017 COUNT* $$/INTER
046018,001501: 00,2017 36242 SHORTT CAF SIX # SCALAR SHORT SHIFTS COME HERE. THE SHIFT
046019,001502: 00,2020 70020 MASK CYR # COUNT-1 IS NOW IN BITS 2-3 OF CYR. THE
046020,001503: 00,2021 54021 TS SR # ROUNDING BIT IS IN BIT1 AT THIS POINT.
046021,001504:
046022,001505: 00,2022 10020 CCS CYR # SEE IF RIGHT OR LEFT SHIFT DESIRED.
046023,001506: 00,2023 12101 TCF TSSL # SHIFT LEFT.
046024,001507:
046025,001508: 00,2024 00024 SRDDV DEC 20 B-14 # MPTEMP SETTING FOR SR BEFORE DDV.
046026,001509:
046027,001510: 00,2025 50021 TSSR INDEX SR # GET SHIFTING BIT.
046028,001511: 00,2026 34736 CAF BIT14
046029,001512: 00,2027 54135 TS MPTEMP
046030,001513:
046031,001514: 00,2030 10020 CCS CYR # SEE IF A ROUND IS DESIRED.
046032,001515: 00,2031 02050 RIGHTR TC MPACSRND # YES - SHIFT RIGHT AND ROUND.
046033,001516: 00,2032 16060 TCF NEWMODE # SET MODE TO DP (C(A) = 0).
046034,001517: 00,2033 30135 MPACSHR CA MPTEMP # DO A TRIPLE PRECISION SHIFT RIGHT.
046035,001518: 00,2034 00006 EXTEND
046036,001519: 00,2035 70156 MP MPAC +2
046037,001520: 00,2036 54156 +3 TS MPAC +2 # (EXIT FROM SQRT AND ABVAL).
046038,001521: 00,2037 30135 CA MPTEMP
046039,001522: 00,2040 00006 EXTEND
046040,001523: 00,2041 70154 MP MPAC # SHIFT MAJOR PART INTO A,L AND PLACE IN
![]() |
Page 1048 |
046042,001525: 00,2042 52155 DXCH MPAC # MPAC,+1.
046043,001526: 00,2043 30135 CA MPTEMP
046044,001527: 00,2044 00006 EXTEND
046045,001528: 00,2045 70001 MP L # ORIGINAL C(MPAC +1).
046046,001529: 00,2046 20156 DAS MPAC +1 # GUARANTEED NO OVERFLOW.
046047,001530: 00,2047 16061 TCF DANZIG
046048,001531:
046049,001532: # MPAC SHIFT RIGHT AND ROUND SUBROUTINES.
046050,001533:
046051,001534: 00,2050 30156 MPACSRND CA MPAC +2 # WE HAVE TO DO ALL THREE MULTIPLIES SINCE
046052,001535: 00,2051 00006 EXTEND # MPAC +1 AND MPAC +2 MIGHT HAVE SIGN
046053,001536: 00,2052 70135 MP MPTEMP # DISAGREEMENT WITH A SHIFT RIGHT OF 1.
046054,001537: 00,2053 56155 XCH MPAC +1
046055,001538: 00,2054 00006 EXTEND
046056,001539: 00,2055 70135 MP MPTEMP
046057,001540: 00,2056 56155 XCH MPAC +1 # TRIAL MINOR PART.
046058,001541: 00,2057 60001 AD L
046059,001542:
046060,001543: 00,2060 60000 VSHR2 DOUBLE # (FINISH VECTOR COMPONENT SHIFT RIGHT
046061,001544: 00,2061 54156 TS MPAC +2 # AND ROUND.
046062,001545: 00,2062 12064 TCF +2
046063,001546: 00,2063 26155 ADS MPAC +1 # GUARANTEED NO OVERFLOW.
046064,001547:
046065,001548: 00,2064 34755 CAF ZERO
046066,001549: 00,2065 54156 TS MPAC +2
046067,001550: 00,2066 56154 XCH MPAC # SETTING TO ZERO SO FOLLOWING DAS WORKS.
046068,001551: 00,2067 00006 EXTEND
046069,001552: 00,2070 70135 MP MPTEMP
046070,001553: 00,2071 20155 DAS MPAC # AGAIN NO OVERFLOW.
046071,001554: 00,2072 00002 TC Q
046072,001555:
046073,001556: 00,2073 30135 VSHRRND CA MPTEMP # ENTRY TO SHIFT RIGHT AND ROUND MPAC WHEN
046074,001557: 00,2074 00006 EXTEND # MPAC CONTAINS A VECTOR COMPONENT.
046075,001558: 00,2075 70155 MP MPAC +1
046076,001559: 00,2076 54155 TS MPAC +1
046077,001560: 00,2077 56001 XCH L
046078,001561: 00,2100 12060 TCF VSHR2 # GO ADD ONE IF NECESSARY AND FINISH.
046079,001562:
![]() |
Page 1049 |
046081,001564: # ROUTINE FOR SHORT SCALAR SHIFT LEFT (AND MAYBE ROUND).
046082,001565:
046083,001566: 00,2101 30021 TSSL CA SR # GET SHIFT COUNT FOR SR.
046084,001567: 00,2102 54135 +1 TS MPTEMP
046085,001568:
046086,001569: 00,2103 00006 +2 EXTEND # ENTRY HERE FROM SL FOR SCALARS.
046087,001570: 00,2104 30156 DCA MPAC +1 # SHIFTING LEFT ONE PLACE AT A TIME IS
046088,001571: 00,2105 20156 DAS MPAC +1 # FASTER THAN DOING THE WHOLE SHIFT WITH
046089,001572: 00,2106 60154 AD MPAC # MULTIPLIES ASSUMING THAT FREQUENCY OF
046090,001573: 00,2107 60154 AD MPAC # SHIFT COUNTS GOES DOWN RAPIDLY AS A
046091,001574: 00,2110 54154 TS MPAC # FUNCTION OF THEIR MAGNITUDE.
046092,001575: 00,2111 12113 TCF +2
046093,001576: 00,2112 54121 TS OVFIND # OVERFLOW. (LEAVES OVERFLOW-CORRECTED
046094,001577: # RESULT ANYWAY).
046095,001578: 00,2113 10135 CCS MPTEMP # LOOP ON DECREMENTED SHIFT COUNT.
046096,001579: 00,2114 12102 TCF TSSL +1
046097,001580:
046098,001581: 00,2115 10020 CCS CYR # SEE IF ROUND WANTED.
046099,001582: 00,2116 07136 ROUND TC ROUNDSUB # YES - ROUND AND EXIT.
046100,001583: 00,2117 16061 TCF DANZIG # SL LEAVES A ZERO IN CYR FOR NO ROUND.
046101,001584: 00,2120 16061 TCF DANZIG # NO - EXIT IMMEDIATL
046102,001585:
![]() |
Page 1050 |
046104,001587: # VECTOR SHIFTING ROUTINES.
046105,001588:
046106,001589: 00,2121 34757 SHORTV CAF LOW3 # SAVE 3 BIT SHIFT COUNT - 1 WITHOUT
046107,001590: 00,2122 70020 MASK CYR # EDITING CYR.
046108,001591: 00,2123 54135 TS MPTEMP
046109,001592: 00,2124 10020 CCS CYR # SEE IF LEFT OR RIGHT SHIFT.
046110,001593: 00,2125 12145 TCF VSSL # VECTOR SHIFT LEFT.
046111,001594: 00,2126 00176 OCT176 OCT 176 # USED IN PROCESSED SHIFTS WITH - COUNT.
046112,001595:
046113,001596: 00,2127 50135 VSSR INDEX MPTEMP # (ENTRY FROM SR). PICK UP SHIFTING BIT.
046114,001597: 00,2130 34736 CAF BIT14 # MPTEMP CONTAINS THE SHIFT COUNT - 1.
046115,001598: 00,2131 54135 TS MPTEMP
046116,001599: 00,2132 02073 TC VSHRRND # SHIFT X COMPONENT.
046117,001600:
046118,001601: 00,2133 52155 DXCH MPAC # SWAP X AND Y COMPONENTS.
046119,001602: 00,2134 52160 DXCH MPAC +3
046120,001603: 00,2135 52155 DXCH MPAC
046121,001604: 00,2136 02073 TC VSHRRND # SHIFT Y COMPONENT.
046122,001605:
046123,001606: 00,2137 52155 DXCH MPAC # SWAP Y AND Z COMPONENTS.
046124,001607: 00,2140 52162 DXCH MPAC +5
046125,001608: 00,2141 52155 DXCH MPAC
046126,001609: 00,2142 02073 TC VSHRRND # SHIFT Z COMPONENT.
046127,001610:
046128,001611: 00,2143 17420 TCF VROTATEX # RESTORE COMPONENTS TO PROPER PLACES.
046129,001612:
![]() |
Page 1051 |
046131,001614: # VECTOR SHIFT LEFT - DONE ONE PLACE AT A TIME.
046132,001615:
046133,001616: 00,2144 54135 -1 TS MPTEMP # SHIFTING LOOP.
046134,001617:
046135,001618: 00,2145 00006 VSSL EXTEND
046136,001619: 00,2146 30155 DCA MPAC
046137,001620: 00,2147 20155 DAS MPAC
046138,001621: 00,2150 00006 EXTEND
046139,001622: 00,2151 12153 BZF +2
046140,001623: 00,2152 07017 TC OVERFLOW
046141,001624:
046142,001625: 00,2153 00006 EXTEND
046143,001626: 00,2154 30160 DCA MPAC +3
046144,001627: 00,2155 20160 DAS MPAC +3
046145,001628: 00,2156 00006 EXTEND
046146,001629: 00,2157 12161 BZF +2
046147,001630: 00,2160 07014 TC OVERFLWY
046148,001631:
046149,001632: 00,2161 00006 EXTEND
046150,001633: 00,2162 30162 DCA MPAC +5
046151,001634: 00,2163 20162 DAS MPAC +5
046152,001635: 00,2164 00006 EXTEND
046153,001636: 00,2165 12167 BZF +2
046154,001637: 00,2166 07011 TC OVERFLWZ
046155,001638:
046156,001639: 00,2167 10135 CCS MPTEMP # LOOP ON DECREMENTED SHIFT COUNTER.
046157,001640: 00,2170 12144 TCF VSSL -1
046158,001641: 00,2171 16061 TCF DANZIG # EXIT.
046159,001642:
![]() |
Page 1052 |
046161,001644: # TSLC - TRIPLE SHIFT LEFT AND COUNT. SHIFTS MPAC LEFT UNTIL GREATER THAN .5 IN MAGNITUDE, LEAVING
046162,001645: # THE COMPLEMENT OF THE NUMBER OF SHIFTS REQUIRED IN X.
046163,001646:
046164,001647: 00,2172 54135 TSLC2 TS MPTEMP # START BY ZEROING SHIFT COUNT (IN A NOW).
046165,001648: 00,2173 06723 TC BRANCH # EXIT WITH NO SHIFTING IF ARGUMENT ZERO.
046166,001649: 00,2174 12176 TCF +2
046167,001650: 00,2175 12212 TCF ENDTSLC # STORES ZERO SHIFT COUNT IN THIS CASE.
046168,001651:
046169,001652: 00,2176 07257 TC TPAGREE # MAY CAUSE UPSHIFT OF ONE EXTRA PLACE.
046170,001653:
046171,001654: 00,2177 30154 CA MPAC # BEGIN NORMALIZATION LOOP.
046172,001655: 00,2200 12207 TCF TSLCTEST
046173,001656:
046174,001657: 00,2201 24135 TSLCLOOP INCR MPTEMP # INCREMENT SHIFT COUNTER.
046175,001658: 00,2202 00006 EXTEND
046176,001659: 00,2203 30156 DCA MPAC +1
046177,001660: 00,2204 20156 DAS MPAC +1
046178,001661: 00,2205 60154 AD MPAC
046179,001662: 00,2206 26154 ADS MPAC
046180,001663: 00,2207 60000 TSLCTEST DOUBLE # SEE IF (ANOTHER) SHIFT IS REQUIRED.
046181,001664: 00,2210 54000 OVSK
046182,001665: 00,2211 12201 TCF TSLCLOOP # YES - INCREMENT COUNT AND SHIFT AGAIN.
046183,001666:
046184,001667: 00,2212 40135 ENDTSLC CS MPTEMP
046185,001668: 00,2213 16623 TCF STORE1 # STORE SHIFT COUNT AND RETURN TO DANZIG.
046186,001669:
![]() |
Page 1053 |
046188,001671: # THE FOLLOWING ROUTINES PROCESSES THE GENERAL SHIFT INSTRUCTIONS SR, SRR, SL, AND SLR.
046189,001672: # THE GIVEN ADDRESS IS DECODED AS FOLLOWS:
046190,001673: # BITS 1-7 SHIFT COUNT (SUBADDRESS) LESS THAN 125 DECIMAL.
046191,001674: # BIT 8 PSEUDO SIGN BIT (DETECTS CHANGE IN SIGN IN INDEXED SHIFTS).
046192,001675: # BIT 9 0 FOR LEFT SHIFT, AND 1 FOR RIGHT SHIFT.
046193,001676: # BIT 10 1 FOR TERMINAL ROUND ON SCALAR SHIFTS, 0 OTHERWISE.
046194,001677: # BITS 11-13 0.
046195,001678: # BIT 14 1.
046196,001679: # BIT 15 0.
046197,001680: # THE ABOVE ENCODING IS DONE BY THE YUL SYSTEM.
046198,001681:
046199,001682: 00,2214 70116 GENSHIFT MASK ADDRWD # GET SHIFT COUNT, TESTING FOR ZERO.
046200,001683: 00,2215 10000 CCS A # (ARRIVES WITH C(A) = LOW7).
046201,001684: 00,2216 12224 TCF GENSHFT2 # IF NON-ZERO, PROCEED WITH DECREMENTED CT
046202,001685:
046203,001686: 00,2217 34742 CAF BIT10 # ZERO SHIFT COUNT. NO SHIFTS NEEDED BUT
046204,001687: 00,2220 70116 MASK ADDRWD # WE MIGHT HAVE TO ROUND MPAC ON SLR AND
046205,001688: 00,2221 10000 CCS A # SRR (SCALAR ONLY).
046206,001689: 00,2222 07136 TC ROUNDSUB
046207,001690: 00,2223 16061 TCF DANZIG
046208,001691:
046209,001692: 00,2224 54135 GENSHFT2 TS MPTEMP # DECREMENTED SHIFT COUNT TO MPTEMP.
046210,001693: 00,2225 34744 CAF BIT8 # TEST MEANING OF LOW SEVEN BIT COUNT IN
046211,001694: 00,2226 00006 EXTEND # MPTEMP NOW.
046212,001695: 00,2227 70116 MP ADDRWD
046213,001696: 00,2230 76245 MASK LOW2 # JUMPS ON SHIFT DIRECTION (BIT8) AND
046214,001697: 00,2231 50000 INDEX A
046215,001698: 00,2232 12233 TCF +1 # ORIGINAL SHIFT DIRECTION (BIT 9).
046216,001699: 00,2233 12332 TCF RIGHT- # NEGATIVE SHIFT COUNT FOR SL OR SLR.
046217,001700: 00,2234 12342 TCF LEFT # SL OR SLR.
046218,001701: 00,2235 12336 TCF LEFT- # NEGATIVE SHIFT COUNT WITH SR OR SRR.
046219,001702:
![]() |
Page 1054 |
046221,001704: # GENERAL SHIFT RIGHT.
046222,001705:
046223,001706: 00,2236 10163 RIGHT CCS MODE # SEE IF VECTOR OR SCALAR.
046224,001707: 00,2237 12277 TCF GENSCR
046225,001708: 00,2240 12277 TCF GENSCR
046226,001709:
046227,001710: 00,2241 30135 CA MPTEMP # SEE IF SHIFT COUNT LESS THAN 14D.
046228,001711: 00,2242 63733 VRIGHT2 AD NEG12
046229,001712: 00,2243 00006 EXTEND
046230,001713: 00,2244 62127 BZMF VSSR # IF SO, BRANCH AND SHIFT IMMEDIATELY.
046231,001714:
046232,001715: 00,2245 67747 AD NEGONE # IF NOT, REDUCE MPTEMP BY A TOTAL OF 14,
046233,001716: 00,2246 54135 TS MPTEMP # AND DO A SHIFT RIGHT AND ROUND BY 14.
046234,001717: 00,2247 34755 CAF ZERO # THE ROUND AT THIS STAGE MAY INTRODUCE A
046235,001718: 00,2250 54001 TS L # ONE BIT ERROR IN A SHIFT RIGHT 15D.
046236,001719: 00,2251 56154 XCH MPAC
046237,001720: 00,2252 56155 XCH MPAC +1
046238,001721: 00,2253 02272 TC SETROUND # X COMPONENT NOW SHIFTED, SO MAKE UP THE
046239,001722: 00,2254 20155 DAS MPAC # ROUNDING QUANTITY (0 IN A AND 0 OR +-1
046240,001723: # IN L).
046241,001724: 00,2255 56157 XCH MPAC +3 # REPEAT THE ABOVE PROCESS FOR Y AND Z.
046242,001725: 00,2256 56160 XCH MPAC +4
046243,001726: 00,2257 02272 TC SETROUND
046244,001727: 00,2260 20160 DAS MPAC +3 # NO OVERFLOW ON THESE ADDS.
046245,001728:
046246,001729: 00,2261 56161 XCH MPAC +5
046247,001730: 00,2262 56162 XCH MPAC +6
046248,001731: 00,2263 02272 TC SETROUND
046249,001732: 00,2264 20162 DAS MPAC +5
046250,001733:
046251,001734: 00,2265 10135 CCS MPTEMP # SEE IF DONE, DOING FINAL DECREMENT.
046252,001735: 00,2266 54135 TS MPTEMP
046253,001736: 00,2267 12242 TCF VRIGHT2
046254,001737: 00,2270 04604 BIASLO DEC .2974 B-1 # SQRT CONSTANT
046255,001738:
046256,001739: 00,2271 16061 TCF DANZIG
046257,001740:
046258,001741: 00,2272 60000 SETROUND DOUBLE # MAKES UP ROUNDING QUANTITY FROM ARRIVING
046259,001742: 00,2273 54156 TS MPAC +2 # C(A). L IS ZERO INITIALLY.
046260,001743: 00,2274 34755 CAF ZERO
046261,001744: 00,2275 56001 XCH L
046262,001745: 00,2276 00002 TC Q # RETURN AND DO THE DAS, RESETTING L TO 0.
046263,001746:
![]() |
Page 1055 |
046265,001748: # PROCESS SR AND SRR FOR SCALARS.
046266,001749:
046267,001750: 00,2277 30135 GENSCR CA MPTEMP # SEE IF THE ORIGINAL SHIFT COUNT WAS LESS
046268,001751: 00,2300 63733 +1 AD NEG12 # THAN 14D.
046269,001752: 00,2301 00006 EXTEND
046270,001753: 00,2302 62322 BZMF DOSSHFT # DO THE SHIFT IMMEDIATELY IF SO.
046271,001754:
046272,001755: 00,2303 67747 +4 AD NEGONE # IF NOT, DECREMENT SHIFT COUNT BY 14D AND
046273,001756: 00,2304 54135 TS MPTEMP # SHIFT MPAC RIGHT 14 PLACES.
046274,001757: 00,2305 34755 CAF ZERO
046275,001758: 00,2306 56154 XCH MPAC
046276,001759: 00,2307 56155 XCH MPAC +1
046277,001760: 00,2310 54156 TS MPAC +2
046278,001761: 00,2311 10135 CCS MPTEMP # SEE IF FINISHED, DO FINAL DECREMENT.
046279,001762: 00,2312 54135 TS MPTEMP
046280,001763: 00,2313 02300 TC GENSCR +1
046281,001764: 00,2314 22650 SLOPEHI DEC .5884 # SQRT CONSTANT.
046282,001765: 00,2315 34742 CAF BIT10 # FINISHED WITH SHIFT. SEE IF ROUND
046283,001766: 00,2316 70116 MASK ADDRWD # WANTED.
046284,001767: 00,2317 10000 CCS A
046285,001768: 00,2320 07136 TC ROUNDSUB
046286,001769: 00,2321 16061 TCF DANZIG # DO SO AND/OR EXIT.
046287,001770:
046288,001771: 00,2322 50135 DOSSHFT INDEX MPTEMP # PICK UP SHIFTING BIT.
046289,001772: 00,2323 34736 CAF BIT14
046290,001773: 00,2324 54135 TS MPTEMP
046291,001774: 00,2325 34742 CAF BIT10 # SEE IF TERMINAL ROUND DESIRED.
046292,001775: 00,2326 70116 MASK ADDRWD
046293,001776: 00,2327 10000 CCS A
046294,001777: 00,2330 12031 TCF RIGHTR # YES.
046295,001778: 00,2331 12033 TCF MPACSHR # JUST SHIFT RIGHT.
046296,001779:
![]() |
Page 1056 |
046298,001781: # PROCESS THE RIGHT- (SL(R) WITH A NEGATIVE COUNT), LEFT-, AND LEFT OPTIONS.
046299,001782:
046300,001783: 00,2332 40135 RIGHT- CS MPTEMP # GET ABSOLUTE VALUE - 1 OF SHIFT COUNT
046301,001784: 00,2333 62126 AD OCT176 # UNDERSTANDING THAT BIT8 (PSEUDO-SIGN)
046302,001785: 00,2334 54135 TS MPTEMP # WAS 1 INITIALLY.
046303,001786: 00,2335 12236 TCF RIGHT # DO NORMAL SHIFT RIGHT.
046304,001787:
046305,001788: 00,2336 42126 LEFT- CS OCT176 # SAME PROLOGUE TO LEFT FOR INDEXED RIGHT
046306,001789: 00,2337 60135 AD MPTEMP # SHIFTS WHOSE NET SHIFT COUNT IS NEGATIVE
046307,001790: 00,2340 40000 COM
046308,001791: 00,2341 54135 TS MPTEMP
046309,001792:
046310,001793: 00,2342 10163 LEFT CCS MODE # SINCE LEFT SHIFTING IS SONE ONE PLACE AT
046311,001794: 00,2343 12346 TCF GENSCL # A TIME, NO COMPARISON WITH 14 NEED BE
046312,001795: 00,2344 12346 TCF GENSCL # DONE. FOR SCALARS, SEE IF TERMINAL ROUND
046313,001796: 00,2345 12145 TCF VSSL # DESIRED. FOR VECTORS, SHIFT IMMEDIATELY.
046314,001797:
046315,001798: 00,2346 40116 GENSCL CS ADDRWD # PUT ROUNDING BIT (BIT 10 OF ADDRWD) INTO
046316,001799: 00,2347 00006 EXTEND # BIT 15 OF CYR WHERE THE ROUNDING BIT OF
046317,001800: 00,2350 74746 MP BIT6 # A SHORT SHIFT LEFT WOULD BE
046318,001801: 00,2351 54020 TS CYR
046319,001802: 00,2352 12103 TCF TSSL +2 # DO THE SHIFT.
046320,001803:
![]() |
Page 1057 |
046322,001805: # SCALAR DIVISION INSTRUCTIONS, DDV AND BDDV, ARE EXECUTED HERE. AT THIS POINT, THE DIVIDEND IS IN MPAC
046323,001806: # AND THE DIVISOR IN BUF.
046324,001807:
046325,001808: 00,2353 44753 DDV/BDDV CS ONE # INITIALIZATION.
046326,001809: 00,2354 54136 TS DVSIGN # +-1 FOR POSITIVE QUOTIENT - -0 FOR NEG.
046327,001810: 00,2355 54137 TS DVNORMCT # DIVIDEND NORMALIZATION COUNT.
046328,001811: 00,2356 54140 TS MAXDVSW # NEAR-ONE DIVIDE FLAG.
046329,001812:
046330,001813: 00,2357 10130 CCS BUF # FORCE BUF POSITIVE WITH THE MAJOR PART
046331,001814: 00,2360 12516 TCF BUFPOS # NON-ZERO.
046332,001815: 00,2361 12363 TCF +2
046333,001816: 00,2362 12531 TCF BUFNEG
046334,001817:
046335,001818: 00,2363 54156 BUFZERO TS MPAC +2 # ZERO THIS.
046336,001819: 00,2364 07257 TC TPAGREE # FORCE SIGN AGREEMENT BEFORE OVERFLOW
046337,001820:
046338,001821: 00,2365 10154 CCS MPAC # TEST TO SEE IF MPAC NON-ZERO. (TOO BIG)
046339,001822: 00,2366 12414 TCF OVF+ # MAJOR PART OF DIVIDEND IS POSITIVE NON-0
046340,001823: 00,2367 12371 TCF +2
046341,001824: 00,2370 12413 TCF OVF+ -1 # MAJOR PART OF DIVIDEND IS NEG. NON-ZERO
046342,001825:
046343,001826: 00,2371 56131 XCH BUF +1 # SHIFT DIVIDEND AND DIVISOR LEFT 14.
046344,001827: 00,2372 56130 XCH BUF
046345,001828: 00,2373 56155 XCH MPAC +1
046346,001829: 00,2374 56154 XCH MPAC
046347,001830: 00,2375 10130 CCS BUF # TRY AGAIN ON FORMER MINOR PART.
046348,001831: 00,2376 12422 TCF BUF+
046349,001832: 00,2377 12401 TCF +2 # OVERFLOW ON ZERO DIVISOR.
046350,001833: 00,2400 12416 TCF BUF-
046351,001834:
046352,001835: 00,2401 40154 CS MPAC # SIGN OF MPAC DETERMINES SIGN OF RESULT.
046353,001836: 00,2402 00006 SGNDVOVF EXTEND
046354,001837: 00,2403 62405 BZMF +2
046355,001838: 00,2404 24136 INCR DVSIGN # NEGMAX IN MPAC PERHAPS.
046356,001839: 00,2405 34733 DVOVF CAF POSMAX # ON DIVISION OVERFLOW OF ANY SORT, SET
046357,001840: 00,2406 54154 TS MPAC # SET DP MPAC TO +-POSMAX.
046358,001841: 00,2407 02630 TC FINALDV +3
046359,001842: 00,2410 34753 CAF ONE # SET OVERFLOW INDICATOR AND EXIT.
046360,001843: 00,2411 54121 TS OVFIND
046361,001844: 00,2412 06061 TC DANZIG
046362,001845:
046363,001846: 00,2413 24136 -1 INCR DVSIGN
046364,001847: 00,2414 40131 OVF+ CS BUF +1 # LOAD LOWER ORDER PART OF DIVISOR.
046365,001848: 00,2415 12402 TCF SGNDVOVF # GET SIGN OF RESULT.
046366,001849:
046367,001850: 00,2416 00006 BUF- EXTEND # IF BUF IS NEGATIVE, COMPLEMENT IT AND
046368,001851: 00,2417 40131 DCS BUF # MAINTAIN DVSIGN FOR FINAL QUOTIENT SIGN.
046369,001852: 00,2420 52131 DXCH BUF
046370,001853: 00,2421 24136 INCR DVSIGN # NOW -0.
046371,001854:
![]() |
Page 1058 |
046373,001856: 00,2422 10154 BUF+ CCS MPAC # FORCE MPAC POSITIVE, CHECKING FOR ZERO
046374,001857: 00,2423 12437 TCF MPAC+ # DIVIDEND IN THE PROCESS.
046375,001858: 00,2424 12426 TCF +2
046376,001859: 00,2425 12433 TCF MPAC-
046377,001860: 00,2426 10155 CCS MPAC +1
046378,001861: 00,2427 12437 TCF MPAC+
046379,001862: 00,2430 16061 TCF DANZIG # EXIT IMMEDIATELY ON ZERO DIVIDEND.
046380,001863: 00,2431 12433 TCF MPAC-
046381,001864: 00,2432 16061 TCF DANZIG
046382,001865:
046383,001866: 00,2433 00006 MPAC- EXTEND # FORCE MPAC POSITIVE AS BUF IN BUF-.
046384,001867: 00,2434 40155 DCS MPAC
046385,001868: 00,2435 52155 DXCH MPAC
046386,001869: 00,2436 24136 INCR DVSIGN # NOW +1 OR -0.
046387,001870:
![]() |
Page 1059 |
046389,001872: 00,2437 40154 MPAC+ CS MPAC # CHECK FOR DIVISION OVERFLOW. IF THE
046390,001873: 00,2440 67747 AD NEGONE # MAJOR PART OF THE DIVIDEND IS LESS THAN
046391,001874: 00,2441 60130 AD BUF # THE MAJOR PART OF THE DIVISOR BY AT
046392,001875: 00,2442 10000 CCS A # LEAST TWO, WE CAN PROCEED IMMEDIATELY
046393,001876: 00,2443 12505 TCF DVNORM # WITHOUT NORMALIZATION PRODUCING A DVMAX.
046394,001877: 00,2444 60001 -1/2+2 OCT 60001 # USED IN SQRTSUB.
046395,001878:
046396,001879: 00,2445 12446 TCF +1 # IF THE ABOVE DOES NOT HOLD, FORCE SIGN
046397,001880: 00,2446 34736 CAF HALF # AGREEMENT IN NUMERATOR AND DENOMINATOR
046398,001881: 00,2447 60000 DOUBLE # TO FACILITATE OVERFLOW AND NEAR-ONE
046399,001882: 00,2450 60155 AD MPAC +1 # CHECKING.
046400,001883: 00,2451 54155 TS MPAC +1
046401,001884: 00,2452 34755 CAF ZERO
046402,001885: 00,2453 64733 AD POSMAX
046403,001886: 00,2454 26154 ADS MPAC
046404,001887:
046405,001888: 00,2455 34736 CAF HALF # SAME FOR BUF.
046406,001889: 00,2456 60000 DOUBLE
046407,001890: 00,2457 60131 AD BUF +1
046408,001891: 00,2460 54131 TS BUF +1
046409,001892: 00,2461 34755 CAF ZERO
046410,001893: 00,2462 64733 AD POSMAX
046411,001894: 00,2463 26130 ADS BUF
046412,001895:
046413,001896: 00,2464 40154 CS MPAC # CHECK MAGNITUDE OF SIGN-CORRECTED
046414,001897: 00,2465 60130 AD BUF # OPERANDS.
046415,001898: 00,2466 10000 CCS A
046416,001899: 00,2467 12505 TCF DVNORM # DIVIDE OK - WILL NOT BECOME MAXDV CASE.
046417,001900: 00,2470 00133 LBUF2 ADRES BUF2
046418,001901: 00,2471 12405 TCF DVOVF # DIVISOR NOT LESS THAN DIVIDEND - OVF.
046419,001902:
046420,001903: 00,2472 54140 TS MAXDVSW # IF THE MAJOR PARTS OF THE DIVIDEND AND
046421,001904: 00,2473 40155 CS MPAC +1 # DIVISOR ARE EQUAL, A SPECIAL APPROXIMA-
046422,001905: 00,2474 60131 AD BUF +1 # TION IS USED (PROVIDED THE DIVISION IS
046423,001906: 00,2475 00006 EXTEND # POSSIBLE, OF COURSE).
046424,001907: 00,2476 62405 BZMF DVOVF
046425,001908: 00,2477 12505 TCF DVNORM # IF NO OVERFLOW.
046426,001909:
![]() |
Page 1060 |
046428,001911: 00,2500 00006 BUFNORM EXTEND # ADD -1 TO AUGMENT SHIFT COUNT AND SHIFT
046429,001912: 00,2501 24137 AUG DVNORMCT # LEFT ONE PLACE.
046430,001913: 00,2502 00006 EXTEND
046431,001914: 00,2503 30131 DCA BUF
046432,001915: 00,2504 20131 DAS BUF
046433,001916:
046434,001917: 00,2505 30130 DVNORM CA BUF # SEE IF DIVISOR NORMALIZED YET.
046435,001918: 00,2506 60000 DOUBLE
046436,001919: 00,2507 54000 OVSK
046437,001920: 00,2510 12500 TCF BUFNORM # NO - SHIFT LEFT ONE AND TRY AGAIN.
046438,001921:
046439,001922: 00,2511 52155 DXCH MPAC # CALL DIVIDEND NORMALIZATION SEQUENCE
046440,001923: 00,2512 50137 INDEX DVNORMCT # PRIOR TO DOING THE DIVIDE.
046441,001924: 00,2513 02565 TC MAXTEST
046442,001925:
046443,001926: 00,2514 54156 TS MPAC +2 # RETURNS WITH DIVISION DONE AND C(A) = 0.
046444,001927: 00,2515 16061 TCF DANZIG
046445,001928:
046446,001929: 00,2516 10000 BUFPOS CCS A
046447,001930: 00,2517 12422 TCF BUF+ # TO BUF+ IF BUF IS GREATER THAN +1.
046448,001931:
046449,001932: 00,2520 40131 CS BUF +1 # IF BUF IS +1, FORCING SIGN AGREEMENT
046450,001933: 00,2521 00006 EXTEND # MAY CAUSE BUF TO BECOME ZERO.
046451,001934: 00,2522 62422 BZMF BUF+ # BRANCH IF SIGNS AGREE.
046452,001935:
046453,001936: 00,2523 34736 CA HALF # SIGNS DISAGREE. FORCE AGREEMENT.
046454,001937: 00,2524 60000 +6 DOUBLE
046455,001938: 00,2525 26131 ADS BUF +1
046456,001939: 00,2526 34755 CA ZERO
046457,001940: 00,2527 54130 TS BUF
046458,001941: 00,2530 12363 TCF BUFZERO
046459,001942:
046460,001943: 00,2531 10000 BUFNEG CCS A
046461,001944: 00,2532 12416 TCF BUF- # TO BUF- IF BUF IS LESS THAN -1.
046462,001945:
046463,001946: 00,2533 30131 CA BUF +1 # IF BUF IS -1, FORCING SIGN AGREEMENT
046464,001947: 00,2534 00006 EXTEND # MAY CAUSE BUF TO BECOME ZERO.
046465,001948: 00,2535 62416 BZMF BUF- # BRANCH IF SIGNS AGREE.
046466,001949:
046467,001950: 00,2536 44736 CS HALF # SIGNS DISAGREE. FORCE AGREEMENT.
046468,001951: 00,2537 12524 TCF BUFPOS +6
046469,001952:
![]() |
Page 1061 |
046471,001954: # THE FOLLOWING ARE PROLOGUES TO SHIFT THE DIVIDEND ARRIVING IN A AND L BEFORE THE DIVIDE.
046472,001955:
046473,001956: 00,2540 22021 -21D LXCH SR # SPECIAL PROLOGUE FOR UNIT WHEN THE
046474,001957: 00,2541 00006 EXTEND # LENGTH OF THE ARGUMENT WAS NOT LESS THAN
046475,001958: 00,2542 74736 MP HALF # .5. IN THIS CASE, EACH COMPONENT MUST BE
046476,001959: 00,2543 56001 XCH L # SHIFTED RIGHT ONE TO PRODUCE A HALF-UNIT
046477,001960: 00,2544 60021 AD SR # VECTOR.
046478,001961: 00,2545 56001 XCH L
046479,001962: 00,2546 12571 TCF GENDDV +1 # WITH DP DIVIDEND IN A,L.
046480,001963:
046481,001964: 00,2547 20001 DDOUBL # PROLOGUE WHICH NORMALIZES THE DIVIDEND
046482,001965: 00,2550 20001 DDOUBL # WHEN IT IS KNOWN THAT NO DIVISION
046483,001966: 00,2551 20001 DDOUBL # OVERFLOW WILL OCCUR.
046484,001967: 00,2552 20001 DDOUBL
046485,001968: 00,2553 20001 DDOUBL
046486,001969: 00,2554 20001 DDOUBL
046487,001970: 00,2555 20001 DDOUBL
046488,001971: 00,2556 20001 DDOUBL
046489,001972: 00,2557 20001 DDOUBL
046490,001973: 00,2560 20001 DDOUBL
046491,001974: 00,2561 20001 DDOUBL
046492,001975: 00,2562 20001 DDOUBL
046493,001976: 00,2563 20001 DDOUBL
046494,001977: 00,2564 52155 DXCH MPAC
046495,001978:
046496,001979: 00,2565 10140 MAXTEST CCS MAXDVSW # 0 IF MAJORS MIGHT BE =, -1 OTHERWISE.
046497,001980: 00,2566 06552 BIASHI DEC .4192 B-1 # SQRT CONSTANTS
046498,001981:
046499,001982: 00,2567 12642 TCF MAXDV # CHECK TO SEE IF THEY ARE NOW EQUAL.
046500,001983:
![]() |
Page 1062 |
046502,001985: # THE FOLLOWING IS A GENERAL PURPOSE DOUBLE PRECISION DIVISION ROUTINE. IT DIVIDES MPAC BY BUF AND LEAVES
046503,001986: # THE RESULT IN MPAC. THE FOLLOWING CONDITIONS MUST BE SATISFIED:
046504,001987:
046505,001988: # 1. THE DIVISOR (BUF) MUST BE POSITIVE AND NOT LESS THAN .5.
046506,001989:
046507,001990: # 2. THE DIVIDEND (MPAC) MUST BE POSITIVE WITH THE MAJOR PART OF MPAC STRICTLY LESS THAN THAT OF BUF
046508,001991: # (A SPECIAL APPROXIMATION, MAXDV, IS USED WHEN THE MAJOR PARTS ARE EQUAL).
046509,001992:
046510,001993: # UNDERSTANDING THAT A/B = Q + S(R/B) WHERE S = 2(-14) AND Q AND R ARE QUOTIENT AND REMAINDER, RESPEC-
046511,001994: # TIVELY, THE FOLLOWING APPROXIMATION IS OBTAINED BY MULTIPLYING ABOVE AND BELOW BY C - SD AND NEGLECTING TERMS OF
046512,001995: # ORDER S-SQUARED (POSSIBLY INTRODUCING ERROR INTO THE LOW TWO BITS OF THE RESULT). SIGN AGREEMENT IS UNNECESSARY.
046513,001996:
046514,001997: # A + SB . (R - QD) A + SB
046515,001998: # ------ = Q + S(------) WHERE Q AND R ARE QUOTIENT AND REMAINDER OF ------ RESPECTIVELY.
046516,001999: # C + SD ( C ) C
046517,002000:
046518,002001: 00,2570 52155 GENDDV DXCH MPAC # WE NEED A AND B ONLY FOR FIRST DV.
046519,002002: 00,2571 00006 +1 EXTEND # (SPECIAL UNIT PROLOGUE ENTERS HERE).
046520,002003: 00,2572 10130 DV BUF # A NOW CONTAINS Q AND L, R.
046521,002004: 00,2573 52155 DXCH MPAC
046522,002005:
046523,002006: 00,2574 40154 CS MPAC # FORM DIVIDEND FOR MINOR PART OF RESULT.
046524,002007: 00,2575 00006 EXTEND
046525,002008: 00,2576 70131 MP BUF +1
046526,002009: 00,2577 60155 AD MPAC +1 # OVERFLOW AT THIS POINT IS POSITIVE SINCE
046527,002010: 00,2600 54000 OVSK # R IS POSITIVE IN EVERY CASE.
046528,002011: 00,2601 12606 TCF +5
046529,002012:
046530,002013: 00,2602 00006 EXTEND # OVERFLOW CAN BE REMOVED BY SUBTRACTING C
046531,002014: 00,2603 60130 SU BUF # (BUF) ONCE SINCE R IS ALWAYS LESS THAN C
046532,002015: 00,2604 24154 INCR MPAC # IN THIS CASE. INCR COMPENSATES SUBTRACT.
046533,002016: 00,2605 12610 TCF +DOWN # (SINCE C(A) IS STILL POSITIVE).
046534,002017:
046535,002018: 00,2606 00006 +5 EXTEND # C(A) CAN BE MADE LESS THAN C IN MAGNI-
046536,002019: 00,2607 62620 BZMF -UP # TUDE BY DIMINISHING IT BY C (SINCE C IS
046537,002020: # NOT LESS THAN .5) UNLESS C(A) = 0.
046538,002021:
![]() |
Page 1063 |
046540,002023: 00,2610 00006 +DOWN EXTEND
046541,002024: 00,2611 60130 SU BUF # IF POSITIVE, REDUCE ONLY IF NECESSARY
046542,002025: 00,2612 00006 EXTEND # SINCE THE COMPENSATING INCR MIGHT CAUSE
046543,002026: 00,2613 12616 BZF +3 # OVERFLOW.
046544,002027: 00,2614 00006 EXTEND # DONT SUBTRACT UNLESS RESULT IS POSITIVE
046545,002028: 00,2615 62624 BZMF ENDMAXDV # OR ZERO.
046546,002029:
046547,002030: 00,2616 24154 +3 INCR MPAC # KEEP SUBTRACT HERE AND COMPENSATE.
046548,002031: 00,2617 12625 TCF FINALDV
046549,002032:
046550,002033: 00,2620 00006 -UP EXTEND # IF ZERO, SET MINOR PART OF RESULT TO
046551,002034: 00,2621 12630 BZF FINALDV +3 # ZERO.
046552,002035:
046553,002036: 00,2622 00006 EXTEND # IF NEGATIVE, ADD C TO A, SUBTRACTING ONE
046554,002037: 00,2623 26154 DIM MPAC # TO COMPENSATE. DIM IS OK HERE SINCE THE
046555,002038: 00,2624 60130 ENDMAXDV AD BUF # MAJOR PART NEVER GOES NEGATIVE.
046556,002039:
![]() |
Page 1064 |
046558,002041: 00,2625 22007 FINALDV ZL # DO DV TO OBTAIN MINOR PART OF RESULT.
046559,002042: 00,2626 00006 EXTEND
046560,002043: 00,2627 10130 DV BUF
046561,002044: 00,2630 54155 +3 TS MPAC +1
046562,002045:
046563,002046: 00,2631 10136 CCS DVSIGN # LEAVE RESULT POSITIVE UNLESS C(DVSIGN)=
046564,002047: 00,2632 00002 TC Q # -0.
046565,002048: 00,2633 00002 TC Q
046566,002049: 00,2634 00002 TC Q
046567,002050:
046568,002051: 00,2635 00006 EXTEND
046569,002052: 00,2636 40155 DCS MPAC
046570,002053: 00,2637 52155 DXCH MPAC
046571,002054: 00,2640 34755 CAF ZERO # SO WE ALWAYS RETURN WITH C(A) = 0.
046572,002055: 00,2641 00002 TC Q
046573,002056:
![]() |
Page 1065 |
046575,002058: # IF THE MAJOR PARTS OF THE DIVISOR AND DIVIDEND ARE EQUAL, BUT THE MINOR PARTS ARE SUCH THAT THE
046576,002059: # DIVIDEND IS STRICTLY LESS THAN THE DIVISOR IN MAGNITUDE, THE FOLLOWING APPROXIMATION IS USED. THE ASSUMPTIONS
046577,002060: # ARE THE SAME AS THE GENERAL ROUTINE WITH THE ADDITION THAT SIGN AGREEMENT IS NECESSARY (B, C, & D POSITIVE).
046578,002061:
046579,002062: # C + SB . (C + B - D)
046580,002063: # ------ = 37777 + S(---------)
046581,002064: # C + SD ( C )
046582,002065:
046583,002066: # THE DIVISION MAY BE PERFORMED IMMEDIATELY SINCE B IS STRICTLY LESS THAN D AND C IS NOT LESS THAN .5.
046584,002067:
046585,002068: 00,2642 40154 MAXDV CS MPAC # SEE IF MAXDV CASE STILL HOLDS AFTER
046586,002069: 00,2643 60130 AD BUF # NORMALIZATION.
046587,002070: 00,2644 00006 EXTEND
046588,002071: 00,2645 12647 BZF +2
046589,002072: 00,2646 12570 TCF GENDDV # MPAC NOW LESS THAN BUF - DIVIDE AS USUAL
046590,002073:
046591,002074: 00,2647 34733 +2 CAF POSMAX # SET MAJOR PART OF RESULT.
046592,002075: 00,2650 54154 TS MPAC
046593,002076:
046594,002077: 00,2651 40131 CS BUF +1 # FORM DIVIDEND OF MINOR PART OF RESULT.
046595,002078: 00,2652 60155 AD MPAC +1
046596,002079: 00,2653 12624 TCF ENDMAXDV # GO ADD C AND DO DIVIDE, ATTACHING SIGN
046597,002080: # BEFORE EXITING.
046598,002081:
![]() |
Page 1066 |
046600,002083: # VECTOR DIVIDED BY SCALAR, V/SC, IS EXECUTED HERE. THE VECTOR IS NOW IN MPAC WITH SCALAR IN BUF.
046601,002084:
046602,002085: 00,2654 44753 V/SC2 CS ONE # INITIALIZE DIVIDEND NORMALIZATION COUNT
046603,002086: 00,2655 54137 TS DVNORMCT # AND DIVISION SIGN REGISTER.
046604,002087: 00,2656 54127 TS VBUF +5
046605,002088:
046606,002089: 00,2657 03010 TC VECAGREE # FORCE SIGN AGREEMENT IN VECTOR
046607,002090:
046608,002091: 00,2660 52131 DXCH BUF
046609,002092: 00,2661 07544 TC ALSIGNAG # SIGN AGREE BUF
046610,002093: 00,2662 52131 DXCH BUF
046611,002094: 00,2663 10130 CCS BUF # FORCE DIVISOR POSITIVE WITH MAJOR PART
046612,002095: 00,2664 12721 TCF /BUF+ # NON-ZERO (IF POSSIBLE).
046613,002096: 00,2665 12667 TCF +2
046614,002097: 00,2666 12715 TCF /BUF-
046615,002098:
046616,002099: 00,2667 56131 XCH BUF +1 # SHIFT VECTOR AND SCALAR LEFT 14.
046617,002100: 00,2670 56130 XCH BUF
046618,002101: 00,2671 56155 XCH MPAC +1
046619,002102: 00,2672 56154 XCH MPAC
046620,002103: 00,2673 00006 EXTEND # CHECK FOR OVERFLOW IN EACH CASE.
046621,002104: 00,2674 12676 BZF +2
046622,002105: 00,2675 12405 TCF DVOVF
046623,002106:
046624,002107: 00,2676 56160 XCH MPAC +4
046625,002108: 00,2677 56157 XCH MPAC +3
046626,002109: 00,2700 00006 EXTEND
046627,002110: 00,2701 12703 BZF +2
046628,002111: 00,2702 12405 TCF DVOVF
046629,002112:
046630,002113: 00,2703 56162 XCH MPAC +6
046631,002114: 00,2704 56161 XCH MPAC +5
046632,002115: 00,2705 00006 EXTEND
046633,002116: 00,2706 12710 BZF +2
046634,002117: 00,2707 12405 TCF DVOVF
046635,002118:
046636,002119: 00,2710 10130 CCS BUF
046637,002120: 00,2711 12721 TCF /BUF+
046638,002121: 00,2712 12405 TCF DVOVF # ZERO DIVISOR - OVERFLOW.
046639,002122: 00,2713 12715 TCF /BUF-
046640,002123: 00,2714 12405 TCF DVOVF
046641,002124:
046642,002125: 00,2715 00006 /BUF- EXTEND # ON NEGATIVE, COMPLEMENT BUF AND MAINTAIN
046643,002126: 00,2716 40131 DCS BUF # DVSIGN IN VBUF +5.
046644,002127: 00,2717 52131 DXCH BUF
046645,002128: 00,2720 24127 INCR VBUF +5
046646,002129:
![]() |
Page 1067 |
046648,002131: 00,2721 00006 /BUF+ EXTEND
046649,002132: 00,2722 30131 DCA BUF # LEAVE ABS(ORIG DIVISOR) IN BUF2
046650,002133: 00,2723 52134 DXCH BUF2 # FOR OVERFLOW TESTING
046651,002134: 00,2724 12732 TCF /NORM # NORMALIZE DIVISOR IN BUF.
046652,002135:
046653,002136: 00,2725 00006 /NORM2 EXTEND # IF LESS THAN .5, AUGMENT DVNORMCT AND
046654,002137: 00,2726 24137 AUG DVNORMCT # DOUBLE DIVISOR.
046655,002138: 00,2727 00006 EXTEND
046656,002139: 00,2730 30131 DCA BUF
046657,002140: 00,2731 20131 DAS BUF
046658,002141:
046659,002142: 00,2732 30130 /NORM CA BUF # SEE IF DIVISOR NORMALIZED.
046660,002143: 00,2733 60000 DOUBLE
046661,002144: 00,2734 54000 OVSK
046662,002145: 00,2735 12725 TCF /NORM2 # DOUBLE AND TRY AGAIN IF NOT.
046663,002146:
046664,002147: 00,2736 02750 TC V/SCDV # DO X COMPONENT DIVIDE.
046665,002148: 00,2737 52160 DXCH MPAC +3 # SUPPLY ARGUMENTS IN USUAL SEQUENCE.
046666,002149: 00,2740 52155 DXCH MPAC
046667,002150: 00,2741 52160 DXCH MPAC +3
046668,002151:
046669,002152: 00,2742 02750 TC V/SCDV # Y COMPONENT.
046670,002153: 00,2743 52162 DXCH MPAC +5
046671,002154: 00,2744 52155 DXCH MPAC
046672,002155: 00,2745 52162 DXCH MPAC +5
046673,002156:
046674,002157: 00,2746 02750 TC V/SCDV # Z COMPONENT.
046675,002158: 00,2747 17420 TCF VROTATEX # GO RE-ARRANGE COMPONENTS BEFORE EXIT.
046676,002159:
![]() |
Page 1068 |
046678,002161: # SUBROUTINE USED BY V/SC TO DIVIDE VECTOR COMPONENT IN MPAC,+1 BY THE SCALAR GIVEN IN BUF.
046679,002162:
046680,002163: 00,2750 30127 V/SCDV CA VBUF +5 # REFLECTS SIGN OF SCALAR.
046681,002164: 00,2751 54136 TS DVSIGN
046682,002165:
046683,002166: 00,2752 10154 CCS MPAC # FORCE MPAC POSITIVE, EXITING ON ZERO.
046684,002167: 00,2753 12767 TCF /MPAC+
046685,002168: 00,2754 12756 TCF +2
046686,002169: 00,2755 12763 TCF /MPAC-
046687,002170:
046688,002171: 00,2756 10155 CCS MPAC +1
046689,002172: 00,2757 12767 TCF /MPAC+
046690,002173: 00,2760 00002 TC Q
046691,002174: 00,2761 12763 TCF /MPAC-
046692,002175: 00,2762 00002 TC Q
046693,002176:
046694,002177: 00,2763 00006 /MPAC- EXTEND # USUAL COMPLEMENTING AND SETTING OF SIGN.
046695,002178: 00,2764 40155 DCS MPAC
046696,002179: 00,2765 52155 DXCH MPAC
046697,002180: 00,2766 24136 INCR DVSIGN
046698,002181:
046699,002182: 00,2767 44753 /MPAC+ CS ONE # INITIALIZE NEAR-ONE SWITCH.
046700,002183: 00,2770 54140 TS MAXDVSW
046701,002184:
046702,002185: 00,2771 40154 CS MPAC # CHECK POSSIBLE OVERFLOW.
046703,002186: 00,2772 60133 AD BUF2 # UNNORMALIZED INPUT DIVISOR.
046704,002187: 00,2773 10000 CCS A
046705,002188: 00,2774 13004 TCF DDVCALL # NOT NEAR-ONE
046706,002189: 00,2775 12777 TCF +2 # +0 IS JUST POSSIBLE
046707,002190: 00,2776 12405 TCF DVOVF # NO HOPE
046708,002191: 00,2777 54140 TS MAXDVSW # SIGNAL POSSIBLE NEAR-ONE CASE
046709,002192: 00,3000 40155 CS MPAC +1 # SEE IF DIVISION CAN BE DONE
046710,002193: 00,3001 60134 AD BUF2 +1
046711,002194: 00,3002 00006 EXTEND
046712,002195: 00,3003 62405 BZMF DVOVF
046713,002196:
046714,002197: 00,3004 52155 DDVCALL DXCH MPAC # CALL PRE-DIVIDE NORMALIZATION.
046715,002198: 00,3005 50137 INDEX DVNORMCT
046716,002199: 00,3006 12565 TCF MAXTEST
046717,002200:
![]() |
Page 1069 |
046719,002202: 00,3007 32506 SLOPELO DEC .8324
046720,002203:
046721,002204: 00,3010 56002 VECAGREE XCH Q # SAVE Q IN A
046722,002205: 00,3011 52155 DXCH MPAC
046723,002206: 00,3012 07544 TC ALSIGNAG # SIGNAGREE MPAC
046724,002207: 00,3013 52155 DXCH MPAC
046725,002208: 00,3014 52160 DXCH MPAC +3
046726,002209: 00,3015 07544 TC ALSIGNAG # SIGN AGREE MPAC +3
046727,002210: 00,3016 52160 DXCH MPAC +3
046728,002211: 00,3017 52162 DXCH MPAC +5
046729,002212: 00,3020 07544 TC ALSIGNAG # SIGNAGREE MPAC +5
046730,002213: 00,3021 52162 DXCH MPAC +5
046731,002214: 00,3022 00000 TC A
046732,002215:
![]() |
Page 1070 |
046734,002217: # THE FOLLOWING ROUTINE EXECUTES THE UNIT INSTRUCTION, WHICH TAKES THE UNIT OF THE VECTOR IN MPAC.
046735,002218:
046736,002219: 00,3023 03010 UNIT TC VECAGREE # FORCE SIGN AGREEMENT IN VECTOR
046737,002220: 00,3024 07532 TC MPACVBUF # SAVE ARGUMENT IN VBUF
046738,002221: 00,3025 34755 CAF ZERO # MUST SENSE OVERFLOW IN FOLLOWING DOT.
046739,002222: 00,3026 56121 XCH OVFIND
046740,002223: 00,3027 54141 TS TEM1
046741,002224: 00,3030 03317 TC VSQSUB # DOT MPAC WITH ITSELF.
046742,002225: 00,3031 30141 CA TEM1
046743,002226: 00,3032 56121 XCH OVFIND
046744,002227: 00,3033 00006 EXTEND
046745,002228: 00,3034 13036 BZF +2
046746,002229: 00,3035 12405 TCF DVOVF
046747,002230: 00,3036 00006 EXTEND
046748,002231: 00,3037 30155 DCA MPAC # LEAVE THE SQUARE OF THE LENGTH OF THE
046749,002232: 00,3040 50120 INDEX FIXLOC # ARGUMENT IN LVSQUARE.
046750,002233: 00,3041 52043 DXCH LVSQUARE
046751,002234:
046752,002235: 00,3042 03343 TC SQRTSUB # GO TAKE THE NORMALIZED SQUARE ROOT.
046753,002236:
046754,002237: 00,3043 10154 CCS MPAC # CHECK FOR UNIT OVERFLOW.
046755,002238: 00,3044 13051 TCF +5 # MPAC IS NOT LESS THAN .5 UNLESS
046756,002239: 00,3045 54001 TS L
046757,002240: 00,3046 50120 INDEX FIXLOC
046758,002241: 00,3047 52045 DXCH LV
046759,002242: 00,3050 12405 TCF DVOVF # INPUT TO SQRTSUB WAS 0.
046760,002243:
046761,002244: 00,3051 44317 CS FOURTEEN # SEE IF THE INPUT WAS SO SMALL THE THE
046762,002245: 00,3052 60135 AD MPTEMP # FIRST TWO REGISTERS OF THE SQUARE WERE 0
046763,002246: 00,3053 10000 CCS A
046764,002247: 00,3054 40000 COM # IF SO, SAVE THE NEGATIVE OF THE SHIFT
046765,002248: 00,3055 13133 TCF SMALL # COUNT -15D.
046766,002249:
046767,002250: 00,3056 13065 TCF LARGE # (THIS IS USUALLY THE CASE.)
046768,002251:
046769,002252: 00,3057 44761 CS THIRTEEN # IF THE SHIFT COUNT WAS EXACTLY 14, SET
046770,002253: 00,3060 54135 TS MPTEMP # THE PRE-DIVIDE NORM COUNT TO -13D.
046771,002254:
046772,002255: 00,3061 30154 CA MPAC # SHIFT THE LENGTH RIGHT 14 BEFORE STORING
046773,002256: 00,3062 54001 SMALL2 TS L # (SMALL EXITS TO THIS POINT).
046774,002257: 00,3063 34755 CAF ZERO
046775,002258: 00,3064 13112 TCF LARGE2 # GO TO STORE LENGTH AND PROCEED.
046776,002259:
046777,002260: 00,3065 10135 LARGE CCS MPTEMP # MOST ALL CASES COME HERE.
046778,002261: 00,3066 13074 TCF LARGE3 # SEE IF NO NORMALIZATION WAS REQUIRED BY
046779,002262:
046780,002263: 00,3067 42024 CS SRDDV # SQRT, AND IF SO, SET UP FOR A SHIFT
046781,002264: 00,3070 54135 TS MPTEMP # RIGHT 1 BEFORE DIVIDING TO PRODUCE
046782,002265: 00,3071 00006 EXTEND # THE DESIRED HALF UNIT VECTOR.
046783,002266: 00,3072 30155 DCA MPAC
![]() |
Page 1071 |
046785,002268: 00,3073 13112 TCF LARGE2
046786,002269:
![]() |
Page 1072 |
046788,002271: 00,3074 40000 LARGE3 COM # LEAVE NEGATIVE OF SHIFT COUNT-1 FOR
046789,002272: 00,3075 54135 TS MPTEMP # PREDIVIDE LEFT SHIFT.
046790,002273:
046791,002274: 00,3076 40000 COM # PICK UP REQUIRED SHIFTING BIT TO UNNORM-
046792,002275: 00,3077 50000 INDEX A # ALIZE THE SQRT RESULT.
046793,002276: 00,3100 34736 CAF BIT14
046794,002277: 00,3101 54130 TS BUF
046795,002278: 00,3102 00006 EXTEND
046796,002279: 00,3103 70155 MP MPAC +1
046797,002280: 00,3104 56130 XCH BUF
046798,002281: 00,3105 00006 EXTEND # (UNNORMALIZE THE SQRT FOR LV).
046799,002282: 00,3106 70154 MP MPAC
046800,002283: 00,3107 56001 XCH L
046801,002284: 00,3110 60130 AD BUF
046802,002285: 00,3111 56001 XCH L
046803,002286:
046804,002287: 00,3112 50120 LARGE2 INDEX FIXLOC
046805,002288: 00,3113 52045 DXCH LV # LENGTH NOW STORED IN WORK AREA.
046806,002289:
046807,002290: 00,3114 44753 CS ONE
046808,002291: 00,3115 54140 TS MAXDVSW # NO MAXDV CASES IN UNIT.
046809,002292:
046810,002293: 00,3116 52123 DXCH VBUF # PREPARE X COMPONENT FOR DIVIDE, SETTING
046811,002294: 00,3117 52155 DXCH MPAC # LENGTH OF VECTOR AS DIVISOR IN BUF.
046812,002295: 00,3120 52131 DXCH BUF
046813,002296: 00,3121 03151 TC UNITDV
046814,002297:
046815,002298: 00,3122 52125 DXCH VBUF +2 # DO Y AND Z IN USUAL FASHION SO WE CAN
046816,002299: 00,3123 52155 DXCH MPAC # EXIT THROUGH VROTATEX.
046817,002300: 00,3124 52160 DXCH MPAC +3
046818,002301: 00,3125 03151 TC UNITDV
046819,002302:
046820,002303: 00,3126 52127 DXCH VBUF +4
046821,002304: 00,3127 52155 DXCH MPAC
046822,002305: 00,3130 52162 DXCH MPAC +5
046823,002306: 00,3131 03151 TC UNITDV
046824,002307: 00,3132 17420 TCF VROTATEX # AND EXIT.
046825,002308:
![]() |
Page 1073 |
046827,002310: # IF THE LENGTH OF THE ARGUMENT VECTOR WAS LESS THAN 2(-28), EACH COMPONENT MUST BE SHIFTED LEFT AT LEAST
046828,002311: # 14 PLACES BEFORE THE DIVIDE. NOTE THAT IN THIS CASE, THE MAJOR PART OF EACH COMPONENT IS ZERO.
046829,002312:
046830,002313: 00,3133 54135 SMALL TS MPTEMP # NEGATIVE OF PRE-DIVIDE SHIFT COUNT.
046831,002314:
046832,002315: 00,3134 34755 CAF ZERO # SHIFT EACH COMPONENT LEFT 14.
046833,002316: 00,3135 56123 XCH VBUF +1
046834,002317: 00,3136 56122 XCH VBUF
046835,002318: 00,3137 56125 XCH VBUF +3
046836,002319: 00,3140 56124 XCH VBUF +2
046837,002320: 00,3141 56127 XCH VBUF +5
046838,002321: 00,3142 56126 XCH VBUF +4
046839,002322:
046840,002323: 00,3143 40135 CS MPTEMP
046841,002324: 00,3144 50000 INDEX A
046842,002325: 00,3145 34736 CAF BIT14
046843,002326: 00,3146 00006 EXTEND
046844,002327: 00,3147 70154 MP MPAC
046845,002328: 00,3150 13062 TCF SMALL2
046846,002329:
046847,002330: 00,3151 4761 THIRTEEN = OCT15
046848,002331: 00,3151 4317 FOURTEEN = OCT16
046849,002332: 00,3151 4317 OCT16 = R1D1
046850,002333:
![]() |
Page 1074 |
046852,002335: # THE FOLLOWING ROUTINE SETS UP THE CALL TO THE DIVIDE ROUTINES.
046853,002336:
046854,002337: 00,3151 10154 UNITDV CCS MPAC # FORCE MPAC POSITIVE IF POSSIBLE, SETTING
046855,002338: 00,3152 13170 TCF UMPAC+ # DVSIGN ACCORDING TO THE SIGN OF MPAC
046856,002339: 00,3153 13155 TCF +2 # SINCE THE DIVISOR IS ALWAYS POSITIVE
046857,002340: 00,3154 13162 TCF UMPAC- # HERE.
046858,002341:
046859,002342: 00,3155 10155 CCS MPAC +1
046860,002343: 00,3156 13170 TCF UMPAC+
046861,002344: 00,3157 00002 TC Q # EXIT IMMEDIATELY ON ZERO.
046862,002345: 00,3160 13162 TCF UMPAC-
046863,002346: 00,3161 00002 TC Q
046864,002347:
046865,002348: 00,3162 44755 UMPAC- CS ZERO # IF NEGATIVE, SET -0 IN DVSIGN FOR FINAL
046866,002349: 00,3163 54136 TS DVSIGN # COMPLEMENT.
046867,002350: 00,3164 00006 EXTEND
046868,002351: 00,3165 40155 DCS MPAC # PICK UP ABSOLUTE VALUE OF ARG AND JUMP.
046869,002352: 00,3166 50135 INDEX MPTEMP
046870,002353: 00,3167 12564 TCF MAXTEST -1
046871,002354:
046872,002355: 00,3170 54136 UMPAC+ TS DVSIGN # SET DVSIGN FOR POSITIVE QUOTIENT.
046873,002356: 00,3171 52155 DXCH MPAC
046874,002357: 00,3172 50135 INDEX MPTEMP
046875,002358: 00,3173 12564 TCF MAXTEST -1
046876,002359:
![]() |
Page 1075 |
046878,002361: # MISCELLANEOUS UNARY OPERATIONS.
046879,002362:
046880,002363: 00,3174 03300 DSQ TC DSQSUB # SQUARE THE DP CONTENTS OF MPAC.
046881,002364: 00,3175 16061 TCF DANZIG
046882,002365:
046883,002366: 00,3176 10163 ABVALABS CCS MODE # ABVAL OR ABS INSTRUCTION.
046884,002367: 00,3177 13226 TCF ABS # DO ABS ON SCALAR.
046885,002368: 00,3200 13226 TCF ABS
046886,002369:
046887,002370: 00,3201 03317 ABVAL TC VSQSUB # DOT MPAC WITH ITSELF.
046888,002371: 00,3202 22163 LXCH MODE # MODE IS NOW DP (L ZERO AFTER DAS).
046889,002372:
046890,002373: 00,3203 00006 EXTEND # STORE SQUARE OF LENGTH IN WORK AREA.
046891,002374: 00,3204 30155 DCA MPAC
046892,002375: 00,3205 50120 INDEX FIXLOC
046893,002376: 00,3206 52043 DXCH LVSQUARE
046894,002377:
![]() |
Page 1076 |
046896,002379: # PROGRAM DESCRIPTION - SUBROUTINE SQRT
046897,002380:
046898,002381: # FUNCTIONAL DESCRIPTION - DOUBLE PRECISION SQUARE ROOT ROUTINE
046899,002382: # THIS PROGRAM TAKES THE SQUARE ROOT OF THE 27 OR 28 MOST SIGNIFICANT BITS IN THE TRIPLE PRECISION SET OF
046900,002383: # NUMBERS - MPAC, MPAC+1, AND MPAC+2. THE ROOT IS RETURNED DOUBLE PRECISION IN MPAC AND MPAC+1.
046901,002384:
046902,002385: # WARNING - THIS SUBROUTINE USES A TRIPLE PRECISION INPUT. THE PROGRAMMER MUST ASSURE THE CONTENTS OF MPAC+2
046903,002386: # ESPECIALLY IF THE CONTENTS OF MPAC IS SMALL OR ZERO. FOR DETAILS SEE STG MEMO NO.949.
046904,002387:
046905,002388: # CALLING SEQUENCE - IN INTERPRETIVE MODE I.E., FOLLOWING TC INTPRET,SQRT NO ADDRESS IS ALLOWED
046906,002389: # INPUT SCALING THE BINARY POINT IS ASSUMED TO THE RIGHT OF BIT 15. THE ANSWER IS RETURNED WITH THE SAME SCALING
046907,002390:
046908,002391: # SUBROUTINES - GENSCR, MPACSHR, SQRTSUB, ABORT
046909,002392:
046910,002393: # ABORT EXIT MODE - ABORTS ON NEGATIVE INPUT -1.2X10E-4 (77775 OCTAL) OR LESS.
046911,002394: # DISPLAYS ERROR CODE 1302
046912,002395: # TC ABORT
046913,002396: # OCT 1302
046914,002397:
046915,002398: # DEBRIS - LOCATIONS BUF, MPTEMP, ADDRWD ARE USED
046916,002399:
046917,002400: 00,3207 03343 SQRT TC SQRTSUB # TAKE THE SQUARE ROOT OF MPAC.
046918,002401: 00,3210 10135 CCS MPTEMP # RETURNED NORMALIZED SQUARE ROOT. SEE IF
046919,002402: 00,3211 13213 TCF +2 # ANY UN-NORMALIZATION REQUIRED AND EXIT
046920,002403: 00,3212 16061 TCF DANZIG # IF NOT.
046921,002404:
046922,002405: 00,3213 63733 AD NEG12 # A RIGHT SHIFT OF MORE THAN 13 COULD BE
046923,002406: 00,3214 00006 EXTEND # REQUIRED IF INPUT WAS ZERO IN MPAC,+1.
046924,002407: 00,3215 63221 BZMF SQRTSHFT # GOES HERE IN MOST CASES.
046925,002408: 00,3216 22007 ZL # IF A LONG SHIFT IS REQUIRED, GO TO
046926,002409: 00,3217 22116 LXCH ADDRWD # GENERAL RIGHT SHIFT ROUTINES.
046927,002410: 00,3220 12303 TCF GENSCR +4 # ADDRWD WAS ZERO TO PREVENT ROUND.
046928,002411:
046929,002412: 00,3221 50135 SQRTSHFT INDEX MPTEMP # SELECT SHIFTING BIT AND EXIT THROUGH
046930,002413: 00,3222 34735 CAF BIT15 # SHIFT ROUTINES.
046931,002414: 00,3223 54135 TS MPTEMP
046932,002415: 00,3224 34755 CAF ZERO # TO ZERO MPAC +2 IN THE PROCESS.
046933,002416: 00,3225 12036 TCF MPACSHR +3
046934,002417:
046935,002418: 00,3226 06723 ABS TC BRANCH # TEST SIGN OF MPAC AND COMPLEMENT IF
046936,002419: 00,3227 16061 TCF DANZIG
046937,002420: 00,3230 16061 TCF DANZIG
046938,002421: 00,3231 17670 TCF COMP
046939,002422:
![]() |
Page 1077 |
046941,002424: 00,3232 44751 VDEF CS FOUR # VECTOR DEFINE - ESSENTIALLY TREATS
046942,002425: 00,3233 26166 ADS PUSHLOC # SCALAR IN MPAC AS X COMPONENT, PUSHES UP
046943,002426: 00,3234 00006 EXTEND # FOR Y AND THEN AGAIN FOR Z.
046944,002427: 00,3235 50000 INDEX A
046945,002428: 00,3236 30003 DCA 2
046946,002429: 00,3237 52160 DXCH MPAC +3
046947,002430: 00,3240 00006 EXTEND
046948,002431: 00,3241 50166 INDEX PUSHLOC
046949,002432: 00,3242 30001 DCA 0
046950,002433: 00,3243 52162 DXCH MPAC +5
046951,002434: 00,3244 16521 TCF VMODE # MODE IS NON VECTOR.
046952,002435:
046953,002436: 00,3245 03317 VSQ TC VSQSUB # DOT MPAC WITH ITSELF.
046954,002437: 00,3246 17332 TCF DMODE # MODE IS NOW DP.
046955,002438:
046956,002439: 00,3247 00006 PUSH EXTEND # PUSH DOWN MPAC LEAVING IT LOADED.
046957,002440: 00,3250 30155 DCA MPAC
046958,002441: 00,3251 50166 INDEX PUSHLOC # PUSH DOWN FIRST TWO REGISTERS IN EACH
046959,002442: 00,3252 52001 DXCH 0
046960,002443:
046961,002444: 00,3253 50163 INDEX MODE # INCREMENT PUSHDOWN POINTER.
046962,002445: 00,3254 36244 CAF NO.WDS
046963,002446: 00,3255 26166 ADS PUSHLOC
046964,002447:
046965,002448: 00,3256 10163 CCS MODE
046966,002449: 00,3257 13272 TCF TPUSH # PUSH DOWN MPAC +2.
046967,002450: 00,3260 16061 TCF DANZIG # DONE FOR DP.
046968,002451:
046969,002452: 00,3261 00006 EXTEND # ON VECTOR, PUSH DOWN Y AND Z COMPONENTS.
046970,002453: 00,3262 30160 DCA MPAC +3
046971,002454: 00,3263 50166 INDEX PUSHLOC
046972,002455: 00,3264 51775 DXCH 0 -4
046973,002456: 00,3265 00006 EXTEND
046974,002457: 00,3266 30162 DCA MPAC +5
046975,002458: 00,3267 50166 INDEX PUSHLOC
046976,002459: 00,3270 51777 DXCH 0 -2
046977,002460: 00,3271 16061 TCF DANZIG
046978,002461:
046979,002462: 00,3272 30156 TPUSH CA MPAC +2
046980,002463: 00,3273 16554 TCF ENDTPUSH +2
046981,002464:
046982,002465: 00,3274 50120 RVQ INDEX FIXLOC # RVQ - RETURN IVA QPRET.
046983,002466: 00,3275 30052 CA QPRET
046984,002467: 00,3276 54117 TS POLISH
046985,002468: 00,3277 16652 TCF GOTO +4 # (ASSUME QPRET POINTS TO FIXED ONLY.)
046986,002469:
![]() |
Page 1078 |
046988,002471: # THE FOLLOWING SUBROUTINES ARE USED IN SQUARING MPAC, IN BOTH THE SCALAR AND VECTOR SENSE. THEY ARE
046989,002472: # SPECIAL CASES OF DMPSUB AND DOTSUB, PUT IN TO SAVE SOME TIME.
046990,002473:
046991,002474: 00,3300 30155 DSQSUB CA MPAC +1 # SQUARES THE SCALAR CONTENTS OF MPAC.
046992,002475: 00,3301 00006 EXTEND
046993,002476: 00,3302 70000 SQUARE
046994,002477: 00,3303 54156 TS MPAC +2
046995,002478: 00,3304 34755 CAF ZERO # FORM 2(CROSS TERM).
046996,002479: 00,3305 56155 XCH MPAC +1
046997,002480: 00,3306 00006 EXTEND
046998,002481: 00,3307 70154 MP MPAC
046999,002482: 00,3310 20001 DDOUBL # AND MAYBE OVERFLOW.
047000,002483: 00,3311 20156 DAS MPAC +1 # AND SET A TO NET OVERFLOW.
047001,002484: 00,3312 56154 XCH MPAC
047002,002485: 00,3313 00006 EXTEND
047003,002486: 00,3314 70000 SQUARE
047004,002487: 00,3315 20155 DAS MPAC
047005,002488: 00,3316 00002 TC Q
047006,002489:
047007,002490: 00,3317 00006 VSQSUB EXTEND # DOTS THE VECTOR IN MPAC WITH ITSELF.
047008,002491: 00,3320 22137 QXCH DOTRET
047009,002492: 00,3321 03300 TC DSQSUB # SQUARE THE X COMPONENT.
047010,002493: 00,3322 52160 DXCH MPAC +3
047011,002494: 00,3323 52155 DXCH MPAC
047012,002495: 00,3324 52131 DXCH BUF # SO WE CAN END IN DOTSUB.
047013,002496: 00,3325 30156 CA MPAC +2
047014,002497: 00,3326 54132 TS BUF +2
047015,002498:
047016,002499: 00,3327 03300 TC DSQSUB # SQUARE Y COMPONENT.
047017,002500: 00,3330 52156 DXCH MPAC +1
047018,002501: 00,3331 20132 DAS BUF +1
047019,002502: 00,3332 60154 AD MPAC
047020,002503: 00,3333 60130 AD BUF
047021,002504: 00,3334 54130 TS BUF
047022,002505: 00,3335 13337 TCF +2
047023,002506: 00,3336 54121 TS OVFIND # IF OVERFLOW.
047024,002507:
047025,002508: 00,3337 52162 DXCH MPAC +5
047026,002509: 00,3340 52155 DXCH MPAC
047027,002510: 00,3341 03300 TC DSQSUB # SQUARE Z COMPONENT.
047028,002511: 00,3342 17205 TCF ENDDOT # END AS IN DOTSUB.
047029,002512:
![]() |
Page 1079 |
047031,002514: # DOUBLE PRECISION SQUARE ROOT ROUTINE. TAKE THE SQUARE ROOT OF THE TRIPLE PRECISION (MPAC +2 USED ONLY
047032,002515: # IN NORMALIZATION) CONTENTS OF MPAC AND LEAVE THE NORMALIZED RESULT IN MPAC (C(MPAC) GREATER THAN OR EQUAL TO
047033,002516: # .5). THE RIGHT SHIFT COUNT (TO UNNORMALIZE) IS LEFT IN MPTEMP.
047034,002517:
047035,002518: 00,3343 34755 SQRTSUB CAF ZERO # START BY ZEROING RIGHT SHIFT COUNT.
047036,002519: 00,3344 54135 TS MPTEMP
047037,002520:
047038,002521: 00,3345 10154 CCS MPAC # CHECK FOR POSITIVE ARGUMENT, SHIFTING
047039,002522: 00,3346 13405 TCF SMPAC+ # FIRST SIGNIFICANT MPAC REGISTER INTO
047040,002523: 00,3347 13351 TCF +2 # MPAC ITSELF.
047041,002524: 00,3350 13373 TCF SQRTNEG # SEE IF MAG OF ARGUMENT LESS THAN 10(-4).
047042,002525:
047043,002526: 00,3351 56156 XCH MPAC +2 # MPAC IS ZERO - SHIFT LEFT 14.
047044,002527: 00,3352 56155 XCH MPAC +1
047045,002528: 00,3353 54154 TS MPAC
047046,002529: 00,3354 34757 CAF SEVEN # AUGMENT RIGHT SHIFT COUNTER.
047047,002530: 00,3355 54135 TS MPTEMP
047048,002531:
047049,002532: 00,3356 10154 CCS MPAC # SEE IF MPAC NOW PNZ.
047050,002533: 00,3357 13405 TCF SMPAC+
047051,002534: 00,3360 13362 TCF +2
047052,002535: 00,3361 13376 TCF ZEROANS # NEGATIVE BUT LESS THAN 10(-4) IN MAG.
047053,002536:
047054,002537: 00,3362 56155 XCH MPAC +1 # ZERO - SHIFT LEFT 14 AGAIN.
047055,002538: 00,3363 54154 TS MPAC
047056,002539: 00,3364 34757 CAF SEVEN # AUGMENT RIGHT SHIFT COUNTER.
047057,002540: 00,3365 26135 ADS MPTEMP
047058,002541:
047059,002542: 00,3366 10154 CCS MPAC
047060,002543: 00,3367 13405 TCF SMPAC+
047061,002544: 00,3370 00002 TC Q # SQRT(0) = 0.
047062,002545: 00,3371 13376 TCF ZEROANS
047063,002546: 00,3372 13453 TCF FIXROOT # DO NOT LEAVE SQRTSUB WITH -0 IN MPAC.
047064,002547:
047065,002548: 00,3373 10000 SQRTNEG CCS A # ARGUMENT IS NEGATIVE, BUT SEE IF SIGN-
047066,002549: 00,3374 13402 TCF SQRTABRT # CORRECTED ARGUMENT IS LESS THAN 10(-4)
047067,002550:
047068,002551: 00,3375 10155 CCS MPAC +1 # IN MAGNITUDE. IF SO, CALL ANSWER ZERO.
047069,002552: 00,3376 34755 ZEROANS CAF ZERO # FORCE ANSWER TO ZERO HERE.
047070,002553: 00,3377 13453 TCF FIXROOT
047071,002554: 00,3400 13402 TCF SQRTABRT
047072,002555: 00,3401 13453 TCF FIXROOT
047073,002556:
047074,002557: 00,3402 52165 SQRTABRT DXCH LOC
047075,002558: 00,3403 05726 TC POODOO1
047076,002559: 00,3404 01302 OCT 1302
047077,002560:
![]() |
Page 1080 |
047079,002562: 00,3405 62444 SMPAC+ AD -1/2+2 # SEE IF ARGUMENT GREATER THAN OR EQUAL TO
047080,002563: 00,3406 00006 EXTEND # .5.
047081,002564: 00,3407 63456 BZMF SRTEST # IF SO, SEE IF LESS THAN .25.
047082,002565:
047083,002566: 00,3410 52155 DXCH MPAC # WE WILL TAKE THE SQUARE ROOT OF MPAC/2.
047084,002567: 00,3411 22021 LXCH SR # SHIFT RIGHT 1 AND GO TO THE SQRT ROUTINE
047085,002568: 00,3412 00006 EXTEND
047086,002569: 00,3413 74736 MP HALF
047087,002570: 00,3414 52155 DXCH MPAC
047088,002571: 00,3415 56021 XCH SR
047089,002572: 00,3416 26155 ADS MPAC +1 # GUARANTEED NO OVERFLOW.
047090,002573:
047091,002574: 00,3417 32314 ARGHI CAF SLOPEHI # ARGUMENT BETWEEN .25 AND .5. GET A
047092,002575: 00,3420 00006 EXTEND # LINEAR APPROXIMATION FOR THIS RANGE.
047093,002576: 00,3421 70154 MP MPAC
047094,002577: 00,3422 62566 AD BIASHI # X0/2 = (MPAC/2)(SLOPEHI) + BIASHI/2.
047095,002578:
047096,002579: 00,3423 54130 +4 TS BUF # X0/2 (ARGLO ENTERS HERE).
047097,002580: 00,3424 30154 CA MPAC # SINGLE-PRECISION THROUGHOUT.
047098,002581: 00,3425 22007 ZL
047099,002582: 00,3426 00006 EXTEND
047100,002583: 00,3427 10130 DV BUF # (MPAC/2)/(X0/2)
047101,002584: 00,3430 00006 EXTEND
047102,002585: 00,3431 74736 MP HALF
047103,002586: 00,3432 26130 ADS BUF # X1 = X0/2 + .5(MPAC/2)/(X0/2).
047104,002587:
047105,002588: 00,3433 00006 EXTEND
047106,002589: 00,3434 74736 MP HALF # FORM UP X1/2.
047107,002590: 00,3435 52155 DXCH MPAC # SAVE AND BRING OUT ARGUMENT.
047108,002591: 00,3436 00006 EXTEND # TAKE DP QUOTIENT WITH X1.
047109,002592: 00,3437 10130 DV BUF
047110,002593: 00,3440 54131 TS BUF +1 # SAVE MAJOR PART OF QUOTIENT.
047111,002594: 00,3441 34755 CAF ZERO # FORM MINOR PART OF QUOTIENT USING
047112,002595: 00,3442 56001 XCH L # (REMAINDER,0).
047113,002596: 00,3443 00006 EXTEND
047114,002597: 00,3444 10130 DV BUF
047115,002598: 00,3445 54001 TS L # IN PREPARATION FOR DAS.
047116,002599: 00,3446 30131 CA BUF +1
047117,002600: 00,3447 20155 DAS MPAC # X2 = X1/2 + (MPAC/2)X1
047118,002601:
047119,002602: 00,3450 00006 EXTEND # OVERFLOWS IF ARG. NEAR POSMAX.
047120,002603: 00,3451 13455 BZF TCQBNK00
047121,002604: 00,3452 34733 CAF POSMAX
047122,002605: 00,3453 54154 FIXROOT TS MPAC
047123,002606: 00,3454 54155 TS MPAC +1
047124,002607: 00,3455 00002 TCQBNK00 TC Q # RETURN TO CALLER TO UNNORMALIZE, ETC.
047125,002608:
![]() |
Page 1081 |
047127,002610: 00,3456 64737 SRTEST AD QUARTER # ARGUMENT WAS LESS THAN .5, SEE IF LESS
047128,002611: 00,3457 00006 EXTEND # THAN .25.
047129,002612: 00,3460 63502 BZMF SQRTNORM # IF SO, BEGIN NORMALIZATION.
047130,002613:
047131,002614: 00,3461 52155 DXCH MPAC # IF BETWEEN .5 AND .25, SHIFT RIGHT 1 AND
047132,002615: 00,3462 22021 LXCH SR # START AT ARGLO.
047133,002616: 00,3463 00006 EXTEND
047134,002617: 00,3464 74736 MP HALF
047135,002618: 00,3465 52155 DXCH MPAC
047136,002619: 00,3466 56021 XCH SR
047137,002620: 00,3467 26155 ADS MPAC +1 # NO OVERFLOW.
047138,002621:
047139,002622: 00,3470 33007 ARGLO CAF SLOPELO # (NORMALIZED) ARGUMENT BETWEEN .125 AND
047140,002623: 00,3471 00006 EXTEND # .25
047141,002624: 00,3472 70154 MP MPAC
047142,002625: 00,3473 62270 AD BIASLO
047143,002626: 00,3474 13423 TCF ARGHI +4 # BEGIN SQUARE ROOT.
047144,002627:
047145,002628: 00,3475 00006 SQRTNM2 EXTEND # SHIFT LEFT 2 AND INCREMENT RIGHT SHIFT
047146,002629: 00,3476 30156 DCA MPAC +1 # COUNT (FOR TERMINAL UNNORMALIZATION).
047147,002630: 00,3477 20156 DAS MPAC +1
047148,002631: 00,3500 60154 AD MPAC
047149,002632: 00,3501 26154 ADS MPAC # (NO OVERFLOW).
047150,002633:
047151,002634: 00,3502 24135 SQRTNORM INCR MPTEMP # FIRST TIME THROUGH, JUST SHIFT LEFT 1
047152,002635: 00,3503 00006 EXTEND # (PUTS IN EFFECTIVE RIGHT SHIFT SINCE
047153,002636: 00,3504 30156 DCA MPAC +1 # WE WANT MPAC/2).
047154,002637: 00,3505 20156 DAS MPAC +1
047155,002638: 00,3506 60154 AD MPAC
047156,002639: 00,3507 26154 ADS MPAC # (AGAIN NO OVERFLOW).
047157,002640: 00,3510 60000 DOUBLE
047158,002641: 00,3511 54022 TS CYL
047159,002642:
047160,002643: 00,3512 10022 NORMTEST CCS CYL # SEE IF ARGUMENT NOW NORMALIZED AT
047161,002644: 00,3513 10022 CCS CYL # GREATER THAN .125.
047162,002645: 00,3514 13475 TCF SQRTNM2 # NO - SHIFT LEFT 2 MORE AND TRY AGAIN.
047163,002646: 00,3515 13417 TCF ARGHI # YES - NOW BETWEEN .5 AND .25.
047164,002647: 00,3516 13470 TCF ARGLO # ARGUMENT NOW BETWEEN .25 AND .125.
047165,002648:
![]() |
Page 1082 |
047167,002650: # TRIGONOMETRIC FUNCTION PACKAGE.
047168,002651: # THE FOLLOWING TRIGONOMETRIC FUNCTIONS ARE AVAILABLE AS INTERPRETIVE OPERATIONS:
047169,002652: # 1. SIN COMPUTES (1/2)SINE(2 PI MPAC).
047170,002653: # 2. COS COMPUTES (1/2)COSINE(2 PI MPAC).
047171,002654: # 3. ASIN COMPUTES (1/2PI)ARCSINE(2 MPAC).
047172,002655: # 4. ACOS COMPUTES (1/2PI)ARCCOSINE(2 MPAC).
047173,002656:
047174,002657: # SIN-ASIN AND COS-ACOS ARE MUTUALLY INVERSE, IE SIN(ASIN(X)) = X.
047175,002658:
047176,002659: 00,3517 06723 COSINE TC BRANCH # FINDS COSINE USING THE IDENTITY
047177,002660: 00,3520 13523 TCF +3 # COS(X) = SIN(PI/2 - ABS(X)).
047178,002661: 00,3521 13526 TCF PRESINE
047179,002662: 00,3522 13526 TCF PRESINE
047180,002663:
047181,002664: 00,3523 00006 +3 EXTEND
047182,002665: 00,3524 40155 DCS MPAC
047183,002666: 00,3525 52155 DXCH MPAC
047184,002667:
047185,002668: 00,3526 34737 PRESINE CAF QUARTER # PI/2 SCALED.
047186,002669: 00,3527 26154 ADS MPAC
047187,002670:
047188,002671: 00,3530 52155 SINE DXCH MPAC # DOUBLE ARGUMENT.
047189,002672: 00,3531 20001 DDOUBL
047190,002673: 00,3532 54000 OVSK # SEE IF OVERFLOW PRESENT.
047191,002674: 00,3533 13536 TCF +3 # IF NOT, ARGUMENT OK AS IS.
047192,002675:
047193,002676: 00,3534 00006 EXTEND # IF SO, WE LOST (OR GAINED) PI, SO
047194,002677: 00,3535 40001 DCOM # COMPLEMENT MPAC USING THE IDENTITY
047195,002678: # SIN(X-(+)PI) = SIN(-X).
047196,002679: 00,3536 52155 +3 DXCH MPAC
047197,002680: 00,3537 30154 CA MPAC # SEE IF ARGUMENT GREATER THAN .5 IN
047198,002681: 00,3540 60000 DOUBLE # MAGNITUDE. IF SO, REDUCE IT TO LESS THAN
047199,002682: 00,3541 54001 TS L # .5 (+-PI/2 SCALED) AS FOLLOWS:
047200,002683: 00,3542 13553 TCF SN1
047201,002684:
047202,002685: 00,3543 50000 INDEX A # IF POSITIVE, FORM PI - X, IF NEGATIVE
047203,002686: 00,3544 34735 CAF NEG1/2 +1 # USE -PI -X.
047204,002687: 00,3545 60000 DOUBLE
047205,002688: 00,3546 00006 EXTEND
047206,002689: 00,3547 60154 SU MPAC # GUARANTEED NO OVERFLOW.
047207,002690: 00,3550 54154 TS MPAC
047208,002691: 00,3551 40155 CS MPAC +1
047209,002692: 00,3552 54155 TS MPAC +1
047210,002693:
![]() |
Page 1083 |
047212,002695: 00,3553 00006 SN1 EXTEND # SET UP TO EVALUATE HASTINGS POLYNOMIAL
047213,002696: 00,3554 30155 DCA MPAC
047214,002697: 00,3555 52134 DXCH BUF2
047215,002698: 00,3556 03300 TC DSQSUB # SQUARE MPAC.
047216,002699:
047217,002700: 00,3557 07222 TC POLY # EVALUATE FOURTH ORDER POLYNOMIAL.
047218,002701: 00,3560 00003 DEC 3 B-14
047219,002702: 00,3561 14441 37325 2DEC +.3926990796
047220,002703:
047221,002704: 00,3563 53250 60764 2DEC -.6459637111
047222,002705:
047223,002706: 00,3565 12146 21276 2DEC +.318758717
047224,002707:
047225,002708: 00,3567 75466 71471 2DEC -.074780249
047226,002709:
047227,002710: 00,3571 00236 32757 2DEC +.009694988
047228,002711:
047229,002712: 00,3573 32470 CAF LBUF2 # MULTIPLY BY ARGUMENT AND SHIFT LEFT 2.
047230,002713: 00,3574 07106 TC DMPSUB -1
047231,002714:
047232,002715: 00,3575 00006 EXTEND
047233,002716: 00,3576 30156 DCA MPAC +1
047234,002717: 00,3577 20156 DAS MPAC +1
047235,002718: 00,3600 60154 AD MPAC
047236,002719: 00,3601 26154 ADS MPAC # NEITHER SHIFT OVERFLOWS.
047237,002720: 00,3602 00006 EXTEND
047238,002721: 00,3603 30156 DCA MPAC +1
047239,002722: 00,3604 20156 DAS MPAC +1
047240,002723: 00,3605 60154 AD MPAC
047241,002724: 00,3606 26154 ADS MPAC
047242,002725: 00,3607 16061 TCF DANZIG
047243,002726:
![]() |
Page 1084 |
047245,002728: # ARCSIN/ARCCOS ROUTINE.
047246,002729:
047247,002730: 00,3610 33631 ARCSIN CAF LASINEX # COMPUTE ARCSIN BY USING THE IDENTITY
047248,002731: 00,3611 13613 TCF +2 # ARCSIN(X) = PI/2 - ARCCOS(X).
047249,002732:
047250,002733: 00,3612 33713 ARCCOS CAF LDANZIG # (EXITS IMMEDIATELY).
047251,002734: 00,3613 54136 TS ESCAPE
047252,002735: 00,3614 06723 TC BRANCH # TEST SIGN OF INPUT.
047253,002736: 00,3615 13625 TCF ACOSST # START IMMEDIATELY IF POSITIVE.
047254,002737: 00,3616 13731 TCF ACOSZERO # ARCCOS(0) = PI/2 = .25.
047255,002738: 00,3617 00006 EXTEND # IF NEGATIVE, USE THE IDENTITY
047256,002739: 00,3620 40155 DCS MPAC # ARCCOS(X) = PI - ARCCOS(-X), FORCING
047257,002740: 00,3621 52155 DXCH MPAC # ARGUMENT POSITIVE.
047258,002741: 00,3622 33734 CAF TCSUBTR # SET EXIT TO DO ABOVE BEFORE
047259,002742: 00,3623 56136 XCH ESCAPE # ARCSIN/ARCCOS CONSIDERATIONS.
047260,002743: 00,3624 54137 TS ESCAPE2
047261,002744:
047262,002745: 00,3625 44736 ACOSST CS HALF # TEST MAGNITUDE OF INPUT.
047263,002746: 00,3626 60154 AD MPAC
047264,002747: 00,3627 10000 CCS A
047265,002748: 00,3630 13721 TCF ACOSOVF # THIS IS PROBABLY AN OVERFLOW CASE.
047266,002749:
047267,002750: 00,3631 13707 LASINEX TCF ASINEX
047268,002751:
047269,002752: 00,3632 13642 TCF ACOSST2 # NO OVERFLOW - PROCEED.
047270,002753:
047271,002754: 00,3633 10155 CCS MPAC +1 # IF MAJOR PART IS .5, CALL ANSWER 0
047272,002755: 00,3634 34755 CAF ZERO # UNLESS MINOR PART NEGATIVE.
047273,002756: 00,3635 13637 TCF ACOS=0
047274,002757:
047275,002758: 00,3636 13642 TCF ACOSST2
047276,002759:
047277,002760: 00,3637 54155 ACOS=0 TS MPAC +1
047278,002761: 00,3640 54154 TS MPAC
047279,002762: 00,3641 00136 TC ESCAPE
047280,002763:
047281,002764: 00,3642 00006 ACOSST2 EXTEND # NOW THAT ARGUMENT IS IN PROPER RANGE,
047282,002765: 00,3643 40155 DCS MPAC # BEGIN COMPUTATION. USE HASTINGS
047283,002766: 00,3644 64736 AD HALF # APPROXIMATION ARCCOS(X) = SQRT(1-X)P(X)
047284,002767: 00,3645 52155 DXCH MPAC # IN A SCALED VERSION WHERE P(X) IS A
047285,002768: 00,3646 52134 DXCH BUF2 # SEVENTH ORDER POLYNOMIAL.
047286,002769:
047287,002770: 00,3647 03343 TC SQRTSUB # RETURNS WITH NORMALIZED SQUARE ROOT.
047288,002771:
047289,002772: 00,3650 10135 CCS MPTEMP # SEE IF UN-NORMALIZATION REQUIRED.
047290,002773: 00,3651 13714 TCF ACOSSHR # IF SO.
047291,002774:
![]() |
Page 1085 |
047293,002776: 00,3652 52155 ACOS3 DXCH MPAC # SET UP FOR POLYNOMIAL EVALUATION.
047294,002777: 00,3653 52134 DXCH BUF2
047295,002778: 00,3654 52155 DXCH MPAC
047296,002779:
047297,002780: 00,3655 07222 TC POLY
047298,002781: 00,3656 00006 DEC 6 B-14
047299,002782: 00,3657 13240 23630 2DEC +.353553385 # COEFFICIENTS ARE C 2(+I)/PISQRT(2) WHERE
047300,002783:
047301,002784: 00,3661 74721 47775 2DEC* -.0483017006 B+1* # I
047302,002785:
047303,002786: 00,3663 02440 20237 2DEC* +.0200273085 B+2* # WHERE C STANDS FOR ORIGINAL COEFFS.
047304,002787:
047305,002788: 00,3665 75067 70742 2DEC* -.0112931863 B+3*
047306,002789:
047307,002790: 00,3667 03436 26756 2DEC* +.00695311612 B+4*
047308,002791:
047309,002792: 00,3671 74037 57640 2DEC* -.00384617957 B+5*
047310,002793:
047311,002794: 00,3673 03046 07143 2DEC* +.001501297736 B+6*
047312,002795:
047313,002796: 00,3675 76654 42244 2DEC* -.000284160334 B+7*
047314,002797:
047315,002798:
047316,002799: 00,3677 32470 CAF LBUF2 # DO FINAL MULTIPLY AND GO TO ANY
047317,002800: 00,3700 07106 TC DMPSUB -1 # EPILOGUE SEQUENCES.
047318,002801: 00,3701 00136 TC ESCAPE
047319,002802:
047320,002803: 00,3702 00006 SUBTR EXTEND # EPILOGUE FOR NEGATIVE INPUTS TO ARCCOS.
047321,002804: 00,3703 40155 DCS MPAC
047322,002805: 00,3704 64736 AD HALF # FORMS PI - ARCCOS(-X) = ARCCOS(X).
047323,002806: 00,3705 52155 DXCH MPAC
047324,002807: 00,3706 00137 TC ESCAPE2 # GO TO POSSIBLE ARCSIN EPILOGUE.
047325,002808:
047326,002809: 00,3707 00006 ASINEX EXTEND
047327,002810: 00,3710 40155 DCS MPAC # ARCSIN EPILOGUE - GET ARCSIN(X)
047328,002811: 00,3711 64737 AD QUARTER # = PI/2 - ARCCOS(X).
047329,002812: 00,3712 52155 DXCH MPAC
047330,002813: 00,3713 16061 LDANZIG TCF DANZIG
047331,002814:
![]() |
Page 1086 |
047333,002816: 00,3714 50000 ACOSSHR INDEX A # THE SHIFT RIGHT IS LESS THAN 14 SINCE
047334,002817: 00,3715 34736 CAF BIT14 # THE INPUT WAS NON-ZERO DP.
047335,002818: 00,3716 54135 TS MPTEMP
047336,002819: 00,3717 02073 TC VSHRRND # DP SHIFT RIGHT AND ROUND.
047337,002820: 00,3720 13652 TCF ACOS3 # PROCEED.
047338,002821:
047339,002822: 00,3721 00006 ACOSOVF EXTEND # IF MAJOR PART WAS ONLY 1 MORE THAN .5,
047340,002823: 00,3722 13637 BZF ACOS=0 # CALL ANSWER ZERO.
047341,002824:
047342,002825: 00,3723 00006 ACOSABRT EXTEND # IF OVERFLOW, CALL ANSWER ZERO BUT
047343,002826: 00,3724 30165 DCA LOC # SOUND AN ALARM.
047344,002827: 00,3725 05732 TC ALARM1
047345,002828: 00,3726 01301 OCT 1301
047346,002829:
047347,002830: 00,3727 34755 CAF ZERO
047348,002831: 00,3730 13637 TCF ACOS=0
047349,002832:
047350,002833: 00,3731 34737 ACOSZERO CAF QUARTER # ACOS(0) = PI/2.
047351,002834: 00,3732 13640 TCF ACOS=0 +1 # SET MPAC AND EXIT VIA ESCAPE.
047352,002835:
047353,002836: 00,3733 77763 NEG12 DEC -12 B-14
047354,002837: 00,3734 13702 TCSUBTR TCF SUBTR
047355,002838:
![]() |
Page 1087 |
047357,002840: # THE FOLLOWING INSTRUCTIONS ARE AVAILABLE FOR SETTING, MODIFYING, AND BRANCHING ON INDEX REGISTERS:
047358,002841: # 1. AXT ADDRESS TO INDEX TRUE.
047359,002842: # 1. AXC ADDRESS TO INDEX COMPLEMENTED.
047360,002843: # 3. LXA LOAD INDEX FROM ERASABLE.
047361,002844: # 4. LXC LOAD INDEX COMPLEMENTED FROM ERASABLE.
047362,002845: # 5. SXA STORE INDEX IN ERASABLE.
047363,002846: # 6. XCHX EXCHANGE INDEX REIGSTER WITH ERASABLE.
047364,002847:
047365,002848: # 7. INCR INCREMENT INDEX REGISTER.
047366,002849: # 8. XAD ERASABLE ADD TO INDEX REGISTER.
047367,002850: # 9. XSU ERASABLE SUBTRACT FROM INDEX REGISTER.
047368,002851: # 10. TIX BRANCH ON INDEX REGISTER AND DECREMENT.
047369,002852:
047370,002853: 01,2344 BANK 01
047371,002854:
047372,002855: 01,2344 COUNT* $$/INTER
047373,002856: 01,2344 02441 AXT TC TAGSUB # SELECT APPROPRIATE INDEX REGISTER.
047374,002857: 01,2345 30117 CA POLISH
047375,002858: 01,2346 50130 XSTORE INDEX INDEXLOC # CONTAINS C(FIXLOC) OR C(FIXLOC)+1.
047376,002859: 01,2347 54046 TS X1
047377,002860: 01,2350 16061 TCF DANZIG
047378,002861:
047379,002862: 01,2351 02441 AXC TC TAGSUB
047380,002863: 01,2352 40117 CS POLISH
047381,002864: 01,2353 02346 TC XSTORE
047382,002865:
047383,002866: 01,2354 02427 LXA TC 15ADRERS # LOAD INDEX REGISTER FROM ERASABLE.
047384,002867: 01,2355 50117 INDEX POLISH
047385,002868: 01,2356 30000 CA 0
047386,002869: 01,2357 12346 TCF XSTORE
047387,002870:
047388,002871: 01,2360 02427 LXC TC 15ADRERS # LOAD NDX REG FROM ERASABLE COMPLEMENTED.
047389,002872: 01,2361 50117 INDEX POLISH
047390,002873: 01,2362 40000 CS 0
047391,002874: 01,2363 12346 TCF XSTORE
047392,002875:
047393,002876: 01,2364 02427 SXA TC 15ADRERS # STORE INDEX REGISTER IN ERASABLE.
047394,002877: 01,2365 50130 INDEX INDEXLOC
047395,002878: 01,2366 30046 CA X1
047396,002879: 01,2367 50117 MSTORE1 INDEX POLISH
047397,002880: 01,2370 54000 TS 0
047398,002881: 01,2371 16061 TCF DANZIG
047399,002882:
![]() |
Page 1088 |
047401,002884: 01,2372 02427 XCHX TC 15ADRERS # EXCHANGE INDEX REGISTER WITH ERASABLE.
047402,002885: 01,2373 50117 INDEX POLISH
047403,002886: 01,2374 30000 CA 0
047404,002887: 01,2375 50130 INDEX INDEXLOC
047405,002888: 01,2376 56046 XCH X1
047406,002889: 01,2377 12367 TCF MSTORE1
047407,002890:
047408,002891: 01,2400 02427 XAD TC 15ADRERS # ADD ERASABLE TO INDEX REGISTER.
047409,002892: 01,2401 50117 INDEX POLISH
047410,002893: 01,2402 30000 CA 0
047411,002894: 01,2403 50130 XAD2 INDEX INDEXLOC
047412,002895: 01,2404 26046 ADS X1 # IGNORING OVERFLOWS.
047413,002896: 01,2405 16061 TCF DANZIG
047414,002897:
047415,002898: 01,2406 02441 INCR TC TAGSUB # INCREMENT INDEX REGISTER.
047416,002899: 01,2407 30117 CA POLISH
047417,002900: 01,2410 12403 TCF XAD2
047418,002901:
047419,002902: 01,2411 02427 XSU TC 15ADRERS # SUBTRACT ERASABLE FROM INDEX REGISTER.
047420,002903: 01,2412 50117 INDEX POLISH
047421,002904: 01,2413 40000 CS 0
047422,002905: 01,2414 12403 TCF XAD2
047423,002906:
047424,002907: 01,2415 02441 TIX TC TAGSUB # BRANCH AND DECREMENT ON INDEX.
047425,002908: 01,2416 50130 INDEX INDEXLOC
047426,002909: 01,2417 40050 CS S1
047427,002910: 01,2420 50130 INDEX INDEXLOC
047428,002911: 01,2421 60046 AD X1
047429,002912: 01,2422 00006 EXTEND # NO OPERATION IF DECREMENTED INDEX IS
047430,002913: 01,2423 66061 BZMF DANZIG # NEGATIVE OR ZERO.
047431,002914:
047432,002915: 01,2424 50130 DOTIXBR INDEX INDEXLOC
047433,002916: 01,2425 56046 XCH X1 # IGNORING OVERFLOWS.
047434,002917:
047435,002918: 01,2426 16646 TCF GOTO # DO THE BRANCH USING THE CADR IN POLISH.
047436,002919:
![]() |
Page 1089 |
047438,002921: # SUBROUTINE TO CONVERT AN ERASABLE ADDRESS (11 BITS) TO AN EBANK SETTING AND SUBADDRESS.
047439,002922:
047440,002923: 01,2427 40117 15ADRERS CS POLISH
047441,002924: 01,2430 64772 AD DEC45
047442,002925: 01,2431 10000 CCS A # DOES THE ADDRESS POINT TO THE WORK AREA?
047443,002926: 01,2432 30120 CA FIXLOC # YES. ADD FIXLOC. EBANK OK AS IS.
047444,002927: 01,2433 12440 TCF +5
047445,002928:
047446,002929: 01,2434 35007 CA OCT1400 # NO. SET EBANK & MAKE UP SUBADDRESS.
047447,002930: 01,2435 56117 XCH POLISH
047448,002931: 01,2436 54003 TS EBANK
047449,002932: 01,2437 74357 MASK LOW8
047450,002933: 01,2440 26117 +5 ADS POLISH # FALL INTO TAGSUB, AND RETURN VIA Q.
047451,002934:
047452,002935: # SUBROUTINE WHICH SETS THE ADDRESS OF THE SPECIFIED INDEX IN INDEXLOC. (ACTUALLY, THE ADDRESS -38D.)
047453,002936:
047454,002937: 01,2441 30120 TAGSUB CA FIXLOC
047455,002938: 01,2442 54130 TS INDEXLOC
047456,002939:
047457,002940: 01,2443 10020 CCS CYR # BIT 15 SPECIFIES INDEX.
047458,002941: 01,2444 24130 INCR INDEXLOC # 0 MEANS USE X2.
047459,002942: 01,2445 00002 TC Q
047460,002943: 01,2446 00002 TC Q # 1 FOR X1.
047461,002944:
![]() |
Page 1090 |
047463,002946: # MISCELLANEOUS OPERATION CODES WITH DIRECT ADDRESSES. INCLUDED HERE ARE:
047464,002947: # 1. ITA STORE QPRET (RETURN ADDRESS) IN ERASABLE.
047465,002948: # 2. CALL CALL A SUBROUTINE, LEAVING RETURN IN QPRET.
047466,002949: # 3. RTB RETURN TO BASIC LANGUAGE AT THE GIVEN ADDRESS.
047467,002950: # 4. BHIZ BRANCH IF THE HIGH ORDER OF MPAC IS ZERO (SINGLE PRECISION).
047468,002951: # 5. BOV BRANCH ON OVERFLOW.
047469,002952: # 6. GOTO SIMPLE SEQUENCE CHANGE.
047470,002953:
047471,002954: 01,2447 10020 RTB/BHIZ CCS CYR
047472,002955: 01,2450 30117 RTB CA POLISH
047473,002956: 01,2451 04621 TC SWCALL -1 # SO A "TC Q" FROM ROUTINE LEADS TO DANZIG
047474,002957:
047475,002958: 01,2452 10154 BHIZ CCS MPAC
047476,002959: 01,2453 16061 TCF DANZIG
047477,002960: 01,2454 16646 TCF GOTO
047478,002961: 01,2455 16061 TCF DANZIG
047479,002962: 01,2456 16646 TCF GOTO
047480,002963:
047481,002964: 01,2457 10121 BOV(B) CCS OVFIND # BRANCH ON OVERFLOW TO BASIC OR INTERP.
047482,002965: 01,2460 12462 TCF +2
047483,002966: 01,2461 16061 TCF DANZIG
047484,002967: 01,2462 54121 TS OVFIND
047485,002968: 01,2463 10020 CCS CYR
047486,002969: 01,2464 12450 TCF RTB # IF BASIC.
047487,002970: 01,2465 00360 B5TOBB OCT 360
047488,002971: 01,2466 16646 TCF GOTO
047489,002972:
![]() |
Page 1091 |
047491,002974: 01,2467 10020 BZE/GOTO CCS CYR # SEE WHICH OP-CODE IS DESIRED.
047492,002975: 01,2470 06723 TC BRANCH # DO BZE.
047493,002976: 01,2471 16061 TCF DANZIG
047494,002977: 01,2472 16646 TCF GOTO # DO GOTO.
047495,002978: 01,2473 16061 TCF DANZIG
047496,002979:
047497,002980: 01,2474 10020 BPL/BMN CCS CYR
047498,002981: 01,2475 12503 TCF BPL
047499,002982: 01,2476 12000 5B10 DEC 5 B+10 # SHIFTS OP CODE IN SWITCH INSTRUCTION ADR
047500,002983:
047501,002984: 01,2477 06723 TC BRANCH # DO BMN.
047502,002985: 01,2500 16061 TCF DANZIG
047503,002986: 01,2501 16061 TCF DANZIG
047504,002987: 01,2502 16646 TCF GOTO # ONLY IF NNZ.
047505,002988:
047506,002989: 01,2503 06723 BPL TC BRANCH
047507,002990: 01,2504 16646 TCF GOTO # IF POSITIVE OR ZERO.
047508,002991: 01,2505 16646 TCF GOTO
047509,002992: 01,2506 16061 TCF DANZIG
047510,002993:
047511,002994: 01,2507 10020 CALL/ITA CCS CYR
047512,002995: 01,2510 16640 TCF CALL
047513,002996:
047514,002997: 01,2511 05705 TC CCSHOLE
047515,002998: 01,2512 02427 TC 15ADRERS # STORE QPRET. (TAGSUB AFTER 15ADRERS IS
047516,002999: 01,2513 50120 INDEX FIXLOC # SLOW IN THIS CASE, BUT SAVES STORAGE.)
047517,003000: 01,2514 30052 CA QPRET
047518,003001: 01,2515 12367 TCF MSTORE1
047519,003002:
![]() |
Page 1092 |
047521,003004: # THE FOLLOWING OPERATIONS ARE AVAILABLE FOR ALTERING AND TESTING INTERPRETIVE SWITCHES:
047522,003005:
047523,003006: # 00 BONSET SET A SWITCH AND DO A GOTO IF IT WAS ON.
047524,003007: # 01 SETGO SET A SWITCH AND DO A GOTO.
047525,003008: # 02 BOFSET SET A SWITCH AND DOA GOTO IF IT WAS OFF
047526,003009: # 03 SET SET A SWITCH.
047527,003010:
047528,003011: # 04 BONINV INVERT A SWITCH AND BRANCH IF IT WAS ON.
047529,003012: # 05 INVGO INVERT A SWITCH AND DO A GOTO.
047530,003013: # 06 BOFINV INVERT A SWITCH AND BRANCH IF IT WAS OFF
047531,003014: # 07 INVERT INVERT A SWITCH.
047532,003015:
047533,003016: # 10 BONCLR CLEAR A SWITCH AND BRANCH IF IT WAS ON.
047534,003017: # 11 CLRGO CLEAR A SWITCH AND DO A GOTO.
047535,003018: # 12 BOFCLR CLEAR A SWITCH AND BRANCH IF IT WAS OFF.
047536,003019: # 13 CLEAR CLEAR A SWITCH.
047537,003020:
047538,003021: # 14 BON BRANCH IF A SWITCH WAS ON.
047539,003022: # 16 BOFF BRANCH IF A SWITCH WAS OFF.
047540,003023:
047541,003024: # THE ADDRESS SUPPLIED WITH THE SWITCH INSTRUCTION IS INTERPRETED AS FOLLOWS:
047542,003025:
047543,003026: # BITS 1-4 SWITCH BIT NUMBER (1-15).
047544,003027: # BITS 5-8 SWITCH OPERATION NUMBER.
047545,003028: # BITS 9- SWITCH WORD NUMBER (UP TO 64 SWITCH WORDS).
047546,003029:
047547,003030: # THE ADDRESS ITSELF IS MADE UP BY THE YUL SYSTEM ASSEMBLER. THE BRANCH INSTRUCTIONS REQUIRE TWO
047548,003031: # ADDRESSES, THE SECOND TAKEN AS THE DIRECT (OR INDIRECT IF IN ERASABLE) ADDRESS OF THE BRANCH.
047549,003032:
047550,003033: 01,2516 34762 SWITCHES CAF LOW4 # LEAVE THE SWITCH BIT IN SWBIT.
047551,003034: 01,2517 70117 MASK POLISH
047552,003035: 01,2520 50000 INDEX A
047553,003036: 01,2521 34735 CAF BIT15 # (NUMBER FROM LEFT TO RIGHT.)
047554,003037: 01,2522 54131 TS SWBIT
047555,003038:
047556,003039: 01,2523 34745 CAF BIT7 # LEAVE THE SWITCH NUMBER IN SWWORD.
047557,003040: 01,2524 00006 EXTEND
047558,003041: 01,2525 70117 MP POLISH
047559,003042: 01,2526 54130 TS SWWORD
047560,003043:
047561,003044: 01,2527 00004 INHINT # DURING SWITCH CHANGE SO RUPT CAN USE TOO
047562,003045: 01,2530 50000 INDEX A # LEAVE THE SWITCH WORD ITSELF IN L.
047563,003046: 01,2531 30074 CA STATE
047564,003047: 01,2532 54002 TS Q # Q WILL BE USED AS A CHANNEL.
![]() |
Page 1093 |
047566,003049: 01,2533 34741 CAF BIT11
047567,003050: 01,2534 00006 EXTEND # DISPATCH SWITCH BIT OPERATION AS IN BITS
047568,003051: 01,2535 70117 MP POLISH # 7-8 OF POLISH.
047569,003052: 01,2536 72576 MASK B3TOB4 # GETS 4X2-BIT CODE.
047570,003053: 01,2537 50000 INDEX A
047571,003054: 01,2540 12541 TCF +1
047572,003055:
047573,003056: 01,2541 30131 +1 CA SWBIT # 00 - SET SWITCH IN QUESTION.
047574,003057: 01,2542 00006 EXTEND
047575,003058: 01,2543 04002 ROR QCHAN
047576,003059: 01,2544 12553 TCF SWSTORE
047577,003060:
047578,003061: 01,2545 30131 +5 CA SWBIT # 01 - INVERT SWITCH.
047579,003062: 01,2546 00006 EXTEND
047580,003063: 01,2547 06002 RXOR QCHAN
047581,003064: 01,2550 12553 TCF SWSTORE
047582,003065:
047583,003066: 01,2551 40131 +9D CS SWBIT # 10 - CLEAR.
047584,003067: 01,2552 70002 MASK Q
047585,003068: 01,2553 50130 SWSTORE INDEX SWWORD
047586,003069: 01,2554 54074 TS STATE # NEW SWITCH WORD.
047587,003070:
![]() |
Page 1094 |
047589,003072: 01,2555 00003 +13D RELINT # 11 - NOOP.
047590,003073: 01,2556 34737 CAF BIT13
047591,003074: 01,2557 00006 EXTEND # DISPATCH SEQUCE CHANGING OR BRANCHING
047592,003075: 01,2560 70117 MP POLISH # CODE.
047593,003076: 01,2561 72576 MASK B3TOB4
047594,003077: 01,2562 50000 INDEX A
047595,003078: 01,2563 12564 TCF +1 # ORIGINALLY STORED IN BITS 5-6.
047596,003079:
047597,003080: 01,2564 40002 +1 CS Q # 00 - BRANCH IF ON.
047598,003081: 01,2565 70131 TEST MASK SWBIT
047599,003082: 01,2566 10000 CCS A
047600,003083: 01,2567 12577 TCF SWSKIP
047601,003084:
047602,003085: 01,2570 16715 +5 TCF SWBRANCH # 01 - GO TO.
047603,003086:
047604,003087: 01,2571 12577 TCF SWSKIP # HERE ONLY ON BIT 15.
047605,003088:
047606,003089: 01,2572 05705 TC CCSHOLE
047607,003090: 01,2573 05705 TC CCSHOLE
047608,003091:
047609,003092: 01,2574 30002 +9D CA Q # 10 - BRANCH IF OFF.
047610,003093: 01,2575 12565 TCF TEST
047611,003094:
047612,003095: 01,2576 00014 B3TOB4 OCT 0014
047613,003096: 01,2577 24164 SWSKIP INCR LOC
047614,003097:
047615,003098: 01,2600 01,2516 SW/ EQUALS SWITCHES
047616,003099:
047617,003100: 01,2600 16061 +13D TCF DANZIG # 11 - NOOP.
047618,003101:
End of include-file INTERPRETER.agc. Parent file is MAIN.agc