|
|
1.1 root 1: /*
2: * Copyright(C) 2006 Cameron Rich
3: *
4: * This library is free software; you can redistribute it and/or modify
5: * it under the terms of the GNU Lesser General Public License as published by
6: * the Free Software Foundation; either version 2.1 of the License, or
7: * (at your option) any later version.
8: *
9: * This library is distributed in the hope that it will be useful,
10: * but WITHOUT ANY WARRANTY; without even the implied warranty of
11: * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12: * GNU Lesser General Public License for more details.
13: *
14: * You should have received a copy of the GNU Lesser General Public License
15: * along with this library; if not, write to the Free Software
16: * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17: */
18:
19: #ifndef BIGINT_IMPL_HEADER
20: #define BIGINT_IMPL_HEADER
21:
22: /* Maintain a number of precomputed variables when doing reduction */
23: #define BIGINT_M_OFFSET 0 /**< Normal modulo offset. */
24: #ifdef CONFIG_BIGINT_CRT
25: #define BIGINT_P_OFFSET 1 /**< p modulo offset. */
26: #define BIGINT_Q_OFFSET 2 /**< q module offset. */
27: #define BIGINT_NUM_MODS 3 /**< The number of modulus constants used. */
28: #else
29: #define BIGINT_NUM_MODS 1
30: #endif
31:
32: /* Architecture specific functions for big ints */
33: #ifdef WIN32
34: #define COMP_RADIX 4294967296i64
35: #define COMP_BIG_MSB 0x8000000000000000i64
36: #else
37: #define COMP_RADIX 4294967296ULL /**< Max component + 1 */
38: #define COMP_BIG_MSB 0x8000000000000000ULL /**< (Max dbl comp + 1)/ 2 */
39: #endif
40: #define COMP_BIT_SIZE 32 /**< Number of bits in a component. */
41: #define COMP_BYTE_SIZE 4 /**< Number of bytes in a component. */
42: #define COMP_NUM_NIBBLES 8 /**< Used For diagnostics only. */
43:
44: typedef uint32_t comp; /**< A single precision component. */
45: typedef uint64_t long_comp; /**< A double precision component. */
46: typedef int64_t slong_comp; /**< A signed double precision component. */
47:
48: /**
49: * @struct _bigint
50: * @brief A big integer basic object
51: */
52: struct _bigint
53: {
54: struct _bigint* next; /**< The next bigint in the cache. */
55: short size; /**< The number of components in this bigint. */
56: short max_comps; /**< The heapsize allocated for this bigint */
57: int refs; /**< An internal reference count. */
58: comp* comps; /**< A ptr to the actual component data */
59: };
60:
61: typedef struct _bigint bigint; /**< An alias for _bigint */
62:
63: /**
64: * Maintains the state of the cache, and a number of variables used in
65: * reduction.
66: */
67: typedef struct /**< A big integer "session" context. */
68: {
69: bigint *active_list; /**< Bigints currently used. */
70: bigint *free_list; /**< Bigints not used. */
71: bigint *bi_radix; /**< The radix used. */
72: bigint *bi_mod[BIGINT_NUM_MODS]; /**< modulus */
73:
74: #if defined(CONFIG_BIGINT_MONTGOMERY)
75: bigint *bi_RR_mod_m[BIGINT_NUM_MODS]; /**< R^2 mod m */
76: bigint *bi_R_mod_m[BIGINT_NUM_MODS]; /**< R mod m */
77: comp N0_dash[BIGINT_NUM_MODS];
78: #elif defined(CONFIG_BIGINT_BARRETT)
79: bigint *bi_mu[BIGINT_NUM_MODS]; /**< Storage for mu */
80: #endif
81: bigint *bi_normalised_mod[BIGINT_NUM_MODS]; /**< Normalised mod storage. */
82: bigint **g; /**< Used by sliding-window. */
83: int window; /**< The size of the sliding window */
84: int active_count; /**< Number of active bigints. */
85: int free_count; /**< Number of free bigints. */
86:
87: #ifdef CONFIG_BIGINT_MONTGOMERY
88: uint8_t use_classical; /**< Use classical reduction. */
89: #endif
90: uint8_t mod_offset; /**< The mod offset we are using */
91: } BI_CTX;
92:
93: #ifndef WIN32
94: #define max(a,b) ((a)>(b)?(a):(b)) /**< Find the maximum of 2 numbers. */
95: #define min(a,b) ((a)<(b)?(a):(b)) /**< Find the minimum of 2 numbers. */
96: #endif
97:
98: #define PERMANENT 0x7FFF55AA /**< A magic number for permanents. */
99:
100: #define V1 v->comps[v->size-1] /**< v1 for division */
101: #define V2 v->comps[v->size-2] /**< v2 for division */
102: #define U(j) tmp_u->comps[tmp_u->size-j-1] /**< uj for division */
103: #define Q(j) quotient->comps[quotient->size-j-1] /**< qj for division */
104:
105: #endif
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.