IDA_new/gf-complete/examples/gf_example_2.c

   1 /*
   2  * GF-Complete: A Comprehensive Open Source Library for Galois Field Arithmetic
   3  * James S. Plank, Ethan L. Miller, Kevin M. Greenan,
   4  * Benjamin A. Arnold, John A. Burnum, Adam W. Disney, Allen C. McBride.
   5  *
   6  * gf_example_2.c
   7  *
   8  * Demonstrates using the procedures for examples in GF(2^w) for w <= 32.
   9  */
  10
  11 #include <stdio.h>
  12 #include <getopt.h>
  13 #include <stdint.h>
  14 #include <string.h>
  15 #include <stdlib.h>
  16 #include <time.h>
  17
  18 #include "gf_complete.h"
  19 #include "gf_rand.h"
  20
  21 void usage(char *s)
  22 {
  23   fprintf(stderr, "usage: gf_example_2 w - w must be between 1 and 32\n");
  24   exit(1);
  25 }
  26
  27 int main(int argc, char **argv)
  28 {
  29   uint32_t a, b, c;
  30   uint8_t *r1, *r2;
  31   uint16_t *r16 = NULL;
  32   uint32_t *r32 = NULL;
  33   int w, i;
  34   gf_t gf;
  35
  36   if (argc != 2) usage(NULL);
  37   w = atoi(argv[1]);
  38   if (w <= 0 || w > 32) usage("Bad w");
  39
  40   /* Get two random numbers in a and b */
  41
  42   MOA_Seed(time(0));
  43   a = MOA_Random_W(w, 0);
  44   b = MOA_Random_W(w, 0);
  45
  46   /* Create the proper instance of the gf_t object using defaults: */
  47
  48   gf_init_easy(&gf, w);
  49
  50   /* And multiply a and b using the galois field: */
  51
  52   c = gf.multiply.w32(&gf, a, b);
  53   printf("%u * %u = %u\n", a, b, c);
  54
  55   /* Divide the product by a and b */
  56
  57   printf("%u / %u = %u\n", c, a, gf.divide.w32(&gf, c, a));
  58   printf("%u / %u = %u\n", c, b, gf.divide.w32(&gf, c, b));
  59
  60   /* If w is 4, 8, 16 or 32, do a very small region operation */
  61
  62   if (w == 4 || w == 8 || w == 16 || w == 32) {
  63     r1 = (uint8_t *) malloc(16);
  64     r2 = (uint8_t *) malloc(16);
  65
  66     if (w == 4 || w == 8) {
  67       r1[0] = b;
  68       for (i = 1; i < 16; i++) r1[i] = MOA_Random_W(8, 1);
  69     } else if (w == 16) {
  70       r16 = (uint16_t *) r1;
  71       r16[0] = b;
  72       for (i = 1; i < 8; i++) r16[i] = MOA_Random_W(16, 1);
  73     } else {
  74       r32 = (uint32_t *) r1;
  75       r32[0] = b;
  76       for (i = 1; i < 4; i++) r32[i] = MOA_Random_W(32, 1);
  77     }
  78
  79     gf.multiply_region.w32(&gf, r1, r2, a, 16, 0);
  80
  81     printf("\nmultiply_region by 0x%x (%u)\n\n", a, a);
  82     printf("R1 (the source):  ");
  83     if (w == 4) {
  84       for (i = 0; i < 16; i++) printf(" %x %x", r1[i] >> 4, r1[i] & 0xf);
  85     } else if (w == 8) {
  86       for (i = 0; i < 16; i++) printf(" %02x", r1[i]);
  87     } else if (w == 16) {
  88       for (i = 0; i < 8; i++) printf(" %04x", r16[i]);
  89     } else if (w == 32) {
  90       for (i = 0; i < 4; i++) printf(" %08x", r32[i]);
  91     }
  92     printf("\nR2 (the product): ");
  93     if (w == 4) {
  94       for (i = 0; i < 16; i++) printf(" %x %x", r2[i] >> 4, r2[i] & 0xf);
  95     } else if (w == 8) {
  96       for (i = 0; i < 16; i++) printf(" %02x", r2[i]);
  97     } else if (w == 16) {
  98       r16 = (uint16_t *) r2;
  99       for (i = 0; i < 8; i++) printf(" %04x", r16[i]);
 100     } else if (w == 32) {
 101       r32 = (uint32_t *) r2;
 102       for (i = 0; i < 4; i++) printf(" %08x", r32[i]);
 103     }
 104     printf("\n");
 105   }
 106   exit(0);
 107 }