fix miller-rabin

Doug Coleman 2009-05-06 14:10:29 -05:00
parent 15bc3625e8
commit c0156d462e
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