| OLD | NEW |
| (Empty) |
| 1 /* crypto/bn/bn_lib.c */ | |
| 2 /* Copyright (C) 1995-1998 Eric Young (eay@cryptsoft.com) | |
| 3 * All rights reserved. | |
| 4 * | |
| 5 * This package is an SSL implementation written | |
| 6 * by Eric Young (eay@cryptsoft.com). | |
| 7 * The implementation was written so as to conform with Netscapes SSL. | |
| 8 * | |
| 9 * This library is free for commercial and non-commercial use as long as | |
| 10 * the following conditions are aheared to. The following conditions | |
| 11 * apply to all code found in this distribution, be it the RC4, RSA, | |
| 12 * lhash, DES, etc., code; not just the SSL code. The SSL documentation | |
| 13 * included with this distribution is covered by the same copyright terms | |
| 14 * except that the holder is Tim Hudson (tjh@cryptsoft.com). | |
| 15 * | |
| 16 * Copyright remains Eric Young's, and as such any Copyright notices in | |
| 17 * the code are not to be removed. | |
| 18 * If this package is used in a product, Eric Young should be given attribution | |
| 19 * as the author of the parts of the library used. | |
| 20 * This can be in the form of a textual message at program startup or | |
| 21 * in documentation (online or textual) provided with the package. | |
| 22 * | |
| 23 * Redistribution and use in source and binary forms, with or without | |
| 24 * modification, are permitted provided that the following conditions | |
| 25 * are met: | |
| 26 * 1. Redistributions of source code must retain the copyright | |
| 27 * notice, this list of conditions and the following disclaimer. | |
| 28 * 2. Redistributions in binary form must reproduce the above copyright | |
| 29 * notice, this list of conditions and the following disclaimer in the | |
| 30 * documentation and/or other materials provided with the distribution. | |
| 31 * 3. All advertising materials mentioning features or use of this software | |
| 32 * must display the following acknowledgement: | |
| 33 * "This product includes cryptographic software written by | |
| 34 * Eric Young (eay@cryptsoft.com)" | |
| 35 * The word 'cryptographic' can be left out if the rouines from the library | |
| 36 * being used are not cryptographic related :-). | |
| 37 * 4. If you include any Windows specific code (or a derivative thereof) from | |
| 38 * the apps directory (application code) you must include an acknowledgement: | |
| 39 * "This product includes software written by Tim Hudson (tjh@cryptsoft.com)" | |
| 40 * | |
| 41 * THIS SOFTWARE IS PROVIDED BY ERIC YOUNG ``AS IS'' AND | |
| 42 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE | |
| 43 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE | |
| 44 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE | |
| 45 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL | |
| 46 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS | |
| 47 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) | |
| 48 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT | |
| 49 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY | |
| 50 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF | |
| 51 * SUCH DAMAGE. | |
| 52 * | |
| 53 * The licence and distribution terms for any publically available version or | |
| 54 * derivative of this code cannot be changed. i.e. this code cannot simply be | |
| 55 * copied and put under another distribution licence | |
| 56 * [including the GNU Public Licence.] | |
| 57 */ | |
| 58 | |
| 59 #ifndef BN_DEBUG | |
| 60 # undef NDEBUG /* avoid conflicting definitions */ | |
| 61 # define NDEBUG | |
| 62 #endif | |
| 63 | |
| 64 #include <assert.h> | |
| 65 #include <limits.h> | |
| 66 #include <stdio.h> | |
| 67 #include "cryptlib.h" | |
| 68 #include "bn_lcl.h" | |
| 69 | |
| 70 const char BN_version[]="Big Number" OPENSSL_VERSION_PTEXT; | |
| 71 | |
| 72 /* This stuff appears to be completely unused, so is deprecated */ | |
| 73 #ifndef OPENSSL_NO_DEPRECATED | |
| 74 /* For a 32 bit machine | |
| 75 * 2 - 4 == 128 | |
| 76 * 3 - 8 == 256 | |
| 77 * 4 - 16 == 512 | |
| 78 * 5 - 32 == 1024 | |
| 79 * 6 - 64 == 2048 | |
| 80 * 7 - 128 == 4096 | |
| 81 * 8 - 256 == 8192 | |
| 82 */ | |
| 83 static int bn_limit_bits=0; | |
| 84 static int bn_limit_num=8; /* (1<<bn_limit_bits) */ | |
| 85 static int bn_limit_bits_low=0; | |
| 86 static int bn_limit_num_low=8; /* (1<<bn_limit_bits_low) */ | |
| 87 static int bn_limit_bits_high=0; | |
| 88 static int bn_limit_num_high=8; /* (1<<bn_limit_bits_high) */ | |
| 89 static int bn_limit_bits_mont=0; | |
| 90 static int bn_limit_num_mont=8; /* (1<<bn_limit_bits_mont) */ | |
| 91 | |
| 92 void BN_set_params(int mult, int high, int low, int mont) | |
| 93 { | |
| 94 if (mult >= 0) | |
| 95 { | |
| 96 if (mult > (int)(sizeof(int)*8)-1) | |
| 97 mult=sizeof(int)*8-1; | |
| 98 bn_limit_bits=mult; | |
| 99 bn_limit_num=1<<mult; | |
| 100 } | |
| 101 if (high >= 0) | |
| 102 { | |
| 103 if (high > (int)(sizeof(int)*8)-1) | |
| 104 high=sizeof(int)*8-1; | |
| 105 bn_limit_bits_high=high; | |
| 106 bn_limit_num_high=1<<high; | |
| 107 } | |
| 108 if (low >= 0) | |
| 109 { | |
| 110 if (low > (int)(sizeof(int)*8)-1) | |
| 111 low=sizeof(int)*8-1; | |
| 112 bn_limit_bits_low=low; | |
| 113 bn_limit_num_low=1<<low; | |
| 114 } | |
| 115 if (mont >= 0) | |
| 116 { | |
| 117 if (mont > (int)(sizeof(int)*8)-1) | |
| 118 mont=sizeof(int)*8-1; | |
| 119 bn_limit_bits_mont=mont; | |
| 120 bn_limit_num_mont=1<<mont; | |
| 121 } | |
| 122 } | |
| 123 | |
| 124 int BN_get_params(int which) | |
| 125 { | |
| 126 if (which == 0) return(bn_limit_bits); | |
| 127 else if (which == 1) return(bn_limit_bits_high); | |
| 128 else if (which == 2) return(bn_limit_bits_low); | |
| 129 else if (which == 3) return(bn_limit_bits_mont); | |
| 130 else return(0); | |
| 131 } | |
| 132 #endif | |
| 133 | |
| 134 const BIGNUM *BN_value_one(void) | |
| 135 { | |
| 136 static const BN_ULONG data_one=1L; | |
| 137 static const BIGNUM const_one={(BN_ULONG *)&data_one,1,1,0,BN_FLG_STATIC
_DATA}; | |
| 138 | |
| 139 return(&const_one); | |
| 140 } | |
| 141 | |
| 142 int BN_num_bits_word(BN_ULONG l) | |
| 143 { | |
| 144 static const unsigned char bits[256]={ | |
| 145 0,1,2,2,3,3,3,3,4,4,4,4,4,4,4,4, | |
| 146 5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5, | |
| 147 6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6, | |
| 148 6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6, | |
| 149 7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7, | |
| 150 7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7, | |
| 151 7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7, | |
| 152 7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7, | |
| 153 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8, | |
| 154 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8, | |
| 155 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8, | |
| 156 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8, | |
| 157 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8, | |
| 158 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8, | |
| 159 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8, | |
| 160 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8, | |
| 161 }; | |
| 162 | |
| 163 #if defined(SIXTY_FOUR_BIT_LONG) | |
| 164 if (l & 0xffffffff00000000L) | |
| 165 { | |
| 166 if (l & 0xffff000000000000L) | |
| 167 { | |
| 168 if (l & 0xff00000000000000L) | |
| 169 { | |
| 170 return(bits[(int)(l>>56)]+56); | |
| 171 } | |
| 172 else return(bits[(int)(l>>48)]+48); | |
| 173 } | |
| 174 else | |
| 175 { | |
| 176 if (l & 0x0000ff0000000000L) | |
| 177 { | |
| 178 return(bits[(int)(l>>40)]+40); | |
| 179 } | |
| 180 else return(bits[(int)(l>>32)]+32); | |
| 181 } | |
| 182 } | |
| 183 else | |
| 184 #else | |
| 185 #ifdef SIXTY_FOUR_BIT | |
| 186 if (l & 0xffffffff00000000LL) | |
| 187 { | |
| 188 if (l & 0xffff000000000000LL) | |
| 189 { | |
| 190 if (l & 0xff00000000000000LL) | |
| 191 { | |
| 192 return(bits[(int)(l>>56)]+56); | |
| 193 } | |
| 194 else return(bits[(int)(l>>48)]+48); | |
| 195 } | |
| 196 else | |
| 197 { | |
| 198 if (l & 0x0000ff0000000000LL) | |
| 199 { | |
| 200 return(bits[(int)(l>>40)]+40); | |
| 201 } | |
| 202 else return(bits[(int)(l>>32)]+32); | |
| 203 } | |
| 204 } | |
| 205 else | |
| 206 #endif | |
| 207 #endif | |
| 208 { | |
| 209 #if defined(THIRTY_TWO_BIT) || defined(SIXTY_FOUR_BIT) || defined(SIXTY_FOUR_BIT
_LONG) | |
| 210 if (l & 0xffff0000L) | |
| 211 { | |
| 212 if (l & 0xff000000L) | |
| 213 return(bits[(int)(l>>24L)]+24); | |
| 214 else return(bits[(int)(l>>16L)]+16); | |
| 215 } | |
| 216 else | |
| 217 #endif | |
| 218 { | |
| 219 #if defined(THIRTY_TWO_BIT) || defined(SIXTY_FOUR_BIT) || defined(SIXTY_FOUR_BIT
_LONG) | |
| 220 if (l & 0xff00L) | |
| 221 return(bits[(int)(l>>8)]+8); | |
| 222 else | |
| 223 #endif | |
| 224 return(bits[(int)(l )] ); | |
| 225 } | |
| 226 } | |
| 227 } | |
| 228 | |
| 229 int BN_num_bits(const BIGNUM *a) | |
| 230 { | |
| 231 int i = a->top - 1; | |
| 232 bn_check_top(a); | |
| 233 | |
| 234 if (BN_is_zero(a)) return 0; | |
| 235 return ((i*BN_BITS2) + BN_num_bits_word(a->d[i])); | |
| 236 } | |
| 237 | |
| 238 void BN_clear_free(BIGNUM *a) | |
| 239 { | |
| 240 int i; | |
| 241 | |
| 242 if (a == NULL) return; | |
| 243 bn_check_top(a); | |
| 244 if (a->d != NULL) | |
| 245 { | |
| 246 OPENSSL_cleanse(a->d,a->dmax*sizeof(a->d[0])); | |
| 247 if (!(BN_get_flags(a,BN_FLG_STATIC_DATA))) | |
| 248 OPENSSL_free(a->d); | |
| 249 } | |
| 250 i=BN_get_flags(a,BN_FLG_MALLOCED); | |
| 251 OPENSSL_cleanse(a,sizeof(BIGNUM)); | |
| 252 if (i) | |
| 253 OPENSSL_free(a); | |
| 254 } | |
| 255 | |
| 256 void BN_free(BIGNUM *a) | |
| 257 { | |
| 258 if (a == NULL) return; | |
| 259 bn_check_top(a); | |
| 260 if ((a->d != NULL) && !(BN_get_flags(a,BN_FLG_STATIC_DATA))) | |
| 261 OPENSSL_free(a->d); | |
| 262 if (a->flags & BN_FLG_MALLOCED) | |
| 263 OPENSSL_free(a); | |
| 264 else | |
| 265 { | |
| 266 #ifndef OPENSSL_NO_DEPRECATED | |
| 267 a->flags|=BN_FLG_FREE; | |
| 268 #endif | |
| 269 a->d = NULL; | |
| 270 } | |
| 271 } | |
| 272 | |
| 273 void BN_init(BIGNUM *a) | |
| 274 { | |
| 275 memset(a,0,sizeof(BIGNUM)); | |
| 276 bn_check_top(a); | |
| 277 } | |
| 278 | |
| 279 BIGNUM *BN_new(void) | |
| 280 { | |
| 281 BIGNUM *ret; | |
| 282 | |
| 283 if ((ret=(BIGNUM *)OPENSSL_malloc(sizeof(BIGNUM))) == NULL) | |
| 284 { | |
| 285 BNerr(BN_F_BN_NEW,ERR_R_MALLOC_FAILURE); | |
| 286 return(NULL); | |
| 287 } | |
| 288 ret->flags=BN_FLG_MALLOCED; | |
| 289 ret->top=0; | |
| 290 ret->neg=0; | |
| 291 ret->dmax=0; | |
| 292 ret->d=NULL; | |
| 293 bn_check_top(ret); | |
| 294 return(ret); | |
| 295 } | |
| 296 | |
| 297 /* This is used both by bn_expand2() and bn_dup_expand() */ | |
| 298 /* The caller MUST check that words > b->dmax before calling this */ | |
| 299 static BN_ULONG *bn_expand_internal(const BIGNUM *b, int words) | |
| 300 { | |
| 301 BN_ULONG *A,*a = NULL; | |
| 302 const BN_ULONG *B; | |
| 303 int i; | |
| 304 | |
| 305 bn_check_top(b); | |
| 306 | |
| 307 if (words > (INT_MAX/(4*BN_BITS2))) | |
| 308 { | |
| 309 BNerr(BN_F_BN_EXPAND_INTERNAL,BN_R_BIGNUM_TOO_LONG); | |
| 310 return NULL; | |
| 311 } | |
| 312 if (BN_get_flags(b,BN_FLG_STATIC_DATA)) | |
| 313 { | |
| 314 BNerr(BN_F_BN_EXPAND_INTERNAL,BN_R_EXPAND_ON_STATIC_BIGNUM_DATA)
; | |
| 315 return(NULL); | |
| 316 } | |
| 317 a=A=(BN_ULONG *)OPENSSL_malloc(sizeof(BN_ULONG)*words); | |
| 318 if (A == NULL) | |
| 319 { | |
| 320 BNerr(BN_F_BN_EXPAND_INTERNAL,ERR_R_MALLOC_FAILURE); | |
| 321 return(NULL); | |
| 322 } | |
| 323 #if 1 | |
| 324 B=b->d; | |
| 325 /* Check if the previous number needs to be copied */ | |
| 326 if (B != NULL) | |
| 327 { | |
| 328 for (i=b->top>>2; i>0; i--,A+=4,B+=4) | |
| 329 { | |
| 330 /* | |
| 331 * The fact that the loop is unrolled | |
| 332 * 4-wise is a tribute to Intel. It's | |
| 333 * the one that doesn't have enough | |
| 334 * registers to accomodate more data. | |
| 335 * I'd unroll it 8-wise otherwise:-) | |
| 336 * | |
| 337 * <appro@fy.chalmers.se> | |
| 338 */ | |
| 339 BN_ULONG a0,a1,a2,a3; | |
| 340 a0=B[0]; a1=B[1]; a2=B[2]; a3=B[3]; | |
| 341 A[0]=a0; A[1]=a1; A[2]=a2; A[3]=a3; | |
| 342 } | |
| 343 switch (b->top&3) | |
| 344 { | |
| 345 case 3: A[2]=B[2]; | |
| 346 case 2: A[1]=B[1]; | |
| 347 case 1: A[0]=B[0]; | |
| 348 case 0: /* workaround for ultrix cc: without 'case 0', the optim
izer does | |
| 349 * the switch table by doing a=top&3; a--; goto jump_tab
le[a]; | |
| 350 * which fails for top== 0 */ | |
| 351 ; | |
| 352 } | |
| 353 } | |
| 354 | |
| 355 #else | |
| 356 memset(A,0,sizeof(BN_ULONG)*words); | |
| 357 memcpy(A,b->d,sizeof(b->d[0])*b->top); | |
| 358 #endif | |
| 359 | |
| 360 return(a); | |
| 361 } | |
| 362 | |
| 363 /* This is an internal function that can be used instead of bn_expand2() | |
| 364 * when there is a need to copy BIGNUMs instead of only expanding the | |
| 365 * data part, while still expanding them. | |
| 366 * Especially useful when needing to expand BIGNUMs that are declared | |
| 367 * 'const' and should therefore not be changed. | |
| 368 * The reason to use this instead of a BN_dup() followed by a bn_expand2() | |
| 369 * is memory allocation overhead. A BN_dup() followed by a bn_expand2() | |
| 370 * will allocate new memory for the BIGNUM data twice, and free it once, | |
| 371 * while bn_dup_expand() makes sure allocation is made only once. | |
| 372 */ | |
| 373 | |
| 374 #ifndef OPENSSL_NO_DEPRECATED | |
| 375 BIGNUM *bn_dup_expand(const BIGNUM *b, int words) | |
| 376 { | |
| 377 BIGNUM *r = NULL; | |
| 378 | |
| 379 bn_check_top(b); | |
| 380 | |
| 381 /* This function does not work if | |
| 382 * words <= b->dmax && top < words | |
| 383 * because BN_dup() does not preserve 'dmax'! | |
| 384 * (But bn_dup_expand() is not used anywhere yet.) | |
| 385 */ | |
| 386 | |
| 387 if (words > b->dmax) | |
| 388 { | |
| 389 BN_ULONG *a = bn_expand_internal(b, words); | |
| 390 | |
| 391 if (a) | |
| 392 { | |
| 393 r = BN_new(); | |
| 394 if (r) | |
| 395 { | |
| 396 r->top = b->top; | |
| 397 r->dmax = words; | |
| 398 r->neg = b->neg; | |
| 399 r->d = a; | |
| 400 } | |
| 401 else | |
| 402 { | |
| 403 /* r == NULL, BN_new failure */ | |
| 404 OPENSSL_free(a); | |
| 405 } | |
| 406 } | |
| 407 /* If a == NULL, there was an error in allocation in | |
| 408 bn_expand_internal(), and NULL should be returned */ | |
| 409 } | |
| 410 else | |
| 411 { | |
| 412 r = BN_dup(b); | |
| 413 } | |
| 414 | |
| 415 bn_check_top(r); | |
| 416 return r; | |
| 417 } | |
| 418 #endif | |
| 419 | |
| 420 /* This is an internal function that should not be used in applications. | |
| 421 * It ensures that 'b' has enough room for a 'words' word number | |
| 422 * and initialises any unused part of b->d with leading zeros. | |
| 423 * It is mostly used by the various BIGNUM routines. If there is an error, | |
| 424 * NULL is returned. If not, 'b' is returned. */ | |
| 425 | |
| 426 BIGNUM *bn_expand2(BIGNUM *b, int words) | |
| 427 { | |
| 428 bn_check_top(b); | |
| 429 | |
| 430 if (words > b->dmax) | |
| 431 { | |
| 432 BN_ULONG *a = bn_expand_internal(b, words); | |
| 433 if(!a) return NULL; | |
| 434 if(b->d) OPENSSL_free(b->d); | |
| 435 b->d=a; | |
| 436 b->dmax=words; | |
| 437 } | |
| 438 | |
| 439 /* None of this should be necessary because of what b->top means! */ | |
| 440 #if 0 | |
| 441 /* NB: bn_wexpand() calls this only if the BIGNUM really has to grow */ | |
| 442 if (b->top < b->dmax) | |
| 443 { | |
| 444 int i; | |
| 445 BN_ULONG *A = &(b->d[b->top]); | |
| 446 for (i=(b->dmax - b->top)>>3; i>0; i--,A+=8) | |
| 447 { | |
| 448 A[0]=0; A[1]=0; A[2]=0; A[3]=0; | |
| 449 A[4]=0; A[5]=0; A[6]=0; A[7]=0; | |
| 450 } | |
| 451 for (i=(b->dmax - b->top)&7; i>0; i--,A++) | |
| 452 A[0]=0; | |
| 453 assert(A == &(b->d[b->dmax])); | |
| 454 } | |
| 455 #endif | |
| 456 bn_check_top(b); | |
| 457 return b; | |
| 458 } | |
| 459 | |
| 460 BIGNUM *BN_dup(const BIGNUM *a) | |
| 461 { | |
| 462 BIGNUM *t; | |
| 463 | |
| 464 if (a == NULL) return NULL; | |
| 465 bn_check_top(a); | |
| 466 | |
| 467 t = BN_new(); | |
| 468 if (t == NULL) return NULL; | |
| 469 if(!BN_copy(t, a)) | |
| 470 { | |
| 471 BN_free(t); | |
| 472 return NULL; | |
| 473 } | |
| 474 bn_check_top(t); | |
| 475 return t; | |
| 476 } | |
| 477 | |
| 478 BIGNUM *BN_copy(BIGNUM *a, const BIGNUM *b) | |
| 479 { | |
| 480 int i; | |
| 481 BN_ULONG *A; | |
| 482 const BN_ULONG *B; | |
| 483 | |
| 484 bn_check_top(b); | |
| 485 | |
| 486 if (a == b) return(a); | |
| 487 if (bn_wexpand(a,b->top) == NULL) return(NULL); | |
| 488 | |
| 489 #if 1 | |
| 490 A=a->d; | |
| 491 B=b->d; | |
| 492 for (i=b->top>>2; i>0; i--,A+=4,B+=4) | |
| 493 { | |
| 494 BN_ULONG a0,a1,a2,a3; | |
| 495 a0=B[0]; a1=B[1]; a2=B[2]; a3=B[3]; | |
| 496 A[0]=a0; A[1]=a1; A[2]=a2; A[3]=a3; | |
| 497 } | |
| 498 switch (b->top&3) | |
| 499 { | |
| 500 case 3: A[2]=B[2]; | |
| 501 case 2: A[1]=B[1]; | |
| 502 case 1: A[0]=B[0]; | |
| 503 case 0: ; /* ultrix cc workaround, see comments in bn_expand_int
ernal */ | |
| 504 } | |
| 505 #else | |
| 506 memcpy(a->d,b->d,sizeof(b->d[0])*b->top); | |
| 507 #endif | |
| 508 | |
| 509 a->top=b->top; | |
| 510 a->neg=b->neg; | |
| 511 bn_check_top(a); | |
| 512 return(a); | |
| 513 } | |
| 514 | |
| 515 void BN_swap(BIGNUM *a, BIGNUM *b) | |
| 516 { | |
| 517 int flags_old_a, flags_old_b; | |
| 518 BN_ULONG *tmp_d; | |
| 519 int tmp_top, tmp_dmax, tmp_neg; | |
| 520 | |
| 521 bn_check_top(a); | |
| 522 bn_check_top(b); | |
| 523 | |
| 524 flags_old_a = a->flags; | |
| 525 flags_old_b = b->flags; | |
| 526 | |
| 527 tmp_d = a->d; | |
| 528 tmp_top = a->top; | |
| 529 tmp_dmax = a->dmax; | |
| 530 tmp_neg = a->neg; | |
| 531 | |
| 532 a->d = b->d; | |
| 533 a->top = b->top; | |
| 534 a->dmax = b->dmax; | |
| 535 a->neg = b->neg; | |
| 536 | |
| 537 b->d = tmp_d; | |
| 538 b->top = tmp_top; | |
| 539 b->dmax = tmp_dmax; | |
| 540 b->neg = tmp_neg; | |
| 541 | |
| 542 a->flags = (flags_old_a & BN_FLG_MALLOCED) | (flags_old_b & BN_FLG_STATI
C_DATA); | |
| 543 b->flags = (flags_old_b & BN_FLG_MALLOCED) | (flags_old_a & BN_FLG_STATI
C_DATA); | |
| 544 bn_check_top(a); | |
| 545 bn_check_top(b); | |
| 546 } | |
| 547 | |
| 548 void BN_clear(BIGNUM *a) | |
| 549 { | |
| 550 bn_check_top(a); | |
| 551 if (a->d != NULL) | |
| 552 memset(a->d,0,a->dmax*sizeof(a->d[0])); | |
| 553 a->top=0; | |
| 554 a->neg=0; | |
| 555 } | |
| 556 | |
| 557 BN_ULONG BN_get_word(const BIGNUM *a) | |
| 558 { | |
| 559 if (a->top > 1) | |
| 560 return BN_MASK2; | |
| 561 else if (a->top == 1) | |
| 562 return a->d[0]; | |
| 563 /* a->top == 0 */ | |
| 564 return 0; | |
| 565 } | |
| 566 | |
| 567 int BN_set_word(BIGNUM *a, BN_ULONG w) | |
| 568 { | |
| 569 bn_check_top(a); | |
| 570 if (bn_expand(a,(int)sizeof(BN_ULONG)*8) == NULL) return(0); | |
| 571 a->neg = 0; | |
| 572 a->d[0] = w; | |
| 573 a->top = (w ? 1 : 0); | |
| 574 bn_check_top(a); | |
| 575 return(1); | |
| 576 } | |
| 577 | |
| 578 BIGNUM *BN_bin2bn(const unsigned char *s, int len, BIGNUM *ret) | |
| 579 { | |
| 580 unsigned int i,m; | |
| 581 unsigned int n; | |
| 582 BN_ULONG l; | |
| 583 BIGNUM *bn = NULL; | |
| 584 | |
| 585 if (ret == NULL) | |
| 586 ret = bn = BN_new(); | |
| 587 if (ret == NULL) return(NULL); | |
| 588 bn_check_top(ret); | |
| 589 l=0; | |
| 590 n=len; | |
| 591 if (n == 0) | |
| 592 { | |
| 593 ret->top=0; | |
| 594 return(ret); | |
| 595 } | |
| 596 i=((n-1)/BN_BYTES)+1; | |
| 597 m=((n-1)%(BN_BYTES)); | |
| 598 if (bn_wexpand(ret, (int)i) == NULL) | |
| 599 { | |
| 600 if (bn) BN_free(bn); | |
| 601 return NULL; | |
| 602 } | |
| 603 ret->top=i; | |
| 604 ret->neg=0; | |
| 605 while (n--) | |
| 606 { | |
| 607 l=(l<<8L)| *(s++); | |
| 608 if (m-- == 0) | |
| 609 { | |
| 610 ret->d[--i]=l; | |
| 611 l=0; | |
| 612 m=BN_BYTES-1; | |
| 613 } | |
| 614 } | |
| 615 /* need to call this due to clear byte at top if avoiding | |
| 616 * having the top bit set (-ve number) */ | |
| 617 bn_correct_top(ret); | |
| 618 return(ret); | |
| 619 } | |
| 620 | |
| 621 /* ignore negative */ | |
| 622 int BN_bn2bin(const BIGNUM *a, unsigned char *to) | |
| 623 { | |
| 624 int n,i; | |
| 625 BN_ULONG l; | |
| 626 | |
| 627 bn_check_top(a); | |
| 628 n=i=BN_num_bytes(a); | |
| 629 while (i--) | |
| 630 { | |
| 631 l=a->d[i/BN_BYTES]; | |
| 632 *(to++)=(unsigned char)(l>>(8*(i%BN_BYTES)))&0xff; | |
| 633 } | |
| 634 return(n); | |
| 635 } | |
| 636 | |
| 637 int BN_ucmp(const BIGNUM *a, const BIGNUM *b) | |
| 638 { | |
| 639 int i; | |
| 640 BN_ULONG t1,t2,*ap,*bp; | |
| 641 | |
| 642 bn_check_top(a); | |
| 643 bn_check_top(b); | |
| 644 | |
| 645 i=a->top-b->top; | |
| 646 if (i != 0) return(i); | |
| 647 ap=a->d; | |
| 648 bp=b->d; | |
| 649 for (i=a->top-1; i>=0; i--) | |
| 650 { | |
| 651 t1= ap[i]; | |
| 652 t2= bp[i]; | |
| 653 if (t1 != t2) | |
| 654 return((t1 > t2) ? 1 : -1); | |
| 655 } | |
| 656 return(0); | |
| 657 } | |
| 658 | |
| 659 int BN_cmp(const BIGNUM *a, const BIGNUM *b) | |
| 660 { | |
| 661 int i; | |
| 662 int gt,lt; | |
| 663 BN_ULONG t1,t2; | |
| 664 | |
| 665 if ((a == NULL) || (b == NULL)) | |
| 666 { | |
| 667 if (a != NULL) | |
| 668 return(-1); | |
| 669 else if (b != NULL) | |
| 670 return(1); | |
| 671 else | |
| 672 return(0); | |
| 673 } | |
| 674 | |
| 675 bn_check_top(a); | |
| 676 bn_check_top(b); | |
| 677 | |
| 678 if (a->neg != b->neg) | |
| 679 { | |
| 680 if (a->neg) | |
| 681 return(-1); | |
| 682 else return(1); | |
| 683 } | |
| 684 if (a->neg == 0) | |
| 685 { gt=1; lt= -1; } | |
| 686 else { gt= -1; lt=1; } | |
| 687 | |
| 688 if (a->top > b->top) return(gt); | |
| 689 if (a->top < b->top) return(lt); | |
| 690 for (i=a->top-1; i>=0; i--) | |
| 691 { | |
| 692 t1=a->d[i]; | |
| 693 t2=b->d[i]; | |
| 694 if (t1 > t2) return(gt); | |
| 695 if (t1 < t2) return(lt); | |
| 696 } | |
| 697 return(0); | |
| 698 } | |
| 699 | |
| 700 int BN_set_bit(BIGNUM *a, int n) | |
| 701 { | |
| 702 int i,j,k; | |
| 703 | |
| 704 if (n < 0) | |
| 705 return 0; | |
| 706 | |
| 707 i=n/BN_BITS2; | |
| 708 j=n%BN_BITS2; | |
| 709 if (a->top <= i) | |
| 710 { | |
| 711 if (bn_wexpand(a,i+1) == NULL) return(0); | |
| 712 for(k=a->top; k<i+1; k++) | |
| 713 a->d[k]=0; | |
| 714 a->top=i+1; | |
| 715 } | |
| 716 | |
| 717 a->d[i]|=(((BN_ULONG)1)<<j); | |
| 718 bn_check_top(a); | |
| 719 return(1); | |
| 720 } | |
| 721 | |
| 722 int BN_clear_bit(BIGNUM *a, int n) | |
| 723 { | |
| 724 int i,j; | |
| 725 | |
| 726 bn_check_top(a); | |
| 727 if (n < 0) return 0; | |
| 728 | |
| 729 i=n/BN_BITS2; | |
| 730 j=n%BN_BITS2; | |
| 731 if (a->top <= i) return(0); | |
| 732 | |
| 733 a->d[i]&=(~(((BN_ULONG)1)<<j)); | |
| 734 bn_correct_top(a); | |
| 735 return(1); | |
| 736 } | |
| 737 | |
| 738 int BN_is_bit_set(const BIGNUM *a, int n) | |
| 739 { | |
| 740 int i,j; | |
| 741 | |
| 742 bn_check_top(a); | |
| 743 if (n < 0) return 0; | |
| 744 i=n/BN_BITS2; | |
| 745 j=n%BN_BITS2; | |
| 746 if (a->top <= i) return 0; | |
| 747 return (int)(((a->d[i])>>j)&((BN_ULONG)1)); | |
| 748 } | |
| 749 | |
| 750 int BN_mask_bits(BIGNUM *a, int n) | |
| 751 { | |
| 752 int b,w; | |
| 753 | |
| 754 bn_check_top(a); | |
| 755 if (n < 0) return 0; | |
| 756 | |
| 757 w=n/BN_BITS2; | |
| 758 b=n%BN_BITS2; | |
| 759 if (w >= a->top) return 0; | |
| 760 if (b == 0) | |
| 761 a->top=w; | |
| 762 else | |
| 763 { | |
| 764 a->top=w+1; | |
| 765 a->d[w]&= ~(BN_MASK2<<b); | |
| 766 } | |
| 767 bn_correct_top(a); | |
| 768 return(1); | |
| 769 } | |
| 770 | |
| 771 void BN_set_negative(BIGNUM *a, int b) | |
| 772 { | |
| 773 if (b && !BN_is_zero(a)) | |
| 774 a->neg = 1; | |
| 775 else | |
| 776 a->neg = 0; | |
| 777 } | |
| 778 | |
| 779 int bn_cmp_words(const BN_ULONG *a, const BN_ULONG *b, int n) | |
| 780 { | |
| 781 int i; | |
| 782 BN_ULONG aa,bb; | |
| 783 | |
| 784 aa=a[n-1]; | |
| 785 bb=b[n-1]; | |
| 786 if (aa != bb) return((aa > bb)?1:-1); | |
| 787 for (i=n-2; i>=0; i--) | |
| 788 { | |
| 789 aa=a[i]; | |
| 790 bb=b[i]; | |
| 791 if (aa != bb) return((aa > bb)?1:-1); | |
| 792 } | |
| 793 return(0); | |
| 794 } | |
| 795 | |
| 796 /* Here follows a specialised variants of bn_cmp_words(). It has the | |
| 797 property of performing the operation on arrays of different sizes. | |
| 798 The sizes of those arrays is expressed through cl, which is the | |
| 799 common length ( basicall, min(len(a),len(b)) ), and dl, which is the | |
| 800 delta between the two lengths, calculated as len(a)-len(b). | |
| 801 All lengths are the number of BN_ULONGs... */ | |
| 802 | |
| 803 int bn_cmp_part_words(const BN_ULONG *a, const BN_ULONG *b, | |
| 804 int cl, int dl) | |
| 805 { | |
| 806 int n,i; | |
| 807 n = cl-1; | |
| 808 | |
| 809 if (dl < 0) | |
| 810 { | |
| 811 for (i=dl; i<0; i++) | |
| 812 { | |
| 813 if (b[n-i] != 0) | |
| 814 return -1; /* a < b */ | |
| 815 } | |
| 816 } | |
| 817 if (dl > 0) | |
| 818 { | |
| 819 for (i=dl; i>0; i--) | |
| 820 { | |
| 821 if (a[n+i] != 0) | |
| 822 return 1; /* a > b */ | |
| 823 } | |
| 824 } | |
| 825 return bn_cmp_words(a,b,cl); | |
| 826 } | |
| OLD | NEW |