shadata.c

   1 #include <math.h>
   2 #include <stdio.h>
   3
   4 static const unsigned primes[64] =
   5 {
   6   2, 3, 5, 7, 11, 13, 17, 19, 23, 29,
   7   31, 37, 41, 43, 47, 53, 59, 61, 67, 71,
   8   73, 79, 83, 89, 97, 101, 103, 107, 109, 113,
   9   127, 131, 137, 139, 149, 151, 157, 163, 167, 173,
  10   179, 181, 191, 193, 197, 199, 211, 223, 227, 229,
  11   233, 239, 241, 251, 257, 263, 269, 271, 277, 281,
  12   283, 293, 307, 311
  13 };
  14
  15 int main(int argc, char **argv)
  16 {
  17   int i;
  18   static const double third = 1.0/3;
  19
  20   printf("SHA-256 constants: \n");
  21   for (i = 0; i < 64; )
  22     {
  23       double root = pow(primes[i++], third);
  24       double fraction = root - floor(root);
  25       double value = floor(ldexp(fraction, 32));
  26
  27       printf("0x%08lxUL, ", (unsigned long) value);
  28       if (!(i % 4))
  29         printf("\n");
  30     }
  31
  32   printf("\nSHA-256 initial values: \n");
  33
  34   for (i = 0; i < 8; )
  35     {
  36       double root = pow(primes[i++], 0.5);
  37       double fraction = root - (floor(root));
  38       double value = floor(ldexp(fraction, 32));
  39
  40       printf("0x%08lxUL, ", (unsigned long) value);
  41       if (!(i % 4))
  42         printf("\n");
  43     }
  44
  45   return 0;
  46 }