| /* |
| * 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; } |
| |
| //made static:extern void pstm_set(pstm_int *a, pstm_digit b); |
| |
| //made static:extern void pstm_zero(pstm_int * a); |
| |
| //bbox: pool unused |
| #define pstm_init(pool, a) \ |
| pstm_init( a) |
| //made static: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) FAST_FUNC; |
| |
| //bbox: pool unused |
| #define pstm_init_copy(pool, a, b, toSqr) \ |
| pstm_init_copy( a, b, toSqr) |
| //made static:extern int32 pstm_init_copy(psPool_t *pool, pstm_int * a, pstm_int * b, |
| //made static: int toSqr); //bbox: was int16 toSqr |
| |
| //made static:extern int pstm_count_bits (pstm_int * a) FAST_FUNC; //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) FAST_FUNC; |
| |
| extern int32 pstm_read_unsigned_bin(pstm_int *a, unsigned char *b, int32 c) FAST_FUNC; |
| |
| extern int32 pstm_unsigned_bin_size(pstm_int *a) FAST_FUNC; |
| |
| extern int32 pstm_copy(pstm_int * a, pstm_int * b); |
| |
| //made static:extern void pstm_exch(pstm_int * a, pstm_int * b); |
| |
| extern void pstm_clear(pstm_int * a) FAST_FUNC; |
| |
| 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) FAST_FUNC; |
| |
| extern int32 pstm_grow(pstm_int * a, int size) FAST_FUNC; //bbox: was int16 size |
| |
| extern void pstm_clamp(pstm_int * a) FAST_FUNC; |
| |
| extern int32 pstm_cmp(pstm_int * a, pstm_int * b) FAST_FUNC; |
| |
| extern int32 pstm_cmp_mag(pstm_int * a, pstm_int * b) FAST_FUNC; |
| |
| //made static:extern void pstm_rshd(pstm_int *a, int x); //bbox: was int16 x |
| |
| //made static: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) |
| //made static:extern int32 pstm_div(psPool_t *pool, pstm_int *a, pstm_int *b, pstm_int *c, |
| //made static: pstm_int *d); |
| |
| //bbox: pool unused |
| #define pstm_div_2d(pool, a, b, c, d) \ |
| pstm_div_2d( a, b, c, d) |
| //made static:extern int32 pstm_div_2d(psPool_t *pool, pstm_int *a, int b, pstm_int *c, |
| //made static: pstm_int *d); //bbox: was int16 b |
| |
| extern int32 pstm_div_2(pstm_int * a, pstm_int * b) FAST_FUNC; |
| |
| extern int32 s_pstm_sub(pstm_int *a, pstm_int *b, pstm_int *c) FAST_FUNC; |
| |
| extern int32 pstm_sub(pstm_int *a, pstm_int *b, pstm_int *c) FAST_FUNC; |
| |
| //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) FAST_FUNC; |
| |
| extern int32 pstm_mul_2(pstm_int * a, pstm_int * b) FAST_FUNC; |
| |
| //bbox: pool unused |
| #define pstm_mod(pool, a, b, c) \ |
| pstm_mod( a, b, c) |
| //made static: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) FAST_FUNC; |
| |
| //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) FAST_FUNC; |
| |
| //made static: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) FAST_FUNC; |
| |
| //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) FAST_FUNC; |
| |
| //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) FAST_FUNC; |
| |
| //made static: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) FAST_FUNC; |
| |
| #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) FAST_FUNC; |
| |
| //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) FAST_FUNC; |
| |
| //made static:extern int32 pstm_cmp_d(pstm_int *a, pstm_digit b); |
| |
| //made static:extern int32 pstm_montgomery_calc_normalization(pstm_int *a, pstm_int *b); |
| |
| //made static: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) FAST_FUNC; |
| |
| #else /* DISABLE_PSTM */ |
| typedef int32 pstm_int; |
| #endif /* !DISABLE_PSTM */ |
| #endif /* _h_PSTMATH */ |
| |