123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286 |
- /*
- * Copyright (C) 2017 Denys Vlasenko
- *
- * Licensed under GPLv2, see file LICENSE in this source tree.
- */
- /* The file is taken almost verbatim from matrixssl-3-7-2b-open/crypto/math/.
- * Changes are flagged with //bbox
- */
- /**
- * @file pstm.h
- * @version 33ef80f (HEAD, tag: MATRIXSSL-3-7-2-OPEN, tag: MATRIXSSL-3-7-2-COMM, origin/master, origin/HEAD, master)
- *
- * multiple-precision integer library.
- */
- /*
- * Copyright (c) 2013-2015 INSIDE Secure Corporation
- * Copyright (c) PeerSec Networks, 2002-2011
- * All Rights Reserved
- *
- * The latest version of this code is available at http://www.matrixssl.org
- *
- * This software is open source; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or
- * (at your option) any later version.
- *
- * This General Public License does NOT permit incorporating this software
- * into proprietary programs. If you are unable to comply with the GPL, a
- * commercial license for this software may be purchased from INSIDE at
- * http://www.insidesecure.com/eng/Company/Locations
- *
- * This program is distributed in WITHOUT ANY WARRANTY; without even the
- * implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
- * See the GNU General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
- * http://www.gnu.org/copyleft/gpl.html
- */
- /******************************************************************************/
- #ifndef _h_PSTMATH
- #define _h_PSTMATH
- #ifndef DISABLE_PSTM
- /* Define this here to avoid including circular limits.h on some platforms */
- #ifndef CHAR_BIT
- #define CHAR_BIT 8
- #endif
- /******************************************************************************/
- /*
- If native 64 bit integers are not supported, we do not support 32x32->64
- in hardware, so we must set the 16 bit flag to produce 16x16->32 products.
- */
- #ifndef HAVE_NATIVE_INT64
- #define PSTM_16BIT
- #endif /* ! HAVE_NATIVE_INT64 */
- /******************************************************************************/
- /*
- Some default configurations.
- pstm_word should be the largest value the processor can hold as the product
- of a multiplication. Most platforms support a 32x32->64 MAC instruction,
- so 64bits is the default pstm_word size.
- pstm_digit should be half the size of pstm_word
- */
- #ifdef PSTM_8BIT
- /* 8-bit digits, 16-bit word products */
- typedef unsigned char pstm_digit;
- typedef unsigned short pstm_word;
- #define DIGIT_BIT 8
- #elif defined(PSTM_16BIT)
- /* 16-bit digits, 32-bit word products */
- typedef unsigned short pstm_digit;
- typedef unsigned long pstm_word;
- #define DIGIT_BIT 16
- #elif defined(PSTM_64BIT)
- /* 64-bit digits, 128-bit word products */
- #ifndef __GNUC__
- #error "64bit digits requires GCC"
- #endif
- typedef unsigned long pstm_digit;
- typedef unsigned long pstm_word __attribute__ ((mode(TI)));
- #define DIGIT_BIT 64
- #else
- /* This is the default case, 32-bit digits, 64-bit word products */
- typedef uint32 pstm_digit;
- typedef uint64 pstm_word;
- #define DIGIT_BIT 32
- #define PSTM_32BIT
- #endif /* digit and word size */
- #define PSTM_MASK (pstm_digit)(-1)
- #define PSTM_DIGIT_MAX PSTM_MASK
- /******************************************************************************/
- /*
- equalities
- */
- #define PSTM_LT -1 /* less than */
- #define PSTM_EQ 0 /* equal to */
- #define PSTM_GT 1 /* greater than */
- #define PSTM_ZPOS 0 /* positive integer */
- #define PSTM_NEG 1 /* negative */
- #define PSTM_OKAY PS_SUCCESS
- #define PSTM_MEM PS_MEM_FAIL
- /******************************************************************************/
- /*
- Various build options
- */
- #define PSTM_DEFAULT_INIT 64 /* default (64) digits of allocation */
- #define PSTM_MAX_SIZE 4096
- typedef struct {
- int used, alloc, sign; //bbox: was int16
- pstm_digit *dp;
- //bbox psPool_t *pool;
- } pstm_int;
- /******************************************************************************/
- /*
- Operations on large integers
- */
- #define pstm_iszero(a) (((a)->used == 0) ? PS_TRUE : PS_FALSE)
- #define pstm_iseven(a) (((a)->used > 0 && (((a)->dp[0] & 1) == 0)) ? PS_TRUE : PS_FALSE)
- #define pstm_isodd(a) (((a)->used > 0 && (((a)->dp[0] & 1) == 1)) ? PS_TRUE : PS_FALSE)
- #define pstm_abs(a, b) { pstm_copy(a, b); (b)->sign = 0; }
- extern void pstm_set(pstm_int *a, pstm_digit b);
- extern void pstm_zero(pstm_int * a);
- //bbox: pool unused
- #define pstm_init(pool, a) \
- pstm_init( a)
- extern int32 pstm_init(psPool_t *pool, pstm_int * a);
- //bbox: pool unused
- #define pstm_init_size(pool, a, size) \
- pstm_init_size( a, size)
- extern int32 pstm_init_size(psPool_t *pool, pstm_int * a, uint32 size);
- //bbox: pool unused
- #define pstm_init_copy(pool, a, b, toSqr) \
- pstm_init_copy( a, b, toSqr)
- extern int32 pstm_init_copy(psPool_t *pool, pstm_int * a, pstm_int * b,
- int toSqr); //bbox: was int16 toSqr
- extern int pstm_count_bits (pstm_int * a); //bbox: was returning int16
- //bbox: pool unused
- #define pstm_init_for_read_unsigned_bin(pool, a, len) \
- pstm_init_for_read_unsigned_bin( a, len)
- extern int32 pstm_init_for_read_unsigned_bin(psPool_t *pool, pstm_int *a,
- uint32 len);
- extern int32 pstm_read_unsigned_bin(pstm_int *a, unsigned char *b, int32 c);
- extern int32 pstm_unsigned_bin_size(pstm_int *a);
- extern int32 pstm_copy(pstm_int * a, pstm_int * b);
- extern void pstm_exch(pstm_int * a, pstm_int * b);
- extern void pstm_clear(pstm_int * a);
- extern void pstm_clear_multi(pstm_int *mp0, pstm_int *mp1, pstm_int *mp2,
- pstm_int *mp3, pstm_int *mp4, pstm_int *mp5, pstm_int *mp6,
- pstm_int *mp7);
- extern int32 pstm_grow(pstm_int * a, int size); //bbox: was int16 size
- extern void pstm_clamp(pstm_int * a);
- extern int32 pstm_cmp(pstm_int * a, pstm_int * b);
- extern int32 pstm_cmp_mag(pstm_int * a, pstm_int * b);
- extern void pstm_rshd(pstm_int *a, int x); //bbox: was int16 x
- extern int32 pstm_lshd(pstm_int * a, int b); //bbox: was int16 b
- //bbox: pool unused
- #define pstm_div(pool, a, b, c, d) \
- pstm_div( a, b, c, d)
- extern int32 pstm_div(psPool_t *pool, pstm_int *a, pstm_int *b, pstm_int *c,
- pstm_int *d);
- //bbox: pool unused
- #define pstm_div_2d(pool, a, b, c, d) \
- pstm_div_2d( a, b, c, d)
- extern int32 pstm_div_2d(psPool_t *pool, pstm_int *a, int b, pstm_int *c,
- pstm_int *d); //bbox: was int16 b
- extern int32 pstm_div_2(pstm_int * a, pstm_int * b);
- extern int32 s_pstm_sub(pstm_int *a, pstm_int *b, pstm_int *c);
- extern int32 pstm_sub(pstm_int *a, pstm_int *b, pstm_int *c);
- //bbox: pool unused
- #define pstm_sub_d(pool, a, b, c) \
- pstm_sub_d( a, b, c)
- extern int32 pstm_sub_d(psPool_t *pool, pstm_int *a, pstm_digit b, pstm_int *c);
- extern int32 pstm_mul_2(pstm_int * a, pstm_int * b);
- //bbox: pool unused
- #define pstm_mod(pool, a, b, c) \
- pstm_mod( a, b, c)
- extern int32 pstm_mod(psPool_t *pool, pstm_int *a, pstm_int *b, pstm_int *c);
- //bbox: pool unused
- #define pstm_mulmod(pool, a, b, c, d) \
- pstm_mulmod( a, b, c, d)
- extern int32 pstm_mulmod(psPool_t *pool, pstm_int *a, pstm_int *b, pstm_int *c,
- pstm_int *d);
- //bbox: pool unused
- #define pstm_exptmod(pool, G, X, P, Y) \
- pstm_exptmod( G, X, P, Y)
- extern int32 pstm_exptmod(psPool_t *pool, pstm_int *G, pstm_int *X, pstm_int *P,
- pstm_int *Y);
- extern int32 pstm_2expt(pstm_int *a, int b); //bbox: was int16 b
- extern int32 pstm_add(pstm_int *a, pstm_int *b, pstm_int *c);
- //bbox: pool unused
- #define pstm_to_unsigned_bin(pool, a, b) \
- pstm_to_unsigned_bin( a, b)
- extern int32 pstm_to_unsigned_bin(psPool_t *pool, pstm_int *a,
- unsigned char *b);
- //bbox: pool unused
- #define pstm_to_unsigned_bin_nr(pool, a, b) \
- pstm_to_unsigned_bin_nr( a, b)
- extern int32 pstm_to_unsigned_bin_nr(psPool_t *pool, pstm_int *a,
- unsigned char *b);
- extern int32 pstm_montgomery_setup(pstm_int *a, pstm_digit *rho);
- //bbox: pool unused
- #define pstm_montgomery_reduce(pool, a, m, mp, paD, paDlen) \
- pstm_montgomery_reduce( a, m, mp, paD, paDlen)
- extern int32 pstm_montgomery_reduce(psPool_t *pool, pstm_int *a, pstm_int *m,
- pstm_digit mp, pstm_digit *paD, uint32 paDlen);
- #define pstm_mul_comba(pool, A, B, C, paD, paDlen) \
- pstm_mul_comba( A, B, C, paD, paDlen)
- extern int32 pstm_mul_comba(psPool_t *pool, pstm_int *A, pstm_int *B,
- pstm_int *C, pstm_digit *paD, uint32 paDlen);
- //bbox: pool unused
- #define pstm_sqr_comba(pool, A, B, paD, paDlen) \
- pstm_sqr_comba( A, B, paD, paDlen)
- extern int32 pstm_sqr_comba(psPool_t *pool, pstm_int *A, pstm_int *B,
- pstm_digit *paD, uint32 paDlen);
- extern int32 pstm_cmp_d(pstm_int *a, pstm_digit b);
- extern int32 pstm_montgomery_calc_normalization(pstm_int *a, pstm_int *b);
- extern int32 pstm_mul_d(pstm_int *a, pstm_digit b, pstm_int *c);
- //bbox: pool unused
- #define pstm_invmod(pool, a, b, c) \
- pstm_invmod( a, b, c)
- extern int32 pstm_invmod(psPool_t *pool, pstm_int * a, pstm_int * b,
- pstm_int * c);
- #else /* DISABLE_PSTM */
- typedef int32 pstm_int;
- #endif /* !DISABLE_PSTM */
- #endif /* _h_PSTMATH */
|