appendix-b.html

Wiley - Applied Cryptography, Protocols, Algorthms, and Source Code in C

static unsigned char pc1[56] = {
       56, 48, 40, 32, 24, 16,  8,   0, 57, 49, 41, 33, 25, 17,
        9,  1, 58, 50, 42, 34, 26,  18, 10,  2, 59, 51, 43, 35,
       62, 54, 46, 38, 30, 22, 14,   6, 61, 53, 45, 37, 29, 21,
       13,  5, 60, 52, 44, 36, 28,  20, 12,  4, 27, 19, 11,  3 };

static unsigned char totrot[16] = {
       1,2,4,6,8,10,12,14,15,17,19,21,23,25,27,28 };

static unsigned char pc2[48] = {
       13, 16, 10, 23,  0,  4,       2, 27, 14,  5, 20,  9,
       22, 18, 11,  3, 25,  7,      15,  6, 26, 19, 12,  1,
       40, 51, 30, 36, 46, 54,      29, 39, 50, 44, 32, 47,
       43, 48, 38, 55, 33, 52,      45, 41, 49, 35, 28, 31 };
void deskey(key, edf)       /* Thanks to James Gillogly &amp Phil Karn! */
unsigned char *key;
short edf;
{
       register int i, j, l, m, n;
       unsigned char pc1m[56], pcr[56];
       unsigned long kn[32];

       for ( j = 0; j < 56; j++ ) {
              l = pc1[j];
              m = l &amp 07;
              pc1m[j] = (key[l >> 3] &amp bytebit[m]) ? 1 : 0;
              }
       for( i = 0; i < 16; i++ ) {
              if( edf == DE1 ) m = (15 – i) << 1;
              else m = i << 1;
              n = m + 1;
              kn[m] = kn[n] = 0L;
              for( j = 0; j < 28; j++ ) {
                     l = j + totrot[i];
                     if( l < 28 ) pcr[j] = pc1m[l];
                     else pcr[j] = pc1m[l – 28];
                     }
              for( j = 28; j < 56; j++ ) {
                  l = j + totrot[i];
                  if( l < 56 ) pcr[j] = pc1m[l];
                  else pcr[j] = pc1m[l – 28];
                  }
              for( j = 0; j < 24; j++ ) {
                     if( pcr[pc2[j]] ) kn[m] |= bigbyte[j];
                     if( pcr[pc2[j+24]] ) kn[n] |= bigbyte[j];
                     }
              }
       cookey(kn);
       return;
}

static void cookey(raw1)
register unsigned long *raw1;
{
       register unsigned long *cook, *raw0;
       unsigned long dough[32];
       register int i;

       cook = dough;
       for( i = 0; i < 16; i++, raw1++ ) {
              raw0 = raw1++;
              *cook   = (*raw0 &amp 0×00fc0000L) << 6;
              *cook  |= (*raw0 &amp 0×00000fc0L) << 10;
              *cook  |= (*raw1 &amp 0×00fc0000L) >> 10;
              *cook++       |= (*raw1 &amp 0×00000fc0L) >> 6;
              *cook   = (*raw0 &amp 0×0003f000L) << 12;
              *cook  |= (*raw0 &amp 0×0000003fL) << 16;
              *cook  |= (*raw1 &amp 0×0003f000L) >> 4;
              *cook++       |= (*raw1 &amp 0×0000003fL);
              }
       usekey(dough);
       return;
}

void cpkey(into)
register unsigned long *into;
{
       register unsigned long *from, *endp;
       from = KnL, endp = &ampKnL[32];
       while( from < endp ) *into++ = *from++;
       return;
}

void usekey(from)
register unsigned long *from;
{
       register unsigned long *to, *endp;
       to = KnL, endp = &ampKnL[32];
       while( to < endp ) *to++ = *from++;
       return;
}

void des(inblock, outblock)
unsigned char *inblock, *outblock;
{
       unsigned long work[2];

       scrunch(inblock, work);
       desfunc(work, KnL);
       unscrun(work, outblock);
       return;
}

static void scrunch(outof, into)
register unsigned char *outof;
register unsigned long *into;
{
       *into   = (*outof++ &amp 0xffL) << 24;
       *into  |= (*outof++ &amp 0xffL) << 16;
       *into  |= (*outof++ &amp 0xffL) << 8;
       *into++ |= (*outof++ &amp 0xffL);
       *into   = (*outof++ &amp 0xffL) << 24;
       *into  |= (*outof++ &amp 0xffL) << 16;
       *into  |= (*outof++ &amp 0xffL) << 8;
       *into  |= (*outof   &amp 0xffL);
       return;
}

static void unscrun(outof, into)
register unsigned long *outof;
register unsigned char *into;
{
       *into++ = (*outof >> 24) &amp 0xffL;
       *into++ = (*outof >> 16) &amp 0xffL;
       *into++ = (*outof >>  8) &amp 0xffL;
       *into++ =  *outof++           &amp 0xffL;
       *into++ = (*outof >> 24) &amp 0xffL;
       *into++ = (*outof >> 16) &amp 0xffL;
       *into++ = (*outof >>  8) &amp 0xffL;
       *into   =  *outof     &amp 0xffL;
       return;
}
static unsigned long SP1[64] = {
       0×01010400L, 0×00000000L, 0×00010000L, 0×01010404L,
       0×01010004L, 0×00010404L, 0×00000004L, 0×00010000L,
       0×00000400L, 0×01010400L, 0×01010404L, 0×00000400L,
       0×01000404L, 0×01010004L, 0×01000000L, 0×00000004L,
       0×00000404L, 0×01000400L, 0×01000400L, 0×00010400L,
       0×00010400L, 0×01010000L, 0×01010000L, 0×01000404L,
       0×00010004L, 0×01000004L, 0×01000004L, 0×00010004L,
       0×00000000L, 0×00000404L, 0×00010404L, 0×01000000L,
       0×00010000L, 0×01010404L, 0×00000004L, 0×01010000L,
       0×01010400L, 0×01000000L, 0×01000000L, 0×00000400L,
       0×01010004L, 0×00010000L, 0×00010400L, 0×01000004L,
       0×00000400L, 0×00000004L, 0×01000404L, 0×00010404L,
       0×01010404L, 0×00010004L, 0×01010000L, 0×01000404L,
       0×01000004L, 0×00000404L, 0×00010404L, 0×01010400L,
       0×00000404L, 0×01000400L, 0×01000400L, 0×00000000L,
       0×00010004L, 0×00010400L, 0×00000000L, 0×01010004L };

static unsigned long SP2[64] = {
       0x80108020L, 0x80008000L, 0×00008000L, 0×00108020L,
       0×00100000L, 0×00000020L, 0x80100020L, 0x80008020L,
       0x80000020L, 0x80108020L, 0x80108000L, 0x80000000L,
       0x80008000L, 0×00100000L, 0×00000020L, 0x80100020L,
       0×00108000L, 0×00100020L, 0x80008020L, 0×00000000L,
       0x80000000L, 0×00008000L, 0×00108020L, 0x80100000L,
       0×00100020L, 0x80000020L, 0×00000000L, 0×00108000L,
       0×00008020L, 0x80108000L, 0x80100000L, 0×00008020L,
       0×00000000L, 0×00108020L, 0x80100020L, 0×00100000L,
       0x80008020L, 0x80100000L, 0x80108000L, 0×00008000L,
       0x80100000L, 0x80008000L, 0×00000020L, 0x80108020L,
       0×00108020L, 0×00000020L, 0×00008000L, 0x80000000L,
       0×00008020L, 0x80108000L, 0×00100000L, 0x80000020L,
       0×00100020L, 0x80008020L, 0x80000020L, 0×00100020L,
       0×00108000L, 0×00000000L, 0x80008000L, 0×00008020L,
       0x80000000L, 0x80100020L, 0x80108020L, 0×00108000L };

static unsigned long SP3[64] = {
       0×00000208L, 0×08020200L, 0×00000000L, 0×08020008L,
       0×08000200L, 0×00000000L, 0×00020208L, 0×08000200L,
       0×00020008L, 0×08000008L, 0×08000008L, 0×00020000L,
       0×08020208L, 0×00020008L, 0×08020000L, 0×00000208L,
       0×08000000L, 0×00000008L, 0×08020200L, 0×00000200L,
       0×00020200L, 0×08020000L, 0×08020008L, 0×00020208L,
       0×08000208L, 0×00020200L, 0×00020000L, 0×08000208L,
       0×00000008L, 0×08020208L, 0×00000200L, 0×08000000L,
       0×08020200L, 0×08000000L, 0×00020008L, 0×00000208L,
       0×00020000L, 0×08020200L, 0×08000200L, 0×00000000L,
       0×00000200L, 0×00020008L, 0×08020208L, 0×08000200L,
       0×08000008L, 0×00000200L, 0×00000000L, 0×08020008L,
       0×08000208L, 0×00020000L, 0×08000000L, 0×08020208L,
       0×00000008L, 0×00020208L, 0×00020200L, 0×08000008L,
       0×08020000L, 0×08000208L, 0×00000208L, 0×08020000L,
       0×00020208L, 0×00000008L, 0×08020008L, 0×00020200L };

static unsigned long SP4[64] = {
       0×00802001L, 0×00002081L, 0×00002081L, 0×00000080L,
       0×00802080L, 0×00800081L, 0×00800001L, 0×00002001L,
       0×00000000L, 0×00802000L, 0×00802000L, 0×00802081L,
       0×00000081L, 0×00000000L, 0×00800080L, 0×00800001L,
       0×00000001L, 0×00002000L, 0×00800000L, 0×00802001L,
       0×00000080L, 0×00800000L, 0×00002001L, 0×00002080L,
       0×00800081L, 0×00000001L, 0×00002080L, 0×00800080L,
       0×00002000L, 0×00802080L, 0×00802081L, 0×00000081L,
       0×00800080L, 0×00800001L, 0×00802000L, 0×00802081L,
       0×00000081L, 0×00000000L, 0×00000000L, 0×00802000L,
       0×00002080L, 0×00800080L, 0×00800081L, 0×00000001L,
       0×00802001L, 0×00002081L, 0×00002081L, 0×00000080L,
       0×00802081L, 0×00000081L, 0×00000001L, 0×00002000L,
       0×00800001L, 0×00002001L, 0×00802080L, 0×00800081L,
       0×00002001L, 0×00002080L, 0×00800000L, 0×00802001L,
       0×00000080L, 0×00800000L, 0×00002000L, 0×00802080L };

static unsigned long SP5[64] = {
       0×00000100L, 0×02080100L, 0×02080000L, 0x42000100L,
       0×00080000L, 0×00000100L, 0x40000000L, 0×02080000L,
       0x40080100L, 0×00080000L, 0×02000100L, 0x40080100L,
       0x42000100L, 0x42080000L, 0×00080100L, 0x40000000L,
       0×02000000L, 0x40080000L, 0x40080000L, 0×00000000L,
       0x40000100L, 0x42080100L, 0x42080100L, 0×02000100L,
       0x42080000L, 0x40000100L, 0×00000000L, 0x42000000L,
       0×02080100L, 0×02000000L, 0x42000000L, 0×00080100L,
       0×00080000L, 0x42000100L, 0×00000100L, 0×02000000L,
       0x40000000L, 0×02080000L, 0x42000100L, 0x40080100L,
       0×02000100L, 0x40000000L, 0x42080000L, 0×02080100L,
       0x40080100L, 0×00000100L, 0×02000000L, 0x42080000L,
       0x42080100L, 0×00080100L, 0x42000000L, 0x42080100L,
       0×02080000L, 0×00000000L, 0x40080000L, 0x42000000L,
       0×00080100L, 0×02000100L, 0x40000100L, 0×00080000L,
       0×00000000L, 0x40080000L, 0×02080100L, 0x40000100L };

static unsigned long SP6[64] = {
       0x20000010L, 0x20400000L, 0×00004000L, 0x20404010L,
       0x20400000L, 0×00000010L, 0x20404010L, 0×00400000L,
       0x20004000L, 0×00404010L, 0×00400000L, 0x20000010L,
       0×00400010L, 0x20004000L, 0x20000000L, 0×00004010L,
       0×00000000L, 0×00400010L, 0x20004010L, 0×00004000L,
       0×00404000L, 0x20004010L, 0×00000010L, 0x20400010L,
       0x20400010L, 0×00000000L, 0×00404010L, 0x20404000L,
       0×00004010L, 0×00404000L, 0x20404000L, 0x20000000L,
       0x20004000L, 0×00000010L, 0x20400010L, 0×00404000L,
       0x20404010L, 0×00400000L, 0×00004010L, 0x20000010L,
       0×00400000L, 0x20004000L, 0x20000000L, 0×00004010L,
       0x20000010L, 0x20404010L, 0×00404000L, 0x20400000L,
       0×00404010L, 0x20404000L, 0×00000000L, 0x20400010L,
       0×00000010L, 0×00004000L, 0x20400000L, 0×00404010L,
       0×00004000L, 0×00400010L, 0x20004010L, 0×00000000L,
       0x20404000L, 0x20000000L, 0×00400010L, 0x20004010L };

static unsigned long SP7[64] = {
       0×00200000L, 0×04200002L, 0×04000802L, 0×00000000L,
       0×00000800L, 0×04000802L, 0×00200802L, 0×04200800L,
       0×04200802L, 0×00200000L, 0×00000000L, 0×04000002L,
       0×00000002L, 0×04000000L, 0×04200002L, 0×00000802L,
       0×04000800L, 0×00200802L, 0×00200002L, 0×04000800L,
       0×04000002L, 0×04200000L, 0×04200800L, 0×00200002L,
       0×04200000L, 0×00000800L, 0×00000802L, 0×04200802L,
       0×00200800L, 0×00000002L, 0×04000000L, 0×00200800L,
       0×04000000L, 0×00200800L, 0×00200000L, 0×04000802L,
       0×04000802L, 0×04200002L, 0×04200002L, 0×00000002L,
       0×00200002L, 0×04000000L, 0×04000800L, 0×00200000L,
       0×04200800L, 0×00000802L, 0×00200802L, 0×04200800L,
       0×00000802L, 0×04000002L, 0×04200802L, 0×04200000L,
       0×00200800L, 0×00000000L, 0×00000002L, 0×04200802L,
       0×00000000L, 0×00200802L, 0×04200000L, 0×00000800L,
       0×04000002L, 0×04000800L, 0×00000800L, 0×00200002L };

static unsigned long SP8[64] = {
       0x10001040L, 0×00001000L, 0×00040000L, 0x10041040L,
       0x10000000L, 0x10001040L, 0×00000040L, 0x10000000L,
       0×00040040L, 0x10040000L, 0x10041040L, 0×00041000L,
       0x10041000L, 0×00041040L, 0×00001000L, 0×00000040L,
       0x10040000L, 0x10000040L, 0x10001000L, 0×00001040L,
       0×00041000L, 0×00040040L, 0x10040040L, 0x10041000L,
       0×00001040L, 0×00000000L, 0×00000000L, 0x10040040L,
       0x10000040L, 0x10001000L, 0×00041040L, 0×00040000L,
       0×00041040L, 0×00040000L, 0x10041000L, 0×00001000L,
       0×00000040L, 0x10040040L, 0×00001000L, 0×00041040L,
       0x10001000L, 0×00000040L, 0x10000040L, 0x10040000L,
       0x10040040L, 0x10000000L, 0×00040000L, 0x10001040L,
       0×00000000L, 0x10041040L, 0×00040040L, 0x10000040L,
       0x10040000L, 0x10001000L, 0x10001040L, 0×00000000L,
       0x10041040L, 0×00041000L, 0×00041000L, 0×00001040L,
       0×00001040L, 0×00040040L, 0x10000000L, 0x10041000L };

static void desfunc(block, keys)
register unsigned long *block, *keys;
{
       register unsigned long fval, work, right, leftt;
       register int round;

       leftt = block[0];
       right = block[1];
       work = ((leftt >> 4) ^ right) &amp 0×0f0f0f0fL;
       right ^= work;
       leftt ^= (work << 4);
       work = ((leftt >> 16) ^ right) &amp 0×0000ffffL;
       right ^= work;
       leftt ^= (work << 16);
       work = ((right >> 2) ^ leftt) &amp 0x33333333L;
       leftt ^= work;
       right ^= (work << 2);
       work = ((right >> 8) ^ leftt) &amp 0×00ff00ffL;
       leftt ^= work;
       right ^= (work << 8);
       right = ((right << 1) | ((right >> 31) &amp 1L)) &amp 0xffffffffL;
       work = (leftt ^ right) &amp 0xaaaaaaaaL;
       leftt ^= work;
       right ^= work;
       leftt = ((leftt << 1) | ((leftt >> 31) &amp 1L)) &amp 0xffffffffL;

       for( round = 0; round < 8; round++ ) {
              work  = (right << 28) | (right >> 4);
              work ^= *keys++;
              fval  = SP7[ work             &amp 0x3fL];
              fval |= SP5[(work >>  8) &amp 0x3fL];
              fval |= SP3[(work >> 16) &amp 0x3fL];
              fval |= SP1[(work >> 24) &amp 0x3fL];
              work  = right ^ *keys++;
              fval |= SP8[ work             &amp 0x3fL];
              fval |= SP6[(work >>  8) &amp 0x3fL];
              fval |= SP4[(work >> 16) &amp 0x3fL];
              fval |= SP2[(work >> 24) &amp 0x3fL];
              leftt ^= fval;
              work  = (leftt << 28) | (leftt >> 4);
              work ^= *keys++;
              fval  = SP7[ work             &amp 0x3fL];
              fval |= SP5[(work >>  8) &amp 0x3fL];
              fval |= SP3[(work >> 16) &amp 0x3fL];
              fval |= SP1[(work >> 24) &amp 0x3fL];
              work  = leftt ^ *keys++;
              fval |= SP8[ work             &amp 0x3fL];
              fval |= SP6[(work >>  8) &amp 0x3fL];
              fval |= SP4[(work >> 16) &amp 0x3fL];
              fval |= SP2[(work >> 24) &amp 0x3fL];
              right ^= fval;
              }

       right = (right << 31) | (right >> 1);
       work = (leftt ^ right) &amp 0xaaaaaaaaL;
       leftt ^= work;
       right ^= work;
       leftt = (leftt << 31) | (leftt >> 1);
       work = ((leftt >> 8) ^ right) &amp 0×00ff00ffL;
       right ^= work;
       leftt ^= (work << 8);
       work = ((leftt >> 2) ^ right) &amp 0x33333333L;
       right ^= work;
       leftt ^= (work << 2);
       work = ((right >> 16) ^ leftt) &amp 0×0000ffffL;
       leftt ^= work;