sunrpc/des_impl.c

   1 /* Copyright (C) 1992 Eric Young */
   2 /* Collected from libdes and modified for SECURE RPC by Martin Kuck 1994 */
   3 /* This file is distributed under the terms of the GNU Lesser General */
   4 /* Public License, version 2.1 or later - see the file COPYING.LIB for details.*/
   5 /* If you did not receive a copy of the license with this program, please*/
   6 /* see <http://www.gnu.org/licenses/> to obtain a copy.  */
   7 #include <string.h>
   8 #include <stdint.h>
   9 #include "des.h"
  10
  11
  12 static const uint32_t des_SPtrans[8][64] =
  13 {
  14   {                             /* nibble 0 */
  15     0x00820200, 0x00020000, 0x80800000, 0x80820200,
  16     0x00800000, 0x80020200, 0x80020000, 0x80800000,
  17     0x80020200, 0x00820200, 0x00820000, 0x80000200,
  18     0x80800200, 0x00800000, 0x00000000, 0x80020000,
  19     0x00020000, 0x80000000, 0x00800200, 0x00020200,
  20     0x80820200, 0x00820000, 0x80000200, 0x00800200,
  21     0x80000000, 0x00000200, 0x00020200, 0x80820000,
  22     0x00000200, 0x80800200, 0x80820000, 0x00000000,
  23     0x00000000, 0x80820200, 0x00800200, 0x80020000,
  24     0x00820200, 0x00020000, 0x80000200, 0x00800200,
  25     0x80820000, 0x00000200, 0x00020200, 0x80800000,
  26     0x80020200, 0x80000000, 0x80800000, 0x00820000,
  27     0x80820200, 0x00020200, 0x00820000, 0x80800200,
  28     0x00800000, 0x80000200, 0x80020000, 0x00000000,
  29     0x00020000, 0x00800000, 0x80800200, 0x00820200,
  30     0x80000000, 0x80820000, 0x00000200, 0x80020200},
  31
  32   {                             /* nibble 1 */
  33     0x10042004, 0x00000000, 0x00042000, 0x10040000,
  34     0x10000004, 0x00002004, 0x10002000, 0x00042000,
  35     0x00002000, 0x10040004, 0x00000004, 0x10002000,
  36     0x00040004, 0x10042000, 0x10040000, 0x00000004,
  37     0x00040000, 0x10002004, 0x10040004, 0x00002000,
  38     0x00042004, 0x10000000, 0x00000000, 0x00040004,
  39     0x10002004, 0x00042004, 0x10042000, 0x10000004,
  40     0x10000000, 0x00040000, 0x00002004, 0x10042004,
  41     0x00040004, 0x10042000, 0x10002000, 0x00042004,
  42     0x10042004, 0x00040004, 0x10000004, 0x00000000,
  43     0x10000000, 0x00002004, 0x00040000, 0x10040004,
  44     0x00002000, 0x10000000, 0x00042004, 0x10002004,
  45     0x10042000, 0x00002000, 0x00000000, 0x10000004,
  46     0x00000004, 0x10042004, 0x00042000, 0x10040000,
  47     0x10040004, 0x00040000, 0x00002004, 0x10002000,
  48     0x10002004, 0x00000004, 0x10040000, 0x00042000},
  49
  50   {                             /* nibble 2 */
  51     0x41000000, 0x01010040, 0x00000040, 0x41000040,
  52     0x40010000, 0x01000000, 0x41000040, 0x00010040,
  53     0x01000040, 0x00010000, 0x01010000, 0x40000000,
  54     0x41010040, 0x40000040, 0x40000000, 0x41010000,
  55     0x00000000, 0x40010000, 0x01010040, 0x00000040,
  56     0x40000040, 0x41010040, 0x00010000, 0x41000000,
  57     0x41010000, 0x01000040, 0x40010040, 0x01010000,
  58     0x00010040, 0x00000000, 0x01000000, 0x40010040,
  59     0x01010040, 0x00000040, 0x40000000, 0x00010000,
  60     0x40000040, 0x40010000, 0x01010000, 0x41000040,
  61     0x00000000, 0x01010040, 0x00010040, 0x41010000,
  62     0x40010000, 0x01000000, 0x41010040, 0x40000000,
  63     0x40010040, 0x41000000, 0x01000000, 0x41010040,
  64     0x00010000, 0x01000040, 0x41000040, 0x00010040,
  65     0x01000040, 0x00000000, 0x41010000, 0x40000040,
  66     0x41000000, 0x40010040, 0x00000040, 0x01010000},
  67
  68   {                             /* nibble 3 */
  69     0x00100402, 0x04000400, 0x00000002, 0x04100402,
  70     0x00000000, 0x04100000, 0x04000402, 0x00100002,
  71     0x04100400, 0x04000002, 0x04000000, 0x00000402,
  72     0x04000002, 0x00100402, 0x00100000, 0x04000000,
  73     0x04100002, 0x00100400, 0x00000400, 0x00000002,
  74     0x00100400, 0x04000402, 0x04100000, 0x00000400,
  75     0x00000402, 0x00000000, 0x00100002, 0x04100400,
  76     0x04000400, 0x04100002, 0x04100402, 0x00100000,
  77     0x04100002, 0x00000402, 0x00100000, 0x04000002,
  78     0x00100400, 0x04000400, 0x00000002, 0x04100000,
  79     0x04000402, 0x00000000, 0x00000400, 0x00100002,
  80     0x00000000, 0x04100002, 0x04100400, 0x00000400,
  81     0x04000000, 0x04100402, 0x00100402, 0x00100000,
  82     0x04100402, 0x00000002, 0x04000400, 0x00100402,
  83     0x00100002, 0x00100400, 0x04100000, 0x04000402,
  84     0x00000402, 0x04000000, 0x04000002, 0x04100400},
  85
  86   {                             /* nibble 4 */
  87     0x02000000, 0x00004000, 0x00000100, 0x02004108,
  88     0x02004008, 0x02000100, 0x00004108, 0x02004000,
  89     0x00004000, 0x00000008, 0x02000008, 0x00004100,
  90     0x02000108, 0x02004008, 0x02004100, 0x00000000,
  91     0x00004100, 0x02000000, 0x00004008, 0x00000108,
  92     0x02000100, 0x00004108, 0x00000000, 0x02000008,
  93     0x00000008, 0x02000108, 0x02004108, 0x00004008,
  94     0x02004000, 0x00000100, 0x00000108, 0x02004100,
  95     0x02004100, 0x02000108, 0x00004008, 0x02004000,
  96     0x00004000, 0x00000008, 0x02000008, 0x02000100,
  97     0x02000000, 0x00004100, 0x02004108, 0x00000000,
  98     0x00004108, 0x02000000, 0x00000100, 0x00004008,
  99     0x02000108, 0x00000100, 0x00000000, 0x02004108,
 100     0x02004008, 0x02004100, 0x00000108, 0x00004000,
 101     0x00004100, 0x02004008, 0x02000100, 0x00000108,
 102     0x00000008, 0x00004108, 0x02004000, 0x02000008},
 103
 104   {                             /* nibble 5 */
 105     0x20000010, 0x00080010, 0x00000000, 0x20080800,
 106     0x00080010, 0x00000800, 0x20000810, 0x00080000,
 107     0x00000810, 0x20080810, 0x00080800, 0x20000000,
 108     0x20000800, 0x20000010, 0x20080000, 0x00080810,
 109     0x00080000, 0x20000810, 0x20080010, 0x00000000,
 110     0x00000800, 0x00000010, 0x20080800, 0x20080010,
 111     0x20080810, 0x20080000, 0x20000000, 0x00000810,
 112     0x00000010, 0x00080800, 0x00080810, 0x20000800,
 113     0x00000810, 0x20000000, 0x20000800, 0x00080810,
 114     0x20080800, 0x00080010, 0x00000000, 0x20000800,
 115     0x20000000, 0x00000800, 0x20080010, 0x00080000,
 116     0x00080010, 0x20080810, 0x00080800, 0x00000010,
 117     0x20080810, 0x00080800, 0x00080000, 0x20000810,
 118     0x20000010, 0x20080000, 0x00080810, 0x00000000,
 119     0x00000800, 0x20000010, 0x20000810, 0x20080800,
 120     0x20080000, 0x00000810, 0x00000010, 0x20080010},
 121
 122   {                             /* nibble 6 */
 123     0x00001000, 0x00000080, 0x00400080, 0x00400001,
 124     0x00401081, 0x00001001, 0x00001080, 0x00000000,
 125     0x00400000, 0x00400081, 0x00000081, 0x00401000,
 126     0x00000001, 0x00401080, 0x00401000, 0x00000081,
 127     0x00400081, 0x00001000, 0x00001001, 0x00401081,
 128     0x00000000, 0x00400080, 0x00400001, 0x00001080,
 129     0x00401001, 0x00001081, 0x00401080, 0x00000001,
 130     0x00001081, 0x00401001, 0x00000080, 0x00400000,
 131     0x00001081, 0x00401000, 0x00401001, 0x00000081,
 132     0x00001000, 0x00000080, 0x00400000, 0x00401001,
 133     0x00400081, 0x00001081, 0x00001080, 0x00000000,
 134     0x00000080, 0x00400001, 0x00000001, 0x00400080,
 135     0x00000000, 0x00400081, 0x00400080, 0x00001080,
 136     0x00000081, 0x00001000, 0x00401081, 0x00400000,
 137     0x00401080, 0x00000001, 0x00001001, 0x00401081,
 138     0x00400001, 0x00401080, 0x00401000, 0x00001001},
 139
 140   {                             /* nibble 7 */
 141     0x08200020, 0x08208000, 0x00008020, 0x00000000,
 142     0x08008000, 0x00200020, 0x08200000, 0x08208020,
 143     0x00000020, 0x08000000, 0x00208000, 0x00008020,
 144     0x00208020, 0x08008020, 0x08000020, 0x08200000,
 145     0x00008000, 0x00208020, 0x00200020, 0x08008000,
 146     0x08208020, 0x08000020, 0x00000000, 0x00208000,
 147     0x08000000, 0x00200000, 0x08008020, 0x08200020,
 148     0x00200000, 0x00008000, 0x08208000, 0x00000020,
 149     0x00200000, 0x00008000, 0x08000020, 0x08208020,
 150     0x00008020, 0x08000000, 0x00000000, 0x00208000,
 151     0x08200020, 0x08008020, 0x08008000, 0x00200020,
 152     0x08208000, 0x00000020, 0x00200020, 0x08008000,
 153     0x08208020, 0x00200000, 0x08200000, 0x08000020,
 154     0x00208000, 0x00008020, 0x08008020, 0x08200000,
 155     0x00000020, 0x08208000, 0x00208020, 0x00000000,
 156     0x08000000, 0x08200020, 0x00008000, 0x00208020}};
 157
 158 static const uint32_t des_skb[8][64] =
 159 {
 160   {                             /* for C bits (numbered as per FIPS 46) 1 2 3 4 5 6 */
 161     0x00000000, 0x00000010, 0x20000000, 0x20000010,
 162     0x00010000, 0x00010010, 0x20010000, 0x20010010,
 163     0x00000800, 0x00000810, 0x20000800, 0x20000810,
 164     0x00010800, 0x00010810, 0x20010800, 0x20010810,
 165     0x00000020, 0x00000030, 0x20000020, 0x20000030,
 166     0x00010020, 0x00010030, 0x20010020, 0x20010030,
 167     0x00000820, 0x00000830, 0x20000820, 0x20000830,
 168     0x00010820, 0x00010830, 0x20010820, 0x20010830,
 169     0x00080000, 0x00080010, 0x20080000, 0x20080010,
 170     0x00090000, 0x00090010, 0x20090000, 0x20090010,
 171     0x00080800, 0x00080810, 0x20080800, 0x20080810,
 172     0x00090800, 0x00090810, 0x20090800, 0x20090810,
 173     0x00080020, 0x00080030, 0x20080020, 0x20080030,
 174     0x00090020, 0x00090030, 0x20090020, 0x20090030,
 175     0x00080820, 0x00080830, 0x20080820, 0x20080830,
 176     0x00090820, 0x00090830, 0x20090820, 0x20090830},
 177   {                             /* for C bits (numbered as per FIPS 46) 7 8 10 11 12 13 */
 178     0x00000000, 0x02000000, 0x00002000, 0x02002000,
 179     0x00200000, 0x02200000, 0x00202000, 0x02202000,
 180     0x00000004, 0x02000004, 0x00002004, 0x02002004,
 181     0x00200004, 0x02200004, 0x00202004, 0x02202004,
 182     0x00000400, 0x02000400, 0x00002400, 0x02002400,
 183     0x00200400, 0x02200400, 0x00202400, 0x02202400,
 184     0x00000404, 0x02000404, 0x00002404, 0x02002404,
 185     0x00200404, 0x02200404, 0x00202404, 0x02202404,
 186     0x10000000, 0x12000000, 0x10002000, 0x12002000,
 187     0x10200000, 0x12200000, 0x10202000, 0x12202000,
 188     0x10000004, 0x12000004, 0x10002004, 0x12002004,
 189     0x10200004, 0x12200004, 0x10202004, 0x12202004,
 190     0x10000400, 0x12000400, 0x10002400, 0x12002400,
 191     0x10200400, 0x12200400, 0x10202400, 0x12202400,
 192     0x10000404, 0x12000404, 0x10002404, 0x12002404,
 193     0x10200404, 0x12200404, 0x10202404, 0x12202404},
 194   {                             /* for C bits (numbered as per FIPS 46) 14 15 16 17 19 20 */
 195     0x00000000, 0x00000001, 0x00040000, 0x00040001,
 196     0x01000000, 0x01000001, 0x01040000, 0x01040001,
 197     0x00000002, 0x00000003, 0x00040002, 0x00040003,
 198     0x01000002, 0x01000003, 0x01040002, 0x01040003,
 199     0x00000200, 0x00000201, 0x00040200, 0x00040201,
 200     0x01000200, 0x01000201, 0x01040200, 0x01040201,
 201     0x00000202, 0x00000203, 0x00040202, 0x00040203,
 202     0x01000202, 0x01000203, 0x01040202, 0x01040203,
 203     0x08000000, 0x08000001, 0x08040000, 0x08040001,
 204     0x09000000, 0x09000001, 0x09040000, 0x09040001,
 205     0x08000002, 0x08000003, 0x08040002, 0x08040003,
 206     0x09000002, 0x09000003, 0x09040002, 0x09040003,
 207     0x08000200, 0x08000201, 0x08040200, 0x08040201,
 208     0x09000200, 0x09000201, 0x09040200, 0x09040201,
 209     0x08000202, 0x08000203, 0x08040202, 0x08040203,
 210     0x09000202, 0x09000203, 0x09040202, 0x09040203},
 211   {                             /* for C bits (numbered as per FIPS 46) 21 23 24 26 27 28 */
 212     0x00000000, 0x00100000, 0x00000100, 0x00100100,
 213     0x00000008, 0x00100008, 0x00000108, 0x00100108,
 214     0x00001000, 0x00101000, 0x00001100, 0x00101100,
 215     0x00001008, 0x00101008, 0x00001108, 0x00101108,
 216     0x04000000, 0x04100000, 0x04000100, 0x04100100,
 217     0x04000008, 0x04100008, 0x04000108, 0x04100108,
 218     0x04001000, 0x04101000, 0x04001100, 0x04101100,
 219     0x04001008, 0x04101008, 0x04001108, 0x04101108,
 220     0x00020000, 0x00120000, 0x00020100, 0x00120100,
 221     0x00020008, 0x00120008, 0x00020108, 0x00120108,
 222     0x00021000, 0x00121000, 0x00021100, 0x00121100,
 223     0x00021008, 0x00121008, 0x00021108, 0x00121108,
 224     0x04020000, 0x04120000, 0x04020100, 0x04120100,
 225     0x04020008, 0x04120008, 0x04020108, 0x04120108,
 226     0x04021000, 0x04121000, 0x04021100, 0x04121100,
 227     0x04021008, 0x04121008, 0x04021108, 0x04121108},
 228   {                             /* for D bits (numbered as per FIPS 46) 1 2 3 4 5 6 */
 229     0x00000000, 0x10000000, 0x00010000, 0x10010000,
 230     0x00000004, 0x10000004, 0x00010004, 0x10010004,
 231     0x20000000, 0x30000000, 0x20010000, 0x30010000,
 232     0x20000004, 0x30000004, 0x20010004, 0x30010004,
 233     0x00100000, 0x10100000, 0x00110000, 0x10110000,
 234     0x00100004, 0x10100004, 0x00110004, 0x10110004,
 235     0x20100000, 0x30100000, 0x20110000, 0x30110000,
 236     0x20100004, 0x30100004, 0x20110004, 0x30110004,
 237     0x00001000, 0x10001000, 0x00011000, 0x10011000,
 238     0x00001004, 0x10001004, 0x00011004, 0x10011004,
 239     0x20001000, 0x30001000, 0x20011000, 0x30011000,
 240     0x20001004, 0x30001004, 0x20011004, 0x30011004,
 241     0x00101000, 0x10101000, 0x00111000, 0x10111000,
 242     0x00101004, 0x10101004, 0x00111004, 0x10111004,
 243     0x20101000, 0x30101000, 0x20111000, 0x30111000,
 244     0x20101004, 0x30101004, 0x20111004, 0x30111004},
 245   {                             /* for D bits (numbered as per FIPS 46) 8 9 11 12 13 14 */
 246     0x00000000, 0x08000000, 0x00000008, 0x08000008,
 247     0x00000400, 0x08000400, 0x00000408, 0x08000408,
 248     0x00020000, 0x08020000, 0x00020008, 0x08020008,
 249     0x00020400, 0x08020400, 0x00020408, 0x08020408,
 250     0x00000001, 0x08000001, 0x00000009, 0x08000009,
 251     0x00000401, 0x08000401, 0x00000409, 0x08000409,
 252     0x00020001, 0x08020001, 0x00020009, 0x08020009,
 253     0x00020401, 0x08020401, 0x00020409, 0x08020409,
 254     0x02000000, 0x0A000000, 0x02000008, 0x0A000008,
 255     0x02000400, 0x0A000400, 0x02000408, 0x0A000408,
 256     0x02020000, 0x0A020000, 0x02020008, 0x0A020008,
 257     0x02020400, 0x0A020400, 0x02020408, 0x0A020408,
 258     0x02000001, 0x0A000001, 0x02000009, 0x0A000009,
 259     0x02000401, 0x0A000401, 0x02000409, 0x0A000409,
 260     0x02020001, 0x0A020001, 0x02020009, 0x0A020009,
 261     0x02020401, 0x0A020401, 0x02020409, 0x0A020409},
 262   {                             /* for D bits (numbered as per FIPS 46) 16 17 18 19 20 21 */
 263     0x00000000, 0x00000100, 0x00080000, 0x00080100,
 264     0x01000000, 0x01000100, 0x01080000, 0x01080100,
 265     0x00000010, 0x00000110, 0x00080010, 0x00080110,
 266     0x01000010, 0x01000110, 0x01080010, 0x01080110,
 267     0x00200000, 0x00200100, 0x00280000, 0x00280100,
 268     0x01200000, 0x01200100, 0x01280000, 0x01280100,
 269     0x00200010, 0x00200110, 0x00280010, 0x00280110,
 270     0x01200010, 0x01200110, 0x01280010, 0x01280110,
 271     0x00000200, 0x00000300, 0x00080200, 0x00080300,
 272     0x01000200, 0x01000300, 0x01080200, 0x01080300,
 273     0x00000210, 0x00000310, 0x00080210, 0x00080310,
 274     0x01000210, 0x01000310, 0x01080210, 0x01080310,
 275     0x00200200, 0x00200300, 0x00280200, 0x00280300,
 276     0x01200200, 0x01200300, 0x01280200, 0x01280300,
 277     0x00200210, 0x00200310, 0x00280210, 0x00280310,
 278     0x01200210, 0x01200310, 0x01280210, 0x01280310},
 279   {                             /* for D bits (numbered as per FIPS 46) 22 23 24 25 27 28 */
 280     0x00000000, 0x04000000, 0x00040000, 0x04040000,
 281     0x00000002, 0x04000002, 0x00040002, 0x04040002,
 282     0x00002000, 0x04002000, 0x00042000, 0x04042000,
 283     0x00002002, 0x04002002, 0x00042002, 0x04042002,
 284     0x00000020, 0x04000020, 0x00040020, 0x04040020,
 285     0x00000022, 0x04000022, 0x00040022, 0x04040022,
 286     0x00002020, 0x04002020, 0x00042020, 0x04042020,
 287     0x00002022, 0x04002022, 0x00042022, 0x04042022,
 288     0x00000800, 0x04000800, 0x00040800, 0x04040800,
 289     0x00000802, 0x04000802, 0x00040802, 0x04040802,
 290     0x00002800, 0x04002800, 0x00042800, 0x04042800,
 291     0x00002802, 0x04002802, 0x00042802, 0x04042802,
 292     0x00000820, 0x04000820, 0x00040820, 0x04040820,
 293     0x00000822, 0x04000822, 0x00040822, 0x04040822,
 294     0x00002820, 0x04002820, 0x00042820, 0x04042820,
 295     0x00002822, 0x04002822, 0x00042822, 0x04042822},
 296 };
 297
 298 #define c2l(c,l)        (l =((unsigned long)(*((c)++)))    , \
 299                          l|=((unsigned long)(*((c)++)))<< 8, \
 300                          l|=((unsigned long)(*((c)++)))<<16, \
 301                          l|=((unsigned long)(*((c)++)))<<24)
 302
 303 #define l2c(l,c)        (*((c)++)=(unsigned char)(((l)    )&0xff), \
 304                          *((c)++)=(unsigned char)(((l)>> 8)&0xff), \
 305                          *((c)++)=(unsigned char)(((l)>>16)&0xff), \
 306                          *((c)++)=(unsigned char)(((l)>>24)&0xff))
 307
 308 /*
 309  * IP and FP
 310  * The problem is more of a geometric problem that random bit fiddling.
 311  *  0  1  2  3  4  5  6  7      62 54 46 38 30 22 14  6
 312  *  8  9 10 11 12 13 14 15      60 52 44 36 28 20 12  4
 313  * 16 17 18 19 20 21 22 23      58 50 42 34 26 18 10  2
 314  * 24 25 26 27 28 29 30 31  to  56 48 40 32 24 16  8  0
 315  *
 316  * 32 33 34 35 36 37 38 39      63 55 47 39 31 23 15  7
 317  * 40 41 42 43 44 45 46 47      61 53 45 37 29 21 13  5
 318  * 48 49 50 51 52 53 54 55      59 51 43 35 27 19 11  3
 319  * 56 57 58 59 60 61 62 63      57 49 41 33 25 17  9  1
 320  *
 321  * The output has been subject to swaps of the form
 322  * 0 1 -> 3 1 but the odd and even bits have been put into
 323  * 2 3    2 0
 324  * different words.  The main trick is to remember that
 325  * t=((l>>size)^r)&(mask);
 326  * r^=t;
 327  * l^=(t<<size);
 328  * can be used to swap and move bits between words.
 329  *
 330  * So l =  0  1  2  3  r = 16 17 18 19
 331  *         4  5  6  7      20 21 22 23
 332  *         8  9 10 11      24 25 26 27
 333  *        12 13 14 15      28 29 30 31
 334  * becomes (for size == 2 and mask == 0x3333)
 335  * t =   2^16  3^17 -- --   l =  0  1 16 17  r =  2  3 18 19
 336  *       6^20  7^21 -- --        4  5 20 21       6  7 22 23
 337  *      10^24 11^25 -- --        8  9 24 25      10 11 24 25
 338  *      14^28 15^29 -- --       12 13 28 29      14 15 28 29
 339  *
 340  * Thanks for hints from Richard Outerbridge - he told me IP&FP
 341  * could be done in 15 xor, 10 shifts and 5 ands.
 342  * When I finally started to think of the problem in 2D
 343  * I first got ~42 operations without xors.  When I remembered
 344  * how to use xors :-) I got it to its final state.
 345  */
 346
 347 #define PERM_OP(a,b,t,n,m) ((t)=((((a)>>(n))^(b))&(m)),\
 348         (b)^=(t),\
 349         (a)^=((t)<<(n)))
 350
 351 #define HPERM_OP(a,t,n,m) ((t)=((((a)<<(16-(n)))^(a))&(m)),\
 352         (a)=(a)^(t)^(t>>(16-(n))))
 353
 354
 355 #define D_ENCRYPT(L,R,S)        \
 356         u=(R^s[S  ]); \
 357         t=R^s[S+1]; \
 358         t=((t>>4)+(t<<28)); \
 359         L^=     des_SPtrans[1][(t    )&0x3f]| \
 360                 des_SPtrans[3][(t>> 8)&0x3f]| \
 361                 des_SPtrans[5][(t>>16)&0x3f]| \
 362                 des_SPtrans[7][(t>>24)&0x3f]| \
 363                 des_SPtrans[0][(u    )&0x3f]| \
 364                 des_SPtrans[2][(u>> 8)&0x3f]| \
 365                 des_SPtrans[4][(u>>16)&0x3f]| \
 366                 des_SPtrans[6][(u>>24)&0x3f];
 367
 368 #define ITERATIONS 16
 369
 370 static const char shifts2[16] =
 371 {0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0};
 372
 373 static void des_set_key (unsigned char *, unsigned long *) internal_function;
 374 static void des_encrypt (unsigned long *, unsigned long *, int)
 375      internal_function;
 376 int _des_crypt (char *, unsigned, struct desparams *);
 377
 378 static void
 379 internal_function
 380 des_set_key (unsigned char *key, unsigned long *schedule)
 381 {
 382   register unsigned long c, d, t, s;
 383   register unsigned char *in;
 384   register unsigned long *k;
 385   register int i;
 386
 387   k = (unsigned long *) schedule;
 388   in = key;
 389
 390   c2l (in, c);
 391   c2l (in, d);
 392
 393   /* I now do it in 47 simple operations :-)
 394    * Thanks to John Fletcher (john_fletcher@lccmail.ocf.llnl.gov)
 395    * for the inspiration. :-) */
 396   PERM_OP (d, c, t, 4, 0x0f0f0f0f);
 397   HPERM_OP (c, t, -2, 0xcccc0000);
 398   HPERM_OP (d, t, -2, 0xcccc0000);
 399   PERM_OP (d, c, t, 1, 0x55555555);
 400   PERM_OP (c, d, t, 8, 0x00ff00ff);
 401   PERM_OP (d, c, t, 1, 0x55555555);
 402   d = (((d & 0x000000ff) << 16) | (d & 0x0000ff00) |
 403        ((d & 0x00ff0000) >> 16) | ((c & 0xf0000000) >> 4));
 404   c &= 0x0fffffff;
 405
 406   for (i = 0; i < ITERATIONS; i++)
 407     {
 408       if (shifts2[i])
 409         {
 410           c = ((c >> 2) | (c << 26));
 411           d = ((d >> 2) | (d << 26));
 412         }
 413       else
 414         {
 415           c = ((c >> 1) | (c << 27));
 416           d = ((d >> 1) | (d << 27));
 417         }
 418       c &= 0x0fffffff;
 419       d &= 0x0fffffff;
 420       /* could be a few less shifts but I am to lazy at this
 421        * point in time to investigate */
 422       s = des_skb[0][(c) & 0x3f] |
 423         des_skb[1][((c >> 6) & 0x03) | ((c >> 7) & 0x3c)] |
 424         des_skb[2][((c >> 13) & 0x0f) | ((c >> 14) & 0x30)] |
 425         des_skb[3][((c >> 20) & 0x01) | ((c >> 21) & 0x06) | ((c >> 22) & 0x38)];
 426       t = des_skb[4][(d) & 0x3f] |
 427         des_skb[5][((d >> 7) & 0x03) | ((d >> 8) & 0x3c)] |
 428         des_skb[6][(d >> 15) & 0x3f] |
 429         des_skb[7][((d >> 21) & 0x0f) | ((d >> 22) & 0x30)];
 430
 431       /* table contained 0213 4657 */
 432       *(k++) = ((t << 16) | (s & 0x0000ffff)) & 0xffffffff;
 433       s = ((s >> 16) | (t & 0xffff0000));
 434
 435       s = (s << 4) | (s >> 28);
 436       *(k++) = s & 0xffffffff;
 437     }
 438 }
 439
 440
 441 static void
 442 internal_function
 443 des_encrypt (unsigned long *buf, unsigned long *schedule, int encrypt)
 444 {
 445   register unsigned long l, r, t, u;
 446   register int i;
 447   register unsigned long *s;
 448
 449   l = buf[0];
 450   r = buf[1];
 451
 452   /* do IP */
 453   PERM_OP (r, l, t, 4, 0x0f0f0f0f);
 454   PERM_OP (l, r, t, 16, 0x0000ffff);
 455   PERM_OP (r, l, t, 2, 0x33333333);
 456   PERM_OP (l, r, t, 8, 0x00ff00ff);
 457   PERM_OP (r, l, t, 1, 0x55555555);
 458   /* r and l are reversed - remember that :-) - fix
 459    * it in the next step */
 460
 461   /* Things have been modified so that the initial rotate is
 462    * done outside the loop.  This required the
 463    * des_SPtrans values in sp.h to be rotated 1 bit to the right.
 464    * One perl script later and things have a 5% speed up on a sparc2.
 465    * Thanks to Richard Outerbridge <71755.204@CompuServe.COM>
 466    * for pointing this out. */
 467   t = (r << 1) | (r >> 31);
 468   r = (l << 1) | (l >> 31);
 469   l = t;
 470
 471   /* clear the top bits on machines with 8byte longs */
 472   l &= 0xffffffff;
 473   r &= 0xffffffff;
 474
 475   s = (unsigned long *) schedule;
 476   /* I don't know if it is worth the effort of loop unrolling the
 477    * inner loop */
 478   if (encrypt)
 479     {
 480       for (i = 0; i < 32; i += 4)
 481         {
 482           D_ENCRYPT (l, r, i + 0);      /*  1 */
 483           D_ENCRYPT (r, l, i + 2);      /*  2 */
 484         }
 485     }
 486   else
 487     {
 488       for (i = 30; i > 0; i -= 4)
 489         {
 490           D_ENCRYPT (l, r, i - 0);      /* 16 */
 491           D_ENCRYPT (r, l, i - 2);      /* 15 */
 492         }
 493     }
 494   l = (l >> 1) | (l << 31);
 495   r = (r >> 1) | (r << 31);
 496   /* clear the top bits on machines with 8byte longs */
 497   l &= 0xffffffff;
 498   r &= 0xffffffff;
 499
 500   /* swap l and r
 501    * we will not do the swap so just remember they are
 502    * reversed for the rest of the subroutine
 503    * luckily FP fixes this problem :-) */
 504
 505   PERM_OP (r, l, t, 1, 0x55555555);
 506   PERM_OP (l, r, t, 8, 0x00ff00ff);
 507   PERM_OP (r, l, t, 2, 0x33333333);
 508   PERM_OP (l, r, t, 16, 0x0000ffff);
 509   PERM_OP (r, l, t, 4, 0x0f0f0f0f);
 510
 511   buf[0] = l;
 512   buf[1] = r;
 513
 514   l = r = t = u = 0;
 515 }
 516
 517
 518 int
 519 _des_crypt (char *buf, unsigned len, struct desparams *desp)
 520 {
 521   unsigned long schedule[32];
 522   register unsigned long tin0, tin1;
 523   register unsigned long tout0, tout1, xor0, xor1;
 524   register unsigned char *in, *out;
 525   unsigned long tbuf[2];
 526   unsigned char *iv, *oiv;
 527   int cbc_mode;
 528
 529   cbc_mode = (desp->des_mode == CBC) ? 1 : 0;
 530
 531   in = (unsigned char *) buf;
 532   out = (unsigned char *) buf;
 533   oiv = iv = (unsigned char *) desp->des_ivec;
 534
 535   des_set_key (desp->des_key, schedule);
 536
 537   tin0 = tin1 = 0;              /* For GCC */
 538   if (desp->des_dir == ENCRYPT)
 539     {
 540       c2l (iv, tout0);
 541       c2l (iv, tout1);
 542       for (; len > 0; len -= 8)
 543         {
 544           c2l (in, tin0);
 545           c2l (in, tin1);
 546           if (cbc_mode)
 547             {
 548               tin0 ^= tout0;
 549               tin1 ^= tout1;
 550             }
 551           tbuf[0] = tin0;
 552           tbuf[1] = tin1;
 553           des_encrypt (tbuf, schedule, 1);
 554           tout0 = tbuf[0];
 555           tout1 = tbuf[1];
 556           l2c (tout0, out);
 557           l2c (tout1, out);
 558         }
 559       l2c (tout0, oiv);
 560       l2c (tout1, oiv);
 561     }
 562   else
 563     {
 564       c2l (iv, xor0);
 565       c2l (iv, xor1);
 566       for (; len > 0; len -= 8)
 567         {
 568           c2l (in, tin0);
 569           c2l (in, tin1);
 570           tbuf[0] = tin0;
 571           tbuf[1] = tin1;
 572           des_encrypt (tbuf, schedule, 0);
 573           if (cbc_mode)
 574             {
 575               tout0 = tbuf[0] ^ xor0;
 576               tout1 = tbuf[1] ^ xor1;
 577               xor0 = tin0;
 578               xor1 = tin1;
 579             }
 580           else
 581             {
 582               tout0 = tbuf[0];
 583               tout1 = tbuf[1];
 584             }
 585           l2c (tout0, out);
 586           l2c (tout1, out);
 587         }
 588       l2c (tin0, oiv);
 589       l2c (tin1, oiv);
 590     }
 591   tin0 = tin1 = tout0 = tout1 = xor0 = xor1 = 0;
 592   tbuf[0] = tbuf[1] = 0;
 593   __bzero (schedule, sizeof (schedule));
 594
 595   return (1);
 596 }