2015-11-11 19:49:07 -05:00
|
|
|
/* LibTomMath, multiple-precision integer library -- Tom St Denis
|
|
|
|
*
|
|
|
|
* LibTomMath is a library that provides multiple-precision
|
|
|
|
* integer arithmetic as well as number theoretic functionality.
|
|
|
|
*
|
|
|
|
* The library was designed directly after the MPI library by
|
|
|
|
* Michael Fromberger but has been written from scratch with
|
|
|
|
* additional optimizations in place.
|
|
|
|
*
|
2018-12-29 11:56:20 -05:00
|
|
|
* SPDX-License-Identifier: Unlicense
|
2015-11-11 19:49:07 -05:00
|
|
|
*/
|
|
|
|
#ifndef TOMMATH_PRIV_H_
|
|
|
|
#define TOMMATH_PRIV_H_
|
|
|
|
|
2018-05-02 15:43:17 -04:00
|
|
|
#include "tommath.h"
|
2015-11-19 15:34:58 -05:00
|
|
|
#include <ctype.h>
|
2015-11-11 19:49:07 -05:00
|
|
|
|
2017-03-31 18:35:03 -04:00
|
|
|
#ifndef MIN
|
2017-08-30 13:15:27 -04:00
|
|
|
#define MIN(x, y) (((x) < (y)) ? (x) : (y))
|
2017-03-31 18:35:03 -04:00
|
|
|
#endif
|
2015-11-19 15:53:32 -05:00
|
|
|
|
2017-03-31 18:35:03 -04:00
|
|
|
#ifndef MAX
|
2017-08-30 13:15:27 -04:00
|
|
|
#define MAX(x, y) (((x) > (y)) ? (x) : (y))
|
2017-03-31 18:35:03 -04:00
|
|
|
#endif
|
2015-11-19 15:53:32 -05:00
|
|
|
|
2015-11-11 19:49:07 -05:00
|
|
|
#ifdef __cplusplus
|
|
|
|
extern "C" {
|
2015-11-19 15:53:32 -05:00
|
|
|
|
|
|
|
/* C++ compilers don't like assigning void * to mp_digit * */
|
2017-08-29 23:56:38 -04:00
|
|
|
#define OPT_CAST(x) (x *)
|
2015-11-19 15:53:32 -05:00
|
|
|
|
|
|
|
#else
|
|
|
|
|
|
|
|
/* C on the other hand doesn't care */
|
2017-08-29 23:56:38 -04:00
|
|
|
#define OPT_CAST(x)
|
2015-11-19 15:53:32 -05:00
|
|
|
|
|
|
|
#endif
|
|
|
|
|
|
|
|
/* define heap macros */
|
|
|
|
#ifndef XMALLOC
|
2017-08-30 13:15:27 -04:00
|
|
|
/* default to libc stuff */
|
2017-08-29 23:56:38 -04:00
|
|
|
# define XMALLOC malloc
|
|
|
|
# define XFREE free
|
|
|
|
# define XREALLOC realloc
|
|
|
|
# define XCALLOC calloc
|
2015-11-19 15:53:32 -05:00
|
|
|
#else
|
2017-08-30 13:15:27 -04:00
|
|
|
/* prototypes for our heap functions */
|
|
|
|
extern void *XMALLOC(size_t n);
|
|
|
|
extern void *XREALLOC(void *p, size_t n);
|
|
|
|
extern void *XCALLOC(size_t n, size_t s);
|
|
|
|
extern void XFREE(void *p);
|
2015-11-11 19:49:07 -05:00
|
|
|
#endif
|
|
|
|
|
|
|
|
/* lowlevel functions, do not call! */
|
2017-09-20 10:59:43 -04:00
|
|
|
int s_mp_add(const mp_int *a, const mp_int *b, mp_int *c);
|
|
|
|
int s_mp_sub(const mp_int *a, const mp_int *b, mp_int *c);
|
2015-11-11 19:49:07 -05:00
|
|
|
#define s_mp_mul(a, b, c) s_mp_mul_digs(a, b, c, (a)->used + (b)->used + 1)
|
2017-09-20 10:59:43 -04:00
|
|
|
int fast_s_mp_mul_digs(const mp_int *a, const mp_int *b, mp_int *c, int digs);
|
|
|
|
int s_mp_mul_digs(const mp_int *a, const mp_int *b, mp_int *c, int digs);
|
|
|
|
int fast_s_mp_mul_high_digs(const mp_int *a, const mp_int *b, mp_int *c, int digs);
|
|
|
|
int s_mp_mul_high_digs(const mp_int *a, const mp_int *b, mp_int *c, int digs);
|
|
|
|
int fast_s_mp_sqr(const mp_int *a, mp_int *b);
|
|
|
|
int s_mp_sqr(const mp_int *a, mp_int *b);
|
|
|
|
int mp_karatsuba_mul(const mp_int *a, const mp_int *b, mp_int *c);
|
|
|
|
int mp_toom_mul(const mp_int *a, const mp_int *b, mp_int *c);
|
|
|
|
int mp_karatsuba_sqr(const mp_int *a, mp_int *b);
|
|
|
|
int mp_toom_sqr(const mp_int *a, mp_int *b);
|
|
|
|
int fast_mp_invmod(const mp_int *a, const mp_int *b, mp_int *c);
|
|
|
|
int mp_invmod_slow(const mp_int *a, const mp_int *b, mp_int *c);
|
|
|
|
int fast_mp_montgomery_reduce(mp_int *x, const mp_int *n, mp_digit rho);
|
|
|
|
int mp_exptmod_fast(const mp_int *G, const mp_int *X, const mp_int *P, mp_int *Y, int redmode);
|
|
|
|
int s_mp_exptmod(const mp_int *G, const mp_int *X, const mp_int *P, mp_int *Y, int redmode);
|
2015-11-11 19:49:07 -05:00
|
|
|
void bn_reverse(unsigned char *s, int len);
|
|
|
|
|
2018-05-02 21:44:47 -04:00
|
|
|
extern const char *const mp_s_rmap;
|
2017-10-09 08:34:20 -04:00
|
|
|
extern const uint8_t mp_s_rmap_reverse[];
|
|
|
|
extern const size_t mp_s_rmap_reverse_sz;
|
2015-11-11 19:49:07 -05:00
|
|
|
|
|
|
|
/* Fancy macro to set an MPI from another type.
|
|
|
|
* There are several things assumed:
|
|
|
|
* x is the counter and unsigned
|
|
|
|
* a is the pointer to the MPI
|
|
|
|
* b is the original value that should be set in the MPI.
|
|
|
|
*/
|
|
|
|
#define MP_SET_XLONG(func_name, type) \
|
|
|
|
int func_name (mp_int * a, type b) \
|
|
|
|
{ \
|
|
|
|
unsigned int x; \
|
|
|
|
int res; \
|
|
|
|
\
|
|
|
|
mp_zero (a); \
|
|
|
|
\
|
|
|
|
/* set four bits at a time */ \
|
2015-11-13 11:47:30 -05:00
|
|
|
for (x = 0; x < (sizeof(type) * 2u); x++) { \
|
2015-11-11 19:49:07 -05:00
|
|
|
/* shift the number up four bits */ \
|
|
|
|
if ((res = mp_mul_2d (a, 4, a)) != MP_OKAY) { \
|
|
|
|
return res; \
|
|
|
|
} \
|
|
|
|
\
|
|
|
|
/* OR in the top four bits of the source */ \
|
2017-10-15 13:58:35 -04:00
|
|
|
a->dp[0] |= (mp_digit)(b >> ((sizeof(type) * 8u) - 4u)) & 15uL;\
|
2015-11-11 19:49:07 -05:00
|
|
|
\
|
|
|
|
/* shift the source up to the next four bits */ \
|
|
|
|
b <<= 4; \
|
|
|
|
\
|
|
|
|
/* ensure that digits are not clamped off */ \
|
|
|
|
a->used += 1; \
|
|
|
|
} \
|
|
|
|
mp_clamp (a); \
|
|
|
|
return MP_OKAY; \
|
|
|
|
}
|
|
|
|
|
|
|
|
#ifdef __cplusplus
|
2017-08-30 13:15:27 -04:00
|
|
|
}
|
2015-11-11 19:49:07 -05:00
|
|
|
#endif
|
|
|
|
|
|
|
|
#endif
|
|
|
|
|
|
|
|
|
2017-08-28 10:27:26 -04:00
|
|
|
/* ref: $Format:%D$ */
|
|
|
|
/* git commit: $Format:%H$ */
|
|
|
|
/* commit time: $Format:%ai$ */
|