raw
mpi-genesis~            1 /* mpi-internal.h  -  Internal to the Multi Precision Integers
mpi-genesis~            2  *	Copyright (C) 1994, 1996 Free Software Foundation, Inc.
mpi-genesis~            3  *	Copyright (C) 1998, 2000 Free Software Foundation, Inc.
mpi-genesis~            4  *
mpi-genesis~            5  * This file is part of GnuPG.
mpi-genesis~            6  *
mpi-genesis~            7  * GnuPG is free software; you can redistribute it and/or modify
mpi-genesis~            8  * it under the terms of the GNU General Public License as published by
mpi-genesis~            9  * the Free Software Foundation; either version 3 of the License, or
mpi-genesis~           10  * (at your option) any later version.
mpi-genesis~           11  *
mpi-genesis~           12  * GnuPG is distributed in the hope that it will be useful,
mpi-genesis~           13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
mpi-genesis~           14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
mpi-genesis~           15  * GNU General Public License for more details.
mpi-genesis~           16  *
mpi-genesis~           17  * You should have received a copy of the GNU General Public License
mpi-genesis~           18  * along with this program; if not, see <http://www.gnu.org/licenses/>.
mpi-genesis~           19  *
mpi-genesis~           20  * Note: This code is heavily based on the GNU MP Library.
mpi-genesis~           21  *	 Actually it's the same code with only minor changes in the
mpi-genesis~           22  *	 way the data is stored; this is to support the abstraction
mpi-genesis~           23  *	 of an optional secure memory allocation which may be used
mpi-genesis~           24  *	 to avoid revealing of sensitive data due to paging etc.
mpi-genesis~           25  *	 The GNU MP Library itself is published under the LGPL;
mpi-genesis~           26  *	 however I decided to publish this code under the plain GPL.
mpi-genesis~           27  */
mpi-genesis~           28 
mpi-genesis~           29 #ifndef G10_MPI_INTERNAL_H
mpi-genesis~           30 #define G10_MPI_INTERNAL_H
mpi-genesis~           31 
mpi-genesis~           32 #include "mpi.h"
mpi-genesis~           33 #include "mpi-asm-defs.h"
mpi-genesis~           34 
mpi-genesis~           35 #if BYTES_PER_MPI_LIMB == SIZEOF_UNSIGNED_INT
mpi-genesis~           36   typedef unsigned int mpi_limb_t;
mpi-genesis~           37   typedef   signed int mpi_limb_signed_t;
mpi-genesis~           38 #elif BYTES_PER_MPI_LIMB == SIZEOF_UNSIGNED_LONG
mpi-genesis~           39   typedef unsigned long int mpi_limb_t;
mpi-genesis~           40   typedef   signed long int mpi_limb_signed_t;
mpi-genesis~           41 #elif BYTES_PER_MPI_LIMB == SIZEOF_UNSIGNED_LONG_LONG
mpi-genesis~           42   typedef unsigned long long int mpi_limb_t;
mpi-genesis~           43   typedef   signed long long int mpi_limb_signed_t;
mpi-genesis~           44 #elif BYTES_PER_MPI_LIMB == SIZEOF_UNSIGNED_SHORT
mpi-genesis~           45   typedef unsigned short int mpi_limb_t;
mpi-genesis~           46   typedef   signed short int mpi_limb_signed_t;
mpi-genesis~           47 #else
mpi-genesis~           48 #error BYTES_PER_MPI_LIMB does not match any C type
mpi-genesis~           49 #endif
mpi-genesis~           50 #define BITS_PER_MPI_LIMB    (8*BYTES_PER_MPI_LIMB)
mpi-genesis~           51 
mpi-genesis~           52 
mpi-genesis~           53 struct gcry_mpi {
mpi-genesis~           54     int alloced;    /* array size (# of allocated limbs) */
mpi-genesis~           55     int nlimbs;     /* number of valid limbs */
mpi-genesis~           56     unsigned int nbits; /* the real number of valid bits (info only) */
mpi-genesis~           57     int sign;	    /* indicates a negative number */
mpi-genesis~           58     unsigned flags; /* bit 0: array must be allocated in secure memory space */
mpi-genesis~           59 		    /* bit 1: not used */
mpi-genesis~           60 		    /* bit 2: the limb is a pointer to some xmalloced data */
mpi-genesis~           61     mpi_limb_t *d;  /* array with the limbs */
mpi-genesis~           62 };
mpi-genesis~           63 
mpi-genesis~           64 
mpi-genesis~           65 
mpi-genesis~           66 /* If KARATSUBA_THRESHOLD is not already defined, define it to a
mpi-genesis~           67  * value which is good on most machines.  */
mpi-genesis~           68 
mpi-genesis~           69 /* tested 4, 16, 32 and 64, where 16 gave the best performance when
mpi-genesis~           70  * checking a 768 and a 1024 bit ElGamal signature.
mpi-genesis~           71  * (wk 22.12.97) */
mpi-genesis~           72 #ifndef KARATSUBA_THRESHOLD
mpi-genesis~           73 #define KARATSUBA_THRESHOLD 16
mpi-genesis~           74 #endif
mpi-genesis~           75 
mpi-genesis~           76 /* The code can't handle KARATSUBA_THRESHOLD smaller than 2.  */
mpi-genesis~           77 #if KARATSUBA_THRESHOLD < 2
mpi-genesis~           78 #undef KARATSUBA_THRESHOLD
mpi-genesis~           79 #define KARATSUBA_THRESHOLD 2
mpi-genesis~           80 #endif
mpi-genesis~           81 
mpi-genesis~           82 
mpi-genesis~           83 typedef mpi_limb_t *mpi_ptr_t; /* pointer to a limb */
mpi-genesis~           84 typedef int mpi_size_t;        /* (must be a signed type) */
mpi-genesis~           85 
mpi-genesis~           86 #define ABS(x) (x >= 0 ? x : -x)
mpi-genesis~           87 #define MIN(l,o) ((l) < (o) ? (l) : (o))
mpi-genesis~           88 #define MAX(h,i) ((h) > (i) ? (h) : (i))
mpi-genesis~           89 #define RESIZE_IF_NEEDED(a,b) \
mpi-genesis~           90     do {			   \
mpi-genesis~           91 	if( (a)->alloced < (b) )   \
mpi-genesis~           92 	    mpi_resize((a), (b));  \
mpi-genesis~           93     } while(0)
mpi-genesis~           94 
mpi-genesis~           95 /* Copy N limbs from S to D.  */
mpi-genesis~           96 #define MPN_COPY( d, s, n) \
mpi-genesis~           97     do {				\
mpi-genesis~           98 	mpi_size_t _i;			\
mpi-genesis~           99 	for( _i = 0; _i < (n); _i++ )	\
mpi-genesis~          100 	    (d)[_i] = (s)[_i];		\
mpi-genesis~          101     } while(0)
mpi-genesis~          102 
mpi-genesis~          103 #define MPN_COPY_INCR( d, s, n) 	\
mpi-genesis~          104     do {				\
mpi-genesis~          105 	mpi_size_t _i;			\
mpi-genesis~          106 	for( _i = 0; _i < (n); _i++ )	\
mpi-genesis~          107 	    (d)[_i] = (d)[_i];		\
mpi-genesis~          108     } while (0)
mpi-genesis~          109 
mpi-genesis~          110 #define MPN_COPY_DECR( d, s, n ) \
mpi-genesis~          111     do {				\
mpi-genesis~          112 	mpi_size_t _i;			\
mpi-genesis~          113 	for( _i = (n)-1; _i >= 0; _i--) \
mpi-genesis~          114 	   (d)[_i] = (s)[_i];		\
mpi-genesis~          115     } while(0)
mpi-genesis~          116 
mpi-genesis~          117 /* Zero N limbs at D */
mpi-genesis~          118 #define MPN_ZERO(d, n) \
mpi-genesis~          119     do {				  \
mpi-genesis~          120 	int  _i;			  \
mpi-genesis~          121 	for( _i = 0; _i < (n); _i++ )  \
mpi-genesis~          122 	    (d)[_i] = 0;		    \
mpi-genesis~          123     } while (0)
mpi-genesis~          124 
mpi-genesis~          125 #define MPN_NORMALIZE(d, n)  \
mpi-genesis~          126     do {		       \
mpi-genesis~          127 	while( (n) > 0 ) {     \
mpi-genesis~          128 	    if( (d)[(n)-1] ) \
mpi-genesis~          129 		break;	       \
mpi-genesis~          130 	    (n)--;	       \
mpi-genesis~          131 	}		       \
mpi-genesis~          132     } while(0)
mpi-genesis~          133 
mpi-genesis~          134 #define MPN_NORMALIZE_NOT_ZERO(d, n) \
mpi-genesis~          135     do {				    \
mpi-genesis~          136 	for(;;) {			    \
mpi-genesis~          137 	    if( (d)[(n)-1] )		    \
mpi-genesis~          138 		break;			    \
mpi-genesis~          139 	    (n)--;			    \
mpi-genesis~          140 	}				    \
mpi-genesis~          141     } while(0)
mpi-genesis~          142 
mpi-genesis~          143 #define MPN_MUL_N_RECURSE(prodp, up, vp, size, tspace) \
mpi-genesis~          144     do {						\
mpi-genesis~          145 	if( (size) < KARATSUBA_THRESHOLD )		\
mpi-genesis~          146 	    mul_n_basecase (prodp, up, vp, size);	\
mpi-genesis~          147 	else						\
mpi-genesis~          148 	    mul_n (prodp, up, vp, size, tspace);	\
mpi-genesis~          149     } while (0);
mpi-genesis~          150 
mpi-genesis~          151 
mpi-genesis~          152 /* Divide the two-limb number in (NH,,NL) by D, with DI being the largest
mpi-genesis~          153  * limb not larger than (2**(2*BITS_PER_MP_LIMB))/D - (2**BITS_PER_MP_LIMB).
mpi-genesis~          154  * If this would yield overflow, DI should be the largest possible number
mpi-genesis~          155  * (i.e., only ones).  For correct operation, the most significant bit of D
mpi-genesis~          156  * has to be set.  Put the quotient in Q and the remainder in R.
mpi-genesis~          157  */
mpi-genesis~          158 #define UDIV_QRNND_PREINV(q, r, nh, nl, d, di) \
mpi-genesis~          159     do {							    \
mpi-genesis~          160 	mpi_limb_t _q, _ql, _r; 				    \
mpi-genesis~          161 	mpi_limb_t _xh, _xl;					    \
mpi-genesis~          162 	umul_ppmm (_q, _ql, (nh), (di));			    \
mpi-genesis~          163 	_q += (nh);	/* DI is 2**BITS_PER_MPI_LIMB too small */  \
mpi-genesis~          164 	umul_ppmm (_xh, _xl, _q, (d));				    \
mpi-genesis~          165 	sub_ddmmss (_xh, _r, (nh), (nl), _xh, _xl);		    \
mpi-genesis~          166 	if( _xh ) {						    \
mpi-genesis~          167 	    sub_ddmmss (_xh, _r, _xh, _r, 0, (d));		    \
mpi-genesis~          168 	    _q++;						    \
mpi-genesis~          169 	    if( _xh) {						    \
mpi-genesis~          170 		sub_ddmmss (_xh, _r, _xh, _r, 0, (d));		    \
mpi-genesis~          171 		_q++;						    \
mpi-genesis~          172 	    }							    \
mpi-genesis~          173 	}							    \
mpi-genesis~          174 	if( _r >= (d) ) {					    \
mpi-genesis~          175 	    _r -= (d);						    \
mpi-genesis~          176 	    _q++;						    \
mpi-genesis~          177 	}							    \
mpi-genesis~          178 	(r) = _r;						    \
mpi-genesis~          179 	(q) = _q;						    \
mpi-genesis~          180     } while (0)
mpi-genesis~          181 
mpi-genesis~          182 
mpi-genesis~          183 /*-- mpiutil.c --*/
mpi-genesis~          184 #ifdef M_DEBUG
mpi-genesis~          185 #define mpi_alloc_limb_space(n,f)  mpi_debug_alloc_limb_space((n),(f), M_DBGINFO( __LINE__ ) )
mpi-genesis~          186 #define mpi_free_limb_space(n)  mpi_debug_free_limb_space((n),  M_DBGINFO( __LINE__ ) )
mpi-genesis~          187   mpi_ptr_t mpi_debug_alloc_limb_space( unsigned nlimbs, int sec, const char *info  );
mpi-genesis~          188   void mpi_debug_free_limb_space( mpi_ptr_t a, const char *info );
mpi-genesis~          189 #else
mpi-genesis~          190   mpi_ptr_t mpi_alloc_limb_space( unsigned nlimbs, int sec );
mpi-genesis~          191   void mpi_free_limb_space( mpi_ptr_t a );
mpi-genesis~          192 #endif
mpi-genesis~          193 void mpi_assign_limb_space( MPI a, mpi_ptr_t ap, unsigned nlimbs );
mpi-genesis~          194 
mpi-genesis~          195 /*-- mpi-bit.c --*/
mpi-genesis~          196 void mpi_rshift_limbs( MPI a, unsigned int count );
mpi-genesis~          197 void mpi_lshift_limbs( MPI a, unsigned int count );
mpi-genesis~          198 
mpi-genesis~          199 
mpi-genesis~          200 /*-- mpihelp-add.c --*/
mpi-genesis~          201 mpi_limb_t mpihelp_add_1(mpi_ptr_t res_ptr,  mpi_ptr_t s1_ptr,
mpi-genesis~          202 			 mpi_size_t s1_size, mpi_limb_t s2_limb );
mpi-genesis~          203 mpi_limb_t mpihelp_add_n( mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr,
mpi-genesis~          204 			  mpi_ptr_t s2_ptr,  mpi_size_t size);
mpi-genesis~          205 mpi_limb_t mpihelp_add(mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr, mpi_size_t s1_size,
mpi-genesis~          206 		       mpi_ptr_t s2_ptr, mpi_size_t s2_size);
mpi-genesis~          207 
mpi-genesis~          208 /*-- mpihelp-sub.c --*/
mpi-genesis~          209 mpi_limb_t mpihelp_sub_1( mpi_ptr_t res_ptr,  mpi_ptr_t s1_ptr,
mpi-genesis~          210 			  mpi_size_t s1_size, mpi_limb_t s2_limb );
mpi-genesis~          211 mpi_limb_t mpihelp_sub_n( mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr,
mpi-genesis~          212 			  mpi_ptr_t s2_ptr, mpi_size_t size);
mpi-genesis~          213 mpi_limb_t mpihelp_sub(mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr, mpi_size_t s1_size,
mpi-genesis~          214 		       mpi_ptr_t s2_ptr, mpi_size_t s2_size);
mpi-genesis~          215 
mpi-genesis~          216 /*-- mpihelp-cmp.c --*/
mpi-genesis~          217 int mpihelp_cmp( mpi_ptr_t op1_ptr, mpi_ptr_t op2_ptr, mpi_size_t size );
mpi-genesis~          218 
mpi-genesis~          219 /*-- mpihelp-mul.c --*/
mpi-genesis~          220 
mpi-genesis~          221 struct karatsuba_ctx {
mpi-genesis~          222     struct karatsuba_ctx *next;
mpi-genesis~          223     mpi_ptr_t tspace;
mpi-genesis~          224     mpi_size_t tspace_size;
mpi-genesis~          225     mpi_ptr_t tp;
mpi-genesis~          226     mpi_size_t tp_size;
mpi-genesis~          227 };
mpi-genesis~          228 
mpi-genesis~          229 void mpihelp_release_karatsuba_ctx( struct karatsuba_ctx *ctx );
mpi-genesis~          230 
mpi-genesis~          231 mpi_limb_t mpihelp_addmul_1( mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr,
mpi-genesis~          232 			     mpi_size_t s1_size, mpi_limb_t s2_limb);
mpi-genesis~          233 mpi_limb_t mpihelp_submul_1( mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr,
mpi-genesis~          234 			     mpi_size_t s1_size, mpi_limb_t s2_limb);
mpi-genesis~          235 void mpihelp_mul_n( mpi_ptr_t prodp, mpi_ptr_t up, mpi_ptr_t vp,
mpi-genesis~          236 						   mpi_size_t size);
mpi-genesis~          237 mpi_limb_t mpihelp_mul( mpi_ptr_t prodp, mpi_ptr_t up, mpi_size_t usize,
mpi-genesis~          238 					 mpi_ptr_t vp, mpi_size_t vsize);
mpi-genesis~          239 void mpih_sqr_n_basecase( mpi_ptr_t prodp, mpi_ptr_t up, mpi_size_t size );
mpi-genesis~          240 void mpih_sqr_n( mpi_ptr_t prodp, mpi_ptr_t up, mpi_size_t size,
mpi-genesis~          241 						mpi_ptr_t tspace);
mpi-genesis~          242 
mpi-genesis~          243 void mpihelp_mul_karatsuba_case( mpi_ptr_t prodp,
mpi-genesis~          244 				 mpi_ptr_t up, mpi_size_t usize,
mpi-genesis~          245 				 mpi_ptr_t vp, mpi_size_t vsize,
mpi-genesis~          246 				 struct karatsuba_ctx *ctx );
mpi-genesis~          247 
mpi-genesis~          248 
mpi-genesis~          249 /*-- mpihelp-mul_1.c (or xxx/cpu/ *.S) --*/
mpi-genesis~          250 mpi_limb_t mpihelp_mul_1( mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr,
mpi-genesis~          251 			  mpi_size_t s1_size, mpi_limb_t s2_limb);
mpi-genesis~          252 
mpi-genesis~          253 /*-- mpihelp-div.c --*/
mpi-genesis~          254 mpi_limb_t mpihelp_mod_1(mpi_ptr_t dividend_ptr, mpi_size_t dividend_size,
mpi-genesis~          255 						 mpi_limb_t divisor_limb);
mpi-genesis~          256 mpi_limb_t mpihelp_divrem( mpi_ptr_t qp, mpi_size_t qextra_limbs,
mpi-genesis~          257 			   mpi_ptr_t np, mpi_size_t nsize,
mpi-genesis~          258 			   mpi_ptr_t dp, mpi_size_t dsize);
mpi-genesis~          259 mpi_limb_t mpihelp_divmod_1( mpi_ptr_t quot_ptr,
mpi-genesis~          260 			     mpi_ptr_t dividend_ptr, mpi_size_t dividend_size,
mpi-genesis~          261 			     mpi_limb_t divisor_limb);
mpi-genesis~          262 
mpi-genesis~          263 /*-- mpihelp-shift.c --*/
mpi-genesis~          264 mpi_limb_t mpihelp_lshift( mpi_ptr_t wp, mpi_ptr_t up, mpi_size_t usize,
mpi-genesis~          265 							   unsigned cnt);
mpi-genesis~          266 mpi_limb_t mpihelp_rshift( mpi_ptr_t wp, mpi_ptr_t up, mpi_size_t usize,
mpi-genesis~          267 							   unsigned cnt);
mpi-genesis~          268 
mpi-genesis~          269 
mpi-genesis~          270 /* Define stuff for longlong.h.  */
mpi-genesis~          271 #define W_TYPE_SIZE BITS_PER_MPI_LIMB
mpi-genesis~          272   typedef mpi_limb_t   UWtype;
mpi-genesis~          273   typedef unsigned int UHWtype;
mpi-genesis~          274 #if defined (__GNUC__)
mpi-genesis~          275   typedef unsigned int UQItype	  __attribute__ ((mode (QI)));
mpi-genesis~          276   typedef	   int SItype	  __attribute__ ((mode (SI)));
mpi-genesis~          277   typedef unsigned int USItype	  __attribute__ ((mode (SI)));
mpi-genesis~          278   typedef	   int DItype	  __attribute__ ((mode (DI)));
mpi-genesis~          279   typedef unsigned int UDItype	  __attribute__ ((mode (DI)));
mpi-genesis~          280 #else
mpi-genesis~          281   typedef unsigned char UQItype;
mpi-genesis~          282   typedef	   long SItype;
mpi-genesis~          283   typedef unsigned long USItype;
mpi-genesis~          284 #endif
mpi-genesis~          285 
mpi-genesis~          286 #ifdef __GNUC__
mpi-genesis~          287 #include "mpi-inline.h"
mpi-genesis~          288 #endif
mpi-genesis~          289 
mpi-genesis~          290 #endif /*G10_MPI_INTERNAL_H*/