Wed, 31 Dec 2014 06:09:35 +0100
Cloned upstream origin tor-browser at tor-browser-31.3.0esr-4.5-1-build1
revision ID fc1c9ff7c1b2defdbc039f12214767608f46423f for hacking purpose.
michael@0 | 1 | /* |
michael@0 | 2 | * Simple test driver for MPI library |
michael@0 | 3 | * |
michael@0 | 4 | * Test 8: Probabilistic primality tester |
michael@0 | 5 | * |
michael@0 | 6 | * This Source Code Form is subject to the terms of the Mozilla Public |
michael@0 | 7 | * License, v. 2.0. If a copy of the MPL was not distributed with this |
michael@0 | 8 | * file, You can obtain one at http://mozilla.org/MPL/2.0/. */ |
michael@0 | 9 | |
michael@0 | 10 | #include <stdio.h> |
michael@0 | 11 | #include <stdlib.h> |
michael@0 | 12 | #include <string.h> |
michael@0 | 13 | #include <ctype.h> |
michael@0 | 14 | #include <limits.h> |
michael@0 | 15 | #include <time.h> |
michael@0 | 16 | |
michael@0 | 17 | #define MP_IOFUNC 1 |
michael@0 | 18 | #include "mpi.h" |
michael@0 | 19 | |
michael@0 | 20 | #include "mpprime.h" |
michael@0 | 21 | |
michael@0 | 22 | int main(int argc, char *argv[]) |
michael@0 | 23 | { |
michael@0 | 24 | int ix; |
michael@0 | 25 | mp_digit num; |
michael@0 | 26 | mp_int a; |
michael@0 | 27 | |
michael@0 | 28 | srand(time(NULL)); |
michael@0 | 29 | |
michael@0 | 30 | if(argc < 2) { |
michael@0 | 31 | fprintf(stderr, "Usage: %s <a>\n", argv[0]); |
michael@0 | 32 | return 1; |
michael@0 | 33 | } |
michael@0 | 34 | |
michael@0 | 35 | printf("Test 8: Probabilistic primality testing\n\n"); |
michael@0 | 36 | |
michael@0 | 37 | mp_init(&a); |
michael@0 | 38 | |
michael@0 | 39 | mp_read_radix(&a, argv[1], 10); |
michael@0 | 40 | |
michael@0 | 41 | printf("a = "); mp_print(&a, stdout); fputc('\n', stdout); |
michael@0 | 42 | |
michael@0 | 43 | printf("\nChecking for divisibility by small primes ... \n"); |
michael@0 | 44 | num = 170; |
michael@0 | 45 | if(mpp_divis_primes(&a, &num) == MP_YES) { |
michael@0 | 46 | printf("it is not prime\n"); |
michael@0 | 47 | goto CLEANUP; |
michael@0 | 48 | } |
michael@0 | 49 | printf("Passed that test (not divisible by any small primes).\n"); |
michael@0 | 50 | |
michael@0 | 51 | for(ix = 0; ix < 10; ix++) { |
michael@0 | 52 | printf("\nPerforming Rabin-Miller test, iteration %d\n", ix + 1); |
michael@0 | 53 | |
michael@0 | 54 | if(mpp_pprime(&a, 5) == MP_NO) { |
michael@0 | 55 | printf("it is not prime\n"); |
michael@0 | 56 | goto CLEANUP; |
michael@0 | 57 | } |
michael@0 | 58 | } |
michael@0 | 59 | printf("All tests passed; a is probably prime\n"); |
michael@0 | 60 | |
michael@0 | 61 | CLEANUP: |
michael@0 | 62 | mp_clear(&a); |
michael@0 | 63 | |
michael@0 | 64 | return 0; |
michael@0 | 65 | } |