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