Thu, 22 Jan 2015 13:21:57 +0100
Incorporate requested changes from Mozilla in review:
https://bugzilla.mozilla.org/show_bug.cgi?id=1123480#c6
michael@0 | 1 | /* |
michael@0 | 2 | * Simple test driver for MPI library |
michael@0 | 3 | * |
michael@0 | 4 | * Test 7: Random and divisibility tests |
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 | mp_digit num; |
michael@0 | 25 | mp_int a, b; |
michael@0 | 26 | |
michael@0 | 27 | srand(time(NULL)); |
michael@0 | 28 | |
michael@0 | 29 | if(argc < 3) { |
michael@0 | 30 | fprintf(stderr, "Usage: %s <a> <b>\n", argv[0]); |
michael@0 | 31 | return 1; |
michael@0 | 32 | } |
michael@0 | 33 | |
michael@0 | 34 | printf("Test 7: Random & divisibility tests\n\n"); |
michael@0 | 35 | |
michael@0 | 36 | mp_init(&a); |
michael@0 | 37 | mp_init(&b); |
michael@0 | 38 | |
michael@0 | 39 | mp_read_radix(&a, argv[1], 10); |
michael@0 | 40 | mp_read_radix(&b, argv[2], 10); |
michael@0 | 41 | |
michael@0 | 42 | printf("a = "); mp_print(&a, stdout); fputc('\n', stdout); |
michael@0 | 43 | printf("b = "); mp_print(&b, stdout); fputc('\n', stdout); |
michael@0 | 44 | |
michael@0 | 45 | if(mpp_divis(&a, &b) == MP_YES) |
michael@0 | 46 | printf("a is divisible by b\n"); |
michael@0 | 47 | else |
michael@0 | 48 | printf("a is not divisible by b\n"); |
michael@0 | 49 | |
michael@0 | 50 | if(mpp_divis(&b, &a) == MP_YES) |
michael@0 | 51 | printf("b is divisible by a\n"); |
michael@0 | 52 | else |
michael@0 | 53 | printf("b is not divisible by a\n"); |
michael@0 | 54 | |
michael@0 | 55 | printf("\nb = mpp_random()\n"); |
michael@0 | 56 | mpp_random(&b); |
michael@0 | 57 | printf("b = "); mp_print(&b, stdout); fputc('\n', stdout); |
michael@0 | 58 | mpp_random(&b); |
michael@0 | 59 | printf("b = "); mp_print(&b, stdout); fputc('\n', stdout); |
michael@0 | 60 | mpp_random(&b); |
michael@0 | 61 | printf("b = "); mp_print(&b, stdout); fputc('\n', stdout); |
michael@0 | 62 | |
michael@0 | 63 | printf("\nTesting a for divisibility by first 170 primes\n"); |
michael@0 | 64 | num = 170; |
michael@0 | 65 | if(mpp_divis_primes(&a, &num) == MP_YES) |
michael@0 | 66 | printf("It is divisible by at least one of them\n"); |
michael@0 | 67 | else |
michael@0 | 68 | printf("It is not divisible by any of them\n"); |
michael@0 | 69 | |
michael@0 | 70 | mp_clear(&b); |
michael@0 | 71 | mp_clear(&a); |
michael@0 | 72 | |
michael@0 | 73 | return 0; |
michael@0 | 74 | } |