userapps/opensource/sshd/libtommath/etc/mont.c

   1 /* tests the montgomery routines */
   2 #include <tommath.h>
   3
   4 int main(void)
   5 {
   6    mp_int modulus, R, p, pp;
   7    mp_digit mp;
   8    long x, y;
   9
  10    srand(time(NULL));
  11    mp_init_multi(&modulus, &R, &p, &pp, NULL);
  12
  13    /* loop through various sizes */
  14    for (x = 4; x < 256; x++) {
  15        printf("DIGITS == %3ld...", x); fflush(stdout);
  16
  17        /* make up the odd modulus */
  18        mp_rand(&modulus, x);
  19        modulus.dp[0] |= 1;
  20
  21        /* now find the R value */
  22        mp_montgomery_calc_normalization(&R, &modulus);
  23        mp_montgomery_setup(&modulus, &mp);
  24
  25        /* now run through a bunch tests */
  26        for (y = 0; y < 1000; y++) {
  27            mp_rand(&p, x/2);        /* p = random */
  28            mp_mul(&p, &R, &pp);     /* pp = R * p */
  29            mp_montgomery_reduce(&pp, &modulus, mp);
  30
  31            /* should be equal to p */
  32            if (mp_cmp(&pp, &p) != MP_EQ) {
  33               printf("FAILURE!\n");
  34               exit(-1);
  35            }
  36        }
  37        printf("PASSED\n");
  38     }
  39
  40     return 0;
  41 }
  42
  43
  44
  45
  46