crypto/desrfb.c

   1 /*
   2  * This is D3DES (V5.09) by Richard Outerbridge with the double and
   3  * triple-length support removed for use in VNC.  Also the bytebit[] array
   4  * has been reversed so that the most significant bit in each byte of the
   5  * key is ignored, not the least significant.
   6  *
   7  * These changes are:
   8  *  Copyright (C) 1999 AT&T Laboratories Cambridge.  All Rights Reserved.
   9  *
  10  * This software is distributed in the hope that it will be useful,
  11  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
  13  */
  14
  15 /* D3DES (V5.09) -
  16  *
  17  * A portable, public domain, version of the Data Encryption Standard.
  18  *
  19  * Written with Symantec's THINK (Lightspeed) C by Richard Outerbridge.
  20  * Thanks to: Dan Hoey for his excellent Initial and Inverse permutation
  21  * code;  Jim Gillogly & Phil Karn for the DES key schedule code; Dennis
  22  * Ferguson, Eric Young and Dana How for comparing notes; and Ray Lau,
  23  * for humouring me on.
  24  *
  25  * Copyright (c) 1988,1989,1990,1991,1992 by Richard Outerbridge.
  26  * (GEnie : OUTER; CIS : [71755,204]) Graven Imagery, 1992.
  27  */
  28
  29 #include "crypto/desrfb.h"
  30
  31 static void scrunch(unsigned char *, unsigned long *);
  32 static void unscrun(unsigned long *, unsigned char *);
  33 static void desfunc(unsigned long *, unsigned long *);
  34 static void cookey(unsigned long *);
  35
  36 static unsigned long KnL[32] = { 0L };
  37
  38 static const unsigned short bytebit[8]  = {
  39         01, 02, 04, 010, 020, 040, 0100, 0200 };
  40
  41 static const unsigned long bigbyte[24] = {
  42         0x800000L,      0x400000L,      0x200000L,      0x100000L,
  43         0x80000L,       0x40000L,       0x20000L,       0x10000L,
  44         0x8000L,        0x4000L,        0x2000L,        0x1000L,
  45         0x800L,         0x400L,         0x200L,         0x100L,
  46         0x80L,          0x40L,          0x20L,          0x10L,
  47         0x8L,           0x4L,           0x2L,           0x1L    };
  48
  49 /* Use the key schedule specified in the Standard (ANSI X3.92-1981). */
  50
  51 static const unsigned char pc1[56] = {
  52         56, 48, 40, 32, 24, 16,  8,      0, 57, 49, 41, 33, 25, 17,
  53          9,  1, 58, 50, 42, 34, 26,     18, 10,  2, 59, 51, 43, 35,
  54         62, 54, 46, 38, 30, 22, 14,      6, 61, 53, 45, 37, 29, 21,
  55         13,  5, 60, 52, 44, 36, 28,     20, 12,  4, 27, 19, 11,  3 };
  56
  57 static const unsigned char totrot[16] = {
  58         1,2,4,6,8,10,12,14,15,17,19,21,23,25,27,28 };
  59
  60 static const unsigned char pc2[48] = {
  61         13, 16, 10, 23,  0,  4,  2, 27, 14,  5, 20,  9,
  62         22, 18, 11,  3, 25,  7, 15,  6, 26, 19, 12,  1,
  63         40, 51, 30, 36, 46, 54, 29, 39, 50, 44, 32, 47,
  64         43, 48, 38, 55, 33, 52, 45, 41, 49, 35, 28, 31 };
  65
  66 /* Thanks to James Gillogly & Phil Karn! */
  67 void deskey(unsigned char *key, int edf)
  68 {
  69         register int i, j, l, m, n;
  70         unsigned char pc1m[56], pcr[56];
  71         unsigned long kn[32];
  72
  73         for ( j = 0; j < 56; j++ ) {
  74                 l = pc1[j];
  75                 m = l & 07;
  76                 pc1m[j] = (key[l >> 3] & bytebit[m]) ? 1 : 0;
  77                 }
  78         for( i = 0; i < 16; i++ ) {
  79                 if( edf == DE1 ) m = (15 - i) << 1;
  80                 else m = i << 1;
  81                 n = m + 1;
  82                 kn[m] = kn[n] = 0L;
  83                 for( j = 0; j < 28; j++ ) {
  84                         l = j + totrot[i];
  85                         if( l < 28 ) pcr[j] = pc1m[l];
  86                         else pcr[j] = pc1m[l - 28];
  87                         }
  88                 for( j = 28; j < 56; j++ ) {
  89                     l = j + totrot[i];
  90                     if( l < 56 ) pcr[j] = pc1m[l];
  91                     else pcr[j] = pc1m[l - 28];
  92                     }
  93                 for( j = 0; j < 24; j++ ) {
  94                         if( pcr[pc2[j]] ) kn[m] |= bigbyte[j];
  95                         if( pcr[pc2[j+24]] ) kn[n] |= bigbyte[j];
  96                         }
  97                 }
  98         cookey(kn);
  99         return;
 100         }
 101
 102 static void cookey(register unsigned long *raw1)
 103 {
 104         register unsigned long *cook, *raw0;
 105         unsigned long dough[32];
 106         register int i;
 107
 108         cook = dough;
 109         for( i = 0; i < 16; i++, raw1++ ) {
 110                 raw0 = raw1++;
 111                 *cook    = (*raw0 & 0x00fc0000L) << 6;
 112                 *cook   |= (*raw0 & 0x00000fc0L) << 10;
 113                 *cook   |= (*raw1 & 0x00fc0000L) >> 10;
 114                 *cook++ |= (*raw1 & 0x00000fc0L) >> 6;
 115                 *cook    = (*raw0 & 0x0003f000L) << 12;
 116                 *cook   |= (*raw0 & 0x0000003fL) << 16;
 117                 *cook   |= (*raw1 & 0x0003f000L) >> 4;
 118                 *cook++ |= (*raw1 & 0x0000003fL);
 119                 }
 120         usekey(dough);
 121         return;
 122         }
 123
 124 void usekey(register unsigned long *from)
 125 {
 126         register unsigned long *to, *endp;
 127
 128         to = KnL, endp = &KnL[32];
 129         while( to < endp ) *to++ = *from++;
 130         return;
 131         }
 132
 133 void des(unsigned char *inblock, unsigned char *outblock)
 134 {
 135         unsigned long work[2];
 136
 137         scrunch(inblock, work);
 138         desfunc(work, KnL);
 139         unscrun(work, outblock);
 140         return;
 141         }
 142
 143 static void scrunch(register unsigned char *outof, register unsigned long *into)
 144 {
 145         *into    = (*outof++ & 0xffL) << 24;
 146         *into   |= (*outof++ & 0xffL) << 16;
 147         *into   |= (*outof++ & 0xffL) << 8;
 148         *into++ |= (*outof++ & 0xffL);
 149         *into    = (*outof++ & 0xffL) << 24;
 150         *into   |= (*outof++ & 0xffL) << 16;
 151         *into   |= (*outof++ & 0xffL) << 8;
 152         *into   |= (*outof   & 0xffL);
 153         return;
 154         }
 155
 156 static void unscrun(register unsigned long *outof, register unsigned char *into)
 157 {
 158         *into++ = (unsigned char)((*outof >> 24) & 0xffL);
 159         *into++ = (unsigned char)((*outof >> 16) & 0xffL);
 160         *into++ = (unsigned char)((*outof >>  8) & 0xffL);
 161         *into++ = (unsigned char)(*outof++       & 0xffL);
 162         *into++ = (unsigned char)((*outof >> 24) & 0xffL);
 163         *into++ = (unsigned char)((*outof >> 16) & 0xffL);
 164         *into++ = (unsigned char)((*outof >>  8) & 0xffL);
 165         *into   =  (unsigned char)(*outof        & 0xffL);
 166         return;
 167         }
 168
 169 static const unsigned long SP1[64] = {
 170         0x01010400L, 0x00000000L, 0x00010000L, 0x01010404L,
 171         0x01010004L, 0x00010404L, 0x00000004L, 0x00010000L,
 172         0x00000400L, 0x01010400L, 0x01010404L, 0x00000400L,
 173         0x01000404L, 0x01010004L, 0x01000000L, 0x00000004L,
 174         0x00000404L, 0x01000400L, 0x01000400L, 0x00010400L,
 175         0x00010400L, 0x01010000L, 0x01010000L, 0x01000404L,
 176         0x00010004L, 0x01000004L, 0x01000004L, 0x00010004L,
 177         0x00000000L, 0x00000404L, 0x00010404L, 0x01000000L,
 178         0x00010000L, 0x01010404L, 0x00000004L, 0x01010000L,
 179         0x01010400L, 0x01000000L, 0x01000000L, 0x00000400L,
 180         0x01010004L, 0x00010000L, 0x00010400L, 0x01000004L,
 181         0x00000400L, 0x00000004L, 0x01000404L, 0x00010404L,
 182         0x01010404L, 0x00010004L, 0x01010000L, 0x01000404L,
 183         0x01000004L, 0x00000404L, 0x00010404L, 0x01010400L,
 184         0x00000404L, 0x01000400L, 0x01000400L, 0x00000000L,
 185         0x00010004L, 0x00010400L, 0x00000000L, 0x01010004L };
 186
 187 static const unsigned long SP2[64] = {
 188         0x80108020L, 0x80008000L, 0x00008000L, 0x00108020L,
 189         0x00100000L, 0x00000020L, 0x80100020L, 0x80008020L,
 190         0x80000020L, 0x80108020L, 0x80108000L, 0x80000000L,
 191         0x80008000L, 0x00100000L, 0x00000020L, 0x80100020L,
 192         0x00108000L, 0x00100020L, 0x80008020L, 0x00000000L,
 193         0x80000000L, 0x00008000L, 0x00108020L, 0x80100000L,
 194         0x00100020L, 0x80000020L, 0x00000000L, 0x00108000L,
 195         0x00008020L, 0x80108000L, 0x80100000L, 0x00008020L,
 196         0x00000000L, 0x00108020L, 0x80100020L, 0x00100000L,
 197         0x80008020L, 0x80100000L, 0x80108000L, 0x00008000L,
 198         0x80100000L, 0x80008000L, 0x00000020L, 0x80108020L,
 199         0x00108020L, 0x00000020L, 0x00008000L, 0x80000000L,
 200         0x00008020L, 0x80108000L, 0x00100000L, 0x80000020L,
 201         0x00100020L, 0x80008020L, 0x80000020L, 0x00100020L,
 202         0x00108000L, 0x00000000L, 0x80008000L, 0x00008020L,
 203         0x80000000L, 0x80100020L, 0x80108020L, 0x00108000L };
 204
 205 static const unsigned long SP3[64] = {
 206         0x00000208L, 0x08020200L, 0x00000000L, 0x08020008L,
 207         0x08000200L, 0x00000000L, 0x00020208L, 0x08000200L,
 208         0x00020008L, 0x08000008L, 0x08000008L, 0x00020000L,
 209         0x08020208L, 0x00020008L, 0x08020000L, 0x00000208L,
 210         0x08000000L, 0x00000008L, 0x08020200L, 0x00000200L,
 211         0x00020200L, 0x08020000L, 0x08020008L, 0x00020208L,
 212         0x08000208L, 0x00020200L, 0x00020000L, 0x08000208L,
 213         0x00000008L, 0x08020208L, 0x00000200L, 0x08000000L,
 214         0x08020200L, 0x08000000L, 0x00020008L, 0x00000208L,
 215         0x00020000L, 0x08020200L, 0x08000200L, 0x00000000L,
 216         0x00000200L, 0x00020008L, 0x08020208L, 0x08000200L,
 217         0x08000008L, 0x00000200L, 0x00000000L, 0x08020008L,
 218         0x08000208L, 0x00020000L, 0x08000000L, 0x08020208L,
 219         0x00000008L, 0x00020208L, 0x00020200L, 0x08000008L,
 220         0x08020000L, 0x08000208L, 0x00000208L, 0x08020000L,
 221         0x00020208L, 0x00000008L, 0x08020008L, 0x00020200L };
 222
 223 static const unsigned long SP4[64] = {
 224         0x00802001L, 0x00002081L, 0x00002081L, 0x00000080L,
 225         0x00802080L, 0x00800081L, 0x00800001L, 0x00002001L,
 226         0x00000000L, 0x00802000L, 0x00802000L, 0x00802081L,
 227         0x00000081L, 0x00000000L, 0x00800080L, 0x00800001L,
 228         0x00000001L, 0x00002000L, 0x00800000L, 0x00802001L,
 229         0x00000080L, 0x00800000L, 0x00002001L, 0x00002080L,
 230         0x00800081L, 0x00000001L, 0x00002080L, 0x00800080L,
 231         0x00002000L, 0x00802080L, 0x00802081L, 0x00000081L,
 232         0x00800080L, 0x00800001L, 0x00802000L, 0x00802081L,
 233         0x00000081L, 0x00000000L, 0x00000000L, 0x00802000L,
 234         0x00002080L, 0x00800080L, 0x00800081L, 0x00000001L,
 235         0x00802001L, 0x00002081L, 0x00002081L, 0x00000080L,
 236         0x00802081L, 0x00000081L, 0x00000001L, 0x00002000L,
 237         0x00800001L, 0x00002001L, 0x00802080L, 0x00800081L,
 238         0x00002001L, 0x00002080L, 0x00800000L, 0x00802001L,
 239         0x00000080L, 0x00800000L, 0x00002000L, 0x00802080L };
 240
 241 static const unsigned long SP5[64] = {
 242         0x00000100L, 0x02080100L, 0x02080000L, 0x42000100L,
 243         0x00080000L, 0x00000100L, 0x40000000L, 0x02080000L,
 244         0x40080100L, 0x00080000L, 0x02000100L, 0x40080100L,
 245         0x42000100L, 0x42080000L, 0x00080100L, 0x40000000L,
 246         0x02000000L, 0x40080000L, 0x40080000L, 0x00000000L,
 247         0x40000100L, 0x42080100L, 0x42080100L, 0x02000100L,
 248         0x42080000L, 0x40000100L, 0x00000000L, 0x42000000L,
 249         0x02080100L, 0x02000000L, 0x42000000L, 0x00080100L,
 250         0x00080000L, 0x42000100L, 0x00000100L, 0x02000000L,
 251         0x40000000L, 0x02080000L, 0x42000100L, 0x40080100L,
 252         0x02000100L, 0x40000000L, 0x42080000L, 0x02080100L,
 253         0x40080100L, 0x00000100L, 0x02000000L, 0x42080000L,
 254         0x42080100L, 0x00080100L, 0x42000000L, 0x42080100L,
 255         0x02080000L, 0x00000000L, 0x40080000L, 0x42000000L,
 256         0x00080100L, 0x02000100L, 0x40000100L, 0x00080000L,
 257         0x00000000L, 0x40080000L, 0x02080100L, 0x40000100L };
 258
 259 static const unsigned long SP6[64] = {
 260         0x20000010L, 0x20400000L, 0x00004000L, 0x20404010L,
 261         0x20400000L, 0x00000010L, 0x20404010L, 0x00400000L,
 262         0x20004000L, 0x00404010L, 0x00400000L, 0x20000010L,
 263         0x00400010L, 0x20004000L, 0x20000000L, 0x00004010L,
 264         0x00000000L, 0x00400010L, 0x20004010L, 0x00004000L,
 265         0x00404000L, 0x20004010L, 0x00000010L, 0x20400010L,
 266         0x20400010L, 0x00000000L, 0x00404010L, 0x20404000L,
 267         0x00004010L, 0x00404000L, 0x20404000L, 0x20000000L,
 268         0x20004000L, 0x00000010L, 0x20400010L, 0x00404000L,
 269         0x20404010L, 0x00400000L, 0x00004010L, 0x20000010L,
 270         0x00400000L, 0x20004000L, 0x20000000L, 0x00004010L,
 271         0x20000010L, 0x20404010L, 0x00404000L, 0x20400000L,
 272         0x00404010L, 0x20404000L, 0x00000000L, 0x20400010L,
 273         0x00000010L, 0x00004000L, 0x20400000L, 0x00404010L,
 274         0x00004000L, 0x00400010L, 0x20004010L, 0x00000000L,
 275         0x20404000L, 0x20000000L, 0x00400010L, 0x20004010L };
 276
 277 static const unsigned long SP7[64] = {
 278         0x00200000L, 0x04200002L, 0x04000802L, 0x00000000L,
 279         0x00000800L, 0x04000802L, 0x00200802L, 0x04200800L,
 280         0x04200802L, 0x00200000L, 0x00000000L, 0x04000002L,
 281         0x00000002L, 0x04000000L, 0x04200002L, 0x00000802L,
 282         0x04000800L, 0x00200802L, 0x00200002L, 0x04000800L,
 283         0x04000002L, 0x04200000L, 0x04200800L, 0x00200002L,
 284         0x04200000L, 0x00000800L, 0x00000802L, 0x04200802L,
 285         0x00200800L, 0x00000002L, 0x04000000L, 0x00200800L,
 286         0x04000000L, 0x00200800L, 0x00200000L, 0x04000802L,
 287         0x04000802L, 0x04200002L, 0x04200002L, 0x00000002L,
 288         0x00200002L, 0x04000000L, 0x04000800L, 0x00200000L,
 289         0x04200800L, 0x00000802L, 0x00200802L, 0x04200800L,
 290         0x00000802L, 0x04000002L, 0x04200802L, 0x04200000L,
 291         0x00200800L, 0x00000000L, 0x00000002L, 0x04200802L,
 292         0x00000000L, 0x00200802L, 0x04200000L, 0x00000800L,
 293         0x04000002L, 0x04000800L, 0x00000800L, 0x00200002L };
 294
 295 static const unsigned long SP8[64] = {
 296         0x10001040L, 0x00001000L, 0x00040000L, 0x10041040L,
 297         0x10000000L, 0x10001040L, 0x00000040L, 0x10000000L,
 298         0x00040040L, 0x10040000L, 0x10041040L, 0x00041000L,
 299         0x10041000L, 0x00041040L, 0x00001000L, 0x00000040L,
 300         0x10040000L, 0x10000040L, 0x10001000L, 0x00001040L,
 301         0x00041000L, 0x00040040L, 0x10040040L, 0x10041000L,
 302         0x00001040L, 0x00000000L, 0x00000000L, 0x10040040L,
 303         0x10000040L, 0x10001000L, 0x00041040L, 0x00040000L,
 304         0x00041040L, 0x00040000L, 0x10041000L, 0x00001000L,
 305         0x00000040L, 0x10040040L, 0x00001000L, 0x00041040L,
 306         0x10001000L, 0x00000040L, 0x10000040L, 0x10040000L,
 307         0x10040040L, 0x10000000L, 0x00040000L, 0x10001040L,
 308         0x00000000L, 0x10041040L, 0x00040040L, 0x10000040L,
 309         0x10040000L, 0x10001000L, 0x10001040L, 0x00000000L,
 310         0x10041040L, 0x00041000L, 0x00041000L, 0x00001040L,
 311         0x00001040L, 0x00040040L, 0x10000000L, 0x10041000L };
 312
 313 static void desfunc(register unsigned long *block, register unsigned long *keys)
 314 {
 315         register unsigned long fval, work, right, leftt;
 316         register int round;
 317
 318         leftt = block[0];
 319         right = block[1];
 320         work = ((leftt >> 4) ^ right) & 0x0f0f0f0fL;
 321         right ^= work;
 322         leftt ^= (work << 4);
 323         work = ((leftt >> 16) ^ right) & 0x0000ffffL;
 324         right ^= work;
 325         leftt ^= (work << 16);
 326         work = ((right >> 2) ^ leftt) & 0x33333333L;
 327         leftt ^= work;
 328         right ^= (work << 2);
 329         work = ((right >> 8) ^ leftt) & 0x00ff00ffL;
 330         leftt ^= work;
 331         right ^= (work << 8);
 332         right = ((right << 1) | ((right >> 31) & 1L)) & 0xffffffffL;
 333         work = (leftt ^ right) & 0xaaaaaaaaL;
 334         leftt ^= work;
 335         right ^= work;
 336         leftt = ((leftt << 1) | ((leftt >> 31) & 1L)) & 0xffffffffL;
 337
 338         for( round = 0; round < 8; round++ ) {
 339                 work  = (right << 28) | (right >> 4);
 340                 work ^= *keys++;
 341                 fval  = SP7[ work                & 0x3fL];
 342                 fval |= SP5[(work >>  8) & 0x3fL];
 343                 fval |= SP3[(work >> 16) & 0x3fL];
 344                 fval |= SP1[(work >> 24) & 0x3fL];
 345                 work  = right ^ *keys++;
 346                 fval |= SP8[ work                & 0x3fL];
 347                 fval |= SP6[(work >>  8) & 0x3fL];
 348                 fval |= SP4[(work >> 16) & 0x3fL];
 349                 fval |= SP2[(work >> 24) & 0x3fL];
 350                 leftt ^= fval;
 351                 work  = (leftt << 28) | (leftt >> 4);
 352                 work ^= *keys++;
 353                 fval  = SP7[ work                & 0x3fL];
 354                 fval |= SP5[(work >>  8) & 0x3fL];
 355                 fval |= SP3[(work >> 16) & 0x3fL];
 356                 fval |= SP1[(work >> 24) & 0x3fL];
 357                 work  = leftt ^ *keys++;
 358                 fval |= SP8[ work                & 0x3fL];
 359                 fval |= SP6[(work >>  8) & 0x3fL];
 360                 fval |= SP4[(work >> 16) & 0x3fL];
 361                 fval |= SP2[(work >> 24) & 0x3fL];
 362                 right ^= fval;
 363                 }
 364
 365         right = (right << 31) | (right >> 1);
 366         work = (leftt ^ right) & 0xaaaaaaaaL;
 367         leftt ^= work;
 368         right ^= work;
 369         leftt = (leftt << 31) | (leftt >> 1);
 370         work = ((leftt >> 8) ^ right) & 0x00ff00ffL;
 371         right ^= work;
 372         leftt ^= (work << 8);
 373         work = ((leftt >> 2) ^ right) & 0x33333333L;
 374         right ^= work;
 375         leftt ^= (work << 2);
 376         work = ((right >> 16) ^ leftt) & 0x0000ffffL;
 377         leftt ^= work;
 378         right ^= (work << 16);
 379         work = ((right >> 4) ^ leftt) & 0x0f0f0f0fL;
 380         leftt ^= work;
 381         right ^= (work << 4);
 382         *block++ = right;
 383         *block = leftt;
 384         return;
 385         }
 386
 387 /* Validation sets:
 388  *
 389  * Single-length key, single-length plaintext -
 390  * Key    : 0123 4567 89ab cdef
 391  * Plain  : 0123 4567 89ab cde7
 392  * Cipher : c957 4425 6a5e d31d
 393  *
 394  * Double-length key, single-length plaintext -
 395  * Key    : 0123 4567 89ab cdef fedc ba98 7654 3210
 396  * Plain  : 0123 4567 89ab cde7
 397  * Cipher : 7f1d 0a77 826b 8aff
 398  *
 399  * Double-length key, double-length plaintext -
 400  * Key    : 0123 4567 89ab cdef fedc ba98 7654 3210
 401  * Plain  : 0123 4567 89ab cdef 0123 4567 89ab cdff
 402  * Cipher : 27a0 8440 406a df60 278f 47cf 42d6 15d7
 403  *
 404  * Triple-length key, single-length plaintext -
 405  * Key    : 0123 4567 89ab cdef fedc ba98 7654 3210 89ab cdef 0123 4567
 406  * Plain  : 0123 4567 89ab cde7
 407  * Cipher : de0b 7c06 ae5e 0ed5
 408  *
 409  * Triple-length key, double-length plaintext -
 410  * Key    : 0123 4567 89ab cdef fedc ba98 7654 3210 89ab cdef 0123 4567
 411  * Plain  : 0123 4567 89ab cdef 0123 4567 89ab cdff
 412  * Cipher : ad0d 1b30 ac17 cf07 0ed1 1c63 81e4 4de5
 413  *
 414  * d3des V5.0a rwo 9208.07 18:44 Graven Imagery
 415  **********************************************************************/