Home » Archimedes archive » Acorn User » AU 1994-02.adf » !StarInfo_StarInfo » Hill/SQRtest

Hill/SQRtest

This website contains an archive of files for the Acorn Electron, BBC Micro, Acorn Archimedes, Commodore 16 and Commodore 64 computers, which Dominic Ford has rescued from his private collection of floppy disks and cassettes.

Some of these files were originally commercial releases in the 1980s and 1990s, but they are now widely available online. I assume that copyright over them is no longer being asserted. If you own the copyright and would like files to be removed, please contact me.

Tape/disk: Home » Archimedes archive » Acorn User » AU 1994-02.adf » !StarInfo_StarInfo
Filename: Hill/SQRtest
Read OK:
File size: 0DDF bytes
Load address: 0000
Exec address: 0000
File contents
   10acc%=12
   20DIM code% 2^acc%+4096
   30
   40ex=2:ey=3:dx=4:dy=5:tab=6:sx=7:sy=8:ra=9
   50FOR pass%=0 TO 2 STEP 2
   60P%=code%
   70[OPT pass%
   80
   90.components
  100  STMFD r13!,{r14}
  110  MOV   sx,#160*2^acc%
  120  MOV   sy,#128*2^acc%
  130  MOV   ex,r0
  140  MOV   ey,r1
  150  SUB   dx,ex,sx
  160  SUB   dy,ey,sy
  170  BL    sqr
  180  MOV   ra,r0
  190  CMP   ra,#0
  200  BEQ   zero
  210  STR   ra,radius
  220  MOV   r0,dx,ASL #8
  230  MOV   r1,ra,ASR #(acc%-8)
  240  BL    div
  250  STR   r0,cos
  260  MOV   r0,dy,ASL #8
  270  MOV   r1,ra,ASR #(acc%-8)
  280  BL    div
  290  STR   r0,sin
  300  B     done
  310. zero
  320  MOV   r0,#0
  330  STR   r0,sin
  340  STR   r0,cos
  350.done
  360  LDMFD r13!,{r15}
  370
  380.sqr
  390  STMFD r13!,{r2-r8,r14}
  400  ADR   tab,sqr_tab
  410  :
  420  CMP   dx,#0
  430  RSBLT dx,dx,#0
  440  CMP   dy,#0
  450  RSBLT dy,dy,#0
  460  MOVEQ r0,dx
  470  BEQ   end
  480  CMP   dx,#0
  490  MOVEQ r0,dy
  500  BEQ   end
  510  CMP   dx,dy
  520  BGT   acute
  530  BLT   obtuse
  540  LDR   r1,fortyfive
  550  MOV   r0,dx,ASR #acc%-(20-acc%)
  560  MUL   r0,r1,r0
  570  MOV   r0,r0,ASR #20-acc%
  580  B     end
  590.obtuse
  600  MOV   dy,dy,ASR #acc%-(20-acc%)
  610  CMP   dy,#0
  620  MOVEQ r0,dx
  630  BEQ   end
  640  MOV   r0,dx,ASL #20-acc%
  650  MOV   r1,dy
  660  BL    div
  670  MOV   r0,r0,ASR #2
  680  LDR   r1,[tab,r0,ASL #2]
  690  MUL   r0,r1,dy
  700  MOV   r0,r0,ASR #20-acc%
  710  B     end
  720.acute
  730  MOV   dx,dx,ASR #acc%-(20-acc%)
  740  CMP   dx,#0
  750  MOVEQ r0,dy
  760  BEQ   end
  770  MOV   r0,dy,ASL #20-acc%
  780  MOV   r1,dx
  790  BL    div
  800  MOV   r0,r0,ASR #2
  810  LDR   r1,[tab,r0,ASL #2]
  820  MUL   r0,r1,dx
  830  MOV   r0,r0,ASR #20-acc%
  840  B     end
  850.null
  860  MOV   r0,#0
  870.end
  880  LDMFD r13!,{r2-r8,r15}^
  890
  900.fortyfive EQUD SQR(2)*2^acc%
  910.radius    EQUD 0
  920.sin       EQUD 0
  930.cos       EQUD 0
  940.sqr_tab :]:P%+=2^acc%+4:[OPT pass%:
  950
  960.div
  970  STMFD r13!,{r2,r3,r14}
  980  MOV   r3,#0
  990  STR   r3,neg
 1000  CMP   r0,#0
 1010  STRLT r0,neg
 1020  RSBLT r0,r0,#0
 1030  BEQ   abort
 1040  CMP   r1,#0
 1050  BEQ   error
 1060  CMP   r0,r1
 1070  MOVMI r0,#0
 1080  BMI   abort
 1090  MOV   r2,#1
 1100  CMP   r1,#0
 1110.raiseloop
 1120  BMI   raisedone
 1130  CMP   r1,r0
 1140  BHI   nearlydone
 1150  MOVS  r1,r1,LSL#1
 1160  MOV   r2,r2,LSL#1
 1170  B     raiseloop
 1180.nearlydone
 1190  MOV   r1,r1,LSR#1
 1200  MOV   r2,r2,LSR#1
 1210.raisedone
 1220  CMP   r0,r1
 1230  SUBCS r0,r0,r1
 1240  ADDCS r3,r3,r2
 1250  MOV   r1,r1,LSR#1
 1260  MOVS  r2,r2,LSR#1
 1270  BCC   raisedone
 1280  MOV   r0,r3
 1290  LDR   r3,neg
 1300  CMP   r3,#0
 1310  RSBNE r0,r0,#0
 1320.abort
 1330  LDMFD R13!,{r2,r3,r15}
 1340.error
 1350  ADR   r0,divbyzero
 1360  SWI   "OS_GenerateError"
 1370  LDMFD r13!,{r2,r3,r15}
 1380
 1390.divbyzero
 1400  EQUD 18
 1410  EQUS "Division by zero in code."
 1420  EQUB 0
 1430  ALIGN
 1440
 1450.neg EQUD 0
 1460]
 1470NEXT
 1480
 1490fact%=2^acc%
 1500FOR off%=0 TO fact% STEP 4
 1510  y=off%/fact%
 1520  sqr_tab!(off%)=SQR(y^2+1)*fact%
 1530NEXT
 1540
 1550MODE 13
 1560*pointer
 1570MOUSE RECTANGLE -50,-50,1380,1124
 1580MOUSETO 640+100,512+100
 1590GCOL 3,63
 1600REPEAT
 1610MOUSE mx%,my%,mb%
 1620LINE 640,512,mx%,my%
 1630WAIT
 1640LINE 640,512,mx%,my%
 1650dex=(mx%/ 4)-160
 1660dey=(my%/ 4)-128
 1670r=SQR(dex^2+dey^2)
 1680IF r>0 THEN
 1690vx=dex/r
 1700vy=dey/r
 1710ELSE
 1720vx=0
 1730vy=0
 1740ENDIF
 1750PRINTTAB(0,1) "r = "+FNfix((r))+"      "
 1760PRINTTAB(0,2)"vx = "+FNfix(vx)+"   vy = "+FNfix(vy)+"        "
 1770A%=(mx%/4)*2^acc%
 1780B%=(my%/4)*2^acc%
 1790CALL components
 1800test%=INT((!radius)/2^acc%)
 1810IF test% > 500 OR test%<0 VDU 7
 1820PRINTTAB(0,4) "MC r  = "+FNfix(((!radius)/2^acc%))+"        "
 1830PRINTTAB(0,5) "MC vx = "+FNfix((!cos)/2^acc%)+"   MC vy = "+FNfix((!sin)/2^acc%)+"        "
 1840UNTIL0
 1850
 1860DEF FNfix(n)
 1870t$=STR$(n)
 1880REM =LEFT$(t$,INSTR(t$,".")+4)
 1890=t$

acc%=12
� code% 2^acc%+4096

(,ex=2:ey=3:dx=4:dy=5:tab=6:sx=7:sy=8:ra=9
2� pass%=0 � 2 � 2
<P%=code%
F[OPT pass%
P
Z.components
d  STMFD r13!,{r14}
n  MOV   sx,#160*2^acc%
x  MOV   sy,#128*2^acc%
�  MOV   ex,r0
�  MOV   ey,r1
�  SUB   dx,ex,sx
�  SUB   dy,ey,sy
�  BL    sqr
�  MOV   ra,r0
�  CMP   ra,#0
�  BEQ   zero
�  STR   ra,radius
�  MOV   r0,dx,ASL #8
�  MOV   r1,ra,ASR #(acc%-8)
�  BL    div
�  STR   r0,cos
  MOV   r0,dy,ASL #8
  MOV   r1,ra,ASR #(acc%-8)
  BL    div
"  STR   r0,sin
,  B     done
6
. zero
@  MOV   r0,#0
J  STR   r0,sin
T  STR   r0,cos
^	.done
h  LDMFD r13!,{r15}
r
|.sqr
�  STMFD r13!,{r2-r8,r14}
�  ADR   tab,sqr_tab
�  :
�  CMP   dx,#0
�  RSBLT dx,dx,#0
�  CMP   dy,#0
�  RSBLT dy,dy,#0
�  �Q r0,dx
�  BEQ   end
�  CMP   dx,#0
�  �Q r0,dy
�  BEQ   end
�  CMP   dx,dy
  BGT   acute
  BLT   obtuse
  LDR   r1,fortyfive
&%  MOV   r0,dx,ASR #acc%-(20-acc%)
0  MUL   r0,r1,r0
:  MOV   r0,r0,ASR #20-acc%
D  B     end
N.obtuse
X%  MOV   dy,dy,ASR #acc%-(20-acc%)
b  CMP   dy,#0
l  �Q r0,dx
v  BEQ   end
�  MOV   r0,dx,ASL #20-acc%
�  MOV   r1,dy
�  BL    div
�  MOV   r0,r0,ASR #2
�  LDR   r1,[tab,r0,ASL #2]
�  MUL   r0,r1,dy
�  MOV   r0,r0,ASR #20-acc%
�  B     end
�
.acute
�%  MOV   dx,dx,ASR #acc%-(20-acc%)
�  CMP   dx,#0
�  �Q r0,dy
�  BEQ   end
  MOV   r0,dy,ASL #20-acc%
  MOV   r1,dx
  BL    div
   MOV   r0,r0,ASR #2
*  LDR   r1,[tab,r0,ASL #2]
4  MUL   r0,r1,dx
>  MOV   r0,r0,ASR #20-acc%
H  B     end
R	.null
\  MOV   r0,#0
f.end
p  LDMFD r13!,{r2-r8,r15}^
z
�.fortyfive EQUD �(2)*2^acc%
�.radius    EQUD 0
�.sin       EQUD 0
�.cos       EQUD 0
�(.sqr_tab :]:P%+=2^acc%+4:[OPT pass%:
�
�.div
�  STMFD r13!,{r2,r3,r14}
�  MOV   r3,#0
�  STR   r3,neg
�  CMP   r0,#0
�  STRLT r0,neg
�  RSBLT r0,r0,#0
  BEQ   abort
  CMP   r1,#0
  BEQ   error
$  CMP   r0,r1
.  MOVMI r0,#0
8  BMI   abort
B  MOV   r2,#1
L  CMP   r1,#0
V.raiseloop
`  BMI   raisedone
j  CMP   r1,r0
t  BHI   nearlydone
~  MOVS  r1,r1,LSL#1
�  MOV   r2,r2,LSL#1
�  B     raiseloop
�.nearlydone
�  MOV   r1,r1,LSR#1
�  MOV   r2,r2,LSR#1
�.raisedone
�  CMP   r0,r1
�  SUBCS r0,r0,r1
�  ADDCS r3,r3,r2
�  MOV   r1,r1,LSR#1
�  MOVS  r2,r2,LSR#1
�  BCC   raisedone
  MOV   r0,r3

  LDR   r3,neg
  CMP   r3,#0
  RSBNE r0,r0,#0
(
.abort
2  LDMFD R13!,{r2,r3,r15}
<
.error
F  ADR   r0,divbyzero
P  SWI   "OS_GenerateError"
Z  LDMFD r13!,{r2,r3,r15}
d
n.divbyzero
x
  EQUD 18
�&  EQUS "Division by zero in code."
�  EQUB 0
�  ALIGN
�
�.neg EQUD 0
�]
��
�
�fact%=2^acc%
�� off%=0 � fact% � 4
�  y=off%/fact%
�#  sqr_tab!(off%)=�(y^2+1)*fact%
��

� 13
*pointer
"ȗ ȓ -50,-50,1380,1124
,ȗ� 640+100,512+100
6
� 3,63
@�
Jȗ mx%,my%,mb%
T� 640,512,mx%,my%
^Ȗ
h� 640,512,mx%,my%
rdex=(mx%/ 4)-160
|dey=(my%/ 4)-128
�r=�(dex^2+dey^2)
�� r>0 �
�vx=dex/r
�vy=dey/r
��
�vx=0
�vy=0
��
�$�0,1) "r = "+�fix((r))+"      "
�9�0,2)"vx = "+�fix(vx)+"   vy = "+�fix(vy)+"        "
�A%=(mx%/4)*2^acc%
�B%=(my%/4)*2^acc%
�� components
test%=�((!radius)/2^acc%)
� test% > 500 � test%<0 � 7
9�0,4) "MC r  = "+�fix(((!radius)/2^acc%))+"        "
&V�0,5) "MC vx = "+�fix((!cos)/2^acc%)+"   MC vy = "+�fix((!sin)/2^acc%)+"        "
0�0
:
D
� �fix(n)
Nt$=�(n)
X � =LEFT$(t$,INSTR(t$,".")+4)
b=t$
�
00000000  0d 00 0a 0b 61 63 63 25  3d 31 32 0d 00 14 17 de  |....acc%=12.....|
00000010  20 63 6f 64 65 25 20 32  5e 61 63 63 25 2b 34 30  | code% 2^acc%+40|
00000020  39 36 0d 00 1e 04 0d 00  28 2c 65 78 3d 32 3a 65  |96......(,ex=2:e|
00000030  79 3d 33 3a 64 78 3d 34  3a 64 79 3d 35 3a 74 61  |y=3:dx=4:dy=5:ta|
00000040  62 3d 36 3a 73 78 3d 37  3a 73 79 3d 38 3a 72 61  |b=6:sx=7:sy=8:ra|
00000050  3d 39 0d 00 32 15 e3 20  70 61 73 73 25 3d 30 20  |=9..2.. pass%=0 |
00000060  b8 20 32 20 88 20 32 0d  00 3c 0c 50 25 3d 63 6f  |. 2 . 2..<.P%=co|
00000070  64 65 25 0d 00 46 0e 5b  4f 50 54 20 70 61 73 73  |de%..F.[OPT pass|
00000080  25 0d 00 50 04 0d 00 5a  0f 2e 63 6f 6d 70 6f 6e  |%..P...Z..compon|
00000090  65 6e 74 73 0d 00 64 16  20 20 53 54 4d 46 44 20  |ents..d.  STMFD |
000000a0  72 31 33 21 2c 7b 72 31  34 7d 0d 00 6e 1a 20 20  |r13!,{r14}..n.  |
000000b0  4d 4f 56 20 20 20 73 78  2c 23 31 36 30 2a 32 5e  |MOV   sx,#160*2^|
000000c0  61 63 63 25 0d 00 78 1a  20 20 4d 4f 56 20 20 20  |acc%..x.  MOV   |
000000d0  73 79 2c 23 31 32 38 2a  32 5e 61 63 63 25 0d 00  |sy,#128*2^acc%..|
000000e0  82 11 20 20 4d 4f 56 20  20 20 65 78 2c 72 30 0d  |..  MOV   ex,r0.|
000000f0  00 8c 11 20 20 4d 4f 56  20 20 20 65 79 2c 72 31  |...  MOV   ey,r1|
00000100  0d 00 96 14 20 20 53 55  42 20 20 20 64 78 2c 65  |....  SUB   dx,e|
00000110  78 2c 73 78 0d 00 a0 14  20 20 53 55 42 20 20 20  |x,sx....  SUB   |
00000120  64 79 2c 65 79 2c 73 79  0d 00 aa 0f 20 20 42 4c  |dy,ey,sy....  BL|
00000130  20 20 20 20 73 71 72 0d  00 b4 11 20 20 4d 4f 56  |    sqr....  MOV|
00000140  20 20 20 72 61 2c 72 30  0d 00 be 11 20 20 43 4d  |   ra,r0....  CM|
00000150  50 20 20 20 72 61 2c 23  30 0d 00 c8 10 20 20 42  |P   ra,#0....  B|
00000160  45 51 20 20 20 7a 65 72  6f 0d 00 d2 15 20 20 53  |EQ   zero....  S|
00000170  54 52 20 20 20 72 61 2c  72 61 64 69 75 73 0d 00  |TR   ra,radius..|
00000180  dc 18 20 20 4d 4f 56 20  20 20 72 30 2c 64 78 2c  |..  MOV   r0,dx,|
00000190  41 53 4c 20 23 38 0d 00  e6 1f 20 20 4d 4f 56 20  |ASL #8....  MOV |
000001a0  20 20 72 31 2c 72 61 2c  41 53 52 20 23 28 61 63  |  r1,ra,ASR #(ac|
000001b0  63 25 2d 38 29 0d 00 f0  0f 20 20 42 4c 20 20 20  |c%-8)....  BL   |
000001c0  20 64 69 76 0d 00 fa 12  20 20 53 54 52 20 20 20  | div....  STR   |
000001d0  72 30 2c 63 6f 73 0d 01  04 18 20 20 4d 4f 56 20  |r0,cos....  MOV |
000001e0  20 20 72 30 2c 64 79 2c  41 53 4c 20 23 38 0d 01  |  r0,dy,ASL #8..|
000001f0  0e 1f 20 20 4d 4f 56 20  20 20 72 31 2c 72 61 2c  |..  MOV   r1,ra,|
00000200  41 53 52 20 23 28 61 63  63 25 2d 38 29 0d 01 18  |ASR #(acc%-8)...|
00000210  0f 20 20 42 4c 20 20 20  20 64 69 76 0d 01 22 12  |.  BL    div..".|
00000220  20 20 53 54 52 20 20 20  72 30 2c 73 69 6e 0d 01  |  STR   r0,sin..|
00000230  2c 10 20 20 42 20 20 20  20 20 64 6f 6e 65 0d 01  |,.  B     done..|
00000240  36 0a 2e 20 7a 65 72 6f  0d 01 40 11 20 20 4d 4f  |6.. zero..@.  MO|
00000250  56 20 20 20 72 30 2c 23  30 0d 01 4a 12 20 20 53  |V   r0,#0..J.  S|
00000260  54 52 20 20 20 72 30 2c  73 69 6e 0d 01 54 12 20  |TR   r0,sin..T. |
00000270  20 53 54 52 20 20 20 72  30 2c 63 6f 73 0d 01 5e  | STR   r0,cos..^|
00000280  09 2e 64 6f 6e 65 0d 01  68 16 20 20 4c 44 4d 46  |..done..h.  LDMF|
00000290  44 20 72 31 33 21 2c 7b  72 31 35 7d 0d 01 72 04  |D r13!,{r15}..r.|
000002a0  0d 01 7c 08 2e 73 71 72  0d 01 86 1c 20 20 53 54  |..|..sqr....  ST|
000002b0  4d 46 44 20 72 31 33 21  2c 7b 72 32 2d 72 38 2c  |MFD r13!,{r2-r8,|
000002c0  72 31 34 7d 0d 01 90 17  20 20 41 44 52 20 20 20  |r14}....  ADR   |
000002d0  74 61 62 2c 73 71 72 5f  74 61 62 0d 01 9a 07 20  |tab,sqr_tab.... |
000002e0  20 3a 0d 01 a4 11 20 20  43 4d 50 20 20 20 64 78  | :....  CMP   dx|
000002f0  2c 23 30 0d 01 ae 14 20  20 52 53 42 4c 54 20 64  |,#0....  RSBLT d|
00000300  78 2c 64 78 2c 23 30 0d  01 b8 11 20 20 43 4d 50  |x,dx,#0....  CMP|
00000310  20 20 20 64 79 2c 23 30  0d 01 c2 14 20 20 52 53  |   dy,#0....  RS|
00000320  42 4c 54 20 64 79 2c 64  79 2c 23 30 0d 01 cc 0e  |BLT dy,dy,#0....|
00000330  20 20 ec 51 20 72 30 2c  64 78 0d 01 d6 0f 20 20  |  .Q r0,dx....  |
00000340  42 45 51 20 20 20 65 6e  64 0d 01 e0 11 20 20 43  |BEQ   end....  C|
00000350  4d 50 20 20 20 64 78 2c  23 30 0d 01 ea 0e 20 20  |MP   dx,#0....  |
00000360  ec 51 20 72 30 2c 64 79  0d 01 f4 0f 20 20 42 45  |.Q r0,dy....  BE|
00000370  51 20 20 20 65 6e 64 0d  01 fe 11 20 20 43 4d 50  |Q   end....  CMP|
00000380  20 20 20 64 78 2c 64 79  0d 02 08 11 20 20 42 47  |   dx,dy....  BG|
00000390  54 20 20 20 61 63 75 74  65 0d 02 12 12 20 20 42  |T   acute....  B|
000003a0  4c 54 20 20 20 6f 62 74  75 73 65 0d 02 1c 18 20  |LT   obtuse.... |
000003b0  20 4c 44 52 20 20 20 72  31 2c 66 6f 72 74 79 66  | LDR   r1,fortyf|
000003c0  69 76 65 0d 02 26 25 20  20 4d 4f 56 20 20 20 72  |ive..&%  MOV   r|
000003d0  30 2c 64 78 2c 41 53 52  20 23 61 63 63 25 2d 28  |0,dx,ASR #acc%-(|
000003e0  32 30 2d 61 63 63 25 29  0d 02 30 14 20 20 4d 55  |20-acc%)..0.  MU|
000003f0  4c 20 20 20 72 30 2c 72  31 2c 72 30 0d 02 3a 1e  |L   r0,r1,r0..:.|
00000400  20 20 4d 4f 56 20 20 20  72 30 2c 72 30 2c 41 53  |  MOV   r0,r0,AS|
00000410  52 20 23 32 30 2d 61 63  63 25 0d 02 44 0f 20 20  |R #20-acc%..D.  |
00000420  42 20 20 20 20 20 65 6e  64 0d 02 4e 0b 2e 6f 62  |B     end..N..ob|
00000430  74 75 73 65 0d 02 58 25  20 20 4d 4f 56 20 20 20  |tuse..X%  MOV   |
00000440  64 79 2c 64 79 2c 41 53  52 20 23 61 63 63 25 2d  |dy,dy,ASR #acc%-|
00000450  28 32 30 2d 61 63 63 25  29 0d 02 62 11 20 20 43  |(20-acc%)..b.  C|
00000460  4d 50 20 20 20 64 79 2c  23 30 0d 02 6c 0e 20 20  |MP   dy,#0..l.  |
00000470  ec 51 20 72 30 2c 64 78  0d 02 76 0f 20 20 42 45  |.Q r0,dx..v.  BE|
00000480  51 20 20 20 65 6e 64 0d  02 80 1e 20 20 4d 4f 56  |Q   end....  MOV|
00000490  20 20 20 72 30 2c 64 78  2c 41 53 4c 20 23 32 30  |   r0,dx,ASL #20|
000004a0  2d 61 63 63 25 0d 02 8a  11 20 20 4d 4f 56 20 20  |-acc%....  MOV  |
000004b0  20 72 31 2c 64 79 0d 02  94 0f 20 20 42 4c 20 20  | r1,dy....  BL  |
000004c0  20 20 64 69 76 0d 02 9e  18 20 20 4d 4f 56 20 20  |  div....  MOV  |
000004d0  20 72 30 2c 72 30 2c 41  53 52 20 23 32 0d 02 a8  | r0,r0,ASR #2...|
000004e0  1e 20 20 4c 44 52 20 20  20 72 31 2c 5b 74 61 62  |.  LDR   r1,[tab|
000004f0  2c 72 30 2c 41 53 4c 20  23 32 5d 0d 02 b2 14 20  |,r0,ASL #2].... |
00000500  20 4d 55 4c 20 20 20 72  30 2c 72 31 2c 64 79 0d  | MUL   r0,r1,dy.|
00000510  02 bc 1e 20 20 4d 4f 56  20 20 20 72 30 2c 72 30  |...  MOV   r0,r0|
00000520  2c 41 53 52 20 23 32 30  2d 61 63 63 25 0d 02 c6  |,ASR #20-acc%...|
00000530  0f 20 20 42 20 20 20 20  20 65 6e 64 0d 02 d0 0a  |.  B     end....|
00000540  2e 61 63 75 74 65 0d 02  da 25 20 20 4d 4f 56 20  |.acute...%  MOV |
00000550  20 20 64 78 2c 64 78 2c  41 53 52 20 23 61 63 63  |  dx,dx,ASR #acc|
00000560  25 2d 28 32 30 2d 61 63  63 25 29 0d 02 e4 11 20  |%-(20-acc%).... |
00000570  20 43 4d 50 20 20 20 64  78 2c 23 30 0d 02 ee 0e  | CMP   dx,#0....|
00000580  20 20 ec 51 20 72 30 2c  64 79 0d 02 f8 0f 20 20  |  .Q r0,dy....  |
00000590  42 45 51 20 20 20 65 6e  64 0d 03 02 1e 20 20 4d  |BEQ   end....  M|
000005a0  4f 56 20 20 20 72 30 2c  64 79 2c 41 53 4c 20 23  |OV   r0,dy,ASL #|
000005b0  32 30 2d 61 63 63 25 0d  03 0c 11 20 20 4d 4f 56  |20-acc%....  MOV|
000005c0  20 20 20 72 31 2c 64 78  0d 03 16 0f 20 20 42 4c  |   r1,dx....  BL|
000005d0  20 20 20 20 64 69 76 0d  03 20 18 20 20 4d 4f 56  |    div.. .  MOV|
000005e0  20 20 20 72 30 2c 72 30  2c 41 53 52 20 23 32 0d  |   r0,r0,ASR #2.|
000005f0  03 2a 1e 20 20 4c 44 52  20 20 20 72 31 2c 5b 74  |.*.  LDR   r1,[t|
00000600  61 62 2c 72 30 2c 41 53  4c 20 23 32 5d 0d 03 34  |ab,r0,ASL #2]..4|
00000610  14 20 20 4d 55 4c 20 20  20 72 30 2c 72 31 2c 64  |.  MUL   r0,r1,d|
00000620  78 0d 03 3e 1e 20 20 4d  4f 56 20 20 20 72 30 2c  |x..>.  MOV   r0,|
00000630  72 30 2c 41 53 52 20 23  32 30 2d 61 63 63 25 0d  |r0,ASR #20-acc%.|
00000640  03 48 0f 20 20 42 20 20  20 20 20 65 6e 64 0d 03  |.H.  B     end..|
00000650  52 09 2e 6e 75 6c 6c 0d  03 5c 11 20 20 4d 4f 56  |R..null..\.  MOV|
00000660  20 20 20 72 30 2c 23 30  0d 03 66 08 2e 65 6e 64  |   r0,#0..f..end|
00000670  0d 03 70 1d 20 20 4c 44  4d 46 44 20 72 31 33 21  |..p.  LDMFD r13!|
00000680  2c 7b 72 32 2d 72 38 2c  72 31 35 7d 5e 0d 03 7a  |,{r2-r8,r15}^..z|
00000690  04 0d 03 84 1f 2e 66 6f  72 74 79 66 69 76 65 20  |......fortyfive |
000006a0  45 51 55 44 20 b6 28 32  29 2a 32 5e 61 63 63 25  |EQUD .(2)*2^acc%|
000006b0  0d 03 8e 15 2e 72 61 64  69 75 73 20 20 20 20 45  |.....radius    E|
000006c0  51 55 44 20 30 0d 03 98  15 2e 73 69 6e 20 20 20  |QUD 0.....sin   |
000006d0  20 20 20 20 45 51 55 44  20 30 0d 03 a2 15 2e 63  |    EQUD 0.....c|
000006e0  6f 73 20 20 20 20 20 20  20 45 51 55 44 20 30 0d  |os       EQUD 0.|
000006f0  03 ac 28 2e 73 71 72 5f  74 61 62 20 3a 5d 3a 50  |..(.sqr_tab :]:P|
00000700  25 2b 3d 32 5e 61 63 63  25 2b 34 3a 5b 4f 50 54  |%+=2^acc%+4:[OPT|
00000710  20 70 61 73 73 25 3a 0d  03 b6 04 0d 03 c0 08 2e  | pass%:.........|
00000720  64 69 76 0d 03 ca 1c 20  20 53 54 4d 46 44 20 72  |div....  STMFD r|
00000730  31 33 21 2c 7b 72 32 2c  72 33 2c 72 31 34 7d 0d  |13!,{r2,r3,r14}.|
00000740  03 d4 11 20 20 4d 4f 56  20 20 20 72 33 2c 23 30  |...  MOV   r3,#0|
00000750  0d 03 de 12 20 20 53 54  52 20 20 20 72 33 2c 6e  |....  STR   r3,n|
00000760  65 67 0d 03 e8 11 20 20  43 4d 50 20 20 20 72 30  |eg....  CMP   r0|
00000770  2c 23 30 0d 03 f2 12 20  20 53 54 52 4c 54 20 72  |,#0....  STRLT r|
00000780  30 2c 6e 65 67 0d 03 fc  14 20 20 52 53 42 4c 54  |0,neg....  RSBLT|
00000790  20 72 30 2c 72 30 2c 23  30 0d 04 06 11 20 20 42  | r0,r0,#0....  B|
000007a0  45 51 20 20 20 61 62 6f  72 74 0d 04 10 11 20 20  |EQ   abort....  |
000007b0  43 4d 50 20 20 20 72 31  2c 23 30 0d 04 1a 11 20  |CMP   r1,#0.... |
000007c0  20 42 45 51 20 20 20 65  72 72 6f 72 0d 04 24 11  | BEQ   error..$.|
000007d0  20 20 43 4d 50 20 20 20  72 30 2c 72 31 0d 04 2e  |  CMP   r0,r1...|
000007e0  11 20 20 4d 4f 56 4d 49  20 72 30 2c 23 30 0d 04  |.  MOVMI r0,#0..|
000007f0  38 11 20 20 42 4d 49 20  20 20 61 62 6f 72 74 0d  |8.  BMI   abort.|
00000800  04 42 11 20 20 4d 4f 56  20 20 20 72 32 2c 23 31  |.B.  MOV   r2,#1|
00000810  0d 04 4c 11 20 20 43 4d  50 20 20 20 72 31 2c 23  |..L.  CMP   r1,#|
00000820  30 0d 04 56 0e 2e 72 61  69 73 65 6c 6f 6f 70 0d  |0..V..raiseloop.|
00000830  04 60 15 20 20 42 4d 49  20 20 20 72 61 69 73 65  |.`.  BMI   raise|
00000840  64 6f 6e 65 0d 04 6a 11  20 20 43 4d 50 20 20 20  |done..j.  CMP   |
00000850  72 31 2c 72 30 0d 04 74  16 20 20 42 48 49 20 20  |r1,r0..t.  BHI  |
00000860  20 6e 65 61 72 6c 79 64  6f 6e 65 0d 04 7e 17 20  | nearlydone..~. |
00000870  20 4d 4f 56 53 20 20 72  31 2c 72 31 2c 4c 53 4c  | MOVS  r1,r1,LSL|
00000880  23 31 0d 04 88 17 20 20  4d 4f 56 20 20 20 72 32  |#1....  MOV   r2|
00000890  2c 72 32 2c 4c 53 4c 23  31 0d 04 92 15 20 20 42  |,r2,LSL#1....  B|
000008a0  20 20 20 20 20 72 61 69  73 65 6c 6f 6f 70 0d 04  |     raiseloop..|
000008b0  9c 0f 2e 6e 65 61 72 6c  79 64 6f 6e 65 0d 04 a6  |...nearlydone...|
000008c0  17 20 20 4d 4f 56 20 20  20 72 31 2c 72 31 2c 4c  |.  MOV   r1,r1,L|
000008d0  53 52 23 31 0d 04 b0 17  20 20 4d 4f 56 20 20 20  |SR#1....  MOV   |
000008e0  72 32 2c 72 32 2c 4c 53  52 23 31 0d 04 ba 0e 2e  |r2,r2,LSR#1.....|
000008f0  72 61 69 73 65 64 6f 6e  65 0d 04 c4 11 20 20 43  |raisedone....  C|
00000900  4d 50 20 20 20 72 30 2c  72 31 0d 04 ce 14 20 20  |MP   r0,r1....  |
00000910  53 55 42 43 53 20 72 30  2c 72 30 2c 72 31 0d 04  |SUBCS r0,r0,r1..|
00000920  d8 14 20 20 41 44 44 43  53 20 72 33 2c 72 33 2c  |..  ADDCS r3,r3,|
00000930  72 32 0d 04 e2 17 20 20  4d 4f 56 20 20 20 72 31  |r2....  MOV   r1|
00000940  2c 72 31 2c 4c 53 52 23  31 0d 04 ec 17 20 20 4d  |,r1,LSR#1....  M|
00000950  4f 56 53 20 20 72 32 2c  72 32 2c 4c 53 52 23 31  |OVS  r2,r2,LSR#1|
00000960  0d 04 f6 15 20 20 42 43  43 20 20 20 72 61 69 73  |....  BCC   rais|
00000970  65 64 6f 6e 65 0d 05 00  11 20 20 4d 4f 56 20 20  |edone....  MOV  |
00000980  20 72 30 2c 72 33 0d 05  0a 12 20 20 4c 44 52 20  | r0,r3....  LDR |
00000990  20 20 72 33 2c 6e 65 67  0d 05 14 11 20 20 43 4d  |  r3,neg....  CM|
000009a0  50 20 20 20 72 33 2c 23  30 0d 05 1e 14 20 20 52  |P   r3,#0....  R|
000009b0  53 42 4e 45 20 72 30 2c  72 30 2c 23 30 0d 05 28  |SBNE r0,r0,#0..(|
000009c0  0a 2e 61 62 6f 72 74 0d  05 32 1c 20 20 4c 44 4d  |..abort..2.  LDM|
000009d0  46 44 20 52 31 33 21 2c  7b 72 32 2c 72 33 2c 72  |FD R13!,{r2,r3,r|
000009e0  31 35 7d 0d 05 3c 0a 2e  65 72 72 6f 72 0d 05 46  |15}..<..error..F|
000009f0  18 20 20 41 44 52 20 20  20 72 30 2c 64 69 76 62  |.  ADR   r0,divb|
00000a00  79 7a 65 72 6f 0d 05 50  1e 20 20 53 57 49 20 20  |yzero..P.  SWI  |
00000a10  20 22 4f 53 5f 47 65 6e  65 72 61 74 65 45 72 72  | "OS_GenerateErr|
00000a20  6f 72 22 0d 05 5a 1c 20  20 4c 44 4d 46 44 20 72  |or"..Z.  LDMFD r|
00000a30  31 33 21 2c 7b 72 32 2c  72 33 2c 72 31 35 7d 0d  |13!,{r2,r3,r15}.|
00000a40  05 64 04 0d 05 6e 0e 2e  64 69 76 62 79 7a 65 72  |.d...n..divbyzer|
00000a50  6f 0d 05 78 0d 20 20 45  51 55 44 20 31 38 0d 05  |o..x.  EQUD 18..|
00000a60  82 26 20 20 45 51 55 53  20 22 44 69 76 69 73 69  |.&  EQUS "Divisi|
00000a70  6f 6e 20 62 79 20 7a 65  72 6f 20 69 6e 20 63 6f  |on by zero in co|
00000a80  64 65 2e 22 0d 05 8c 0c  20 20 45 51 55 42 20 30  |de."....  EQUB 0|
00000a90  0d 05 96 0b 20 20 41 4c  49 47 4e 0d 05 a0 04 0d  |....  ALIGN.....|
00000aa0  05 aa 0f 2e 6e 65 67 20  45 51 55 44 20 30 0d 05  |....neg EQUD 0..|
00000ab0  b4 05 5d 0d 05 be 05 ed  0d 05 c8 04 0d 05 d2 10  |..].............|
00000ac0  66 61 63 74 25 3d 32 5e  61 63 63 25 0d 05 dc 18  |fact%=2^acc%....|
00000ad0  e3 20 6f 66 66 25 3d 30  20 b8 20 66 61 63 74 25  |. off%=0 . fact%|
00000ae0  20 88 20 34 0d 05 e6 12  20 20 79 3d 6f 66 66 25  | . 4....  y=off%|
00000af0  2f 66 61 63 74 25 0d 05  f0 23 20 20 73 71 72 5f  |/fact%...#  sqr_|
00000b00  74 61 62 21 28 6f 66 66  25 29 3d b6 28 79 5e 32  |tab!(off%)=.(y^2|
00000b10  2b 31 29 2a 66 61 63 74  25 0d 05 fa 05 ed 0d 06  |+1)*fact%.......|
00000b20  04 04 0d 06 0e 08 eb 20  31 33 0d 06 18 0c 2a 70  |....... 13....*p|
00000b30  6f 69 6e 74 65 72 0d 06  22 1b c8 97 20 c8 93 20  |ointer.."... .. |
00000b40  2d 35 30 2c 2d 35 30 2c  31 33 38 30 2c 31 31 32  |-50,-50,1380,112|
00000b50  34 0d 06 2c 17 c8 97 b8  20 36 34 30 2b 31 30 30  |4..,.... 640+100|
00000b60  2c 35 31 32 2b 31 30 30  0d 06 36 0a e6 20 33 2c  |,512+100..6.. 3,|
00000b70  36 33 0d 06 40 05 f5 0d  06 4a 12 c8 97 20 6d 78  |63..@....J... mx|
00000b80  25 2c 6d 79 25 2c 6d 62  25 0d 06 54 15 86 20 36  |%,my%,mb%..T.. 6|
00000b90  34 30 2c 35 31 32 2c 6d  78 25 2c 6d 79 25 0d 06  |40,512,mx%,my%..|
00000ba0  5e 06 c8 96 0d 06 68 15  86 20 36 34 30 2c 35 31  |^.....h.. 640,51|
00000bb0  32 2c 6d 78 25 2c 6d 79  25 0d 06 72 14 64 65 78  |2,mx%,my%..r.dex|
00000bc0  3d 28 6d 78 25 2f 20 34  29 2d 31 36 30 0d 06 7c  |=(mx%/ 4)-160..||
00000bd0  14 64 65 79 3d 28 6d 79  25 2f 20 34 29 2d 31 32  |.dey=(my%/ 4)-12|
00000be0  38 0d 06 86 14 72 3d b6  28 64 65 78 5e 32 2b 64  |8....r=.(dex^2+d|
00000bf0  65 79 5e 32 29 0d 06 90  0b e7 20 72 3e 30 20 8c  |ey^2)..... r>0 .|
00000c00  0d 06 9a 0c 76 78 3d 64  65 78 2f 72 0d 06 a4 0c  |....vx=dex/r....|
00000c10  76 79 3d 64 65 79 2f 72  0d 06 ae 05 cc 0d 06 b8  |vy=dey/r........|
00000c20  08 76 78 3d 30 0d 06 c2  08 76 79 3d 30 0d 06 cc  |.vx=0....vy=0...|
00000c30  05 cd 0d 06 d6 24 f1 8a  30 2c 31 29 20 22 72 20  |.....$..0,1) "r |
00000c40  3d 20 22 2b a4 66 69 78  28 28 72 29 29 2b 22 20  |= "+.fix((r))+" |
00000c50  20 20 20 20 20 22 0d 06  e0 39 f1 8a 30 2c 32 29  |     "...9..0,2)|
00000c60  22 76 78 20 3d 20 22 2b  a4 66 69 78 28 76 78 29  |"vx = "+.fix(vx)|
00000c70  2b 22 20 20 20 76 79 20  3d 20 22 2b a4 66 69 78  |+"   vy = "+.fix|
00000c80  28 76 79 29 2b 22 20 20  20 20 20 20 20 20 22 0d  |(vy)+"        ".|
00000c90  06 ea 15 41 25 3d 28 6d  78 25 2f 34 29 2a 32 5e  |...A%=(mx%/4)*2^|
00000ca0  61 63 63 25 0d 06 f4 15  42 25 3d 28 6d 79 25 2f  |acc%....B%=(my%/|
00000cb0  34 29 2a 32 5e 61 63 63  25 0d 06 fe 10 d6 20 63  |4)*2^acc%..... c|
00000cc0  6f 6d 70 6f 6e 65 6e 74  73 0d 07 08 1d 74 65 73  |omponents....tes|
00000cd0  74 25 3d a8 28 28 21 72  61 64 69 75 73 29 2f 32  |t%=.((!radius)/2|
00000ce0  5e 61 63 63 25 29 0d 07  12 1f e7 20 74 65 73 74  |^acc%)..... test|
00000cf0  25 20 3e 20 35 30 30 20  84 20 74 65 73 74 25 3c  |% > 500 . test%<|
00000d00  30 20 ef 20 37 0d 07 1c  39 f1 8a 30 2c 34 29 20  |0 . 7...9..0,4) |
00000d10  22 4d 43 20 72 20 20 3d  20 22 2b a4 66 69 78 28  |"MC r  = "+.fix(|
00000d20  28 28 21 72 61 64 69 75  73 29 2f 32 5e 61 63 63  |((!radius)/2^acc|
00000d30  25 29 29 2b 22 20 20 20  20 20 20 20 20 22 0d 07  |%))+"        "..|
00000d40  26 56 f1 8a 30 2c 35 29  20 22 4d 43 20 76 78 20  |&V..0,5) "MC vx |
00000d50  3d 20 22 2b a4 66 69 78  28 28 21 63 6f 73 29 2f  |= "+.fix((!cos)/|
00000d60  32 5e 61 63 63 25 29 2b  22 20 20 20 4d 43 20 76  |2^acc%)+"   MC v|
00000d70  79 20 3d 20 22 2b a4 66  69 78 28 28 21 73 69 6e  |y = "+.fix((!sin|
00000d80  29 2f 32 5e 61 63 63 25  29 2b 22 20 20 20 20 20  |)/2^acc%)+"     |
00000d90  20 20 20 22 0d 07 30 06  fd 30 0d 07 3a 04 0d 07  |   "..0..0..:...|
00000da0  44 0d dd 20 a4 66 69 78  28 6e 29 0d 07 4e 0b 74  |D.. .fix(n)..N.t|
00000db0  24 3d c3 28 6e 29 0d 07  58 20 f4 20 3d 4c 45 46  |$=.(n)..X . =LEF|
00000dc0  54 24 28 74 24 2c 49 4e  53 54 52 28 74 24 2c 22  |T$(t$,INSTR(t$,"|
00000dd0  2e 22 29 2b 34 29 0d 07  62 07 3d 74 24 0d ff     |.")+4)..b.=t$..|
00000ddf