Home » Archimedes archive » Archimedes World » AW-1997-02.adf » !AcornAns_AcornAns » Prime/Primes2
Prime/Primes2
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-1997-02.adf » !AcornAns_AcornAns |
| Filename: | Prime/Primes2 |
| Read OK: | ✔ |
| File size: | 0CB4 bytes |
| Load address: | 0000 |
| Exec address: | 0000 |
Duplicates
There is 1 duplicate copy of this file in the archive:
- Archimedes archive » Archimedes World » AW-1997-02.adf » !AcornAns_AcornAns » Prime/Primes2
- Archimedes archive » Archimedes World » AW-1997-01.adf » !AcornAns_AcornAns » Prime/Primes2
File contents
10maxn=2^32
20np%=23+1.115*2*SQR(maxn)/LN(maxn) :REM approx upper bound on #primes<=SQR(maxn)
30DIM multiple% 4*np%, p% 4*np%
40:
50PROCass
60:
70REM PROCprimes:
80REM Count all primes less than or equal to second parameter, and
90REM display all those which are also greater than or equal to the
100REM first parameter.
110:
120REM To see how quick this code is at finding primes, set both
130REM parameters equal, so that no time is wasted displaying them.
140:
150TIME=0
160PROCprimes(65536*2, 65536*4)
170PRINT "Time: ";TIME/100
180:
190END
200:
210:
220DEF PROCass
230DIM code% 1024
240FOR pass%=0 TO 2 STEP 2
250P%=code%
260[OPT pass%
270.amultiple% EQUD multiple%
280.ap% EQUD p%
290.genprimes%
300MOV 6, 1
310MOV 10, 2
320MOV 1, #0 ;anp (actual number of primes found)
330CMP 0, #1
340MOVLS 0, 1
350MOVLSS PC, 14
360STMFD 13!, {14}
370MOV 2, #2
380MOV 3, #3
390MOV 4, #5
400LDR 5, ap%
410STMIA 5, {2,3,4}
420MOV 11, #2
430CMP 0, #5
440BHS gp1%
450MOV 2, #0
460SUB 7, 0, #1
470MOV 7, 7, LSR #1
480.gp2%
490LDR 3, [5, 2, LSL #2]
500CMP 3, 6
510BLHS outputprime%
520ADD 1, 1, #1
530ADD 2, 2, #1
540CMP 2, 7
550BLS gp2%
560MOV 0, 1
570LDMFD 13!, {PC}^
580.gp1%
590MOV 2, #0
600MOV 7, #2
610.gp3%
620LDR 3, [5, 2, LSL #2]
630CMP 3, 6
640BLHS outputprime%
650ADD 1, 1, #1
660ADD 2, 2, #1
670CMP 2, 7
680BLS gp3%
690:
700MOV 3, #5
710MOV 7, #25
720MOV 8, #2
730MOV 9, #2
740LDR 12, amultiple%
750:
760.gpmainloop%
770AND 14, 3, #255
780CMP 14, #2
790BHS gp8%
800SWI "OS_ReadEscapeState"
810LDMCSFD 13!, {PC}^
820.gp8%
830ADD 3, 3, 9
840CMP 3, 0
850MOVHI 0, 1
860LDMHIFD 13!, {PC}^
870RSB 9, 9, #6
880CMP 8, 11
890CMPLS 7, 3
900BHI gp4%
910STR 7, [12, 8, LSL #2]
920ADD 8, 8, #1
930CMP 8, 11
940LDRLS 14, [5, 8, LSL #2]
950MULLS 7, 14, 14
960.gp4%
970MOV 2, #2
980CMP 2, 8
990BHS gpprime%
1000.gp6%
1010LDR 14, [12, 2, LSL #2]
1020LDR 4, [5, 2, LSL #2]
1030.gp7%
1040CMP 14, 3
1050ADDLO 14, 14, 4, LSL #1
1060BLO gp7%
1070STR 14, [12, 2, LSL #2]
1080MVNEQ 2, #1
1090ADD 2, 2, #1
1100.gp5%
1110CMP 2, 8
1120BLO gp6% ;NOTprime on loop exit iff R2=&ffffffff iff NE
1130BNE gpmainloop%
1140.gpprime%
1150CMP 3, 6
1160BLHS outputprime%
1170ADD 1, 1, #1
1180CMP 3, 10
1190ADDLS 11, 11, #1
1200STRLS 3, [5, 11, LSL #2]
1210B gpmainloop%
1220:
1230.buf% EQUD 0:EQUD 0:EQUD 0
1240.bufend%
1250.outputprime% ;you might want to alter this routine to
1260 ;dump prime into an array instead of print it!
1270STMFD 13!, {0,1,2}
1280MOV 0, 3
1290ADR 1, buf%
1300MOV 2, #bufend%-buf%
1310SWI "OS_ConvertCardinal4"
1320SWI "OS_Write0"
1330SWI "OS_NewLine"
1340LDMFD 13!, {0,1,2}
1350MOV PC, 14
1360]
1370NEXT
1380ENDPROC
1390:
1400DEF PROCprimes(nlow%, n%)
1410LOCAL anp%, v
1420A%=n%
1430B%=nlow%
1440IF n%<0 v=n%+(2^32) ELSE v=n%
1450C%=INT SQR(v)
1460anp%=USR genprimes%
1470PRINT'"Found ";anp%;" prime numbers less than or equal to ";n%'
1480ENDPROC
maxn=2^32
Nnp%=23+1.115*2*�(maxn)/�(maxn) :� approx upper bound on #primes<=SQR(maxn)
� multiple% 4*np%, p% 4*np%
(:
2�ass
<:
F� PROCprimes:
PB� Count all primes less than or equal to second parameter, and
ZC� display all those which are also greater than or equal to the
d� first parameter.
n:
x?� To see how quick this code is at finding primes, set both
�B� parameters equal, so that no time is wasted displaying them.
�:
��=0
� �primes(65536*2, 65536*4)
�� "Time: ";�/100
�:
��
�:
�:
�
� �ass
�� code% 1024
�� pass%=0 � 2 � 2
� P%=code%
[OPT pass%
&.amultiple% EQUD multiple%
.ap% EQUD p%
".genprimes%
,MOV 6, 1
6MOV 10, 2
@<MOV 1, #0 ;anp (actual number of primes found)
JCMP 0, #1
TMOVLS 0, 1
^MOVLSS PC, 14
hSTMFD 13!, {14}
rMOV 2, #2
|MOV 3, #3
�MOV 4, #5
�LDR 5, ap%
�STMIA 5, {2,3,4}
�MOV 11, #2
�CMP 0, #5
�BHS gp1%
�MOV 2, #0
�SUB 7, 0, #1
�MOV 7, 7, LSR #1
� .gp2%
� LDR 3, [5, 2, LSL #2]
�CMP 3, 6
�BLHS outputprime%
ADD 1, 1, #1
ADD 2, 2, #1
CMP 2, 7
&BLS gp2%
0MOV 0, 1
:LDMFD 13!, {PC}^
D .gp1%
NMOV 2, #0
XMOV 7, #2
b .gp3%
l LDR 3, [5, 2, LSL #2]
vCMP 3, 6
�BLHS outputprime%
�ADD 1, 1, #1
�ADD 2, 2, #1
�CMP 2, 7
�BLS gp3%
�:
�MOV 3, #5
�MOV 7, #25
�MOV 8, #2
�MOV 9, #2
� LDR 12, amultiple%
�:
�.gpmainloop%
� 14, 3, #255
CMP 14, #2
BHS gp8%
"SWI "OS_ReadEscapeState"
*LDMCSFD 13!, {PC}^
4 .gp8%
>ADD 3, 3, 9
HCMP 3, 0
RMOVHI 0, 1
\LDMHIFD 13!, {PC}^
fRSB 9, 9, #6
pCMP 8, 11
zCMPLS 7, 3
�BHI gp4%
� STR 7, [12, 8, LSL #2]
�ADD 8, 8, #1
�CMP 8, 11
� LDRLS 14, [5, 8, LSL #2]
�MULLS 7, 14, 14
� .gp4%
�MOV 2, #2
�CMP 2, 8
�BHS gpprime%
� .gp6%
�!LDR 14, [12, 2, LSL #2]
� LDR 4, [5, 2, LSL #2]
.gp7%
CMP 14, 3
ADDLO 14, 14, 4, LSL #1
$BLO gp7%
.!STR 14, [12, 2, LSL #2]
8MVNEQ 2, #1
BADD 2, 2, #1
L .gp5%
VCMP 2, 8
`DBLO gp6% ;�prime on loop exit iff R2=&ffffffff iff NE
jBNE gpmainloop%
t
.gpprime%
~CMP 3, 6
�BLHS outputprime%
�ADD 1, 1, #1
�CMP 3, 10
�ADDLS 11, 11, #1
� STRLS 3, [5, 11, LSL #2]
�B gpmainloop%
�:
�".buf% EQUD 0:EQUD 0:EQUD 0
� .bufend%
�@.outputprime% ;you might want to alter this routine to
�F ;dump prime into an array instead of print it!
�STMFD 13!, {0,1,2}
MOV 0, 3
ADR 1, buf%
MOV 2, #bufend%-buf%
#SWI "OS_ConvertCardinal4"
(SWI "OS_Write0"
2SWI "OS_NewLine"
<LDMFD 13!, {0,1,2}
FMOV PC, 14
P]
Z�
d�
n:
x� �primes(nlow%, n%)
�
� anp%, v
� A%=n%
� B%=nlow%
� � n%<0 v=n%+(2^32) � v=n%
�
C%=� �(v)
�anp%=� genprimes%
�?�'"Found ";anp%;" prime numbers less than or equal to ";n%'
��
� 00000000 0d 00 0a 0d 6d 61 78 6e 3d 32 5e 33 32 0d 00 14 |....maxn=2^32...|
00000010 4e 6e 70 25 3d 32 33 2b 31 2e 31 31 35 2a 32 2a |Nnp%=23+1.115*2*|
00000020 b6 28 6d 61 78 6e 29 2f aa 28 6d 61 78 6e 29 20 |.(maxn)/.(maxn) |
00000030 3a f4 20 61 70 70 72 6f 78 20 75 70 70 65 72 20 |:. approx upper |
00000040 62 6f 75 6e 64 20 6f 6e 20 23 70 72 69 6d 65 73 |bound on #primes|
00000050 3c 3d 53 51 52 28 6d 61 78 6e 29 0d 00 1e 1f de |<=SQR(maxn).....|
00000060 20 6d 75 6c 74 69 70 6c 65 25 20 34 2a 6e 70 25 | multiple% 4*np%|
00000070 2c 20 70 25 20 34 2a 6e 70 25 0d 00 28 05 3a 0d |, p% 4*np%..(.:.|
00000080 00 32 08 f2 61 73 73 0d 00 3c 05 3a 0d 00 46 11 |.2..ass..<.:..F.|
00000090 f4 20 50 52 4f 43 70 72 69 6d 65 73 3a 0d 00 50 |. PROCprimes:..P|
000000a0 42 f4 20 43 6f 75 6e 74 20 61 6c 6c 20 70 72 69 |B. Count all pri|
000000b0 6d 65 73 20 6c 65 73 73 20 74 68 61 6e 20 6f 72 |mes less than or|
000000c0 20 65 71 75 61 6c 20 74 6f 20 73 65 63 6f 6e 64 | equal to second|
000000d0 20 70 61 72 61 6d 65 74 65 72 2c 20 61 6e 64 0d | parameter, and.|
000000e0 00 5a 43 f4 20 64 69 73 70 6c 61 79 20 61 6c 6c |.ZC. display all|
000000f0 20 74 68 6f 73 65 20 77 68 69 63 68 20 61 72 65 | those which are|
00000100 20 61 6c 73 6f 20 67 72 65 61 74 65 72 20 74 68 | also greater th|
00000110 61 6e 20 6f 72 20 65 71 75 61 6c 20 74 6f 20 74 |an or equal to t|
00000120 68 65 0d 00 64 16 f4 20 66 69 72 73 74 20 70 61 |he..d.. first pa|
00000130 72 61 6d 65 74 65 72 2e 0d 00 6e 05 3a 0d 00 78 |rameter...n.:..x|
00000140 3f f4 20 54 6f 20 73 65 65 20 68 6f 77 20 71 75 |?. To see how qu|
00000150 69 63 6b 20 74 68 69 73 20 63 6f 64 65 20 69 73 |ick this code is|
00000160 20 61 74 20 66 69 6e 64 69 6e 67 20 70 72 69 6d | at finding prim|
00000170 65 73 2c 20 73 65 74 20 62 6f 74 68 0d 00 82 42 |es, set both...B|
00000180 f4 20 70 61 72 61 6d 65 74 65 72 73 20 65 71 75 |. parameters equ|
00000190 61 6c 2c 20 73 6f 20 74 68 61 74 20 6e 6f 20 74 |al, so that no t|
000001a0 69 6d 65 20 69 73 20 77 61 73 74 65 64 20 64 69 |ime is wasted di|
000001b0 73 70 6c 61 79 69 6e 67 20 74 68 65 6d 2e 0d 00 |splaying them...|
000001c0 8c 05 3a 0d 00 96 07 d1 3d 30 0d 00 a0 1d f2 70 |..:.....=0.....p|
000001d0 72 69 6d 65 73 28 36 35 35 33 36 2a 32 2c 20 36 |rimes(65536*2, 6|
000001e0 35 35 33 36 2a 34 29 0d 00 aa 14 f1 20 22 54 69 |5536*4)..... "Ti|
000001f0 6d 65 3a 20 22 3b 91 2f 31 30 30 0d 00 b4 05 3a |me: ";./100....:|
00000200 0d 00 be 05 e0 0d 00 c8 05 3a 0d 00 d2 05 3a 0d |.........:....:.|
00000210 00 dc 0a dd 20 f2 61 73 73 0d 00 e6 10 de 20 63 |.... .ass..... c|
00000220 6f 64 65 25 20 31 30 32 34 0d 00 f0 15 e3 20 70 |ode% 1024..... p|
00000230 61 73 73 25 3d 30 20 b8 20 32 20 88 20 32 0d 00 |ass%=0 . 2 . 2..|
00000240 fa 0c 50 25 3d 63 6f 64 65 25 0d 01 04 0e 5b 4f |..P%=code%....[O|
00000250 50 54 20 70 61 73 73 25 0d 01 0e 26 2e 61 6d 75 |PT pass%...&.amu|
00000260 6c 74 69 70 6c 65 25 20 20 20 20 20 20 20 20 20 |ltiple% |
00000270 45 51 55 44 20 6d 75 6c 74 69 70 6c 65 25 0d 01 |EQUD multiple%..|
00000280 18 15 2e 61 70 25 20 20 20 20 20 20 45 51 55 44 |...ap% EQUD|
00000290 20 70 25 0d 01 22 0f 2e 67 65 6e 70 72 69 6d 65 | p%.."..genprime|
000002a0 73 25 0d 01 2c 12 4d 4f 56 20 20 20 20 20 20 20 |s%..,.MOV |
000002b0 36 2c 20 31 0d 01 36 13 4d 4f 56 20 20 20 20 20 |6, 1..6.MOV |
000002c0 20 20 31 30 2c 20 32 0d 01 40 3c 4d 4f 56 20 20 | 10, 2..@<MOV |
000002d0 20 20 20 20 20 31 2c 20 23 30 20 20 20 20 20 3b | 1, #0 ;|
000002e0 61 6e 70 20 28 61 63 74 75 61 6c 20 6e 75 6d 62 |anp (actual numb|
000002f0 65 72 20 6f 66 20 70 72 69 6d 65 73 20 66 6f 75 |er of primes fou|
00000300 6e 64 29 0d 01 4a 13 43 4d 50 20 20 20 20 20 20 |nd)..J.CMP |
00000310 20 30 2c 20 23 31 0d 01 54 12 4d 4f 56 4c 53 20 | 0, #1..T.MOVLS |
00000320 20 20 20 20 30 2c 20 31 0d 01 5e 14 4d 4f 56 4c | 0, 1..^.MOVL|
00000330 53 53 20 20 20 20 50 43 2c 20 31 34 0d 01 68 17 |SS PC, 14..h.|
00000340 53 54 4d 46 44 20 20 20 20 20 31 33 21 2c 20 7b |STMFD 13!, {|
00000350 31 34 7d 0d 01 72 13 4d 4f 56 20 20 20 20 20 20 |14}..r.MOV |
00000360 20 32 2c 20 23 32 0d 01 7c 13 4d 4f 56 20 20 20 | 2, #2..|.MOV |
00000370 20 20 20 20 33 2c 20 23 33 0d 01 86 13 4d 4f 56 | 3, #3....MOV|
00000380 20 20 20 20 20 20 20 34 2c 20 23 35 0d 01 90 14 | 4, #5....|
00000390 4c 44 52 20 20 20 20 20 20 20 35 2c 20 61 70 25 |LDR 5, ap%|
000003a0 0d 01 9a 18 53 54 4d 49 41 20 20 20 20 20 35 2c |....STMIA 5,|
000003b0 20 7b 32 2c 33 2c 34 7d 0d 01 a4 14 4d 4f 56 20 | {2,3,4}....MOV |
000003c0 20 20 20 20 20 20 31 31 2c 20 23 32 0d 01 ae 13 | 11, #2....|
000003d0 43 4d 50 20 20 20 20 20 20 20 30 2c 20 23 35 0d |CMP 0, #5.|
000003e0 01 b8 12 42 48 53 20 20 20 20 20 20 20 67 70 31 |...BHS gp1|
000003f0 25 0d 01 c2 13 4d 4f 56 20 20 20 20 20 20 20 32 |%....MOV 2|
00000400 2c 20 23 30 0d 01 cc 16 53 55 42 20 20 20 20 20 |, #0....SUB |
00000410 20 20 37 2c 20 30 2c 20 23 31 0d 01 d6 1a 4d 4f | 7, 0, #1....MO|
00000420 56 20 20 20 20 20 20 20 37 2c 20 37 2c 20 4c 53 |V 7, 7, LS|
00000430 52 20 23 31 0d 01 e0 09 2e 67 70 32 25 0d 01 ea |R #1.....gp2%...|
00000440 1f 4c 44 52 20 20 20 20 20 20 20 33 2c 20 5b 35 |.LDR 3, [5|
00000450 2c 20 32 2c 20 4c 53 4c 20 23 32 5d 0d 01 f4 12 |, 2, LSL #2]....|
00000460 43 4d 50 20 20 20 20 20 20 20 33 2c 20 36 0d 01 |CMP 3, 6..|
00000470 fe 1a 42 4c 48 53 20 20 20 20 20 20 6f 75 74 70 |..BLHS outp|
00000480 75 74 70 72 69 6d 65 25 0d 02 08 16 41 44 44 20 |utprime%....ADD |
00000490 20 20 20 20 20 20 31 2c 20 31 2c 20 23 31 0d 02 | 1, 1, #1..|
000004a0 12 16 41 44 44 20 20 20 20 20 20 20 32 2c 20 32 |..ADD 2, 2|
000004b0 2c 20 23 31 0d 02 1c 12 43 4d 50 20 20 20 20 20 |, #1....CMP |
000004c0 20 20 32 2c 20 37 0d 02 26 12 42 4c 53 20 20 20 | 2, 7..&.BLS |
000004d0 20 20 20 20 67 70 32 25 0d 02 30 12 4d 4f 56 20 | gp2%..0.MOV |
000004e0 20 20 20 20 20 20 30 2c 20 31 0d 02 3a 18 4c 44 | 0, 1..:.LD|
000004f0 4d 46 44 20 20 20 20 20 31 33 21 2c 20 7b 50 43 |MFD 13!, {PC|
00000500 7d 5e 0d 02 44 09 2e 67 70 31 25 0d 02 4e 13 4d |}^..D..gp1%..N.M|
00000510 4f 56 20 20 20 20 20 20 20 32 2c 20 23 30 0d 02 |OV 2, #0..|
00000520 58 13 4d 4f 56 20 20 20 20 20 20 20 37 2c 20 23 |X.MOV 7, #|
00000530 32 0d 02 62 09 2e 67 70 33 25 0d 02 6c 1f 4c 44 |2..b..gp3%..l.LD|
00000540 52 20 20 20 20 20 20 20 33 2c 20 5b 35 2c 20 32 |R 3, [5, 2|
00000550 2c 20 4c 53 4c 20 23 32 5d 0d 02 76 12 43 4d 50 |, LSL #2]..v.CMP|
00000560 20 20 20 20 20 20 20 33 2c 20 36 0d 02 80 1a 42 | 3, 6....B|
00000570 4c 48 53 20 20 20 20 20 20 6f 75 74 70 75 74 70 |LHS outputp|
00000580 72 69 6d 65 25 0d 02 8a 16 41 44 44 20 20 20 20 |rime%....ADD |
00000590 20 20 20 31 2c 20 31 2c 20 23 31 0d 02 94 16 41 | 1, 1, #1....A|
000005a0 44 44 20 20 20 20 20 20 20 32 2c 20 32 2c 20 23 |DD 2, 2, #|
000005b0 31 0d 02 9e 12 43 4d 50 20 20 20 20 20 20 20 32 |1....CMP 2|
000005c0 2c 20 37 0d 02 a8 12 42 4c 53 20 20 20 20 20 20 |, 7....BLS |
000005d0 20 67 70 33 25 0d 02 b2 05 3a 0d 02 bc 13 4d 4f | gp3%....:....MO|
000005e0 56 20 20 20 20 20 20 20 33 2c 20 23 35 0d 02 c6 |V 3, #5...|
000005f0 14 4d 4f 56 20 20 20 20 20 20 20 37 2c 20 23 32 |.MOV 7, #2|
00000600 35 0d 02 d0 13 4d 4f 56 20 20 20 20 20 20 20 38 |5....MOV 8|
00000610 2c 20 23 32 0d 02 da 13 4d 4f 56 20 20 20 20 20 |, #2....MOV |
00000620 20 20 39 2c 20 23 32 0d 02 e4 1c 4c 44 52 20 20 | 9, #2....LDR |
00000630 20 20 20 20 20 31 32 2c 20 61 6d 75 6c 74 69 70 | 12, amultip|
00000640 6c 65 25 0d 02 ee 05 3a 0d 02 f8 10 2e 67 70 6d |le%....:.....gpm|
00000650 61 69 6e 6c 6f 6f 70 25 0d 03 02 17 80 20 20 20 |ainloop%..... |
00000660 20 20 20 20 31 34 2c 20 33 2c 20 23 32 35 35 0d | 14, 3, #255.|
00000670 03 0c 14 43 4d 50 20 20 20 20 20 20 20 31 34 2c |...CMP 14,|
00000680 20 23 32 0d 03 16 12 42 48 53 20 20 20 20 20 20 | #2....BHS |
00000690 20 67 70 38 25 0d 03 20 22 53 57 49 20 20 20 20 | gp8%.. "SWI |
000006a0 20 20 20 22 4f 53 5f 52 65 61 64 45 73 63 61 70 | "OS_ReadEscap|
000006b0 65 53 74 61 74 65 22 0d 03 2a 18 4c 44 4d 43 53 |eState"..*.LDMCS|
000006c0 46 44 20 20 20 31 33 21 2c 20 7b 50 43 7d 5e 0d |FD 13!, {PC}^.|
000006d0 03 34 09 2e 67 70 38 25 0d 03 3e 15 41 44 44 20 |.4..gp8%..>.ADD |
000006e0 20 20 20 20 20 20 33 2c 20 33 2c 20 39 0d 03 48 | 3, 3, 9..H|
000006f0 12 43 4d 50 20 20 20 20 20 20 20 33 2c 20 30 0d |.CMP 3, 0.|
00000700 03 52 12 4d 4f 56 48 49 20 20 20 20 20 30 2c 20 |.R.MOVHI 0, |
00000710 31 0d 03 5c 18 4c 44 4d 48 49 46 44 20 20 20 31 |1..\.LDMHIFD 1|
00000720 33 21 2c 20 7b 50 43 7d 5e 0d 03 66 16 52 53 42 |3!, {PC}^..f.RSB|
00000730 20 20 20 20 20 20 20 39 2c 20 39 2c 20 23 36 0d | 9, 9, #6.|
00000740 03 70 13 43 4d 50 20 20 20 20 20 20 20 38 2c 20 |.p.CMP 8, |
00000750 31 31 0d 03 7a 12 43 4d 50 4c 53 20 20 20 20 20 |11..z.CMPLS |
00000760 37 2c 20 33 0d 03 84 12 42 48 49 20 20 20 20 20 |7, 3....BHI |
00000770 20 20 67 70 34 25 0d 03 8e 20 53 54 52 20 20 20 | gp4%... STR |
00000780 20 20 20 20 37 2c 20 5b 31 32 2c 20 38 2c 20 4c | 7, [12, 8, L|
00000790 53 4c 20 23 32 5d 0d 03 98 16 41 44 44 20 20 20 |SL #2]....ADD |
000007a0 20 20 20 20 38 2c 20 38 2c 20 23 31 0d 03 a2 13 | 8, 8, #1....|
000007b0 43 4d 50 20 20 20 20 20 20 20 38 2c 20 31 31 0d |CMP 8, 11.|
000007c0 03 ac 20 4c 44 52 4c 53 20 20 20 20 20 31 34 2c |.. LDRLS 14,|
000007d0 20 5b 35 2c 20 38 2c 20 4c 53 4c 20 23 32 5d 0d | [5, 8, LSL #2].|
000007e0 03 b6 17 4d 55 4c 4c 53 20 20 20 20 20 37 2c 20 |...MULLS 7, |
000007f0 31 34 2c 20 31 34 0d 03 c0 09 2e 67 70 34 25 0d |14, 14.....gp4%.|
00000800 03 ca 13 4d 4f 56 20 20 20 20 20 20 20 32 2c 20 |...MOV 2, |
00000810 23 32 0d 03 d4 12 43 4d 50 20 20 20 20 20 20 20 |#2....CMP |
00000820 32 2c 20 38 0d 03 de 16 42 48 53 20 20 20 20 20 |2, 8....BHS |
00000830 20 20 67 70 70 72 69 6d 65 25 0d 03 e8 09 2e 67 | gpprime%.....g|
00000840 70 36 25 0d 03 f2 21 4c 44 52 20 20 20 20 20 20 |p6%...!LDR |
00000850 20 31 34 2c 20 5b 31 32 2c 20 32 2c 20 4c 53 4c | 14, [12, 2, LSL|
00000860 20 23 32 5d 0d 03 fc 1f 4c 44 52 20 20 20 20 20 | #2]....LDR |
00000870 20 20 34 2c 20 5b 35 2c 20 32 2c 20 4c 53 4c 20 | 4, [5, 2, LSL |
00000880 23 32 5d 0d 04 06 09 2e 67 70 37 25 0d 04 10 13 |#2].....gp7%....|
00000890 43 4d 50 20 20 20 20 20 20 20 31 34 2c 20 33 0d |CMP 14, 3.|
000008a0 04 1a 1f 41 44 44 4c 4f 20 20 20 20 20 31 34 2c |...ADDLO 14,|
000008b0 20 31 34 2c 20 34 2c 20 4c 53 4c 20 23 31 0d 04 | 14, 4, LSL #1..|
000008c0 24 12 42 4c 4f 20 20 20 20 20 20 20 67 70 37 25 |$.BLO gp7%|
000008d0 0d 04 2e 21 53 54 52 20 20 20 20 20 20 20 31 34 |...!STR 14|
000008e0 2c 20 5b 31 32 2c 20 32 2c 20 4c 53 4c 20 23 32 |, [12, 2, LSL #2|
000008f0 5d 0d 04 38 13 4d 56 4e 45 51 20 20 20 20 20 32 |]..8.MVNEQ 2|
00000900 2c 20 23 31 0d 04 42 16 41 44 44 20 20 20 20 20 |, #1..B.ADD |
00000910 20 20 32 2c 20 32 2c 20 23 31 0d 04 4c 09 2e 67 | 2, 2, #1..L..g|
00000920 70 35 25 0d 04 56 12 43 4d 50 20 20 20 20 20 20 |p5%..V.CMP |
00000930 20 32 2c 20 38 0d 04 60 44 42 4c 4f 20 20 20 20 | 2, 8..`DBLO |
00000940 20 20 20 67 70 36 25 20 20 20 20 20 20 3b ac 70 | gp6% ;.p|
00000950 72 69 6d 65 20 6f 6e 20 6c 6f 6f 70 20 65 78 69 |rime on loop exi|
00000960 74 20 69 66 66 20 52 32 3d 26 66 66 66 66 66 66 |t iff R2=&ffffff|
00000970 66 66 20 69 66 66 20 4e 45 0d 04 6a 19 42 4e 45 |ff iff NE..j.BNE|
00000980 20 20 20 20 20 20 20 67 70 6d 61 69 6e 6c 6f 6f | gpmainloo|
00000990 70 25 0d 04 74 0d 2e 67 70 70 72 69 6d 65 25 0d |p%..t..gpprime%.|
000009a0 04 7e 12 43 4d 50 20 20 20 20 20 20 20 33 2c 20 |.~.CMP 3, |
000009b0 36 0d 04 88 1a 42 4c 48 53 20 20 20 20 20 20 6f |6....BLHS o|
000009c0 75 74 70 75 74 70 72 69 6d 65 25 0d 04 92 16 41 |utputprime%....A|
000009d0 44 44 20 20 20 20 20 20 20 31 2c 20 31 2c 20 23 |DD 1, 1, #|
000009e0 31 0d 04 9c 13 43 4d 50 20 20 20 20 20 20 20 33 |1....CMP 3|
000009f0 2c 20 31 30 0d 04 a6 18 41 44 44 4c 53 20 20 20 |, 10....ADDLS |
00000a00 20 20 31 31 2c 20 31 31 2c 20 23 31 0d 04 b0 20 | 11, 11, #1... |
00000a10 53 54 52 4c 53 20 20 20 20 20 33 2c 20 5b 35 2c |STRLS 3, [5,|
00000a20 20 31 31 2c 20 4c 53 4c 20 23 32 5d 0d 04 ba 19 | 11, LSL #2]....|
00000a30 42 20 20 20 20 20 20 20 20 20 67 70 6d 61 69 6e |B gpmain|
00000a40 6c 6f 6f 70 25 0d 04 c4 05 3a 0d 04 ce 22 2e 62 |loop%....:...".b|
00000a50 75 66 25 20 20 20 20 20 45 51 55 44 20 30 3a 45 |uf% EQUD 0:E|
00000a60 51 55 44 20 30 3a 45 51 55 44 20 30 0d 04 d8 0c |QUD 0:EQUD 0....|
00000a70 2e 62 75 66 65 6e 64 25 0d 04 e2 40 2e 6f 75 74 |.bufend%...@.out|
00000a80 70 75 74 70 72 69 6d 65 25 20 20 20 20 20 20 20 |putprime% |
00000a90 3b 79 6f 75 20 6d 69 67 68 74 20 77 61 6e 74 20 |;you might want |
00000aa0 74 6f 20 61 6c 74 65 72 20 74 68 69 73 20 72 6f |to alter this ro|
00000ab0 75 74 69 6e 65 20 74 6f 0d 04 ec 46 20 20 20 20 |utine to...F |
00000ac0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | |
00000ad0 3b 64 75 6d 70 20 70 72 69 6d 65 20 69 6e 74 6f |;dump prime into|
00000ae0 20 61 6e 20 61 72 72 61 79 20 69 6e 73 74 65 61 | an array instea|
00000af0 64 20 6f 66 20 70 72 69 6e 74 20 69 74 21 0d 04 |d of print it!..|
00000b00 f6 1a 53 54 4d 46 44 20 20 20 20 20 31 33 21 2c |..STMFD 13!,|
00000b10 20 7b 30 2c 31 2c 32 7d 0d 05 00 12 4d 4f 56 20 | {0,1,2}....MOV |
00000b20 20 20 20 20 20 20 30 2c 20 33 0d 05 0a 15 41 44 | 0, 3....AD|
00000b30 52 20 20 20 20 20 20 20 31 2c 20 62 75 66 25 0d |R 1, buf%.|
00000b40 05 14 1e 4d 4f 56 20 20 20 20 20 20 20 32 2c 20 |...MOV 2, |
00000b50 23 62 75 66 65 6e 64 25 2d 62 75 66 25 0d 05 1e |#bufend%-buf%...|
00000b60 23 53 57 49 20 20 20 20 20 20 20 22 4f 53 5f 43 |#SWI "OS_C|
00000b70 6f 6e 76 65 72 74 43 61 72 64 69 6e 61 6c 34 22 |onvertCardinal4"|
00000b80 0d 05 28 19 53 57 49 20 20 20 20 20 20 20 22 4f |..(.SWI "O|
00000b90 53 5f 57 72 69 74 65 30 22 0d 05 32 1a 53 57 49 |S_Write0"..2.SWI|
00000ba0 20 20 20 20 20 20 20 22 4f 53 5f 4e 65 77 4c 69 | "OS_NewLi|
00000bb0 6e 65 22 0d 05 3c 1a 4c 44 4d 46 44 20 20 20 20 |ne"..<.LDMFD |
00000bc0 20 31 33 21 2c 20 7b 30 2c 31 2c 32 7d 0d 05 46 | 13!, {0,1,2}..F|
00000bd0 14 4d 4f 56 20 20 20 20 20 20 20 50 43 2c 20 31 |.MOV PC, 1|
00000be0 34 0d 05 50 05 5d 0d 05 5a 05 ed 0d 05 64 05 e1 |4..P.]..Z....d..|
00000bf0 0d 05 6e 05 3a 0d 05 78 18 dd 20 f2 70 72 69 6d |..n.:..x.. .prim|
00000c00 65 73 28 6e 6c 6f 77 25 2c 20 6e 25 29 0d 05 82 |es(nlow%, n%)...|
00000c10 0d ea 20 61 6e 70 25 2c 20 76 0d 05 8c 09 41 25 |.. anp%, v....A%|
00000c20 3d 6e 25 0d 05 96 0c 42 25 3d 6e 6c 6f 77 25 0d |=n%....B%=nlow%.|
00000c30 05 a0 1d e7 20 6e 25 3c 30 20 76 3d 6e 25 2b 28 |.... n%<0 v=n%+(|
00000c40 32 5e 33 32 29 20 8b 20 76 3d 6e 25 0d 05 aa 0d |2^32) . v=n%....|
00000c50 43 25 3d a8 20 b6 28 76 29 0d 05 b4 15 61 6e 70 |C%=. .(v)....anp|
00000c60 25 3d ba 20 67 65 6e 70 72 69 6d 65 73 25 0d 05 |%=. genprimes%..|
00000c70 be 3f f1 27 22 46 6f 75 6e 64 20 22 3b 61 6e 70 |.?.'"Found ";anp|
00000c80 25 3b 22 20 70 72 69 6d 65 20 6e 75 6d 62 65 72 |%;" prime number|
00000c90 73 20 6c 65 73 73 20 74 68 61 6e 20 6f 72 20 65 |s less than or e|
00000ca0 71 75 61 6c 20 74 6f 20 22 3b 6e 25 27 0d 05 c8 |qual to ";n%'...|
00000cb0 05 e1 0d ff |....|
00000cb4 
.