2013-08-17 12:20:15 +00:00
|
|
|
/**
|
|
|
|
* Copyright (c) 2013 Tomas Dzetkulic
|
|
|
|
* Copyright (c) 2013 Pavol Rusnak
|
|
|
|
*
|
|
|
|
* Permission is hereby granted, free of charge, to any person obtaining
|
|
|
|
* a copy of this software and associated documentation files (the "Software"),
|
|
|
|
* to deal in the Software without restriction, including without limitation
|
|
|
|
* the rights to use, copy, modify, merge, publish, distribute, sublicense,
|
|
|
|
* and/or sell copies of the Software, and to permit persons to whom the
|
|
|
|
* Software is furnished to do so, subject to the following conditions:
|
|
|
|
*
|
|
|
|
* The above copyright notice and this permission notice shall be included
|
|
|
|
* in all copies or substantial portions of the Software.
|
|
|
|
*
|
|
|
|
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
|
|
|
|
* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
|
|
|
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
|
|
|
|
* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES
|
|
|
|
* OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
|
|
|
|
* ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
|
|
|
|
* OTHER DEALINGS IN THE SOFTWARE.
|
|
|
|
*/
|
|
|
|
|
2013-09-12 01:15:22 +00:00
|
|
|
#ifndef __BIGNUM_H__
|
|
|
|
#define __BIGNUM_H__
|
2013-08-17 12:20:15 +00:00
|
|
|
|
|
|
|
#include <stdint.h>
|
|
|
|
|
2013-09-12 01:15:22 +00:00
|
|
|
// use precomputed Inverse Values of powers of two
|
2013-10-03 15:32:27 +00:00
|
|
|
#ifndef USE_PRECOMPUTED_IV
|
2013-09-12 01:15:22 +00:00
|
|
|
#define USE_PRECOMPUTED_IV 1
|
2013-10-03 15:32:27 +00:00
|
|
|
#endif
|
2013-09-12 01:15:22 +00:00
|
|
|
|
|
|
|
// use precomputed Curve Points (some scalar multiples of curve base point G)
|
2013-10-03 15:32:27 +00:00
|
|
|
#ifndef USE_PRECOMPUTED_CP
|
2013-09-12 01:15:22 +00:00
|
|
|
#define USE_PRECOMPUTED_CP 1
|
2013-10-03 15:32:27 +00:00
|
|
|
#endif
|
2013-09-12 01:15:22 +00:00
|
|
|
|
2013-10-03 15:32:27 +00:00
|
|
|
// use fast inverse method
|
|
|
|
#ifndef USE_INVERSE_FAST
|
|
|
|
#define USE_INVERSE_FAST 1
|
|
|
|
#endif
|
2013-09-12 01:15:22 +00:00
|
|
|
|
|
|
|
// bignum256 are 256 bits stored as 8*30 bit + 1*16 bit
|
|
|
|
// val[0] are lowest 30 bits, val[8] highest 16 bits
|
|
|
|
typedef struct {
|
|
|
|
uint32_t val[9];
|
|
|
|
} bignum256;
|
2013-08-17 12:20:15 +00:00
|
|
|
|
|
|
|
// read 4 big endian bytes into uint32
|
2013-09-11 17:02:22 +00:00
|
|
|
uint32_t read_be(const uint8_t *data);
|
2013-08-17 12:20:15 +00:00
|
|
|
|
|
|
|
// write 4 big endian bytes
|
2013-09-11 17:02:22 +00:00
|
|
|
void write_be(uint8_t *data, uint32_t x);
|
2013-08-17 12:20:15 +00:00
|
|
|
|
2013-09-12 01:15:22 +00:00
|
|
|
void bn_read_be(const uint8_t *in_number, bignum256 *out_number);
|
|
|
|
|
|
|
|
void bn_write_be(const bignum256 *in_number, uint8_t *out_number);
|
|
|
|
|
2013-09-19 00:54:10 +00:00
|
|
|
void bn_zero(bignum256 *a);
|
|
|
|
|
2013-09-12 01:15:22 +00:00
|
|
|
int bn_is_zero(const bignum256 *a);
|
|
|
|
|
|
|
|
int bn_is_less(const bignum256 *a, const bignum256 *b);
|
|
|
|
|
2014-01-04 16:39:37 +00:00
|
|
|
int bn_is_equal(const bignum256 *a, const bignum256 *b);
|
|
|
|
|
2013-09-19 12:52:52 +00:00
|
|
|
int bn_bitlen(const bignum256 *a);
|
2013-09-19 00:54:10 +00:00
|
|
|
|
|
|
|
void bn_lshift(bignum256 *a);
|
|
|
|
|
|
|
|
void bn_rshift(bignum256 *a);
|
|
|
|
|
2013-09-12 17:21:24 +00:00
|
|
|
void bn_mod(bignum256 *x, const bignum256 *prime);
|
2013-09-12 01:15:22 +00:00
|
|
|
|
2013-09-12 17:21:24 +00:00
|
|
|
void bn_multiply(const bignum256 *k, bignum256 *x, const bignum256 *prime);
|
2013-09-12 01:15:22 +00:00
|
|
|
|
2013-09-12 17:21:24 +00:00
|
|
|
void bn_fast_mod(bignum256 *x, const bignum256 *prime);
|
2013-09-12 01:15:22 +00:00
|
|
|
|
2013-09-24 13:54:36 +00:00
|
|
|
void bn_sqrt(bignum256 *x, const bignum256 *prime);
|
|
|
|
|
2013-09-12 17:21:24 +00:00
|
|
|
void bn_inverse(bignum256 *x, const bignum256 *prime);
|
|
|
|
|
2013-09-19 00:54:10 +00:00
|
|
|
void bn_normalize(bignum256 *a);
|
|
|
|
|
2013-09-12 17:21:24 +00:00
|
|
|
void bn_addmod(bignum256 *a, const bignum256 *b, const bignum256 *prime);
|
2013-09-12 01:15:22 +00:00
|
|
|
|
2013-09-24 13:54:36 +00:00
|
|
|
void bn_addmodi(bignum256 *a, uint32_t b, const bignum256 *prime);
|
|
|
|
|
2013-09-12 01:15:22 +00:00
|
|
|
void bn_substract(const bignum256 *a, const bignum256 *b, bignum256 *res);
|
|
|
|
|
2013-09-19 12:52:52 +00:00
|
|
|
void bn_substract_noprime(const bignum256 *a, const bignum256 *b, bignum256 *res);
|
|
|
|
|
2013-09-25 10:39:23 +00:00
|
|
|
void bn_divmod58(bignum256 *a, uint32_t *r);
|
2013-09-19 12:52:52 +00:00
|
|
|
|
2013-10-08 16:54:21 +00:00
|
|
|
#ifndef BN_PRINT
|
|
|
|
#define BN_PRINT 0
|
|
|
|
#endif
|
|
|
|
|
2013-10-08 12:01:20 +00:00
|
|
|
#if BN_PRINT
|
2013-09-24 13:54:36 +00:00
|
|
|
void bn_print(const bignum256 *a);
|
2013-10-08 12:01:20 +00:00
|
|
|
void bn_print_raw(const bignum256 *a);
|
2013-09-24 13:54:36 +00:00
|
|
|
#endif
|
|
|
|
|
2013-08-17 12:20:15 +00:00
|
|
|
#endif
|