9
FILE_LICENCE ( GPL2_OR_LATER_OR_UBDL );
15
/** Element of a big integer */
16
typedef uint32_t bigint_element_t;
19
* Initialise big integer
21
* @v value0 Element 0 of big integer to initialise
22
* @v size Number of elements
24
* @v len Length of raw data
26
static inline __attribute__ (( always_inline )) void
27
bigint_init_raw ( uint32_t *value0, unsigned int size,
28
const void *data, size_t len ) {
29
size_t pad_len = ( sizeof ( bigint_t ( size ) ) - len );
30
uint8_t *value_byte = ( ( void * ) value0 );
31
const uint8_t *data_byte = ( data + len );
33
/* Copy raw data in reverse order, padding with zeros */
35
*(value_byte++) = *(--data_byte);
43
* @v addend0 Element 0 of big integer to add
44
* @v value0 Element 0 of big integer to be added to
45
* @v size Number of elements
47
static inline __attribute__ (( always_inline )) void
48
bigint_add_raw ( const uint32_t *addend0, uint32_t *value0,
50
bigint_t ( size ) __attribute__ (( may_alias )) *value =
51
( ( void * ) value0 );
52
uint32_t *discard_addend;
53
uint32_t *discard_value;
54
uint32_t *discard_end;
55
uint32_t discard_addend_i;
56
uint32_t discard_value_i;
58
__asm__ __volatile__ ( "adds %2, %0, %8, lsl #2\n\t" /* clear CF */
66
: "=l" ( discard_addend ),
67
"=l" ( discard_value ),
69
"=l" ( discard_addend_i ),
70
"=l" ( discard_value_i ),
72
: "0" ( addend0 ), "1" ( value0 ), "l" ( size )
77
* Subtract big integers
79
* @v subtrahend0 Element 0 of big integer to subtract
80
* @v value0 Element 0 of big integer to be subtracted from
81
* @v size Number of elements
83
static inline __attribute__ (( always_inline )) void
84
bigint_subtract_raw ( const uint32_t *subtrahend0, uint32_t *value0,
86
bigint_t ( size ) __attribute__ (( may_alias )) *value =
87
( ( void * ) value0 );
88
uint32_t *discard_subtrahend;
89
uint32_t *discard_value;
90
uint32_t *discard_end;
91
uint32_t discard_subtrahend_i;
92
uint32_t discard_value_i;
94
__asm__ __volatile__ ( "add %2, %0, %8, lsl #2\n\t"
95
"cmp %2, %0\n\t" /* set CF */
100
"stmia %1!, {%4}\n\t"
103
: "=l" ( discard_subtrahend ),
104
"=l" ( discard_value ),
105
"=l" ( discard_end ),
106
"=l" ( discard_subtrahend_i ),
107
"=l" ( discard_value_i ),
109
: "0" ( subtrahend0 ), "1" ( value0 ),
115
* Rotate big integer left
117
* @v value0 Element 0 of big integer
118
* @v size Number of elements
120
static inline __attribute__ (( always_inline )) void
121
bigint_rol_raw ( uint32_t *value0, unsigned int size ) {
122
bigint_t ( size ) __attribute__ (( may_alias )) *value =
123
( ( void * ) value0 );
124
uint32_t *discard_value;
125
uint32_t *discard_end;
126
uint32_t discard_value_i;
128
__asm__ __volatile__ ( "adds %1, %0, %5, lsl #2\n\t" /* clear CF */
132
"stmia %0!, {%2}\n\t"
135
: "=l" ( discard_value ),
136
"=l" ( discard_end ),
137
"=l" ( discard_value_i ),
139
: "0" ( value0 ), "1" ( size )
144
* Rotate big integer right
146
* @v value0 Element 0 of big integer
147
* @v size Number of elements
149
static inline __attribute__ (( always_inline )) void
150
bigint_ror_raw ( uint32_t *value0, unsigned int size ) {
151
bigint_t ( size ) __attribute__ (( may_alias )) *value =
152
( ( void * ) value0 );
153
uint32_t *discard_value;
154
uint32_t *discard_end;
155
uint32_t discard_value_i;
157
__asm__ __volatile__ ( "adds %1, %0, %5, lsl #2\n\t" /* clear CF */
159
"ldmdb %1!, {%2}\n\t"
164
: "=l" ( discard_value ),
165
"=l" ( discard_end ),
166
"=l" ( discard_value_i ),
168
: "0" ( value0 ), "1" ( size )
173
* Test if big integer is equal to zero
175
* @v value0 Element 0 of big integer
176
* @v size Number of elements
177
* @ret is_zero Big integer is equal to zero
179
static inline __attribute__ (( always_inline, pure )) int
180
bigint_is_zero_raw ( const uint32_t *value0, unsigned int size ) {
181
const uint32_t *value = value0;
185
value_i = *(value++);
190
return ( value_i == 0 );
194
* Compare big integers
196
* @v value0 Element 0 of big integer
197
* @v reference0 Element 0 of reference big integer
198
* @v size Number of elements
199
* @ret geq Big integer is greater than or equal to the reference
201
static inline __attribute__ (( always_inline, pure )) int
202
bigint_is_geq_raw ( const uint32_t *value0, const uint32_t *reference0,
203
unsigned int size ) {
204
const uint32_t *value = ( value0 + size );
205
const uint32_t *reference = ( reference0 + size );
207
uint32_t reference_i;
210
value_i = *(--value);
211
reference_i = *(--reference);
212
if ( value_i != reference_i )
216
return ( value_i >= reference_i );
220
* Test if bit is set in big integer
222
* @v value0 Element 0 of big integer
223
* @v size Number of elements
225
* @ret is_set Bit is set
227
static inline __attribute__ (( always_inline )) int
228
bigint_bit_is_set_raw ( const uint32_t *value0, unsigned int size,
230
const bigint_t ( size ) __attribute__ (( may_alias )) *value =
231
( ( const void * ) value0 );
232
unsigned int index = ( bit / ( 8 * sizeof ( value->element[0] ) ) );
233
unsigned int subindex = ( bit % ( 8 * sizeof ( value->element[0] ) ) );
235
return ( value->element[index] & ( 1 << subindex ) );
239
* Find highest bit set in big integer
241
* @v value0 Element 0 of big integer
242
* @v size Number of elements
243
* @ret max_bit Highest bit set + 1 (or 0 if no bits set)
245
static inline __attribute__ (( always_inline )) int
246
bigint_max_set_bit_raw ( const uint32_t *value0, unsigned int size ) {
247
const uint32_t *value = ( value0 + size );
248
int max_bit = ( 8 * sizeof ( bigint_t ( size ) ) );
252
value_i = *(--value);
253
max_bit -= ( 32 - fls ( value_i ) );
264
* @v source0 Element 0 of source big integer
265
* @v source_size Number of elements in source big integer
266
* @v dest0 Element 0 of destination big integer
267
* @v dest_size Number of elements in destination big integer
269
static inline __attribute__ (( always_inline )) void
270
bigint_grow_raw ( const uint32_t *source0, unsigned int source_size,
271
uint32_t *dest0, unsigned int dest_size ) {
272
unsigned int pad_size = ( dest_size - source_size );
274
memcpy ( dest0, source0, sizeof ( bigint_t ( source_size ) ) );
275
memset ( ( dest0 + source_size ), 0, sizeof ( bigint_t ( pad_size ) ) );
281
* @v source0 Element 0 of source big integer
282
* @v source_size Number of elements in source big integer
283
* @v dest0 Element 0 of destination big integer
284
* @v dest_size Number of elements in destination big integer
286
static inline __attribute__ (( always_inline )) void
287
bigint_shrink_raw ( const uint32_t *source0, unsigned int source_size __unused,
288
uint32_t *dest0, unsigned int dest_size ) {
290
memcpy ( dest0, source0, sizeof ( bigint_t ( dest_size ) ) );
294
* Finalise big integer
296
* @v value0 Element 0 of big integer to finalise
297
* @v size Number of elements
298
* @v out Output buffer
299
* @v len Length of output buffer
301
static inline __attribute__ (( always_inline )) void
302
bigint_done_raw ( const uint32_t *value0, unsigned int size __unused,
303
void *out, size_t len ) {
304
const uint8_t *value_byte = ( ( const void * ) value0 );
305
uint8_t *out_byte = ( out + len );
307
/* Copy raw data in reverse order */
309
*(--out_byte) = *(value_byte++);
312
extern void bigint_multiply_raw ( const uint32_t *multiplicand0,
313
const uint32_t *multiplier0,
314
uint32_t *value0, unsigned int size );
316
#endif /* _BITS_BIGINT_H */