fix miller-rabin

db4
Doug Coleman 2009-05-06 14:10:29 -05:00
parent 4cb6ecbc3c
commit 8b4815c01e
1 changed files with 3 additions and 0 deletions

View File

@ -8,6 +8,8 @@ IN: math.miller-rabin
: >odd ( n -- int ) dup even? [ 1 + ] when ; foldable
: >even ( n -- int ) 0 clear-bit ; foldable
TUPLE: positive-even-expected n ;
:: (miller-rabin) ( n trials -- ? )
@ -97,6 +99,7 @@ PRIVATE>
} 1&& ;
: next-safe-prime ( n -- q )
1 - >even 2 /
next-safe-prime-candidate
dup >safe-prime-form
dup miller-rabin