main_gfc128.c

#include <stdlib.h>
#include <stdio.h>

#include "gf2_128.h"
#include "gf_complete.h"

const unsigned long MAX_ROUNDS = 100000000;

int main() {
  unsigned int randreg;
  //FILE* urand = fopen("/dev/urandom", "r");
  //fread( &randreg, sizeof(randreg), 1, urand );
  //fclose(urand);
  randreg = 0x1650338b;
  printf( "randreg = %#.8x\n", randreg );

  gf_t gf;
  gf_init_easy(&gf, 128);

  uint64_t x[2];
  uint64_t y[2];
  uint64_t z[2];

  unsigned long i;
  for( i = 0; i < MAX_ROUNDS; i++ ) {

    x[0] = (uint64_t) rand_r( &randreg );
    x[0] <<= 32;
    x[0] |= (uint64_t) rand_r( &randreg );

    x[1] = (uint64_t) rand_r( &randreg );
    x[1] <<= 32;
    x[1] |= (uint64_t) rand_r( &randreg );

    /* y = (uint64_t) rand_r( &randreg ); */
    /* y <<= 32; */
    /* y |= (uint64_t) rand_r( &randreg ); */

    /* z = (uint64_t) rand_r( &randreg ); */
    /* z <<= 32; */
    /* z |= (uint64_t) rand_r( &randreg ); */

    /* // Distrbutivity test */
    /* if( gf.multiply.w64( &gf, x, y ^ z ) == */
    /*     gf.multiply.w64( &gf, x, y ) ^ gf.multiply.w64( &gf, x, z ) ) { */
    /* } else { */
    /*   printf( "Distributivity fail:\n" ); */
    /*   printf( "  x = 0x%.16llx\n", x ); */
    /*   printf( "  y = 0x%.16llx\n", y ); */
    /*   printf( "  z = 0x%.16llx\n", z ); */
    /* } */

    /* // Associativity test */
    /* if( gf.multiply.w64( &gf, x, gf.multiply.w64( &gf, y, z ) ) == */
    /*     gf.multiply.w64( &gf, gf.multiply.w64( &gf, x, y ), z ) ) { */
    /* } else { */
    /*   printf( "Associativity fail:\n" ); */
    /*   printf( "  x = 0x%.16llx\n", x ); */
    /*   printf( "  y = 0x%.16llx\n", y ); */
    /*   printf( "  z = 0x%.16llx\n", z ); */
    /* } */

    /* // inverse test */
    /* if( x == 0 || gf.inverse.w64( &gf, x ) != 0) { */
    /* } else { */
    /*   printf( "Inversion fail:\n" ); */
    /*   printf( "  x = 0x%.16llx\n", x ); */
    /* } */

    // inverse test
    if( ~(x[0] == 0 && x[1] == 0) ) {
      gf.inverse.w128( &gf, x, y );
      gf.multiply.w128( &gf, x, y, z );
      if( z[0] == 0 && z[1] == 1 ) {
      } else {
        printf( "Inversion fail:\n" );
        printf( "  x = 0x%.16llx%.16llx\n", (uint64_t) (x[1]), (uint64_t) x[0] );
      }
    }

  }

  printf( "%lu Tests completed\n", MAX_ROUNDS );
 
}