Home » Archimedes archive » Archimedes World » AW-1996-09.adf » !AcornAns_AcornAns » Uniform
Uniform
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 » Archimedes World » AW-1996-09.adf » !AcornAns_AcornAns |
| Filename: | Uniform |
| Read OK: | ✔ |
| File size: | 0B44 bytes |
| Load address: | 0000 |
| Exec address: | 0000 |
File contents
10REM Demo of pseudo random number generation for a random variable
20REM with uniform distribution on {0,1,2,...,n%-1}.
30REM Arbitrary n%, set to 640 below.
40:
50PROCreadmonitortype
60MODE mo%
70PROCass
80:
90!seed%=-TIME
100A%=640
110CALL calc_m%
120DIM sum%(A%-1)
130scale=2
140REPEAT
150 r% = USR rnd_uniform%
160 sum%(r%) += 1
170 y%=sum%(r%)*scale
180 IF y%<1024 THEN
190 PLOT 69, r%*2, y%
200 ELSE
210 scale=scale/2
220 CLS
230 FOR i%=0 TO A%-1
240 IF sum%(i%) LINE i%*2, 0, i%*2, sum%(i%)*scale
250 NEXT
260 ENDIF
270UNTIL FALSE
280:
290END
300:
310:
320DEF PROCass
330DIM code% 1024
340FOR pass%=0 TO 2 STEP 2
350P%=code%
360[OPT pass%
370:
380.addr_n% EQUD 0
390.addr_m% EQUD 0
400.seed% EQUD -1:EQUD -1
410:
420.rnd_uniform% ;Purpose - generate a random number
430: ; in the range 0,...,n%-1
440:
450: ;Define m% = 1<<k%, where k% is the
460: ;minimal integer such that m%>=n%.
470:
480: ;Must have already stored n% & m%-1
490: ;in addr_n% & addr_m%, prior to
500: ;calling rnd_uniform%.
510: ;NB calc_m% will compute m% for you
520: ; if needed.
530ADR 0, addr_n%
540LDMIA 0, {1,2,3,4} ;r1=n%, r2=m%-1, r3,r4=33-bit seed
550: ; (b0-b31 in r3, b32 in lsb r4)
560.rnd_uniform_gen%
570MOVS 4, 4, LSR #1 ;generate random 32-bit number
580MOVS 0, 3, RRX
590ADC 4, 4, 4
600EOR 0, 0, 3, LSL #12
610EOR 3, 0, 0, LSR #20
620AND 0, 3, 2 ;mask 32-bit number by m%-1
630CMP 0, 1
640BHS rnd_uniform_gen% ;try again if not quite in range
650:
660: ;r0 = random number, range 0,...n%-1
670ADR 1, seed% ;note above loop takes average of
680STMIA 1, {3,4} ;8 to 19 cycles to execute depending
690MOVS PC, 14 ;on n%/m% (1 >= n%/m% > 0.5)
700:
710:
720.calc_m% ;Expects n% in r0, calcs m% & stores n%,m%-1
730CMP 0, #0
740MOVEQ 0, #1 ;a bodge to deal with degenerate case n%=0
750SUBS 3, 0, #1
760MOV 1, #0
770BEQ calc_m_done% ;n%=1
780MOV 2, #1
790MOVS 3, 3, LSR #1
800ADD 1, 2, 1, LSL #1
810BEQ calc_m_done% ;n%=2
820]
830FOR i%=2 TO 30
840[OPT pass%
850MOVS 3, 3, LSR #1
860ADD 1, 2, 1, LSL #1
870BEQ calc_m_done% ;2^(i%-1) < n% <= 2^i%
880]
890NEXT
900[OPT pass%
910MOVS 3, 3, LSR #1
920ADD 1, 2, 1, LSL #1 ;2^30 < n%
930ADDNE 1, 2, 1, LSL #1 ;2^31 < n% < 2^32
940.calc_m_done%
950ADR 2, addr_n%
960STMIA 2, {0,1}
970MOVS PC, 14
980:
990]
1000NEXT
1010ENDPROC
1020:
1030DEF PROCreadmonitortype
1040SYS "OS_Byte",161,133 TO ,,type%
1050type%=(type% AND 12) DIV 4
1060IF type%=1 mo%=18 ELSE mo%=0
1070ENDPROC
C� Demo of pseudo random number generation for a random variable
4� with uniform distribution on {0,1,2,...,n%-1}.
%� Arbitrary n%, set to 640 below.
(:
2�readmonitortype
< � mo%
F�ass
P:
Z
!seed%=-�
d
A%=640
n
� calc_m%
x� sum%(A%-1)
� scale=2
��
� r% = � rnd_uniform%
� sum%(r%) += 1
� y%=sum%(r%)*scale
� � y%<1024 �
� � 69, r%*2, y%
� �
� scale=scale/2
� �
� � i%=0 � A%-1
�1 � sum%(i%) � i%*2, 0, i%*2, sum%(i%)*scale
� �
�
� �
:
"�
,:
6:
@
� �ass
J� code% 1024
T� pass%=0 � 2 � 2
^ P%=code%
h[OPT pass%
r:
|.addr_n% EQUD 0
�.addr_m% EQUD 0
� .seed% EQUD -1:EQUD -1
�:
�E.rnd_uniform% ;Purpose - generate a random number
�D: ; in the range 0,...,n%-1
�:
�E: ;Define m% = 1<<k%, where k% is the
�D: ;minimal integer such that m%>=n%.
�:
�E: ;Must have already stored n% & m%-1
�A: ;in addr_n% & addr_m%, prior to
�8: ;calling rnd_uniform%.
�E: ;NB calc_m% will compute m% for you
0: ; if needed.
ADR 0, addr_n%
DLDMIA 0, {1,2,3,4} ;r1=n%, r2=m%-1, r3,r4=33-bit seed
&D: ; (b0-b31 in r3, b32 in lsb r4)
0.rnd_uniform_gen%
:@MOVS 4, 4, LSR #1 ;generate random 32-bit number
DMOVS 0, 3, RRX
NADC 4, 4, 4
X � 0, 0, 3, LSL #12
b � 3, 0, 0, LSR #20
l;� 0, 3, 2 ;mask 32-bit number by m%-1
vCMP 0, 1
�BBHS rnd_uniform_gen% ;try again if not quite in range
�:
�F: ;r0 = random number, range 0,...n%-1
�CADR 1, seed% ;note above loop takes average of
�FSTMIA 1, {3,4} ;8 to 19 cycles to execute depending
�>MOVS PC, 14 ;on n%/m% (1 >= n%/m% > 0.5)
�:
�:
�D.calc_m% ;Expects n% in r0, calcs m% & stores n%,m%-1
�CMP 0, #0
�?�Q 0, #1 ;a bodge to deal with degenerate case n%=0
�SUBS 3, 0, #1
�MOV 1, #0
'BEQ calc_m_done% ;n%=1
MOV 2, #1
MOVS 3, 3, LSR #1
ADD 1, 2, 1, LSL #1
*'BEQ calc_m_done% ;n%=2
4]
>� i%=2 � 30
H[OPT pass%
RMOVS 3, 3, LSR #1
\ ADD 1, 2, 1, LSL #1
f8BEQ calc_m_done% ;2^(i%-1) < n% <= 2^i%
p]
z�
�[OPT pass%
�MOVS 3, 3, LSR #1
�,ADD 1, 2, 1, LSL #1 ;2^30 < n%
�3ADDNE 1, 2, 1, LSL #1 ;2^31 < n% < 2^32
�.calc_m_done%
�ADR 2, addr_n%
�STMIA 2, {0,1}
�MOVS PC, 14
�:
�]
��
��
�:
� �readmonitortype
"ș "OS_Byte",161,133 � ,,type%
type%=(type% � 12) � 4
$ � type%=1 mo%=18 � mo%=0
.�
� 00000000 0d 00 0a 43 f4 20 44 65 6d 6f 20 6f 66 20 70 73 |...C. Demo of ps|
00000010 65 75 64 6f 20 72 61 6e 64 6f 6d 20 6e 75 6d 62 |eudo random numb|
00000020 65 72 20 67 65 6e 65 72 61 74 69 6f 6e 20 66 6f |er generation fo|
00000030 72 20 61 20 72 61 6e 64 6f 6d 20 76 61 72 69 61 |r a random varia|
00000040 62 6c 65 0d 00 14 34 f4 20 77 69 74 68 20 75 6e |ble...4. with un|
00000050 69 66 6f 72 6d 20 64 69 73 74 72 69 62 75 74 69 |iform distributi|
00000060 6f 6e 20 6f 6e 20 7b 30 2c 31 2c 32 2c 2e 2e 2e |on on {0,1,2,...|
00000070 2c 6e 25 2d 31 7d 2e 0d 00 1e 25 f4 20 41 72 62 |,n%-1}....%. Arb|
00000080 69 74 72 61 72 79 20 6e 25 2c 20 73 65 74 20 74 |itrary n%, set t|
00000090 6f 20 36 34 30 20 62 65 6c 6f 77 2e 0d 00 28 05 |o 640 below...(.|
000000a0 3a 0d 00 32 14 f2 72 65 61 64 6d 6f 6e 69 74 6f |:..2..readmonito|
000000b0 72 74 79 70 65 0d 00 3c 09 eb 20 6d 6f 25 0d 00 |rtype..<.. mo%..|
000000c0 46 08 f2 61 73 73 0d 00 50 05 3a 0d 00 5a 0d 21 |F..ass..P.:..Z.!|
000000d0 73 65 65 64 25 3d 2d 91 0d 00 64 0a 41 25 3d 36 |seed%=-...d.A%=6|
000000e0 34 30 0d 00 6e 0d d6 20 63 61 6c 63 5f 6d 25 0d |40..n.. calc_m%.|
000000f0 00 78 10 de 20 73 75 6d 25 28 41 25 2d 31 29 0d |.x.. sum%(A%-1).|
00000100 00 82 0b 73 63 61 6c 65 3d 32 0d 00 8c 05 f5 0d |...scale=2......|
00000110 00 96 18 20 72 25 20 3d 20 ba 20 72 6e 64 5f 75 |... r% = . rnd_u|
00000120 6e 69 66 6f 72 6d 25 0d 00 a0 12 20 73 75 6d 25 |niform%.... sum%|
00000130 28 72 25 29 20 2b 3d 20 31 0d 00 aa 16 20 79 25 |(r%) += 1.... y%|
00000140 3d 73 75 6d 25 28 72 25 29 2a 73 63 61 6c 65 0d |=sum%(r%)*scale.|
00000150 00 b4 10 20 e7 20 79 25 3c 31 30 32 34 20 8c 0d |... . y%<1024 ..|
00000160 00 be 14 20 20 f0 20 36 39 2c 20 72 25 2a 32 2c |... . 69, r%*2,|
00000170 20 79 25 0d 00 c8 06 20 cc 0d 00 d2 13 20 20 73 | y%.... ..... s|
00000180 63 61 6c 65 3d 73 63 61 6c 65 2f 32 0d 00 dc 07 |cale=scale/2....|
00000190 20 20 db 0d 00 e6 13 20 20 e3 20 69 25 3d 30 20 | ..... . i%=0 |
000001a0 b8 20 41 25 2d 31 0d 00 f0 31 20 20 20 e7 20 73 |. A%-1...1 . s|
000001b0 75 6d 25 28 69 25 29 20 86 20 69 25 2a 32 2c 20 |um%(i%) . i%*2, |
000001c0 30 2c 20 69 25 2a 32 2c 20 73 75 6d 25 28 69 25 |0, i%*2, sum%(i%|
000001d0 29 2a 73 63 61 6c 65 0d 00 fa 07 20 20 ed 0d 01 |)*scale.... ...|
000001e0 04 06 20 cd 0d 01 0e 07 fd 20 a3 0d 01 18 05 3a |.. ...... .....:|
000001f0 0d 01 22 05 e0 0d 01 2c 05 3a 0d 01 36 05 3a 0d |.."....,.:..6.:.|
00000200 01 40 0a dd 20 f2 61 73 73 0d 01 4a 10 de 20 63 |.@.. .ass..J.. c|
00000210 6f 64 65 25 20 31 30 32 34 0d 01 54 15 e3 20 70 |ode% 1024..T.. p|
00000220 61 73 73 25 3d 30 20 b8 20 32 20 88 20 32 0d 01 |ass%=0 . 2 . 2..|
00000230 5e 0c 50 25 3d 63 6f 64 65 25 0d 01 68 0e 5b 4f |^.P%=code%..h.[O|
00000240 50 54 20 70 61 73 73 25 0d 01 72 05 3a 0d 01 7c |PT pass%..r.:..||
00000250 14 2e 61 64 64 72 5f 6e 25 20 20 45 51 55 44 20 |..addr_n% EQUD |
00000260 30 0d 01 86 14 2e 61 64 64 72 5f 6d 25 20 20 45 |0.....addr_m% E|
00000270 51 55 44 20 30 0d 01 90 1d 2e 73 65 65 64 25 20 |QUD 0.....seed% |
00000280 20 20 20 45 51 55 44 20 2d 31 3a 45 51 55 44 20 | EQUD -1:EQUD |
00000290 2d 31 0d 01 9a 05 3a 0d 01 a4 45 2e 72 6e 64 5f |-1....:...E.rnd_|
000002a0 75 6e 69 66 6f 72 6d 25 20 20 20 20 20 20 20 20 |uniform% |
000002b0 20 20 20 20 20 20 20 20 20 3b 50 75 72 70 6f 73 | ;Purpos|
000002c0 65 20 2d 20 67 65 6e 65 72 61 74 65 20 61 20 72 |e - generate a r|
000002d0 61 6e 64 6f 6d 20 6e 75 6d 62 65 72 0d 01 ae 44 |andom number...D|
000002e0 3a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 |: |
000002f0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 3b 20 | ; |
00000300 20 20 20 20 20 20 20 20 20 69 6e 20 74 68 65 20 | in the |
00000310 72 61 6e 67 65 20 30 2c 2e 2e 2e 2c 6e 25 2d 31 |range 0,...,n%-1|
00000320 0d 01 b8 05 3a 0d 01 c2 45 3a 20 20 20 20 20 20 |....:...E: |
00000330 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
00000340 20 20 20 20 20 20 20 3b 44 65 66 69 6e 65 20 6d | ;Define m|
00000350 25 20 3d 20 31 3c 3c 6b 25 2c 20 77 68 65 72 65 |% = 1<<k%, where|
00000360 20 6b 25 20 69 73 20 74 68 65 0d 01 cc 44 3a 20 | k% is the...D: |
00000370 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
00000380 20 20 20 20 20 20 20 20 20 20 20 20 3b 6d 69 6e | ;min|
00000390 69 6d 61 6c 20 69 6e 74 65 67 65 72 20 73 75 63 |imal integer suc|
000003a0 68 20 74 68 61 74 20 6d 25 3e 3d 6e 25 2e 0d 01 |h that m%>=n%...|
000003b0 d6 05 3a 0d 01 e0 45 3a 20 20 20 20 20 20 20 20 |..:...E: |
000003c0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
000003d0 20 20 20 20 20 3b 4d 75 73 74 20 68 61 76 65 20 | ;Must have |
000003e0 61 6c 72 65 61 64 79 20 73 74 6f 72 65 64 20 6e |already stored n|
000003f0 25 20 26 20 6d 25 2d 31 0d 01 ea 41 3a 20 20 20 |% & m%-1...A: |
00000400 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
00000410 20 20 20 20 20 20 20 20 20 20 3b 69 6e 20 61 64 | ;in ad|
00000420 64 72 5f 6e 25 20 26 20 61 64 64 72 5f 6d 25 2c |dr_n% & addr_m%,|
00000430 20 70 72 69 6f 72 20 74 6f 0d 01 f4 38 3a 20 20 | prior to...8: |
00000440 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
00000450 20 20 20 20 20 20 20 20 20 20 20 3b 63 61 6c 6c | ;call|
00000460 69 6e 67 20 72 6e 64 5f 75 6e 69 66 6f 72 6d 25 |ing rnd_uniform%|
00000470 2e 0d 01 fe 45 3a 20 20 20 20 20 20 20 20 20 20 |....E: |
00000480 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
00000490 20 20 20 3b 4e 42 20 63 61 6c 63 5f 6d 25 20 77 | ;NB calc_m% w|
000004a0 69 6c 6c 20 63 6f 6d 70 75 74 65 20 6d 25 20 66 |ill compute m% f|
000004b0 6f 72 20 79 6f 75 0d 02 08 30 3a 20 20 20 20 20 |or you...0: |
000004c0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
000004d0 20 20 20 20 20 20 20 20 3b 20 20 20 69 66 20 6e | ; if n|
000004e0 65 65 64 65 64 2e 0d 02 12 18 41 44 52 20 20 20 |eeded.....ADR |
000004f0 20 20 20 20 30 2c 20 61 64 64 72 5f 6e 25 0d 02 | 0, addr_n%..|
00000500 1c 44 4c 44 4d 49 41 20 20 20 20 20 30 2c 20 7b |.DLDMIA 0, {|
00000510 31 2c 32 2c 33 2c 34 7d 20 20 20 20 20 20 20 20 |1,2,3,4} |
00000520 3b 72 31 3d 6e 25 2c 20 72 32 3d 6d 25 2d 31 2c |;r1=n%, r2=m%-1,|
00000530 20 72 33 2c 72 34 3d 33 33 2d 62 69 74 20 73 65 | r3,r4=33-bit se|
00000540 65 64 0d 02 26 44 3a 20 20 20 20 20 20 20 20 20 |ed..&D: |
00000550 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
00000560 20 20 20 20 3b 20 20 20 20 28 62 30 2d 62 33 31 | ; (b0-b31|
00000570 20 69 6e 20 72 33 2c 20 62 33 32 20 69 6e 20 6c | in r3, b32 in l|
00000580 73 62 20 72 34 29 0d 02 30 15 2e 72 6e 64 5f 75 |sb r4)..0..rnd_u|
00000590 6e 69 66 6f 72 6d 5f 67 65 6e 25 0d 02 3a 40 4d |niform_gen%..:@M|
000005a0 4f 56 53 20 20 20 20 20 20 34 2c 20 34 2c 20 4c |OVS 4, 4, L|
000005b0 53 52 20 23 31 20 20 20 20 20 20 20 20 3b 67 65 |SR #1 ;ge|
000005c0 6e 65 72 61 74 65 20 72 61 6e 64 6f 6d 20 33 32 |nerate random 32|
000005d0 2d 62 69 74 20 6e 75 6d 62 65 72 0d 02 44 17 4d |-bit number..D.M|
000005e0 4f 56 53 20 20 20 20 20 20 30 2c 20 33 2c 20 52 |OVS 0, 3, R|
000005f0 52 58 0d 02 4e 15 41 44 43 20 20 20 20 20 20 20 |RX..N.ADC |
00000600 34 2c 20 34 2c 20 34 0d 02 58 1c 82 20 20 20 20 |4, 4, 4..X.. |
00000610 20 20 20 30 2c 20 30 2c 20 33 2c 20 4c 53 4c 20 | 0, 0, 3, LSL |
00000620 23 31 32 0d 02 62 1c 82 20 20 20 20 20 20 20 33 |#12..b.. 3|
00000630 2c 20 30 2c 20 30 2c 20 4c 53 52 20 23 32 30 0d |, 0, 0, LSR #20.|
00000640 02 6c 3b 80 20 20 20 20 20 20 20 30 2c 20 33 2c |.l;. 0, 3,|
00000650 20 32 20 20 20 20 20 20 20 20 20 20 20 20 20 3b | 2 ;|
00000660 6d 61 73 6b 20 33 32 2d 62 69 74 20 6e 75 6d 62 |mask 32-bit numb|
00000670 65 72 20 62 79 20 6d 25 2d 31 0d 02 76 12 43 4d |er by m%-1..v.CM|
00000680 50 20 20 20 20 20 20 20 30 2c 20 31 0d 02 80 42 |P 0, 1...B|
00000690 42 48 53 20 20 20 20 20 20 20 72 6e 64 5f 75 6e |BHS rnd_un|
000006a0 69 66 6f 72 6d 5f 67 65 6e 25 20 20 20 20 3b 74 |iform_gen% ;t|
000006b0 72 79 20 61 67 61 69 6e 20 69 66 20 6e 6f 74 20 |ry again if not |
000006c0 71 75 69 74 65 20 69 6e 20 72 61 6e 67 65 0d 02 |quite in range..|
000006d0 8a 05 3a 0d 02 94 46 3a 20 20 20 20 20 20 20 20 |..:...F: |
000006e0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
000006f0 20 20 20 20 20 3b 72 30 20 3d 20 72 61 6e 64 6f | ;r0 = rando|
00000700 6d 20 6e 75 6d 62 65 72 2c 20 72 61 6e 67 65 20 |m number, range |
00000710 30 2c 2e 2e 2e 6e 25 2d 31 0d 02 9e 43 41 44 52 |0,...n%-1...CADR|
00000720 20 20 20 20 20 20 20 31 2c 20 73 65 65 64 25 20 | 1, seed% |
00000730 20 20 20 20 20 20 20 20 20 20 20 3b 6e 6f 74 65 | ;note|
00000740 20 61 62 6f 76 65 20 6c 6f 6f 70 20 74 61 6b 65 | above loop take|
00000750 73 20 61 76 65 72 61 67 65 20 6f 66 0d 02 a8 46 |s average of...F|
00000760 53 54 4d 49 41 20 20 20 20 20 31 2c 20 7b 33 2c |STMIA 1, {3,|
00000770 34 7d 20 20 20 20 20 20 20 20 20 20 20 20 3b 38 |4} ;8|
00000780 20 74 6f 20 31 39 20 63 79 63 6c 65 73 20 74 6f | to 19 cycles to|
00000790 20 65 78 65 63 75 74 65 20 64 65 70 65 6e 64 69 | execute dependi|
000007a0 6e 67 0d 02 b2 3e 4d 4f 56 53 20 20 20 20 20 20 |ng...>MOVS |
000007b0 50 43 2c 20 31 34 20 20 20 20 20 20 20 20 20 20 |PC, 14 |
000007c0 20 20 20 20 3b 6f 6e 20 6e 25 2f 6d 25 20 28 31 | ;on n%/m% (1|
000007d0 20 3e 3d 20 6e 25 2f 6d 25 20 3e 20 30 2e 35 29 | >= n%/m% > 0.5)|
000007e0 0d 02 bc 05 3a 0d 02 c6 05 3a 0d 02 d0 44 2e 63 |....:....:...D.c|
000007f0 61 6c 63 5f 6d 25 20 20 20 20 20 20 20 20 20 20 |alc_m% |
00000800 20 20 3b 45 78 70 65 63 74 73 20 6e 25 20 69 6e | ;Expects n% in|
00000810 20 72 30 2c 20 63 61 6c 63 73 20 6d 25 20 26 20 | r0, calcs m% & |
00000820 73 74 6f 72 65 73 20 6e 25 2c 6d 25 2d 31 0d 02 |stores n%,m%-1..|
00000830 da 13 43 4d 50 20 20 20 20 20 20 20 30 2c 20 23 |..CMP 0, #|
00000840 30 0d 02 e4 3f ec 51 20 20 20 20 20 30 2c 20 23 |0...?.Q 0, #|
00000850 31 20 20 20 20 20 3b 61 20 62 6f 64 67 65 20 74 |1 ;a bodge t|
00000860 6f 20 64 65 61 6c 20 77 69 74 68 20 64 65 67 65 |o deal with dege|
00000870 6e 65 72 61 74 65 20 63 61 73 65 20 6e 25 3d 30 |nerate case n%=0|
00000880 0d 02 ee 16 53 55 42 53 20 20 20 20 20 20 33 2c |....SUBS 3,|
00000890 20 30 2c 20 23 31 0d 02 f8 13 4d 4f 56 20 20 20 | 0, #1....MOV |
000008a0 20 20 20 20 31 2c 20 23 30 0d 03 02 27 42 45 51 | 1, #0...'BEQ|
000008b0 20 20 20 20 20 20 20 63 61 6c 63 5f 6d 5f 64 6f | calc_m_do|
000008c0 6e 65 25 20 20 20 20 20 20 20 20 3b 6e 25 3d 31 |ne% ;n%=1|
000008d0 0d 03 0c 13 4d 4f 56 20 20 20 20 20 20 20 32 2c |....MOV 2,|
000008e0 20 23 31 0d 03 16 1a 4d 4f 56 53 20 20 20 20 20 | #1....MOVS |
000008f0 20 33 2c 20 33 2c 20 4c 53 52 20 23 31 0d 03 20 | 3, 3, LSR #1.. |
00000900 1d 41 44 44 20 20 20 20 20 20 20 31 2c 20 32 2c |.ADD 1, 2,|
00000910 20 31 2c 20 4c 53 4c 20 23 31 0d 03 2a 27 42 45 | 1, LSL #1..*'BE|
00000920 51 20 20 20 20 20 20 20 63 61 6c 63 5f 6d 5f 64 |Q calc_m_d|
00000930 6f 6e 65 25 20 20 20 20 20 20 20 20 3b 6e 25 3d |one% ;n%=|
00000940 32 0d 03 34 05 5d 0d 03 3e 0f e3 20 69 25 3d 32 |2..4.]..>.. i%=2|
00000950 20 b8 20 33 30 0d 03 48 0e 5b 4f 50 54 20 70 61 | . 30..H.[OPT pa|
00000960 73 73 25 0d 03 52 1a 4d 4f 56 53 20 20 20 20 20 |ss%..R.MOVS |
00000970 20 33 2c 20 33 2c 20 4c 53 52 20 23 31 0d 03 5c | 3, 3, LSR #1..\|
00000980 1d 41 44 44 20 20 20 20 20 20 20 31 2c 20 32 2c |.ADD 1, 2,|
00000990 20 31 2c 20 4c 53 4c 20 23 31 0d 03 66 38 42 45 | 1, LSL #1..f8BE|
000009a0 51 20 20 20 20 20 20 20 63 61 6c 63 5f 6d 5f 64 |Q calc_m_d|
000009b0 6f 6e 65 25 20 20 20 20 20 20 20 20 3b 32 5e 28 |one% ;2^(|
000009c0 69 25 2d 31 29 20 3c 20 6e 25 20 3c 3d 20 32 5e |i%-1) < n% <= 2^|
000009d0 69 25 0d 03 70 05 5d 0d 03 7a 05 ed 0d 03 84 0e |i%..p.]..z......|
000009e0 5b 4f 50 54 20 70 61 73 73 25 0d 03 8e 1a 4d 4f |[OPT pass%....MO|
000009f0 56 53 20 20 20 20 20 20 33 2c 20 33 2c 20 4c 53 |VS 3, 3, LS|
00000a00 52 20 23 31 0d 03 98 2c 41 44 44 20 20 20 20 20 |R #1...,ADD |
00000a10 20 20 31 2c 20 32 2c 20 31 2c 20 4c 53 4c 20 23 | 1, 2, 1, LSL #|
00000a20 31 20 20 20 20 20 3b 32 5e 33 30 20 3c 20 6e 25 |1 ;2^30 < n%|
00000a30 0d 03 a2 33 41 44 44 4e 45 20 20 20 20 20 31 2c |...3ADDNE 1,|
00000a40 20 32 2c 20 31 2c 20 4c 53 4c 20 23 31 20 20 20 | 2, 1, LSL #1 |
00000a50 20 20 3b 32 5e 33 31 20 3c 20 6e 25 20 3c 20 32 | ;2^31 < n% < 2|
00000a60 5e 33 32 0d 03 ac 11 2e 63 61 6c 63 5f 6d 5f 64 |^32.....calc_m_d|
00000a70 6f 6e 65 25 0d 03 b6 18 41 44 52 20 20 20 20 20 |one%....ADR |
00000a80 20 20 32 2c 20 61 64 64 72 5f 6e 25 0d 03 c0 16 | 2, addr_n%....|
00000a90 53 54 4d 49 41 20 20 20 20 20 32 2c 20 7b 30 2c |STMIA 2, {0,|
00000aa0 31 7d 0d 03 ca 14 4d 4f 56 53 20 20 20 20 20 20 |1}....MOVS |
00000ab0 50 43 2c 20 31 34 0d 03 d4 05 3a 0d 03 de 05 5d |PC, 14....:....]|
00000ac0 0d 03 e8 05 ed 0d 03 f2 05 e1 0d 03 fc 05 3a 0d |..............:.|
00000ad0 04 06 16 dd 20 f2 72 65 61 64 6d 6f 6e 69 74 6f |.... .readmonito|
00000ae0 72 74 79 70 65 0d 04 10 22 c8 99 20 22 4f 53 5f |rtype...".. "OS_|
00000af0 42 79 74 65 22 2c 31 36 31 2c 31 33 33 20 b8 20 |Byte",161,133 . |
00000b00 2c 2c 74 79 70 65 25 0d 04 1a 1a 74 79 70 65 25 |,,type%....type%|
00000b10 3d 28 74 79 70 65 25 20 80 20 31 32 29 20 81 20 |=(type% . 12) . |
00000b20 34 0d 04 24 1c e7 20 74 79 70 65 25 3d 31 20 6d |4..$.. type%=1 m|
00000b30 6f 25 3d 31 38 20 8b 20 6d 6f 25 3d 30 0d 04 2e |o%=18 . mo%=0...|
00000b40 05 e1 0d ff |....|
00000b44 
.