lh | 9ed821d | 2023-04-07 01:36:19 -0700 | [diff] [blame] | 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 | } |