2 * Copyright(C) 2006 Cameron Rich
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.
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.
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
19 #ifndef BIGINT_IMPL_HEADER
20 #define BIGINT_IMPL_HEADER
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. */
29 #define BIGINT_NUM_MODS 1
32 /* Architecture specific functions for big ints */
34 #define COMP_RADIX 4294967296i64
35 #define COMP_BIG_MSB 0x8000000000000000i64
37 #define COMP_RADIX 4294967296ULL /**< Max component + 1 */
38 #define COMP_BIG_MSB 0x8000000000000000ULL /**< (Max dbl comp + 1)/ 2 */
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. */
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. */
50 * @brief A big integer basic object
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 */
61 typedef struct _bigint bigint
; /**< An alias for _bigint */
64 * Maintains the state of the cache, and a number of variables used in
67 typedef struct /**< A big integer "session" context. */
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 */
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 */
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. */
87 #ifdef CONFIG_BIGINT_MONTGOMERY
88 uint8_t use_classical
; /**< Use classical reduction. */
90 uint8_t mod_offset
; /**< The mod offset we are using */
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. */
98 #define PERMANENT 0x7FFF55AA /**< A magic number for permanents. */
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 */