229d3c467d
function old new delta pstm_div 1472 1522 +50 psRsaEncryptPub 403 413 +10 pstm_2expt 91 96 +5 pstm_clear 68 72 +4 pstm_init 39 42 +3 pstm_unsigned_bin_size 36 37 +1 pstm_montgomery_reduce 398 399 +1 pstm_init_size 45 46 +1 pstm_zero 39 38 -1 pstm_set 35 34 -1 pstm_read_unsigned_bin 112 109 -3 pstm_mulmod 123 120 -3 pstm_mod 116 113 -3 pstm_cmp 57 54 -3 pstm_sub 107 102 -5 pstm_to_unsigned_bin 157 151 -6 pstm_clamp 63 57 -6 pstm_add 116 108 -8 pstm_grow 81 72 -9 pstm_count_bits 57 48 -9 pstm_init_copy 84 72 -12 pstm_cmp_mag 93 78 -15 pstm_sqr_comba 567 551 -16 pstm_montgomery_calc_normalization 158 140 -18 pstm_copy 115 92 -23 pstm_lshd 133 109 -24 pstm_mul_comba 525 500 -25 pstm_mul_d 251 224 -27 s_pstm_sub 256 228 -28 s_pstm_add 370 337 -33 pstm_div_2d 444 409 -35 pstm_mul_2 195 156 -39 pstm_rshd 154 104 -50 pstm_mul_2d 247 186 -61 pstm_exptmod 1524 1463 -61 ------------------------------------------------------------------------------ (add/remove: 0/0 grow/shrink: 8/27 up/down: 75/-524) Total: -449 bytes Signed-off-by: Denys Vlasenko <vda.linux@googlemail.com>
287 lines
9.1 KiB
C
287 lines
9.1 KiB
C
/*
|
|
* 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 */
|
|
|