Home » Archimedes archive » Acorn User » AU 1998-06 A.adf » Regulars » StarInfo/Brobecker/Julia_iim
StarInfo/Brobecker/Julia_iim
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 1998-06 A.adf » Regulars |
| Filename: | StarInfo/Brobecker/Julia_iim |
| Read OK: | ✔ |
| File size: | 15E6 bytes |
| Load address: | 0000 |
| Exec address: | 0000 |
File contents
10REM>>> Julia Set, inverse iterative method. (iim)
20REM Alain BROBECKER (baah/Arm's Tech) on ??-Jun-1995
30s%=5 :REMNb of backward iterations to skip.
40REM==== ASM CODE =================================================
50DIMcode% 1024*4:FORopt%=0TO2STEP2:P%=code%:[opt opt%
60; This routine prints the julia set using iim.
70; Parameter r0=adress of middle of videoram.
80; r1=adress of the middle of the x conversion table.
90; r2=cx=real part of the Julia constant.
100; r3=cy=imaginary part of the Julia constant.
110; r4=adress of the buffer for non processed points.
120.julia_iim0
130 stmfd r13!,{r0-r12,r14}
140 mov r5,#0 ; r5 is the number of points in 'stack'.
150 mov r6,#0 ; r6=x'=real part of current z.
160 mov r7,#0 ; r7=y'=imaginary...
170 b julia_iim1
180; Julia_iim2 is totally similar to julia_iim1, except that we load
190; the coordinates from the point in the table. That' s why it' s
200; placed just above julia_iim1. We could avoid the 'b julia_iim1'
210; just above if we were having julia_iim2 independant from julia_iim1.
220.julia_iim2
230 ldmdb r4!,{r6-r7} ; Load a z in the table.
240 rsb r6,r6,#0 ; r6=-x'.
250 rsb r7,r7,#0 ; r7=-y'.
260 subS r5,r5,#1 ; One point extracted from table.
270 ldmEQfd r13!,{r0-r12,pc} ; Quit if all points are seen.
280; julia_iim1 takes in account the fact that r6 and r7 already contains
290; coordinates of a point to process, so we don' t load them.
300.julia_iim1
310 sub r8,r6,r2 ; r8=a=x'-cx.
320 movS r6,r8
330 rsbMI r6,r6,#0 ; r6=abs(a).
340 subS r7,r7,r3 ; r7=b=y'-cy.
350 rsbMI r7,r7,#0 ; r7=abs(b).
360 movMI r14,#1 ; If b was negative, then r14=1.
370 mul r9,r6,r6 ; r9=a^2.
380 mla r6,r7,r7,r9 ; r6=a^2+b^2.
390 FNsqrt(6,9,10) ; r9=sqrt(a^2+b^2).
400 sub r6,r9,r8 ; r6=-a+sqrt(a^2+b^2).
410 mov r6,r6,lsl #13 ; r6=(-a+sqrt(a^2+b^2))/2.
420 FNsqrt(6,7,10) ; r7=y=sqrt((-a+sqrt(a^2+b^2))/2).
430 add r9,r9,r8 ; r9=a+sqrt(a^2+b^2).
440 mov r9,r9,lsl #13 ; r9=(a+sqrt(a^2+b^2))/2.
450 FNsqrt(9,6,10) ; r6=x=sqrt((a+sqrt(a^2+b^2))/2).
460 cmp r14,#1 ; b was negative?
470 rsbEQ r7,r7,#0 ; r7=y=b/(2*x).
480; Now check if (x;y) is already drawn on screen, and if it' s not the case
490; draw (x;y) and (-x;-y), put (x;y) in the table and continue calculations
500; with (-x;-y). If (x;y) is already drawn on screen, then don' t print
510; points anew, and stop processing those points, so jump to julia_iim2
520; in order to load a new point from table.
530; At the same time we change r14, so that it does no longer contain #1.
540; (It will contain the adress of videoram+offset)
550 mov r8,r6,asr #7 ; r8=int(x*128).
560 add r8,r1,r8,lsl #4 ; r8 points on x in the conversion table.
570 ldmia r8,{r8-r11} ; Offsets & patterns of (x;y) & (-x;-y).
580 mov r12,r7,asr #8 ; r12=int(y*64).
590 add r12,r12,r12,lsl #2 ; r12=5*int(y*64)
600 add r14,r8,r12,lsl #4 ; r14=offset in screen of (x;y).
610 ldr r8,[r14,r0]! ; r8=long containing pixel & change r14.
620 tst r8,r9 ; Is (x;y) already drawn?
630 bNE julia_iim2 ; Then stop processing those points...
640 orr r8,r8,r9 ; Else draw the two points.
650 cmp r5,#s% ; Don' t draw if iteration<s%.
660 strGE r8,[r14]
670 sub r14,r10,r12,lsl #4 ; r14=offset of (-x;-y).
680 ldr r8,[r14,r0]! ; r8=long containing pixel & change r14.
690 orr r8,r8,r11
700 strGE r8,[r14]
710 stmia r4!,{r6-r7} ; Put (x;y) in the table.
720 add r5,r5,#1 ; Well, tell it.
730 b julia_iim1 ; Continue calcs with (x;y).
740]:NEXTopt%
750REM==== BASIC CODE ===============================================
760MODE0:OFF:ORIGIN 641,513:MOUSEON
770DIMv% 8:!v%=148:v%!4=-1:SYS"OS_ReadVduVariables",v%,v%
780v%=!v%+80*128+40 :REMv=middle of videoram.
790DIME% 1024*2*4 :REMBuffer for 'waiting' points.
800b%=1 :REMScreen bank number.
810
820REMCreate x conversion table.
830DIMB% 21*32*4*4
840a%=B%
850FORi%=-40TO40STEP4
860 a%!0=i%:a%!4=1:a%!8=-i%:a%!12=1:a%+=16
870 FORj%=1TO31
880 a%!0=i%:a%!4=1<<j%:a%!8=-i%-4:a%!12=1<<(32-j%):a%+=16
890NEXT,
900B%+=10*32*4*4 :REMB%=middle of conversion table.
910
920REPEAT
930 WAIT:SYS"OS_Byte",&71,b%:b%=b%EOR3:SYS"OS_Byte",&70,b%:CLS
940 MOUSEx,y,z
950 PRINTTAB(33,0);"- Julia iim -":PRINTTAB(75,0);"baah"
960 A%=v%+((b%-1)*80<<8) :REMA% point on workscreen.
970 C%=x<<5
980 D%=y<<5
990 CALLjulia_iim0
1000UNTILz<>0
1010END
1020
1030REM==== MACROS ====================================================
1040REMThis macro puts sqrt(m0) into m1.
1050DEFFNsqrt(m0,m1,m2)
1060LOCALn%
1070[opt opt%
1080 mov m1,#0 ; This is ripped from Jan/BASS.
1090 mov m2,#1<<30 ; |
1100 cmp m0,m2 ; |
1110 subHS m0,m0,m2 ; |
1120 adc m1,m1,m1 ; |
1130]:FORn%=2TO30STEP2:[opt opt%
1140 add m2,m1,#1<<30 ; |
1150 cmp m0,m2,ror #n% ; |
1160 subHS m0,m0,m2,ror #n% ; |
1170 adc m1,m1,m1 ; End of ripped code.
1180]:NEXTn%:[opt opt%
1190 cmp m0,m1 ; Flags=val-root.
1200 addPL m1,m1,#1 ; Round to nearest integer.
1210]:=""
3�>>> Julia Set, inverse iterative method. (iim)
9� Alain BROBECKER (baah/Arm's Tech) on ??-Jun-1995
@s%=5 :�Nb of backward iterations to skip.
(D�==== ASM CODE =================================================
20�code% 1024*4:�opt%=0�2�2:P%=code%:[opt opt%
<2; This routine prints the julia set using iim.
F0; Parameter r0=adress of middle of videoram.
PB; r1=adress of the middle of the x conversion table.
Z6; r2=cx=real part of the Julia constant.
d;; r3=cy=imaginary part of the Julia constant.
nA; r4=adress of the buffer for non processed points.
x.julia_iim0
� stmfd r13!,{r0-r12,r14}
�J mov r5,#0 ; r5 is the number of points in 'stack'.
�A mov r6,#0 ; r6=x'=real part of current z.
�6 mov r7,#0 ; r7=y'=imaginary...
� b julia_iim1
�F; Julia_iim2 is totally similar to julia_iim1, except that we load
�D; the coordinates from the point in the table. That' s why it' s
�E; placed just above julia_iim1. We could avoid the 'b julia_iim1'
�J; just above if we were having julia_iim2 independant from julia_iim1.
�.julia_iim2
�: ldmdb r4!,{r6-r7} ; Load a z in the table.
�+ rsb r6,r6,#0 ; r6=-x'.
�+ rsb r7,r7,#0 ; r7=-y'.
C subS r5,r5,#1 ; One point extracted from table.
@ ldmEQfd r13!,{r0-r12,pc} ; Quit if all points are seen.
J; julia_iim1 takes in account the fact that r6 and r7 already contains
"@; coordinates of a point to process, so we don' t load them.
,.julia_iim1
6/ sub r8,r6,r2 ; r8=a=x'-cx.
@ movS r6,r8
J. rsbMI r6,r6,#0 ; r6=abs(a).
T/ subS r7,r7,r3 ; r7=b=y'-cy.
^. rsbMI r7,r7,#0 ; r7=abs(b).
hB movMI r14,#1 ; If b was negative, then r14=1.
r+ mul r9,r6,r6 ; r9=a^2.
|/ mla r6,r7,r7,r9 ; r6=a^2+b^2.
�4 �sqrt(6,9,10) ; r9=sqrt(a^2+b^2).
�8 sub r6,r9,r8 ; r6=-a+sqrt(a^2+b^2).
�< mov r6,r6,lsl #13 ; r6=(-a+sqrt(a^2+b^2))/2.
�C �sqrt(6,7,10) ; r7=y=sqrt((-a+sqrt(a^2+b^2))/2).
�7 add r9,r9,r8 ; r9=a+sqrt(a^2+b^2).
�; mov r9,r9,lsl #13 ; r9=(a+sqrt(a^2+b^2))/2.
�B �sqrt(9,6,10) ; r6=x=sqrt((a+sqrt(a^2+b^2))/2).
�3 cmp r14,#1 ; b was negative?
�1 rsbEQ r7,r7,#0 ; r7=y=b/(2*x).
�N; Now check if (x;y) is already drawn on screen, and if it' s not the case
�N; draw (x;y) and (-x;-y), put (x;y) in the table and continue calculations
�J; with (-x;-y). If (x;y) is already drawn on screen, then don' t print
�J; points anew, and stop processing those points, so jump to julia_iim2
.; in order to load a new point from table.
K; At the same time we change r14, so that it does no longer contain #1.
5; (It will contain the adress of videoram+offset)
&2 mov r8,r6,asr #7 ; r8=int(x*128).
0K add r8,r1,r8,lsl #4 ; r8 points on x in the conversion table.
:J ldmia r8,{r8-r11} ; Offsets & patterns of (x;y) & (-x;-y).
D2 mov r12,r7,asr #8 ; r12=int(y*64).
N3 add r12,r12,r12,lsl #2 ; r12=5*int(y*64)
XB add r14,r8,r12,lsl #4 ; r14=offset in screen of (x;y).
bJ ldr r8,[r14,r0]! ; r8=long containing pixel & change r14.
l; tst r8,r9 ; Is (x;y) already drawn?
vH bNE julia_iim2 ; Then stop processing those points...
�= orr r8,r8,r9 ; Else draw the two points.
�@ cmp r5,#s% ; Don' t draw if iteration<s%.
� strGE r8,[r14]
�: sub r14,r10,r12,lsl #4 ; r14=offset of (-x;-y).
�J ldr r8,[r14,r0]! ; r8=long containing pixel & change r14.
� orr r8,r8,r11
� strGE r8,[r14]
�; stmia r4!,{r6-r7} ; Put (x;y) in the table.
�2 add r5,r5,#1 ; Well, tell it.
�> b julia_iim1 ; Continue calcs with (x;y).
� ]:�opt%
�D�==== BASIC CODE ===============================================
��0:�:ȑ 641,513:ȗ�
7�v% 8:!v%=148:v%!4=-1:ș"OS_ReadVduVariables",v%,v%
3v%=!v%+80*128+40 :�v=middle of videoram.
8�E% 1024*2*4 :�Buffer for 'waiting' points.
1b%=1 :�Screen bank number.
*
4 �Create x conversion table.
>�B% 21*32*4*4
H a%=B%
R�i%=-40�40�4
\+ a%!0=i%:a%!4=1:a%!8=-i%:a%!12=1:a%+=16
f
�j%=1�31
p; a%!0=i%:a%!4=1<<j%:a%!8=-i%-4:a%!12=1<<(32-j%):a%+=16
z�,
�<B%+=10*32*4*4 :�B%=middle of conversion table.
�
��
�7 Ȗ:ș"OS_Byte",&71,b%:b%=b%�3:ș"OS_Byte",&70,b%:�
� ȗx,y,z
�+ �33,0);"- Julia iim -":�75,0);"baah"
�6 A%=v%+((b%-1)*80<<8) :�A% point on workscreen.
� C%=x<<5
� D%=y<<5
� �julia_iim0
� �z<>0
��
�
E�==== MACROS ====================================================
&�This macro puts sqrt(m0) into m1.
ݤsqrt(m0,m1,m2)
$�n%
.
[opt opt%
8G mov m1,#0 ; This is ripped from Jan/BASS.
B, mov m2,#1<<30 ; |
L, cmp m0,m2 ; |
V, subHS m0,m0,m2 ; |
`, adc m1,m1,m1 ; |
j]:�n%=2�30�2:[opt opt%
t, add m2,m1,#1<<30 ; |
~, cmp m0,m2,ror #n% ; |
�, subHS m0,m0,m2,ror #n% ; |
�= adc m1,m1,m1 ; End of ripped code.
�]:�n%:[opt opt%
�9 cmp m0,m1 ; Flags=val-root.
�C addPL m1,m1,#1 ; Round to nearest integer.
� ]:=""
� 00000000 0d 00 0a 33 f4 3e 3e 3e 20 4a 75 6c 69 61 20 53 |...3.>>> Julia S|
00000010 65 74 2c 20 69 6e 76 65 72 73 65 20 69 74 65 72 |et, inverse iter|
00000020 61 74 69 76 65 20 6d 65 74 68 6f 64 2e 20 28 69 |ative method. (i|
00000030 69 6d 29 0d 00 14 39 f4 20 20 20 20 41 6c 61 69 |im)...9. Alai|
00000040 6e 20 42 52 4f 42 45 43 4b 45 52 20 28 62 61 61 |n BROBECKER (baa|
00000050 68 2f 41 72 6d 27 73 20 54 65 63 68 29 20 6f 6e |h/Arm's Tech) on|
00000060 20 3f 3f 2d 4a 75 6e 2d 31 39 39 35 0d 00 1e 40 | ??-Jun-1995...@|
00000070 73 25 3d 35 20 20 20 20 20 20 20 20 20 20 20 20 |s%=5 |
00000080 20 20 20 20 20 20 20 20 3a f4 4e 62 20 6f 66 20 | :.Nb of |
00000090 62 61 63 6b 77 61 72 64 20 69 74 65 72 61 74 69 |backward iterati|
000000a0 6f 6e 73 20 74 6f 20 73 6b 69 70 2e 0d 00 28 44 |ons to skip...(D|
000000b0 f4 3d 3d 3d 3d 20 41 53 4d 20 43 4f 44 45 20 3d |.==== ASM CODE =|
000000c0 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d |================|
*
000000f0 0d 00 32 30 de 63 6f 64 65 25 20 31 30 32 34 2a |..20.code% 1024*|
00000100 34 3a e3 6f 70 74 25 3d 30 b8 32 88 32 3a 50 25 |4:.opt%=0.2.2:P%|
00000110 3d 63 6f 64 65 25 3a 5b 6f 70 74 20 6f 70 74 25 |=code%:[opt opt%|
00000120 0d 00 3c 32 3b 20 54 68 69 73 20 72 6f 75 74 69 |..<2; This routi|
00000130 6e 65 20 70 72 69 6e 74 73 20 74 68 65 20 6a 75 |ne prints the ju|
00000140 6c 69 61 20 73 65 74 20 75 73 69 6e 67 20 69 69 |lia set using ii|
00000150 6d 2e 0d 00 46 30 3b 20 50 61 72 61 6d 65 74 65 |m...F0; Paramete|
00000160 72 20 72 30 3d 61 64 72 65 73 73 20 6f 66 20 6d |r r0=adress of m|
00000170 69 64 64 6c 65 20 6f 66 20 76 69 64 65 6f 72 61 |iddle of videora|
00000180 6d 2e 0d 00 50 42 3b 20 20 20 20 20 20 20 20 20 |m...PB; |
00000190 20 20 72 31 3d 61 64 72 65 73 73 20 6f 66 20 74 | r1=adress of t|
000001a0 68 65 20 6d 69 64 64 6c 65 20 6f 66 20 74 68 65 |he middle of the|
000001b0 20 78 20 63 6f 6e 76 65 72 73 69 6f 6e 20 74 61 | x conversion ta|
000001c0 62 6c 65 2e 0d 00 5a 36 3b 20 20 20 20 20 20 20 |ble...Z6; |
000001d0 20 20 20 20 72 32 3d 63 78 3d 72 65 61 6c 20 70 | r2=cx=real p|
000001e0 61 72 74 20 6f 66 20 74 68 65 20 4a 75 6c 69 61 |art of the Julia|
000001f0 20 63 6f 6e 73 74 61 6e 74 2e 0d 00 64 3b 3b 20 | constant...d;; |
00000200 20 20 20 20 20 20 20 20 20 20 72 33 3d 63 79 3d | r3=cy=|
00000210 69 6d 61 67 69 6e 61 72 79 20 70 61 72 74 20 6f |imaginary part o|
00000220 66 20 74 68 65 20 4a 75 6c 69 61 20 63 6f 6e 73 |f the Julia cons|
00000230 74 61 6e 74 2e 0d 00 6e 41 3b 20 20 20 20 20 20 |tant...nA; |
00000240 20 20 20 20 20 72 34 3d 61 64 72 65 73 73 20 6f | r4=adress o|
00000250 66 20 74 68 65 20 62 75 66 66 65 72 20 66 6f 72 |f the buffer for|
00000260 20 6e 6f 6e 20 70 72 6f 63 65 73 73 65 64 20 70 | non processed p|
00000270 6f 69 6e 74 73 2e 0d 00 78 0f 2e 6a 75 6c 69 61 |oints...x..julia|
00000280 5f 69 69 6d 30 0d 00 82 1f 20 20 73 74 6d 66 64 |_iim0.... stmfd|
00000290 20 20 20 72 31 33 21 2c 7b 72 30 2d 72 31 32 2c | r13!,{r0-r12,|
000002a0 72 31 34 7d 0d 00 8c 4a 20 20 6d 6f 76 20 20 20 |r14}...J mov |
000002b0 20 20 72 35 2c 23 30 20 20 20 20 20 20 20 20 20 | r5,#0 |
000002c0 20 20 20 20 20 20 3b 20 72 35 20 69 73 20 74 68 | ; r5 is th|
000002d0 65 20 6e 75 6d 62 65 72 20 6f 66 20 70 6f 69 6e |e number of poin|
000002e0 74 73 20 69 6e 20 27 73 74 61 63 6b 27 2e 0d 00 |ts in 'stack'...|
000002f0 96 41 20 20 6d 6f 76 20 20 20 20 20 72 36 2c 23 |.A mov r6,#|
00000300 30 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 |0 |
00000310 3b 20 72 36 3d 78 27 3d 72 65 61 6c 20 70 61 72 |; r6=x'=real par|
00000320 74 20 6f 66 20 63 75 72 72 65 6e 74 20 7a 2e 0d |t of current z..|
00000330 00 a0 36 20 20 6d 6f 76 20 20 20 20 20 72 37 2c |..6 mov r7,|
00000340 23 30 20 20 20 20 20 20 20 20 20 20 20 20 20 20 |#0 |
00000350 20 3b 20 72 37 3d 79 27 3d 69 6d 61 67 69 6e 61 | ; r7=y'=imagina|
00000360 72 79 2e 2e 2e 0d 00 aa 18 20 20 62 20 20 20 20 |ry....... b |
00000370 20 20 20 6a 75 6c 69 61 5f 69 69 6d 31 0d 00 b4 | julia_iim1...|
00000380 46 3b 20 4a 75 6c 69 61 5f 69 69 6d 32 20 69 73 |F; Julia_iim2 is|
00000390 20 74 6f 74 61 6c 6c 79 20 73 69 6d 69 6c 61 72 | totally similar|
000003a0 20 74 6f 20 6a 75 6c 69 61 5f 69 69 6d 31 2c 20 | to julia_iim1, |
000003b0 65 78 63 65 70 74 20 74 68 61 74 20 77 65 20 6c |except that we l|
000003c0 6f 61 64 0d 00 be 44 3b 20 74 68 65 20 63 6f 6f |oad...D; the coo|
000003d0 72 64 69 6e 61 74 65 73 20 66 72 6f 6d 20 74 68 |rdinates from th|
000003e0 65 20 70 6f 69 6e 74 20 69 6e 20 74 68 65 20 74 |e point in the t|
000003f0 61 62 6c 65 2e 20 54 68 61 74 27 20 73 20 77 68 |able. That' s wh|
00000400 79 20 69 74 27 20 73 0d 00 c8 45 3b 20 70 6c 61 |y it' s...E; pla|
00000410 63 65 64 20 6a 75 73 74 20 61 62 6f 76 65 20 6a |ced just above j|
00000420 75 6c 69 61 5f 69 69 6d 31 2e 20 57 65 20 63 6f |ulia_iim1. We co|
00000430 75 6c 64 20 61 76 6f 69 64 20 74 68 65 20 27 62 |uld avoid the 'b|
00000440 20 6a 75 6c 69 61 5f 69 69 6d 31 27 0d 00 d2 4a | julia_iim1'...J|
00000450 3b 20 6a 75 73 74 20 61 62 6f 76 65 20 69 66 20 |; just above if |
00000460 77 65 20 77 65 72 65 20 68 61 76 69 6e 67 20 6a |we were having j|
00000470 75 6c 69 61 5f 69 69 6d 32 20 69 6e 64 65 70 65 |ulia_iim2 indepe|
00000480 6e 64 61 6e 74 20 66 72 6f 6d 20 6a 75 6c 69 61 |ndant from julia|
00000490 5f 69 69 6d 31 2e 0d 00 dc 0f 2e 6a 75 6c 69 61 |_iim1......julia|
000004a0 5f 69 69 6d 32 0d 00 e6 3a 20 20 6c 64 6d 64 62 |_iim2...: ldmdb|
000004b0 20 20 20 72 34 21 2c 7b 72 36 2d 72 37 7d 20 20 | r4!,{r6-r7} |
000004c0 20 20 20 20 20 20 20 3b 20 4c 6f 61 64 20 61 20 | ; Load a |
000004d0 7a 20 69 6e 20 74 68 65 20 74 61 62 6c 65 2e 0d |z in the table..|
000004e0 00 f0 2b 20 20 72 73 62 20 20 20 20 20 72 36 2c |..+ rsb r6,|
000004f0 72 36 2c 23 30 20 20 20 20 20 20 20 20 20 20 20 |r6,#0 |
00000500 20 3b 20 72 36 3d 2d 78 27 2e 0d 00 fa 2b 20 20 | ; r6=-x'....+ |
00000510 72 73 62 20 20 20 20 20 72 37 2c 72 37 2c 23 30 |rsb r7,r7,#0|
00000520 20 20 20 20 20 20 20 20 20 20 20 20 3b 20 72 37 | ; r7|
00000530 3d 2d 79 27 2e 0d 01 04 43 20 20 73 75 62 53 20 |=-y'....C subS |
00000540 20 20 20 72 35 2c 72 35 2c 23 31 20 20 20 20 20 | r5,r5,#1 |
00000550 20 20 20 20 20 20 20 3b 20 4f 6e 65 20 70 6f 69 | ; One poi|
00000560 6e 74 20 65 78 74 72 61 63 74 65 64 20 66 72 6f |nt extracted fro|
00000570 6d 20 74 61 62 6c 65 2e 0d 01 0e 40 20 20 6c 64 |m table....@ ld|
00000580 6d 45 51 66 64 20 72 31 33 21 2c 7b 72 30 2d 72 |mEQfd r13!,{r0-r|
00000590 31 32 2c 70 63 7d 20 20 20 20 3b 20 51 75 69 74 |12,pc} ; Quit|
000005a0 20 69 66 20 61 6c 6c 20 70 6f 69 6e 74 73 20 61 | if all points a|
000005b0 72 65 20 73 65 65 6e 2e 0d 01 18 4a 3b 20 6a 75 |re seen....J; ju|
000005c0 6c 69 61 5f 69 69 6d 31 20 74 61 6b 65 73 20 69 |lia_iim1 takes i|
000005d0 6e 20 61 63 63 6f 75 6e 74 20 74 68 65 20 66 61 |n account the fa|
000005e0 63 74 20 74 68 61 74 20 72 36 20 61 6e 64 20 72 |ct that r6 and r|
000005f0 37 20 61 6c 72 65 61 64 79 20 63 6f 6e 74 61 69 |7 already contai|
00000600 6e 73 0d 01 22 40 3b 20 63 6f 6f 72 64 69 6e 61 |ns.."@; coordina|
00000610 74 65 73 20 6f 66 20 61 20 70 6f 69 6e 74 20 74 |tes of a point t|
00000620 6f 20 70 72 6f 63 65 73 73 2c 20 73 6f 20 77 65 |o process, so we|
00000630 20 64 6f 6e 27 20 74 20 6c 6f 61 64 20 74 68 65 | don' t load the|
00000640 6d 2e 0d 01 2c 0f 2e 6a 75 6c 69 61 5f 69 69 6d |m...,..julia_iim|
00000650 31 0d 01 36 2f 20 20 73 75 62 20 20 20 20 20 72 |1..6/ sub r|
00000660 38 2c 72 36 2c 72 32 20 20 20 20 20 20 20 20 20 |8,r6,r2 |
00000670 20 20 20 3b 20 72 38 3d 61 3d 78 27 2d 63 78 2e | ; r8=a=x'-cx.|
00000680 0d 01 40 13 20 20 6d 6f 76 53 20 20 20 20 72 36 |..@. movS r6|
00000690 2c 72 38 0d 01 4a 2e 20 20 72 73 62 4d 49 20 20 |,r8..J. rsbMI |
000006a0 20 72 36 2c 72 36 2c 23 30 20 20 20 20 20 20 20 | r6,r6,#0 |
000006b0 20 20 20 20 20 3b 20 72 36 3d 61 62 73 28 61 29 | ; r6=abs(a)|
000006c0 2e 0d 01 54 2f 20 20 73 75 62 53 20 20 20 20 72 |...T/ subS r|
000006d0 37 2c 72 37 2c 72 33 20 20 20 20 20 20 20 20 20 |7,r7,r3 |
000006e0 20 20 20 3b 20 72 37 3d 62 3d 79 27 2d 63 79 2e | ; r7=b=y'-cy.|
000006f0 0d 01 5e 2e 20 20 72 73 62 4d 49 20 20 20 72 37 |..^. rsbMI r7|
00000700 2c 72 37 2c 23 30 20 20 20 20 20 20 20 20 20 20 |,r7,#0 |
00000710 20 20 3b 20 72 37 3d 61 62 73 28 62 29 2e 0d 01 | ; r7=abs(b)...|
00000720 68 42 20 20 6d 6f 76 4d 49 20 20 20 72 31 34 2c |hB movMI r14,|
00000730 23 31 20 20 20 20 20 20 20 20 20 20 20 20 20 20 |#1 |
00000740 3b 20 49 66 20 62 20 77 61 73 20 6e 65 67 61 74 |; If b was negat|
00000750 69 76 65 2c 20 74 68 65 6e 20 72 31 34 3d 31 2e |ive, then r14=1.|
00000760 0d 01 72 2b 20 20 6d 75 6c 20 20 20 20 20 72 39 |..r+ mul r9|
00000770 2c 72 36 2c 72 36 20 20 20 20 20 20 20 20 20 20 |,r6,r6 |
00000780 20 20 3b 20 72 39 3d 61 5e 32 2e 0d 01 7c 2f 20 | ; r9=a^2...|/ |
00000790 20 6d 6c 61 20 20 20 20 20 72 36 2c 72 37 2c 72 | mla r6,r7,r|
000007a0 37 2c 72 39 20 20 20 20 20 20 20 20 20 3b 20 72 |7,r9 ; r|
000007b0 36 3d 61 5e 32 2b 62 5e 32 2e 0d 01 86 34 20 20 |6=a^2+b^2....4 |
000007c0 a4 73 71 72 74 28 36 2c 39 2c 31 30 29 20 20 20 |.sqrt(6,9,10) |
000007d0 20 20 20 20 20 20 20 20 20 20 20 3b 20 72 39 3d | ; r9=|
000007e0 73 71 72 74 28 61 5e 32 2b 62 5e 32 29 2e 0d 01 |sqrt(a^2+b^2)...|
000007f0 90 38 20 20 73 75 62 20 20 20 20 20 72 36 2c 72 |.8 sub r6,r|
00000800 39 2c 72 38 20 20 20 20 20 20 20 20 20 20 20 20 |9,r8 |
00000810 3b 20 72 36 3d 2d 61 2b 73 71 72 74 28 61 5e 32 |; r6=-a+sqrt(a^2|
00000820 2b 62 5e 32 29 2e 0d 01 9a 3c 20 20 6d 6f 76 20 |+b^2)....< mov |
00000830 20 20 20 20 72 36 2c 72 36 2c 6c 73 6c 20 23 31 | r6,r6,lsl #1|
00000840 33 20 20 20 20 20 20 20 3b 20 72 36 3d 28 2d 61 |3 ; r6=(-a|
00000850 2b 73 71 72 74 28 61 5e 32 2b 62 5e 32 29 29 2f |+sqrt(a^2+b^2))/|
00000860 32 2e 0d 01 a4 43 20 20 a4 73 71 72 74 28 36 2c |2....C .sqrt(6,|
00000870 37 2c 31 30 29 20 20 20 20 20 20 20 20 20 20 20 |7,10) |
00000880 20 20 20 3b 20 72 37 3d 79 3d 73 71 72 74 28 28 | ; r7=y=sqrt((|
00000890 2d 61 2b 73 71 72 74 28 61 5e 32 2b 62 5e 32 29 |-a+sqrt(a^2+b^2)|
000008a0 29 2f 32 29 2e 0d 01 ae 37 20 20 61 64 64 20 20 |)/2)....7 add |
000008b0 20 20 20 72 39 2c 72 39 2c 72 38 20 20 20 20 20 | r9,r9,r8 |
000008c0 20 20 20 20 20 20 20 3b 20 72 39 3d 61 2b 73 71 | ; r9=a+sq|
000008d0 72 74 28 61 5e 32 2b 62 5e 32 29 2e 0d 01 b8 3b |rt(a^2+b^2)....;|
000008e0 20 20 6d 6f 76 20 20 20 20 20 72 39 2c 72 39 2c | mov r9,r9,|
000008f0 6c 73 6c 20 23 31 33 20 20 20 20 20 20 20 3b 20 |lsl #13 ; |
00000900 72 39 3d 28 61 2b 73 71 72 74 28 61 5e 32 2b 62 |r9=(a+sqrt(a^2+b|
00000910 5e 32 29 29 2f 32 2e 0d 01 c2 42 20 20 a4 73 71 |^2))/2....B .sq|
00000920 72 74 28 39 2c 36 2c 31 30 29 20 20 20 20 20 20 |rt(9,6,10) |
00000930 20 20 20 20 20 20 20 20 3b 20 72 36 3d 78 3d 73 | ; r6=x=s|
00000940 71 72 74 28 28 61 2b 73 71 72 74 28 61 5e 32 2b |qrt((a+sqrt(a^2+|
00000950 62 5e 32 29 29 2f 32 29 2e 0d 01 cc 33 20 20 63 |b^2))/2)....3 c|
00000960 6d 70 20 20 20 20 20 72 31 34 2c 23 31 20 20 20 |mp r14,#1 |
00000970 20 20 20 20 20 20 20 20 20 20 20 3b 20 62 20 77 | ; b w|
00000980 61 73 20 6e 65 67 61 74 69 76 65 3f 0d 01 d6 31 |as negative?...1|
00000990 20 20 72 73 62 45 51 20 20 20 72 37 2c 72 37 2c | rsbEQ r7,r7,|
000009a0 23 30 20 20 20 20 20 20 20 20 20 20 20 20 3b 20 |#0 ; |
000009b0 72 37 3d 79 3d 62 2f 28 32 2a 78 29 2e 0d 01 e0 |r7=y=b/(2*x)....|
000009c0 4e 3b 20 4e 6f 77 20 63 68 65 63 6b 20 69 66 20 |N; Now check if |
000009d0 28 78 3b 79 29 20 69 73 20 61 6c 72 65 61 64 79 |(x;y) is already|
000009e0 20 64 72 61 77 6e 20 6f 6e 20 73 63 72 65 65 6e | drawn on screen|
000009f0 2c 20 61 6e 64 20 69 66 20 69 74 27 20 73 20 6e |, and if it' s n|
00000a00 6f 74 20 74 68 65 20 63 61 73 65 0d 01 ea 4e 3b |ot the case...N;|
00000a10 20 64 72 61 77 20 28 78 3b 79 29 20 61 6e 64 20 | draw (x;y) and |
00000a20 28 2d 78 3b 2d 79 29 2c 20 70 75 74 20 28 78 3b |(-x;-y), put (x;|
00000a30 79 29 20 69 6e 20 74 68 65 20 74 61 62 6c 65 20 |y) in the table |
00000a40 61 6e 64 20 63 6f 6e 74 69 6e 75 65 20 63 61 6c |and continue cal|
00000a50 63 75 6c 61 74 69 6f 6e 73 0d 01 f4 4a 3b 20 77 |culations...J; w|
00000a60 69 74 68 20 28 2d 78 3b 2d 79 29 2e 20 49 66 20 |ith (-x;-y). If |
00000a70 28 78 3b 79 29 20 69 73 20 61 6c 72 65 61 64 79 |(x;y) is already|
00000a80 20 64 72 61 77 6e 20 6f 6e 20 73 63 72 65 65 6e | drawn on screen|
00000a90 2c 20 74 68 65 6e 20 64 6f 6e 27 20 74 20 70 72 |, then don' t pr|
00000aa0 69 6e 74 0d 01 fe 4a 3b 20 70 6f 69 6e 74 73 20 |int...J; points |
00000ab0 61 6e 65 77 2c 20 61 6e 64 20 73 74 6f 70 20 70 |anew, and stop p|
00000ac0 72 6f 63 65 73 73 69 6e 67 20 74 68 6f 73 65 20 |rocessing those |
00000ad0 70 6f 69 6e 74 73 2c 20 73 6f 20 6a 75 6d 70 20 |points, so jump |
00000ae0 74 6f 20 6a 75 6c 69 61 5f 69 69 6d 32 0d 02 08 |to julia_iim2...|
00000af0 2e 3b 20 69 6e 20 6f 72 64 65 72 20 74 6f 20 6c |.; in order to l|
00000b00 6f 61 64 20 61 20 6e 65 77 20 70 6f 69 6e 74 20 |oad a new point |
00000b10 66 72 6f 6d 20 74 61 62 6c 65 2e 0d 02 12 4b 3b |from table....K;|
00000b20 20 41 74 20 74 68 65 20 73 61 6d 65 20 74 69 6d | At the same tim|
00000b30 65 20 77 65 20 63 68 61 6e 67 65 20 72 31 34 2c |e we change r14,|
00000b40 20 73 6f 20 74 68 61 74 20 69 74 20 64 6f 65 73 | so that it does|
00000b50 20 6e 6f 20 6c 6f 6e 67 65 72 20 63 6f 6e 74 61 | no longer conta|
00000b60 69 6e 20 23 31 2e 0d 02 1c 35 3b 20 28 49 74 20 |in #1....5; (It |
00000b70 77 69 6c 6c 20 63 6f 6e 74 61 69 6e 20 74 68 65 |will contain the|
00000b80 20 61 64 72 65 73 73 20 6f 66 20 76 69 64 65 6f | adress of video|
00000b90 72 61 6d 2b 6f 66 66 73 65 74 29 0d 02 26 32 20 |ram+offset)..&2 |
00000ba0 20 6d 6f 76 20 20 20 20 20 72 38 2c 72 36 2c 61 | mov r8,r6,a|
00000bb0 73 72 20 23 37 20 20 20 20 20 20 20 20 3b 20 72 |sr #7 ; r|
00000bc0 38 3d 69 6e 74 28 78 2a 31 32 38 29 2e 0d 02 30 |8=int(x*128)...0|
00000bd0 4b 20 20 61 64 64 20 20 20 20 20 72 38 2c 72 31 |K add r8,r1|
00000be0 2c 72 38 2c 6c 73 6c 20 23 34 20 20 20 20 20 3b |,r8,lsl #4 ;|
00000bf0 20 72 38 20 70 6f 69 6e 74 73 20 6f 6e 20 78 20 | r8 points on x |
00000c00 69 6e 20 74 68 65 20 63 6f 6e 76 65 72 73 69 6f |in the conversio|
00000c10 6e 20 74 61 62 6c 65 2e 0d 02 3a 4a 20 20 6c 64 |n table...:J ld|
00000c20 6d 69 61 20 20 20 72 38 2c 7b 72 38 2d 72 31 31 |mia r8,{r8-r11|
00000c30 7d 20 20 20 20 20 20 20 20 20 3b 20 4f 66 66 73 |} ; Offs|
00000c40 65 74 73 20 26 20 70 61 74 74 65 72 6e 73 20 6f |ets & patterns o|
00000c50 66 20 28 78 3b 79 29 20 26 20 28 2d 78 3b 2d 79 |f (x;y) & (-x;-y|
00000c60 29 2e 0d 02 44 32 20 20 6d 6f 76 20 20 20 20 20 |)...D2 mov |
00000c70 72 31 32 2c 72 37 2c 61 73 72 20 23 38 20 20 20 |r12,r7,asr #8 |
00000c80 20 20 20 20 3b 20 72 31 32 3d 69 6e 74 28 79 2a | ; r12=int(y*|
00000c90 36 34 29 2e 0d 02 4e 33 20 20 61 64 64 20 20 20 |64)...N3 add |
00000ca0 20 20 72 31 32 2c 72 31 32 2c 72 31 32 2c 6c 73 | r12,r12,r12,ls|
00000cb0 6c 20 23 32 20 20 3b 20 72 31 32 3d 35 2a 69 6e |l #2 ; r12=5*in|
00000cc0 74 28 79 2a 36 34 29 0d 02 58 42 20 20 61 64 64 |t(y*64)..XB add|
00000cd0 20 20 20 20 20 72 31 34 2c 72 38 2c 72 31 32 2c | r14,r8,r12,|
00000ce0 6c 73 6c 20 23 34 20 20 20 3b 20 72 31 34 3d 6f |lsl #4 ; r14=o|
00000cf0 66 66 73 65 74 20 69 6e 20 73 63 72 65 65 6e 20 |ffset in screen |
00000d00 6f 66 20 28 78 3b 79 29 2e 0d 02 62 4a 20 20 6c |of (x;y)...bJ l|
00000d10 64 72 20 20 20 20 20 72 38 2c 5b 72 31 34 2c 72 |dr r8,[r14,r|
00000d20 30 5d 21 20 20 20 20 20 20 20 20 3b 20 72 38 3d |0]! ; r8=|
00000d30 6c 6f 6e 67 20 63 6f 6e 74 61 69 6e 69 6e 67 20 |long containing |
00000d40 70 69 78 65 6c 20 26 20 63 68 61 6e 67 65 20 72 |pixel & change r|
00000d50 31 34 2e 0d 02 6c 3b 20 20 74 73 74 20 20 20 20 |14...l; tst |
00000d60 20 72 38 2c 72 39 20 20 20 20 20 20 20 20 20 20 | r8,r9 |
00000d70 20 20 20 20 20 3b 20 49 73 20 28 78 3b 79 29 20 | ; Is (x;y) |
00000d80 61 6c 72 65 61 64 79 20 64 72 61 77 6e 3f 0d 02 |already drawn?..|
00000d90 76 48 20 20 62 4e 45 20 20 20 20 20 6a 75 6c 69 |vH bNE juli|
00000da0 61 5f 69 69 6d 32 20 20 20 20 20 20 20 20 20 20 |a_iim2 |
00000db0 3b 20 54 68 65 6e 20 73 74 6f 70 20 70 72 6f 63 |; Then stop proc|
00000dc0 65 73 73 69 6e 67 20 74 68 6f 73 65 20 70 6f 69 |essing those poi|
00000dd0 6e 74 73 2e 2e 2e 0d 02 80 3d 20 20 6f 72 72 20 |nts......= orr |
00000de0 20 20 20 20 72 38 2c 72 38 2c 72 39 20 20 20 20 | r8,r8,r9 |
00000df0 20 20 20 20 20 20 20 20 3b 20 45 6c 73 65 20 64 | ; Else d|
00000e00 72 61 77 20 74 68 65 20 74 77 6f 20 70 6f 69 6e |raw the two poin|
00000e10 74 73 2e 0d 02 8a 40 20 20 63 6d 70 20 20 20 20 |ts....@ cmp |
00000e20 20 72 35 2c 23 73 25 20 20 20 20 20 20 20 20 20 | r5,#s% |
00000e30 20 20 20 20 20 3b 20 44 6f 6e 27 20 74 20 64 72 | ; Don' t dr|
00000e40 61 77 20 69 66 20 69 74 65 72 61 74 69 6f 6e 3c |aw if iteration<|
00000e50 73 25 2e 0d 02 94 16 20 20 73 74 72 47 45 20 20 |s%..... strGE |
00000e60 20 72 38 2c 5b 72 31 34 5d 0d 02 9e 3a 20 20 73 | r8,[r14]...: s|
00000e70 75 62 20 20 20 20 20 72 31 34 2c 72 31 30 2c 72 |ub r14,r10,r|
00000e80 31 32 2c 6c 73 6c 20 23 34 20 20 3b 20 72 31 34 |12,lsl #4 ; r14|
00000e90 3d 6f 66 66 73 65 74 20 6f 66 20 28 2d 78 3b 2d |=offset of (-x;-|
00000ea0 79 29 2e 0d 02 a8 4a 20 20 6c 64 72 20 20 20 20 |y)....J ldr |
00000eb0 20 72 38 2c 5b 72 31 34 2c 72 30 5d 21 20 20 20 | r8,[r14,r0]! |
00000ec0 20 20 20 20 20 3b 20 72 38 3d 6c 6f 6e 67 20 63 | ; r8=long c|
00000ed0 6f 6e 74 61 69 6e 69 6e 67 20 70 69 78 65 6c 20 |ontaining pixel |
00000ee0 26 20 63 68 61 6e 67 65 20 72 31 34 2e 0d 02 b2 |& change r14....|
00000ef0 17 20 20 6f 72 72 20 20 20 20 20 72 38 2c 72 38 |. orr r8,r8|
00000f00 2c 72 31 31 0d 02 bc 16 20 20 73 74 72 47 45 20 |,r11.... strGE |
00000f10 20 20 72 38 2c 5b 72 31 34 5d 0d 02 c6 3b 20 20 | r8,[r14]...; |
00000f20 73 74 6d 69 61 20 20 20 72 34 21 2c 7b 72 36 2d |stmia r4!,{r6-|
00000f30 72 37 7d 20 20 20 20 20 20 20 20 20 3b 20 50 75 |r7} ; Pu|
00000f40 74 20 28 78 3b 79 29 20 69 6e 20 74 68 65 20 74 |t (x;y) in the t|
00000f50 61 62 6c 65 2e 0d 02 d0 32 20 20 61 64 64 20 20 |able....2 add |
00000f60 20 20 20 72 35 2c 72 35 2c 23 31 20 20 20 20 20 | r5,r5,#1 |
00000f70 20 20 20 20 20 20 20 3b 20 57 65 6c 6c 2c 20 74 | ; Well, t|
00000f80 65 6c 6c 20 69 74 2e 0d 02 da 3e 20 20 62 20 20 |ell it....> b |
00000f90 20 20 20 20 20 6a 75 6c 69 61 5f 69 69 6d 31 20 | julia_iim1 |
00000fa0 20 20 20 20 20 20 20 20 20 3b 20 43 6f 6e 74 69 | ; Conti|
00000fb0 6e 75 65 20 63 61 6c 63 73 20 77 69 74 68 20 28 |nue calcs with (|
00000fc0 78 3b 79 29 2e 0d 02 e4 0b 5d 3a ed 6f 70 74 25 |x;y).....]:.opt%|
00000fd0 0d 02 ee 44 f4 3d 3d 3d 3d 20 42 41 53 49 43 20 |...D.==== BASIC |
00000fe0 43 4f 44 45 20 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d |CODE ===========|
00000ff0 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d |================|
*
00001010 3d 3d 3d 3d 0d 02 f8 17 eb 30 3a 87 3a c8 91 20 |====.....0:.:.. |
00001020 36 34 31 2c 35 31 33 3a c8 97 ee 0d 03 02 37 de |641,513:......7.|
00001030 76 25 20 38 3a 21 76 25 3d 31 34 38 3a 76 25 21 |v% 8:!v%=148:v%!|
00001040 34 3d 2d 31 3a c8 99 22 4f 53 5f 52 65 61 64 56 |4=-1:.."OS_ReadV|
00001050 64 75 56 61 72 69 61 62 6c 65 73 22 2c 76 25 2c |duVariables",v%,|
00001060 76 25 0d 03 0c 33 76 25 3d 21 76 25 2b 38 30 2a |v%...3v%=!v%+80*|
00001070 31 32 38 2b 34 30 20 20 20 20 20 20 20 20 3a f4 |128+40 :.|
00001080 76 3d 6d 69 64 64 6c 65 20 6f 66 20 76 69 64 65 |v=middle of vide|
00001090 6f 72 61 6d 2e 0d 03 16 38 de 45 25 20 31 30 32 |oram....8.E% 102|
000010a0 34 2a 32 2a 34 20 20 20 20 20 20 20 20 20 20 3a |4*2*4 :|
000010b0 f4 42 75 66 66 65 72 20 66 6f 72 20 27 77 61 69 |.Buffer for 'wai|
000010c0 74 69 6e 67 27 20 70 6f 69 6e 74 73 2e 0d 03 20 |ting' points... |
000010d0 31 62 25 3d 31 20 20 20 20 20 20 20 20 20 20 20 |1b%=1 |
000010e0 20 20 20 20 20 20 20 20 20 3a f4 53 63 72 65 65 | :.Scree|
000010f0 6e 20 62 61 6e 6b 20 6e 75 6d 62 65 72 2e 0d 03 |n bank number...|
00001100 2a 04 0d 03 34 1f f4 43 72 65 61 74 65 20 78 20 |*...4..Create x |
00001110 63 6f 6e 76 65 72 73 69 6f 6e 20 74 61 62 6c 65 |conversion table|
00001120 2e 0d 03 3e 11 de 42 25 20 32 31 2a 33 32 2a 34 |...>..B% 21*32*4|
00001130 2a 34 0d 03 48 09 61 25 3d 42 25 0d 03 52 10 e3 |*4..H.a%=B%..R..|
00001140 69 25 3d 2d 34 30 b8 34 30 88 34 0d 03 5c 2b 20 |i%=-40.40.4..\+ |
00001150 61 25 21 30 3d 69 25 3a 61 25 21 34 3d 31 3a 61 |a%!0=i%:a%!4=1:a|
00001160 25 21 38 3d 2d 69 25 3a 61 25 21 31 32 3d 31 3a |%!8=-i%:a%!12=1:|
00001170 61 25 2b 3d 31 36 0d 03 66 0d 20 e3 6a 25 3d 31 |a%+=16..f. .j%=1|
00001180 b8 33 31 0d 03 70 3b 20 20 61 25 21 30 3d 69 25 |.31..p; a%!0=i%|
00001190 3a 61 25 21 34 3d 31 3c 3c 6a 25 3a 61 25 21 38 |:a%!4=1<<j%:a%!8|
000011a0 3d 2d 69 25 2d 34 3a 61 25 21 31 32 3d 31 3c 3c |=-i%-4:a%!12=1<<|
000011b0 28 33 32 2d 6a 25 29 3a 61 25 2b 3d 31 36 0d 03 |(32-j%):a%+=16..|
000011c0 7a 06 ed 2c 0d 03 84 3c 42 25 2b 3d 31 30 2a 33 |z..,...<B%+=10*3|
000011d0 32 2a 34 2a 34 20 20 20 20 20 20 20 20 20 20 20 |2*4*4 |
000011e0 3a f4 42 25 3d 6d 69 64 64 6c 65 20 6f 66 20 63 |:.B%=middle of c|
000011f0 6f 6e 76 65 72 73 69 6f 6e 20 74 61 62 6c 65 2e |onversion table.|
00001200 0d 03 8e 04 0d 03 98 05 f5 0d 03 a2 37 20 c8 96 |............7 ..|
00001210 3a c8 99 22 4f 53 5f 42 79 74 65 22 2c 26 37 31 |:.."OS_Byte",&71|
00001220 2c 62 25 3a 62 25 3d 62 25 82 33 3a c8 99 22 4f |,b%:b%=b%.3:.."O|
00001230 53 5f 42 79 74 65 22 2c 26 37 30 2c 62 25 3a db |S_Byte",&70,b%:.|
00001240 0d 03 ac 0c 20 c8 97 78 2c 79 2c 7a 0d 03 b6 2b |.... ..x,y,z...+|
00001250 20 f1 8a 33 33 2c 30 29 3b 22 2d 20 4a 75 6c 69 | ..33,0);"- Juli|
00001260 61 20 69 69 6d 20 2d 22 3a f1 8a 37 35 2c 30 29 |a iim -":..75,0)|
00001270 3b 22 62 61 61 68 22 0d 03 c0 36 20 41 25 3d 76 |;"baah"...6 A%=v|
00001280 25 2b 28 28 62 25 2d 31 29 2a 38 30 3c 3c 38 29 |%+((b%-1)*80<<8)|
00001290 20 20 20 20 3a f4 41 25 20 70 6f 69 6e 74 20 6f | :.A% point o|
000012a0 6e 20 77 6f 72 6b 73 63 72 65 65 6e 2e 0d 03 ca |n workscreen....|
000012b0 0c 20 43 25 3d 78 3c 3c 35 0d 03 d4 0c 20 44 25 |. C%=x<<5.... D%|
000012c0 3d 79 3c 3c 35 0d 03 de 10 20 d6 6a 75 6c 69 61 |=y<<5.... .julia|
000012d0 5f 69 69 6d 30 0d 03 e8 09 fd 7a 3c 3e 30 0d 03 |_iim0.....z<>0..|
000012e0 f2 05 e0 0d 03 fc 04 0d 04 06 45 f4 3d 3d 3d 3d |..........E.====|
000012f0 20 4d 41 43 52 4f 53 20 3d 3d 3d 3d 3d 3d 3d 3d | MACROS ========|
00001300 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d |================|
*
00001320 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 3d 0d 04 10 26 |============...&|
00001330 f4 54 68 69 73 20 6d 61 63 72 6f 20 70 75 74 73 |.This macro puts|
00001340 20 73 71 72 74 28 6d 30 29 20 69 6e 74 6f 20 6d | sqrt(m0) into m|
00001350 31 2e 0d 04 1a 14 dd a4 73 71 72 74 28 6d 30 2c |1.......sqrt(m0,|
00001360 6d 31 2c 6d 32 29 0d 04 24 07 ea 6e 25 0d 04 2e |m1,m2)..$..n%...|
00001370 0d 5b 6f 70 74 20 6f 70 74 25 0d 04 38 47 20 20 |.[opt opt%..8G |
00001380 6d 6f 76 20 20 20 20 20 20 20 6d 31 2c 23 30 20 |mov m1,#0 |
00001390 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
000013a0 20 20 3b 20 54 68 69 73 20 69 73 20 72 69 70 70 | ; This is ripp|
000013b0 65 64 20 66 72 6f 6d 20 4a 61 6e 2f 42 41 53 53 |ed from Jan/BASS|
000013c0 2e 0d 04 42 2c 20 20 6d 6f 76 20 20 20 20 20 20 |...B, mov |
000013d0 20 6d 32 2c 23 31 3c 3c 33 30 20 20 20 20 20 20 | m2,#1<<30 |
000013e0 20 20 20 20 20 20 20 20 20 3b 20 20 7c 0d 04 4c | ; |..L|
000013f0 2c 20 20 63 6d 70 20 20 20 20 20 20 20 6d 30 2c |, cmp m0,|
00001400 6d 32 20 20 20 20 20 20 20 20 20 20 20 20 20 20 |m2 |
00001410 20 20 20 20 20 3b 20 20 7c 0d 04 56 2c 20 20 73 | ; |..V, s|
00001420 75 62 48 53 20 20 20 20 20 6d 30 2c 6d 30 2c 6d |ubHS m0,m0,m|
00001430 32 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 |2 |
00001440 20 3b 20 20 7c 0d 04 60 2c 20 20 61 64 63 20 20 | ; |..`, adc |
00001450 20 20 20 20 20 6d 31 2c 6d 31 2c 6d 31 20 20 20 | m1,m1,m1 |
00001460 20 20 20 20 20 20 20 20 20 20 20 20 20 3b 20 20 | ; |
00001470 7c 0d 04 6a 1a 5d 3a e3 6e 25 3d 32 b8 33 30 88 ||..j.]:.n%=2.30.|
00001480 32 3a 5b 6f 70 74 20 6f 70 74 25 0d 04 74 2c 20 |2:[opt opt%..t, |
00001490 20 61 64 64 20 20 20 20 20 20 20 6d 32 2c 6d 31 | add m2,m1|
000014a0 2c 23 31 3c 3c 33 30 20 20 20 20 20 20 20 20 20 |,#1<<30 |
000014b0 20 20 20 3b 20 20 7c 0d 04 7e 2c 20 20 63 6d 70 | ; |..~, cmp|
000014c0 20 20 20 20 20 20 20 6d 30 2c 6d 32 2c 72 6f 72 | m0,m2,ror|
000014d0 20 23 6e 25 20 20 20 20 20 20 20 20 20 20 20 3b | #n% ;|
000014e0 20 20 7c 0d 04 88 2c 20 20 73 75 62 48 53 20 20 | |..., subHS |
000014f0 20 20 20 6d 30 2c 6d 30 2c 6d 32 2c 72 6f 72 20 | m0,m0,m2,ror |
00001500 23 6e 25 20 20 20 20 20 20 20 20 3b 20 20 7c 0d |#n% ; |.|
00001510 04 92 3d 20 20 61 64 63 20 20 20 20 20 20 20 6d |..= adc m|
00001520 31 2c 6d 31 2c 6d 31 20 20 20 20 20 20 20 20 20 |1,m1,m1 |
00001530 20 20 20 20 20 20 20 3b 20 45 6e 64 20 6f 66 20 | ; End of |
00001540 72 69 70 70 65 64 20 63 6f 64 65 2e 0d 04 9c 13 |ripped code.....|
00001550 5d 3a ed 6e 25 3a 5b 6f 70 74 20 6f 70 74 25 0d |]:.n%:[opt opt%.|
00001560 04 a6 39 20 20 63 6d 70 20 20 20 20 20 20 20 6d |..9 cmp m|
00001570 30 2c 6d 31 20 20 20 20 20 20 20 20 20 20 20 20 |0,m1 |
00001580 20 20 20 20 20 20 20 3b 20 46 6c 61 67 73 3d 76 | ; Flags=v|
00001590 61 6c 2d 72 6f 6f 74 2e 0d 04 b0 43 20 20 61 64 |al-root....C ad|
000015a0 64 50 4c 20 20 20 20 20 6d 31 2c 6d 31 2c 23 31 |dPL m1,m1,#1|
000015b0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
000015c0 3b 20 52 6f 75 6e 64 20 74 6f 20 6e 65 61 72 65 |; Round to neare|
000015d0 73 74 20 69 6e 74 65 67 65 72 2e 0d 04 ba 09 5d |st integer.....]|
000015e0 3a 3d 22 22 0d ff |:=""..|
000015e6 
.