security/nss/lib/freebl/mpi/tests/mptest-8.c

Wed, 31 Dec 2014 06:09:35 +0100

author
Michael Schloh von Bennewitz <michael@schloh.com>
date
Wed, 31 Dec 2014 06:09:35 +0100
changeset 0
6474c204b198
permissions
-rw-r--r--

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 }

mercurial