File RTD

Directory of image this file is from
This file as a plain text file

[RTD
[PROGRAM TO FIND SQUARE ROOT OF X AND ITS RECIPROCAL
[ENTRY RECIPROCAL OF SQUARE ROOT
[JMP 1K
[X
[RESULT
#1K LDA
0
JMP 1X-1
ADD 1R+1
STC p+111
LDA 2
STC p+110
LDA i2
STC p+106+1
LDA i2
STC p+104+2
ADD p+103+1
[ZERO DIVIDE?
AZE i
JMP 1H
[NEGATIVE ARGUMENT?
APO i
JMP p+6
COM
STC p+75+1
ADD p+74+2
COM
STC p+72+2
[CALCULATE RESULT EXPONENT
LDA
p+70
SCR i1
LZE
ADD 1Q+1 [CONSTANT 1
COM
STA 3
STC p+62 [CLEAR ACCUM
[MODIFY FRACTION IF NECESSARY
LZE
ADD 1M+107 [=-1
STC p+57
ADD 3
STC p+42
[APPROXIMATE RESULT
JMP 1M
p+53
p+71
p+54
JMP 1A
p+71
p+51
p+50
SET i4
-4
[A=X**2
[LOOP
 JMP 1M
p+44
p+43
p+50
[A=A*N
JMP 1M
p+46
p+34
p+44
[A=3-A
JMP 1S
p+37
p+41
p+40
[X=A*X
JMP 1M
p+30
p+35
p+26
[HALVE RESULT
LDA
p+24
ADD 1M+7 [=-1
STC p+22
XSK i4
JMP p-25
[DONE WITH LOOP
SET i3
 0
LDA 3
ADD p+14
STA 3
LDA
p+11+1
STA i3
LDA
p+6+2
STA i3
 JMP [RETURN
  0
0
0
  0
0
0
 2
3000
0000
  0
0
0
  0
0
0
[CONSTANTS FOR APPROXIMATION OF
[  1|ROOT X
 0001
5440
4060
 0002
2103
3727
[SQUARE ROOT OF N,
[FOUND BY MULTIPLYING
[1|SQUARE ROOT OF N  BY N
[ENTRY:
[JMP 1T
[X
[RESULT
#1T LDA
0
JMP 1X-1
ADD 1R+1
STC p+32
ADD 3
STA
p+22
STA
p+23
STC p+23
LDA 2
STC p-25
LDA i2
AZE
JMP p+5
[STORE ZERO IF ZERO
STA 3
STA i3
STA i3
JMP p+13
STC p-35+1
LDA i2
STC p-37+2
[X=1|SQUARE ROOT OF N
JMP 1K
p-41
 0 [ADDR RESULT
[X=RESULT*N
JMP 1M
p-44
 0 [ADDR RESULT
 0 [ADDR RESULT
 JMP [RETURN
[END RTD
[MJM 4-1-70
[MJM 4-15-71



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search