X-Git-Url: https://git.cryptolib.org/?p=avr-crypto-lib.git;a=blobdiff_plain;f=test_src%2Fmain-bigint-test.c;h=19ef467f3f08925c52c8cc2f7dd77e07d4f47f84;hp=b757f28f342e42882b02f871ecbfe0a414dd63d5;hb=4b5da1dc27a791b5c448274a3db09cd035b33493;hpb=fdbda6486df88be6ac2ba948f21e84c396499535 diff --git a/test_src/main-bigint-test.c b/test_src/main-bigint-test.c index b757f28..19ef467 100644 --- a/test_src/main-bigint-test.c +++ b/test_src/main-bigint-test.c @@ -1,7 +1,7 @@ /* main-bigint-test.c */ /* This file is part of the AVR-Crypto-Lib. - Copyright (C) 2008, 2009, 2010 Daniel Otte (daniel.otte@rub.de) + Copyright (C) 2006-2015 Daniel Otte (bg@nerilex.org) This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by @@ -30,17 +30,20 @@ #include "performance_test.h" -char* algo_name = "BigInt"; +char *algo_name = "BigInt"; + +#define MAX(a,b) ((a) > (b) ? (a) : (b)) +#define MIN(a,b) ((a) < (b) ? (a) : (b)) /***************************************************************************** * additional validation-functions * *****************************************************************************/ -void test_echo_bigint(void){ +void test_echo_bigint(void) { bigint_t a; cli_putstr_P(PSTR("\r\necho test\r\n")); - for(;;){ + for (;;) { cli_putstr_P(PSTR("\r\nenter hex number:")); - if(bigint_read_hex_echo(&a)){ + if (bigint_read_hex_echo(&a)) { cli_putstr_P(PSTR("\r\n end echo test")); return; } @@ -54,14 +57,14 @@ void test_echo_bigint(void){ void test_add_bigint(void){ bigint_t a, b, c; cli_putstr_P(PSTR("\r\nadd test\r\n")); - for(;;){ + for (;;) { cli_putstr_P(PSTR("\r\nenter a:")); - if(bigint_read_hex_echo(&a)){ + if (bigint_read_hex_echo(&a)) { cli_putstr_P(PSTR("\r\n end add test")); return; } cli_putstr_P(PSTR("\r\nenter b:")); - if(bigint_read_hex_echo(&b)){ + if (bigint_read_hex_echo(&b)) { free(a.wordv); cli_putstr_P(PSTR("\r\n end add test")); return; @@ -71,9 +74,9 @@ void test_add_bigint(void){ cli_putstr_P(PSTR(" + ")); bigint_print_hex(&b); cli_putstr_P(PSTR(" = ")); - uint8_t *c_b; - c_b = malloc(((a.length_W>b.length_W)?a.length_W:b.length_W)+2); - if(c_b==NULL){ + bigint_word_t *c_b; + c_b = malloc((MAX(a.length_W, b.length_W) + 2) * sizeof(bigint_word_t)); + if(c_b == NULL){ cli_putstr_P(PSTR("\n\rERROR: Out of memory!")); free(a.wordv); free(b.wordv); @@ -93,14 +96,14 @@ void test_add_scale_bigint(void){ bigint_t a, b, c; uint16_t scale; cli_putstr_P(PSTR("\r\nadd-scale test\r\n")); - for(;;){ + for (;;) { cli_putstr_P(PSTR("\r\nenter a:")); - if(bigint_read_hex_echo(&a)){ + if (bigint_read_hex_echo(&a)) { cli_putstr_P(PSTR("\r\n end add-scale test")); return; } cli_putstr_P(PSTR("\r\nenter b:")); - if(bigint_read_hex_echo(&b)){ + if (bigint_read_hex_echo(&b)) { cli_putstr_P(PSTR("\r\n end add-scale test")); return; } @@ -117,8 +120,8 @@ void test_add_scale_bigint(void){ return; } */ - uint8_t *c_b; - c_b = malloc(((a.length_W>(b.length_W+scale))?a.length_W:(b.length_W+scale))+2); + bigint_word_t *c_b; + c_b = malloc((MAX(a.length_W, b.length_W+scale) + 2) * sizeof(bigint_word_t)); if(c_b==NULL){ cli_putstr_P(PSTR("\n\rERROR: Out of memory!")); free(a.wordv); @@ -146,14 +149,14 @@ void test_add_scale_bigint(void){ void test_mul_bigint(void){ bigint_t a, b, c; cli_putstr_P(PSTR("\r\nmul test\r\n")); - for(;;){ + for (;;) { cli_putstr_P(PSTR("\r\nenter a:")); - if(bigint_read_hex_echo(&a)){ + if (bigint_read_hex_echo(&a)) { cli_putstr_P(PSTR("\r\n end mul test")); return; } cli_putstr_P(PSTR("\r\nenter b:")); - if(bigint_read_hex_echo(&b)){ + if (bigint_read_hex_echo(&b)) { free(a.wordv); cli_putstr_P(PSTR("\r\n end mul test")); return; @@ -163,9 +166,9 @@ void test_mul_bigint(void){ cli_putstr_P(PSTR(" * ")); bigint_print_hex(&b); cli_putstr_P(PSTR(" = ")); - uint8_t *c_b; - c_b = malloc((((a.length_W>b.length_W)?a.length_W:b.length_W)+1)*2); - if(c_b==NULL){ + bigint_word_t *c_b; + c_b = malloc((MAX(a.length_W, b.length_W) + 1) * 2 * sizeof(bigint_word_t)); + if (c_b==NULL) { cli_putstr_P(PSTR("\n\rERROR: Out of memory!")); free(a.wordv); free(b.wordv); @@ -181,6 +184,99 @@ void test_mul_bigint(void){ } } +void test_mul_mont_bigint(void){ + bigint_t a, b, c, a_, b_, m_, res; + bigint_length_t s; + cli_putstr_P(PSTR("\r\nmul-mont test ( (a * b) % c )\r\n")); + for (;;) { + cli_putstr_P(PSTR("\r\nenter a:")); + if (bigint_read_hex_echo(&a)) { + cli_putstr_P(PSTR("\r\n end mul test")); + return; + } + cli_putstr_P(PSTR("\r\nenter b:")); + if (bigint_read_hex_echo(&b)) { + free(a.wordv); + cli_putstr_P(PSTR("\r\n end mul test")); + return; + } + cli_putstr_P(PSTR("\r\nenter c:")); + if (bigint_read_hex_echo(&c)) { + free(a.wordv); + free(b.wordv); + cli_putstr_P(PSTR("\r\n end mul test")); + return; + } + s = c.length_W; + cli_putstr_P(PSTR("\r\n (")); + bigint_print_hex(&a); + cli_putstr_P(PSTR(" * ")); + bigint_print_hex(&b); + cli_putstr_P(PSTR(") % ")); + bigint_print_hex(&c); + cli_putstr_P(PSTR(" = ")); + bigint_word_t res_w[s], a_w_[s], b_w_[s], m_w_[s + 1]; + res.wordv = res_w; + a_.wordv = a_w_; + b_.wordv = b_w_; + m_.wordv = m_w_; + bigint_mont_gen_m_(&m_, &c); + bigint_mont_trans(&a_, &a, &c); + bigint_mont_trans(&b_, &b, &c); + bigint_mont_mul(&res, &a_, &b_, &c, &m_); + bigint_mont_red(&res, &res, &c, &m_); + bigint_print_hex(&res); + putchar('\n'); + free(a.wordv); + free(b.wordv); + free(c.wordv); + } +} + +void test_mul_word_bigint(void){ + bigint_t a, b; + bigint_word_t *t; + cli_putstr_P(PSTR("\r\nmul test\r\n")); + for (;;) { + cli_putstr_P(PSTR("\r\nenter a:")); + if (bigint_read_hex_echo(&a)) { + cli_putstr_P(PSTR("\r\n end mul test")); + return; + } + cli_putstr_P(PSTR("\r\nenter b:")); + if (bigint_read_hex_echo(&b)) { + free(a.wordv); + cli_putstr_P(PSTR("\r\n end mul test")); + return; + } + cli_putstr_P(PSTR("\r\n ")); + bigint_print_hex(&a); + cli_putstr_P(PSTR(" * ")); + bigint_print_hex(&b); + cli_putstr_P(PSTR(" = ")); + + if (b.length_W > 1) { + free(a.wordv); + free(b.wordv); + cli_putstr_P(PSTR("\r\n end mul test")); + } + + t = realloc(a.wordv, (a.length_W + 3) * sizeof(bigint_word_t)); + if (t == NULL) { + cli_putstr_P(PSTR("\n\rERROR: Out of memory!")); + free(a.wordv); + free(b.wordv); + continue; + } + a.wordv = t; + bigint_mul_word_u(&a, b.wordv[0]); + bigint_print_hex(&a); + cli_putstr_P(PSTR("\r\n")); + free(a.wordv); + free(b.wordv); + } +} + void test_square_bigint(void){ bigint_t a, c; cli_putstr_P(PSTR("\r\nsquare test\r\n")); @@ -193,9 +289,9 @@ void test_square_bigint(void){ cli_putstr_P(PSTR("\r\n ")); bigint_print_hex(&a); cli_putstr_P(PSTR("**2 = ")); - uint8_t *c_b; - c_b = malloc(a.length_W*2); - if(c_b==NULL){ + bigint_word_t *c_b; + c_b = malloc(a.length_W * 2 * sizeof(bigint_word_t)); + if(c_b == NULL){ cli_putstr_P(PSTR("\n\rERROR: Out of memory!")); free(a.wordv); continue; @@ -212,14 +308,14 @@ void test_square_bigint(void){ void test_reduce_bigint(void){ bigint_t a, b; cli_putstr_P(PSTR("\r\nreduce test\r\n")); - for(;;){ + for (;;) { cli_putstr_P(PSTR("\r\nenter a:")); - if(bigint_read_hex_echo(&a)){ + if (bigint_read_hex_echo(&a)) { cli_putstr_P(PSTR("\r\n end reduce test")); return; } cli_putstr_P(PSTR("\r\nenter b:")); - if(bigint_read_hex_echo(&b)){ + if (bigint_read_hex_echo(&b)) { free(a.wordv); cli_putstr_P(PSTR("\r\n end reduce test")); return; @@ -239,28 +335,28 @@ void test_reduce_bigint(void){ /* d = a**b % c */ void test_expmod_bigint(void){ bigint_t a, b, c, d; - uint8_t *d_b; - cli_putstr_P(PSTR("\r\nreduce test\r\n")); - for(;;){ + bigint_word_t *d_b; + cli_putstr_P(PSTR("\r\nexpnonentiation-modulo test\r\n")); + for (;;) { cli_putstr_P(PSTR("\r\nenter a:")); - if(bigint_read_hex_echo(&a)){ + if (bigint_read_hex_echo(&a)) { cli_putstr_P(PSTR("\r\n end expmod test")); return; } cli_putstr_P(PSTR("\r\nenter b:")); - if(bigint_read_hex_echo(&b)){ + if (bigint_read_hex_echo(&b)) { free(a.wordv); cli_putstr_P(PSTR("\r\n end expmod test")); return; } cli_putstr_P(PSTR("\r\nenter c:")); - if(bigint_read_hex_echo(&c)){ + if (bigint_read_hex_echo(&c)) { free(a.wordv); free(b.wordv); cli_putstr_P(PSTR("\r\n end expmod test")); return; } - d_b = malloc(c.length_W); + d_b = malloc(c.length_W * sizeof(bigint_word_t)); if(d_b==NULL){ cli_putstr_P(PSTR("\n\rERROR: Out of memory!")); free(a.wordv); @@ -276,7 +372,7 @@ void test_expmod_bigint(void){ cli_putstr_P(PSTR(" % ")); bigint_print_hex(&c); cli_putstr_P(PSTR(" = ")); - bigint_expmod_u(&d, &a, &b, &c); + bigint_expmod_u_sam(&d, &a, &b, &c); bigint_print_hex(&d); cli_putstr_P(PSTR("\r\n")); free(a.wordv); @@ -287,24 +383,75 @@ void test_expmod_bigint(void){ } } +/* d = a**b % c */ +void test_expmod_mont_bigint(void){ + bigint_t a, b, c, d; + bigint_word_t *d_b; + cli_putstr_P(PSTR("\r\nexpnonentiation-modulo-montgomory test\r\n")); + for (;;) { + cli_putstr_P(PSTR("\r\nenter a:")); + if (bigint_read_hex_echo(&a)) { + cli_putstr_P(PSTR("\r\n end expmod test")); + return; + } + cli_putstr_P(PSTR("\r\nenter b:")); + if (bigint_read_hex_echo(&b)) { + free(a.wordv); + cli_putstr_P(PSTR("\r\n end expmod test")); + return; + } + cli_putstr_P(PSTR("\r\nenter c:")); + if (bigint_read_hex_echo(&c)) { + free(a.wordv); + free(b.wordv); + cli_putstr_P(PSTR("\r\n end expmod test")); + return; + } + d_b = malloc(c.length_W * sizeof(bigint_word_t)); + if (d_b == NULL) { + cli_putstr_P(PSTR("\n\rERROR: Out of memory!")); + free(a.wordv); + free(b.wordv); + free(c.wordv); + continue; + } + d.wordv = d_b; + cli_putstr_P(PSTR("\r\n ")); + bigint_print_hex(&a); + cli_putstr_P(PSTR("**")); + bigint_print_hex(&b); + cli_putstr_P(PSTR(" % ")); + bigint_print_hex(&c); + cli_putstr_P(PSTR(" = ")); + bigint_expmod_u_mont_sam(&d, &a, &b, &c); + bigint_print_hex(&d); + cli_putstr_P(PSTR("\r\n")); + free(a.wordv); + free(b.wordv); + free(c.wordv); + free(d.wordv); + + } +} + void test_gcdext_bigint(void){ bigint_t a, b, c, d, e; cli_putstr_P(PSTR("\r\ngcdext test\r\n")); - for(;;){ + for (;;) { cli_putstr_P(PSTR("\r\nenter a:")); - if(bigint_read_hex_echo(&a)){ + if (bigint_read_hex_echo(&a)) { cli_putstr_P(PSTR("\r\n end gcdext test")); return; } cli_putstr_P(PSTR("\r\nenter b:")); - if(bigint_read_hex_echo(&b)){ + if (bigint_read_hex_echo(&b)) { free(a.wordv); cli_putstr_P(PSTR("\r\n end gcdext test")); return; } - c.wordv = malloc((a.length_Wb.length_W)?a.length_W:b.length_W); - e.wordv = malloc(1+(a.length_W>b.length_W)?a.length_W:b.length_W); + c.wordv = malloc(MIN(a.length_W, b.length_W) * sizeof(bigint_word_t)); + d.wordv = malloc((MAX(a.length_W, b.length_W) + 1) * sizeof(bigint_word_t)); + e.wordv = malloc((MAX(a.length_W, b.length_W) + 1) * sizeof(bigint_word_t)); cli_putstr_P(PSTR("\r\n gcdext( ")); bigint_print_hex(&a); @@ -330,10 +477,10 @@ void test_gcdext_bigint(void){ void test_simple(void){ bigint_t a, b, c; - uint8_t a_b[1], b_b[1], c_b[2]; - a.wordv=a_b; - b.wordv=b_b; - c.wordv=c_b; + bigint_word_t a_b[1], b_b[1], c_b[2]; + a.wordv = a_b; + b.wordv = b_b; + c.wordv = c_b; a.length_W = 1; b.length_W = 1; a_b[0] = 1; @@ -373,12 +520,12 @@ void test_mul_simple(void){ uint8_t b_b[8] = {0xe6, 0xdd, 0xce, 0x00, 0x44, 0x60, 0xda, 0x0d}; uint8_t c_b[16]; - a.wordv=a_b; - b.wordv=b_b; - c.wordv=c_b; - a.length_W = 8; - b.length_W = 8; - a.info=0x80; + a.wordv = (bigint_word_t*)a_b; + b.wordv = (bigint_word_t*)b_b; + c.wordv = (bigint_word_t*)c_b; + a.length_W = 8 / sizeof(bigint_word_t); + b.length_W = 8 / sizeof(bigint_word_t); + a.info = 0x80; bigint_adjust(&a); bigint_adjust(&b); bigint_mul_s(&c, &a, &b); @@ -414,12 +561,12 @@ const uint8_t square_test_data[] PROGMEM = { void test_square_simple(void){ bigint_t a, c; - uint8_t a_b[11] = {0xe6, 0x70, 0x7d, 0x43, 0x74, 0x07, 0x20, 0x22, 0x6a, 0xb8, 0xf4}; - uint8_t c_b[22]; - a.wordv=a_b; - c.wordv=c_b; - a.length_W = 11; - a.info=0x00; + uint8_t a_b[16] = {0xe6, 0x70, 0x7d, 0x43, 0x74, 0x07, 0x20, 0x22, 0x6a, 0xb8, 0xf4, 0, 0, 0, 0, 0}; + uint8_t c_b[32]; + a.wordv = (bigint_word_t*)a_b; + c.wordv = (bigint_word_t*)c_b; + a.length_W = 16 / sizeof(bigint_word_t); + a.info = 0x00; bigint_adjust(&a); bigint_square(&c, &a); cli_putstr_P(PSTR("\r\n test: ")); @@ -432,18 +579,18 @@ void test_square_simple(void){ void test_reduce_simple(void){ bigint_t a, b, c; - uint8_t a_b[2] = {0x62, 0xA8}; - uint8_t b_b[2] = {0x52, 0x27}; - uint8_t c_b[2]; - a.wordv=a_b; - a.length_W = 2; - a.info=0x00; + uint8_t a_b[4] = {0x62, 0xA8}; + uint8_t b_b[4] = {0x52, 0x27}; + uint8_t c_b[4]; + a.wordv = (bigint_word_t*)a_b; + a.length_W = 4 / sizeof(bigint_word_t); + a.info = 0x00; bigint_adjust(&a); - b.wordv=b_b; - b.length_W = 2; - b.info=0x00; + b.wordv = (bigint_word_t*)b_b; + b.length_W = 4 / sizeof(bigint_word_t); + b.info = 0x00; bigint_adjust(&b); - c.wordv = c_b; + c.wordv = (bigint_word_t*)c_b; bigint_copy(&c, &a); bigint_reduce(&c, &b); cli_putstr_P(PSTR("\r\n test: ")); @@ -461,20 +608,20 @@ void test_reduce_simple(void){ void test_gcdext_simple(void){ bigint_t a, b, c, d, e; - uint8_t a_b[5] = {0x71, 0x07, 0x00, 0x09, 0x16}; - uint8_t b_b[5] = {0x72, 0x7D, 0x57, 0xAC, 0X6F}; - uint8_t c_b[6], d_b[6], e_b[6]; - a.wordv=a_b; - a.length_W = 5; - a.info=0x00; + uint8_t a_b[8] = {0x71, 0x07, 0x00, 0x09, 0x16}; + uint8_t b_b[8] = {0x72, 0x7D, 0x57, 0xAC, 0X6F}; + uint8_t c_b[8], d_b[8], e_b[8]; + a.wordv = (bigint_word_t*)a_b; + a.length_W = 8 / sizeof(bigint_word_t); + a.info = 0x00; bigint_adjust(&a); - b.wordv=b_b; - b.length_W = 5; - b.info=0x00; + b.wordv = (bigint_word_t*)b_b; + b.length_W = 8 / sizeof(bigint_word_t); + b.info = 0x00; bigint_adjust(&b); - c.wordv = c_b; - d.wordv = d_b; - e.wordv = e_b; + c.wordv = (bigint_word_t*)c_b; + d.wordv = (bigint_word_t*)d_b; + e.wordv = (bigint_word_t*)e_b; bigint_gcdext(&c, &d, &e, &a, &b); cli_putstr_P(PSTR("\r\n test: gcd( ")); bigint_print_hex(&a); @@ -499,9 +646,12 @@ const char echo_test_str[] PROGMEM = "echo-test"; const char add_test_str[] PROGMEM = "add-test"; const char add_scale_test_str[] PROGMEM = "add-scale-test"; const char mul_test_str[] PROGMEM = "mul-test"; +const char mul_mont_test_str[] PROGMEM = "mul-mont-test"; +const char mul_word_test_str[] PROGMEM = "mul-word-test"; const char square_test_str[] PROGMEM = "square-test"; const char reduce_test_str[] PROGMEM = "reduce-test"; const char expmod_test_str[] PROGMEM = "expmod-test"; +const char expmod_mont_test_str[] PROGMEM = "expmod-mont-test"; const char gcdext_test_str[] PROGMEM = "gcdext-test"; const char quick_test_str[] PROGMEM = "quick-test"; const char performance_str[] PROGMEM = "performance"; @@ -511,9 +661,12 @@ const cmdlist_entry_t cmdlist[] PROGMEM = { { add_test_str, NULL, test_add_bigint }, { add_scale_test_str, NULL, test_add_scale_bigint }, { mul_test_str, NULL, test_mul_bigint }, + { mul_mont_test_str, NULL, test_mul_mont_bigint }, + { mul_word_test_str, NULL, test_mul_word_bigint }, { square_test_str, NULL, test_square_bigint }, { reduce_test_str, NULL, test_reduce_bigint }, - { expmod_test_str, NULL, test_expmod_bigint }, + { expmod_test_str, NULL, test_expmod_bigint }, + { expmod_mont_test_str, NULL, test_expmod_mont_bigint }, { gcdext_test_str, NULL, test_gcdext_bigint }, { quick_test_str, NULL, test_gcdext_simple }, { echo_test_str, NULL, test_echo_bigint },