blob: 6f84f6e4ea04b8b0b780b0aeb6ffea54fae214f4 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
|
/* { dg-options "-fdata-sections" } */
typedef unsigned char unit;
typedef unit *unitptr;
extern short global_precision;
typedef __SIZE_TYPE__ size_t;
extern void *memcpy (void *dest, const void *src, size_t n);
short mp_compare(const unit* r1, const unit* r2)
{
register short precision;
precision = global_precision;
(r1) = ((r1)+(precision)-1);
(r2) = ((r2)+(precision)-1);
do
{ if (*r1 < *r2)
return(-1);
if (*((r1)--) > *((r2)--))
return(1);
} while (--precision);
}
static unit modulus[((1280+(2*8))/8)];
static unit d_data[((1280+(2*8))/8)*2];
int upton_modmult (unitptr prod, unitptr multiplicand, unitptr multiplier)
{
unitptr d = d_data;
while (mp_compare(d,modulus) > 0)
memcpy((void*)(prod), (const void*)(d), (global_precision));
}
|