Sieve of Eratosthenes: Difference between revisions

Content deleted Content added

Inline

@@ Line 15,187: / Line 15,187: @@
     number on the primes ancillary stack is set.
     Initially all the bits are set except for 0 and 1,
-    which are not prime numbers by definition.            )
+    which are not prime numbers by definition.
+    "eratosthenes" unsets all bits above those specified
+    by it's argument. )
   [ bit ~
@@ Line 15,194: / Line 15,196: @@
   [ bit primes share & 0 != ]   is isprime      ( n --> b )
-  [ dup sqrt times
+  [ dup dup sqrt times
       [ i^ 1+
         dup isprime if
           [ dup 2 **
             [ dup -composite
               over +
               rot 2dup >
               dip unrot until ]
             drop ]
         drop ]
+     drop
-     drop ]                     is eratosthenes ( n -->   )
++ bit 1 -
+     primes take &
+     primes put ]                 is eratosthenes ( n -->   )
 eratosthenes