summaryrefslogtreecommitdiff
path: root/lib/freebl/mpi/tests/mptest-4.c
diff options
context:
space:
mode:
Diffstat (limited to 'lib/freebl/mpi/tests/mptest-4.c')
-rw-r--r--lib/freebl/mpi/tests/mptest-4.c124
1 files changed, 71 insertions, 53 deletions
diff --git a/lib/freebl/mpi/tests/mptest-4.c b/lib/freebl/mpi/tests/mptest-4.c
index 300173977..0f326ac2c 100644
--- a/lib/freebl/mpi/tests/mptest-4.c
+++ b/lib/freebl/mpi/tests/mptest-4.c
@@ -15,79 +15,97 @@
#include "mpi.h"
-int main(int argc, char *argv[])
+int
+main(int argc, char *argv[])
{
- int ix;
- mp_int a, b, c, m;
- mp_digit r;
+ int ix;
+ mp_int a, b, c, m;
+ mp_digit r;
- if(argc < 4) {
- fprintf(stderr, "Usage: %s <a> <b> <m>\n", argv[0]);
- return 1;
- }
+ if (argc < 4) {
+ fprintf(stderr, "Usage: %s <a> <b> <m>\n", argv[0]);
+ return 1;
+ }
- printf("Test 4: Modular arithmetic\n\n");
+ printf("Test 4: Modular arithmetic\n\n");
- mp_init(&a);
- mp_init(&b);
- mp_init(&m);
+ mp_init(&a);
+ mp_init(&b);
+ mp_init(&m);
- mp_read_radix(&a, argv[1], 10);
- mp_read_radix(&b, argv[2], 10);
- mp_read_radix(&m, argv[3], 10);
- printf("a = "); mp_print(&a, stdout); fputc('\n', stdout);
- printf("b = "); mp_print(&b, stdout); fputc('\n', stdout);
- printf("m = "); mp_print(&m, stdout); fputc('\n', stdout);
-
- mp_init(&c);
- printf("\nc = a (mod m)\n");
+ mp_read_radix(&a, argv[1], 10);
+ mp_read_radix(&b, argv[2], 10);
+ mp_read_radix(&m, argv[3], 10);
+ printf("a = ");
+ mp_print(&a, stdout);
+ fputc('\n', stdout);
+ printf("b = ");
+ mp_print(&b, stdout);
+ fputc('\n', stdout);
+ printf("m = ");
+ mp_print(&m, stdout);
+ fputc('\n', stdout);
- mp_mod(&a, &m, &c);
- printf("c = "); mp_print(&c, stdout); fputc('\n', stdout);
+ mp_init(&c);
+ printf("\nc = a (mod m)\n");
- printf("\nc = b (mod m)\n");
+ mp_mod(&a, &m, &c);
+ printf("c = ");
+ mp_print(&c, stdout);
+ fputc('\n', stdout);
- mp_mod(&b, &m, &c);
- printf("c = "); mp_print(&c, stdout); fputc('\n', stdout);
+ printf("\nc = b (mod m)\n");
- printf("\nc = b (mod 1853)\n");
+ mp_mod(&b, &m, &c);
+ printf("c = ");
+ mp_print(&c, stdout);
+ fputc('\n', stdout);
- mp_mod_d(&b, 1853, &r);
- printf("c = %04X\n", r);
+ printf("\nc = b (mod 1853)\n");
- printf("\nc = (a + b) mod m\n");
+ mp_mod_d(&b, 1853, &r);
+ printf("c = %04X\n", r);
- mp_addmod(&a, &b, &m, &c);
- printf("c = "); mp_print(&c, stdout); fputc('\n', stdout);
+ printf("\nc = (a + b) mod m\n");
- printf("\nc = (a - b) mod m\n");
+ mp_addmod(&a, &b, &m, &c);
+ printf("c = ");
+ mp_print(&c, stdout);
+ fputc('\n', stdout);
- mp_submod(&a, &b, &m, &c);
- printf("c = "); mp_print(&c, stdout); fputc('\n', stdout);
+ printf("\nc = (a - b) mod m\n");
- printf("\nc = (a * b) mod m\n");
+ mp_submod(&a, &b, &m, &c);
+ printf("c = ");
+ mp_print(&c, stdout);
+ fputc('\n', stdout);
- mp_mulmod(&a, &b, &m, &c);
- printf("c = "); mp_print(&c, stdout); fputc('\n', stdout);
+ printf("\nc = (a * b) mod m\n");
- printf("\nc = (a ** b) mod m\n");
+ mp_mulmod(&a, &b, &m, &c);
+ printf("c = ");
+ mp_print(&c, stdout);
+ fputc('\n', stdout);
- mp_exptmod(&a, &b, &m, &c);
- printf("c = "); mp_print(&c, stdout); fputc('\n', stdout);
+ printf("\nc = (a ** b) mod m\n");
- printf("\nIn-place modular squaring test:\n");
- for(ix = 0; ix < 5; ix++) {
- printf("a = (a * a) mod m a = ");
- mp_sqrmod(&a, &m, &a);
- mp_print(&a, stdout);
+ mp_exptmod(&a, &b, &m, &c);
+ printf("c = ");
+ mp_print(&c, stdout);
fputc('\n', stdout);
- }
-
- mp_clear(&c);
- mp_clear(&m);
- mp_clear(&b);
- mp_clear(&a);
+ printf("\nIn-place modular squaring test:\n");
+ for (ix = 0; ix < 5; ix++) {
+ printf("a = (a * a) mod m a = ");
+ mp_sqrmod(&a, &m, &a);
+ mp_print(&a, stdout);
+ fputc('\n', stdout);
+ }
+
+ mp_clear(&c);
+ mp_clear(&m);
+ mp_clear(&b);
+ mp_clear(&a);
- return 0;
+ return 0;
}