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