change some parameters in the library to meet the needs of the website httpbin.org
Fork of MiniTLS-GPL by
math/mont/fp_montgomery_calc_normalization.c@5:95f70ebfe61f, 2015-02-06 (annotated)
- Committer:
- shiyilei
- Date:
- Fri Feb 06 06:17:33 2015 +0000
- Revision:
- 5:95f70ebfe61f
- Parent:
- 0:35aa5be3b78d
change some parameters in the library to meet the needs of httpbin.org
Who changed what in which revision?
User | Revision | Line number | New contents of line |
---|---|---|---|
MiniTLS | 0:35aa5be3b78d | 1 | /* TomsFastMath, a fast ISO C bignum library. |
MiniTLS | 0:35aa5be3b78d | 2 | * |
MiniTLS | 0:35aa5be3b78d | 3 | * This project is meant to fill in where LibTomMath |
MiniTLS | 0:35aa5be3b78d | 4 | * falls short. That is speed ;-) |
MiniTLS | 0:35aa5be3b78d | 5 | * |
MiniTLS | 0:35aa5be3b78d | 6 | * This project is public domain and free for all purposes. |
MiniTLS | 0:35aa5be3b78d | 7 | * |
MiniTLS | 0:35aa5be3b78d | 8 | * Tom St Denis, tomstdenis@gmail.com |
MiniTLS | 0:35aa5be3b78d | 9 | */ |
MiniTLS | 0:35aa5be3b78d | 10 | #include <tfm.h> |
MiniTLS | 0:35aa5be3b78d | 11 | |
MiniTLS | 0:35aa5be3b78d | 12 | /* computes a = B**n mod b without division or multiplication useful for |
MiniTLS | 0:35aa5be3b78d | 13 | * normalizing numbers in a Montgomery system. |
MiniTLS | 0:35aa5be3b78d | 14 | */ |
MiniTLS | 0:35aa5be3b78d | 15 | void fp_montgomery_calc_normalization(fp_int *a, fp_int *b) |
MiniTLS | 0:35aa5be3b78d | 16 | { |
MiniTLS | 0:35aa5be3b78d | 17 | int x, bits; |
MiniTLS | 0:35aa5be3b78d | 18 | |
MiniTLS | 0:35aa5be3b78d | 19 | /* how many bits of last digit does b use */ |
MiniTLS | 0:35aa5be3b78d | 20 | bits = fp_count_bits (b) % DIGIT_BIT; |
MiniTLS | 0:35aa5be3b78d | 21 | if (!bits) bits = DIGIT_BIT; |
MiniTLS | 0:35aa5be3b78d | 22 | |
MiniTLS | 0:35aa5be3b78d | 23 | /* compute A = B^(n-1) * 2^(bits-1) */ |
MiniTLS | 0:35aa5be3b78d | 24 | if (b->used > 1) { |
MiniTLS | 0:35aa5be3b78d | 25 | fp_2expt (a, (b->used - 1) * DIGIT_BIT + bits - 1); |
MiniTLS | 0:35aa5be3b78d | 26 | } else { |
MiniTLS | 0:35aa5be3b78d | 27 | fp_set(a, 1); |
MiniTLS | 0:35aa5be3b78d | 28 | bits = 1; |
MiniTLS | 0:35aa5be3b78d | 29 | } |
MiniTLS | 0:35aa5be3b78d | 30 | |
MiniTLS | 0:35aa5be3b78d | 31 | /* now compute C = A * B mod b */ |
MiniTLS | 0:35aa5be3b78d | 32 | for (x = bits - 1; x < (int)DIGIT_BIT; x++) { |
MiniTLS | 0:35aa5be3b78d | 33 | fp_mul_2 (a, a); |
MiniTLS | 0:35aa5be3b78d | 34 | if (fp_cmp_mag (a, b) != FP_LT) { |
MiniTLS | 0:35aa5be3b78d | 35 | s_fp_sub (a, b, a); |
MiniTLS | 0:35aa5be3b78d | 36 | } |
MiniTLS | 0:35aa5be3b78d | 37 | } |
MiniTLS | 0:35aa5be3b78d | 38 | } |
MiniTLS | 0:35aa5be3b78d | 39 | |
MiniTLS | 0:35aa5be3b78d | 40 | |
MiniTLS | 0:35aa5be3b78d | 41 | /* $Source: /cvs/libtom/tomsfastmath/src/mont/fp_montgomery_calc_normalization.c,v $ */ |
MiniTLS | 0:35aa5be3b78d | 42 | /* $Revision: 1.1 $ */ |
MiniTLS | 0:35aa5be3b78d | 43 | /* $Date: 2006/12/31 21:25:53 $ */ |