This is a live mirror of the Perl 5 development currently hosted at https://github.com/perl/perl5
Eliminate empty conditional branch
[perl5.git] / win32 / fcrypt.c
1 /* fcrypt.c */
2 /* Copyright (C) 1993 Eric Young - see README for more details */
3 #include <stdio.h>
4 #include <errno.h>
5
6 /* Eric Young.
7  * This version of crypt has been developed from my MIT compatable
8  * DES library.
9  * The library is available at pub/DES at ftp.psy.uq.oz.au
10  * eay@psych.psy.uq.oz.au
11  */
12
13 typedef unsigned char des_cblock[8];
14
15 typedef struct des_ks_struct
16         {
17         union   {
18                 des_cblock _;
19                 /* make sure things are correct size on machines with
20                  * 8 byte longs */
21                 unsigned long pad[2];
22                 } ks;
23 #define _       ks._
24         } des_key_schedule[16];
25
26 #define DES_KEY_SZ      (sizeof(des_cblock))
27 #define DES_ENCRYPT     1
28 #define DES_DECRYPT     0
29
30 #define ITERATIONS 16
31 #define HALF_ITERATIONS 8
32
33 #define c2l(c,l)        (l =((unsigned long)(*((c)++)))    , \
34                          l|=((unsigned long)(*((c)++)))<< 8, \
35                          l|=((unsigned long)(*((c)++)))<<16, \
36                          l|=((unsigned long)(*((c)++)))<<24)
37
38 #define l2c(l,c)        (*((c)++)=(unsigned char)(((l)    )&0xff), \
39                          *((c)++)=(unsigned char)(((l)>> 8)&0xff), \
40                          *((c)++)=(unsigned char)(((l)>>16)&0xff), \
41                          *((c)++)=(unsigned char)(((l)>>24)&0xff))
42
43 static const unsigned long SPtrans[8][64]={
44 { /* nibble 0 */
45 0x00820200, 0x00020000, 0x80800000, 0x80820200,
46 0x00800000, 0x80020200, 0x80020000, 0x80800000,
47 0x80020200, 0x00820200, 0x00820000, 0x80000200,
48 0x80800200, 0x00800000, 0x00000000, 0x80020000,
49 0x00020000, 0x80000000, 0x00800200, 0x00020200,
50 0x80820200, 0x00820000, 0x80000200, 0x00800200,
51 0x80000000, 0x00000200, 0x00020200, 0x80820000,
52 0x00000200, 0x80800200, 0x80820000, 0x00000000,
53 0x00000000, 0x80820200, 0x00800200, 0x80020000,
54 0x00820200, 0x00020000, 0x80000200, 0x00800200,
55 0x80820000, 0x00000200, 0x00020200, 0x80800000,
56 0x80020200, 0x80000000, 0x80800000, 0x00820000,
57 0x80820200, 0x00020200, 0x00820000, 0x80800200,
58 0x00800000, 0x80000200, 0x80020000, 0x00000000,
59 0x00020000, 0x00800000, 0x80800200, 0x00820200,
60 0x80000000, 0x80820000, 0x00000200, 0x80020200},
61 { /* nibble 1 */
62 0x10042004, 0x00000000, 0x00042000, 0x10040000,
63 0x10000004, 0x00002004, 0x10002000, 0x00042000,
64 0x00002000, 0x10040004, 0x00000004, 0x10002000,
65 0x00040004, 0x10042000, 0x10040000, 0x00000004,
66 0x00040000, 0x10002004, 0x10040004, 0x00002000,
67 0x00042004, 0x10000000, 0x00000000, 0x00040004,
68 0x10002004, 0x00042004, 0x10042000, 0x10000004,
69 0x10000000, 0x00040000, 0x00002004, 0x10042004,
70 0x00040004, 0x10042000, 0x10002000, 0x00042004,
71 0x10042004, 0x00040004, 0x10000004, 0x00000000,
72 0x10000000, 0x00002004, 0x00040000, 0x10040004,
73 0x00002000, 0x10000000, 0x00042004, 0x10002004,
74 0x10042000, 0x00002000, 0x00000000, 0x10000004,
75 0x00000004, 0x10042004, 0x00042000, 0x10040000,
76 0x10040004, 0x00040000, 0x00002004, 0x10002000,
77 0x10002004, 0x00000004, 0x10040000, 0x00042000},
78 { /* nibble 2 */
79 0x41000000, 0x01010040, 0x00000040, 0x41000040,
80 0x40010000, 0x01000000, 0x41000040, 0x00010040,
81 0x01000040, 0x00010000, 0x01010000, 0x40000000,
82 0x41010040, 0x40000040, 0x40000000, 0x41010000,
83 0x00000000, 0x40010000, 0x01010040, 0x00000040,
84 0x40000040, 0x41010040, 0x00010000, 0x41000000,
85 0x41010000, 0x01000040, 0x40010040, 0x01010000,
86 0x00010040, 0x00000000, 0x01000000, 0x40010040,
87 0x01010040, 0x00000040, 0x40000000, 0x00010000,
88 0x40000040, 0x40010000, 0x01010000, 0x41000040,
89 0x00000000, 0x01010040, 0x00010040, 0x41010000,
90 0x40010000, 0x01000000, 0x41010040, 0x40000000,
91 0x40010040, 0x41000000, 0x01000000, 0x41010040,
92 0x00010000, 0x01000040, 0x41000040, 0x00010040,
93 0x01000040, 0x00000000, 0x41010000, 0x40000040,
94 0x41000000, 0x40010040, 0x00000040, 0x01010000},
95 { /* nibble 3 */
96 0x00100402, 0x04000400, 0x00000002, 0x04100402,
97 0x00000000, 0x04100000, 0x04000402, 0x00100002,
98 0x04100400, 0x04000002, 0x04000000, 0x00000402,
99 0x04000002, 0x00100402, 0x00100000, 0x04000000,
100 0x04100002, 0x00100400, 0x00000400, 0x00000002,
101 0x00100400, 0x04000402, 0x04100000, 0x00000400,
102 0x00000402, 0x00000000, 0x00100002, 0x04100400,
103 0x04000400, 0x04100002, 0x04100402, 0x00100000,
104 0x04100002, 0x00000402, 0x00100000, 0x04000002,
105 0x00100400, 0x04000400, 0x00000002, 0x04100000,
106 0x04000402, 0x00000000, 0x00000400, 0x00100002,
107 0x00000000, 0x04100002, 0x04100400, 0x00000400,
108 0x04000000, 0x04100402, 0x00100402, 0x00100000,
109 0x04100402, 0x00000002, 0x04000400, 0x00100402,
110 0x00100002, 0x00100400, 0x04100000, 0x04000402,
111 0x00000402, 0x04000000, 0x04000002, 0x04100400},
112 { /* nibble 4 */
113 0x02000000, 0x00004000, 0x00000100, 0x02004108,
114 0x02004008, 0x02000100, 0x00004108, 0x02004000,
115 0x00004000, 0x00000008, 0x02000008, 0x00004100,
116 0x02000108, 0x02004008, 0x02004100, 0x00000000,
117 0x00004100, 0x02000000, 0x00004008, 0x00000108,
118 0x02000100, 0x00004108, 0x00000000, 0x02000008,
119 0x00000008, 0x02000108, 0x02004108, 0x00004008,
120 0x02004000, 0x00000100, 0x00000108, 0x02004100,
121 0x02004100, 0x02000108, 0x00004008, 0x02004000,
122 0x00004000, 0x00000008, 0x02000008, 0x02000100,
123 0x02000000, 0x00004100, 0x02004108, 0x00000000,
124 0x00004108, 0x02000000, 0x00000100, 0x00004008,
125 0x02000108, 0x00000100, 0x00000000, 0x02004108,
126 0x02004008, 0x02004100, 0x00000108, 0x00004000,
127 0x00004100, 0x02004008, 0x02000100, 0x00000108,
128 0x00000008, 0x00004108, 0x02004000, 0x02000008},
129 { /* nibble 5 */
130 0x20000010, 0x00080010, 0x00000000, 0x20080800,
131 0x00080010, 0x00000800, 0x20000810, 0x00080000,
132 0x00000810, 0x20080810, 0x00080800, 0x20000000,
133 0x20000800, 0x20000010, 0x20080000, 0x00080810,
134 0x00080000, 0x20000810, 0x20080010, 0x00000000,
135 0x00000800, 0x00000010, 0x20080800, 0x20080010,
136 0x20080810, 0x20080000, 0x20000000, 0x00000810,
137 0x00000010, 0x00080800, 0x00080810, 0x20000800,
138 0x00000810, 0x20000000, 0x20000800, 0x00080810,
139 0x20080800, 0x00080010, 0x00000000, 0x20000800,
140 0x20000000, 0x00000800, 0x20080010, 0x00080000,
141 0x00080010, 0x20080810, 0x00080800, 0x00000010,
142 0x20080810, 0x00080800, 0x00080000, 0x20000810,
143 0x20000010, 0x20080000, 0x00080810, 0x00000000,
144 0x00000800, 0x20000010, 0x20000810, 0x20080800,
145 0x20080000, 0x00000810, 0x00000010, 0x20080010},
146 { /* nibble 6 */
147 0x00001000, 0x00000080, 0x00400080, 0x00400001,
148 0x00401081, 0x00001001, 0x00001080, 0x00000000,
149 0x00400000, 0x00400081, 0x00000081, 0x00401000,
150 0x00000001, 0x00401080, 0x00401000, 0x00000081,
151 0x00400081, 0x00001000, 0x00001001, 0x00401081,
152 0x00000000, 0x00400080, 0x00400001, 0x00001080,
153 0x00401001, 0x00001081, 0x00401080, 0x00000001,
154 0x00001081, 0x00401001, 0x00000080, 0x00400000,
155 0x00001081, 0x00401000, 0x00401001, 0x00000081,
156 0x00001000, 0x00000080, 0x00400000, 0x00401001,
157 0x00400081, 0x00001081, 0x00001080, 0x00000000,
158 0x00000080, 0x00400001, 0x00000001, 0x00400080,
159 0x00000000, 0x00400081, 0x00400080, 0x00001080,
160 0x00000081, 0x00001000, 0x00401081, 0x00400000,
161 0x00401080, 0x00000001, 0x00001001, 0x00401081,
162 0x00400001, 0x00401080, 0x00401000, 0x00001001},
163 { /* nibble 7 */
164 0x08200020, 0x08208000, 0x00008020, 0x00000000,
165 0x08008000, 0x00200020, 0x08200000, 0x08208020,
166 0x00000020, 0x08000000, 0x00208000, 0x00008020,
167 0x00208020, 0x08008020, 0x08000020, 0x08200000,
168 0x00008000, 0x00208020, 0x00200020, 0x08008000,
169 0x08208020, 0x08000020, 0x00000000, 0x00208000,
170 0x08000000, 0x00200000, 0x08008020, 0x08200020,
171 0x00200000, 0x00008000, 0x08208000, 0x00000020,
172 0x00200000, 0x00008000, 0x08000020, 0x08208020,
173 0x00008020, 0x08000000, 0x00000000, 0x00208000,
174 0x08200020, 0x08008020, 0x08008000, 0x00200020,
175 0x08208000, 0x00000020, 0x00200020, 0x08008000,
176 0x08208020, 0x00200000, 0x08200000, 0x08000020,
177 0x00208000, 0x00008020, 0x08008020, 0x08200000,
178 0x00000020, 0x08208000, 0x00208020, 0x00000000,
179 0x08000000, 0x08200020, 0x00008000, 0x00208020}
180 };
181 static const unsigned long skb[8][64]={
182 { /* for C bits (numbered as per FIPS 46) 1 2 3 4 5 6 */
183 0x00000000,0x00000010,0x20000000,0x20000010,
184 0x00010000,0x00010010,0x20010000,0x20010010,
185 0x00000800,0x00000810,0x20000800,0x20000810,
186 0x00010800,0x00010810,0x20010800,0x20010810,
187 0x00000020,0x00000030,0x20000020,0x20000030,
188 0x00010020,0x00010030,0x20010020,0x20010030,
189 0x00000820,0x00000830,0x20000820,0x20000830,
190 0x00010820,0x00010830,0x20010820,0x20010830,
191 0x00080000,0x00080010,0x20080000,0x20080010,
192 0x00090000,0x00090010,0x20090000,0x20090010,
193 0x00080800,0x00080810,0x20080800,0x20080810,
194 0x00090800,0x00090810,0x20090800,0x20090810,
195 0x00080020,0x00080030,0x20080020,0x20080030,
196 0x00090020,0x00090030,0x20090020,0x20090030,
197 0x00080820,0x00080830,0x20080820,0x20080830,
198 0x00090820,0x00090830,0x20090820,0x20090830},
199 { /* for C bits (numbered as per FIPS 46) 7 8 10 11 12 13 */
200 0x00000000,0x02000000,0x00002000,0x02002000,
201 0x00200000,0x02200000,0x00202000,0x02202000,
202 0x00000004,0x02000004,0x00002004,0x02002004,
203 0x00200004,0x02200004,0x00202004,0x02202004,
204 0x00000400,0x02000400,0x00002400,0x02002400,
205 0x00200400,0x02200400,0x00202400,0x02202400,
206 0x00000404,0x02000404,0x00002404,0x02002404,
207 0x00200404,0x02200404,0x00202404,0x02202404,
208 0x10000000,0x12000000,0x10002000,0x12002000,
209 0x10200000,0x12200000,0x10202000,0x12202000,
210 0x10000004,0x12000004,0x10002004,0x12002004,
211 0x10200004,0x12200004,0x10202004,0x12202004,
212 0x10000400,0x12000400,0x10002400,0x12002400,
213 0x10200400,0x12200400,0x10202400,0x12202400,
214 0x10000404,0x12000404,0x10002404,0x12002404,
215 0x10200404,0x12200404,0x10202404,0x12202404},
216 { /* for C bits (numbered as per FIPS 46) 14 15 16 17 19 20 */
217 0x00000000,0x00000001,0x00040000,0x00040001,
218 0x01000000,0x01000001,0x01040000,0x01040001,
219 0x00000002,0x00000003,0x00040002,0x00040003,
220 0x01000002,0x01000003,0x01040002,0x01040003,
221 0x00000200,0x00000201,0x00040200,0x00040201,
222 0x01000200,0x01000201,0x01040200,0x01040201,
223 0x00000202,0x00000203,0x00040202,0x00040203,
224 0x01000202,0x01000203,0x01040202,0x01040203,
225 0x08000000,0x08000001,0x08040000,0x08040001,
226 0x09000000,0x09000001,0x09040000,0x09040001,
227 0x08000002,0x08000003,0x08040002,0x08040003,
228 0x09000002,0x09000003,0x09040002,0x09040003,
229 0x08000200,0x08000201,0x08040200,0x08040201,
230 0x09000200,0x09000201,0x09040200,0x09040201,
231 0x08000202,0x08000203,0x08040202,0x08040203,
232 0x09000202,0x09000203,0x09040202,0x09040203},
233 { /* for C bits (numbered as per FIPS 46) 21 23 24 26 27 28 */
234 0x00000000,0x00100000,0x00000100,0x00100100,
235 0x00000008,0x00100008,0x00000108,0x00100108,
236 0x00001000,0x00101000,0x00001100,0x00101100,
237 0x00001008,0x00101008,0x00001108,0x00101108,
238 0x04000000,0x04100000,0x04000100,0x04100100,
239 0x04000008,0x04100008,0x04000108,0x04100108,
240 0x04001000,0x04101000,0x04001100,0x04101100,
241 0x04001008,0x04101008,0x04001108,0x04101108,
242 0x00020000,0x00120000,0x00020100,0x00120100,
243 0x00020008,0x00120008,0x00020108,0x00120108,
244 0x00021000,0x00121000,0x00021100,0x00121100,
245 0x00021008,0x00121008,0x00021108,0x00121108,
246 0x04020000,0x04120000,0x04020100,0x04120100,
247 0x04020008,0x04120008,0x04020108,0x04120108,
248 0x04021000,0x04121000,0x04021100,0x04121100,
249 0x04021008,0x04121008,0x04021108,0x04121108},
250 { /* for D bits (numbered as per FIPS 46) 1 2 3 4 5 6 */
251 0x00000000,0x10000000,0x00010000,0x10010000,
252 0x00000004,0x10000004,0x00010004,0x10010004,
253 0x20000000,0x30000000,0x20010000,0x30010000,
254 0x20000004,0x30000004,0x20010004,0x30010004,
255 0x00100000,0x10100000,0x00110000,0x10110000,
256 0x00100004,0x10100004,0x00110004,0x10110004,
257 0x20100000,0x30100000,0x20110000,0x30110000,
258 0x20100004,0x30100004,0x20110004,0x30110004,
259 0x00001000,0x10001000,0x00011000,0x10011000,
260 0x00001004,0x10001004,0x00011004,0x10011004,
261 0x20001000,0x30001000,0x20011000,0x30011000,
262 0x20001004,0x30001004,0x20011004,0x30011004,
263 0x00101000,0x10101000,0x00111000,0x10111000,
264 0x00101004,0x10101004,0x00111004,0x10111004,
265 0x20101000,0x30101000,0x20111000,0x30111000,
266 0x20101004,0x30101004,0x20111004,0x30111004},
267 { /* for D bits (numbered as per FIPS 46) 8 9 11 12 13 14 */
268 0x00000000,0x08000000,0x00000008,0x08000008,
269 0x00000400,0x08000400,0x00000408,0x08000408,
270 0x00020000,0x08020000,0x00020008,0x08020008,
271 0x00020400,0x08020400,0x00020408,0x08020408,
272 0x00000001,0x08000001,0x00000009,0x08000009,
273 0x00000401,0x08000401,0x00000409,0x08000409,
274 0x00020001,0x08020001,0x00020009,0x08020009,
275 0x00020401,0x08020401,0x00020409,0x08020409,
276 0x02000000,0x0A000000,0x02000008,0x0A000008,
277 0x02000400,0x0A000400,0x02000408,0x0A000408,
278 0x02020000,0x0A020000,0x02020008,0x0A020008,
279 0x02020400,0x0A020400,0x02020408,0x0A020408,
280 0x02000001,0x0A000001,0x02000009,0x0A000009,
281 0x02000401,0x0A000401,0x02000409,0x0A000409,
282 0x02020001,0x0A020001,0x02020009,0x0A020009,
283 0x02020401,0x0A020401,0x02020409,0x0A020409},
284 { /* for D bits (numbered as per FIPS 46) 16 17 18 19 20 21 */
285 0x00000000,0x00000100,0x00080000,0x00080100,
286 0x01000000,0x01000100,0x01080000,0x01080100,
287 0x00000010,0x00000110,0x00080010,0x00080110,
288 0x01000010,0x01000110,0x01080010,0x01080110,
289 0x00200000,0x00200100,0x00280000,0x00280100,
290 0x01200000,0x01200100,0x01280000,0x01280100,
291 0x00200010,0x00200110,0x00280010,0x00280110,
292 0x01200010,0x01200110,0x01280010,0x01280110,
293 0x00000200,0x00000300,0x00080200,0x00080300,
294 0x01000200,0x01000300,0x01080200,0x01080300,
295 0x00000210,0x00000310,0x00080210,0x00080310,
296 0x01000210,0x01000310,0x01080210,0x01080310,
297 0x00200200,0x00200300,0x00280200,0x00280300,
298 0x01200200,0x01200300,0x01280200,0x01280300,
299 0x00200210,0x00200310,0x00280210,0x00280310,
300 0x01200210,0x01200310,0x01280210,0x01280310},
301 { /* for D bits (numbered as per FIPS 46) 22 23 24 25 27 28 */
302 0x00000000,0x04000000,0x00040000,0x04040000,
303 0x00000002,0x04000002,0x00040002,0x04040002,
304 0x00002000,0x04002000,0x00042000,0x04042000,
305 0x00002002,0x04002002,0x00042002,0x04042002,
306 0x00000020,0x04000020,0x00040020,0x04040020,
307 0x00000022,0x04000022,0x00040022,0x04040022,
308 0x00002020,0x04002020,0x00042020,0x04042020,
309 0x00002022,0x04002022,0x00042022,0x04042022,
310 0x00000800,0x04000800,0x00040800,0x04040800,
311 0x00000802,0x04000802,0x00040802,0x04040802,
312 0x00002800,0x04002800,0x00042800,0x04042800,
313 0x00002802,0x04002802,0x00042802,0x04042802,
314 0x00000820,0x04000820,0x00040820,0x04040820,
315 0x00000822,0x04000822,0x00040822,0x04040822,
316 0x00002820,0x04002820,0x00042820,0x04042820,
317 0x00002822,0x04002822,0x00042822,0x04042822}
318 };
319
320 /* See ecb_encrypt.c for a pseudo description of these macros. */
321 #define PERM_OP(a,b,t,n,m) ((t)=((((a)>>(n))^(b))&(m)),\
322         (b)^=(t),\
323         (a)^=((t)<<(n)))
324
325 #define HPERM_OP(a,t,n,m) ((t)=((((a)<<(16-(n)))^(a))&(m)),\
326         (a)=(a)^(t)^(t>>(16-(n))))\
327
328 static const char shifts2[16]={0,0,1,1,1,1,1,1,0,1,1,1,1,1,1,0};
329
330 static int body(
331         unsigned long *out0,
332         unsigned long *out1,
333         des_key_schedule ks,
334         unsigned long Eswap0,
335         unsigned long Eswap1);
336
337 static int
338 des_set_key(des_cblock *key, des_key_schedule schedule)
339         {
340         unsigned long c,d,t,s;
341         unsigned char *in;
342         unsigned long *k;
343         int i;
344
345         k=(unsigned long *)schedule;
346         in=(unsigned char *)key;
347
348         c2l(in,c);
349         c2l(in,d);
350
351         /* I now do it in 47 simple operations :-)
352          * Thanks to John Fletcher (john_fletcher@lccmail.ocf.llnl.gov)
353          * for the inspiration. :-) */
354         PERM_OP (d,c,t,4,0x0f0f0f0f);
355         HPERM_OP(c,t,-2,0xcccc0000);
356         HPERM_OP(d,t,-2,0xcccc0000);
357         PERM_OP (d,c,t,1,0x55555555);
358         PERM_OP (c,d,t,8,0x00ff00ff);
359         PERM_OP (d,c,t,1,0x55555555);
360         d=      (((d&0x000000ff)<<16)| (d&0x0000ff00)     |
361                  ((d&0x00ff0000)>>16)|((c&0xf0000000)>>4));
362         c&=0x0fffffff;
363
364         for (i=0; i<ITERATIONS; i++)
365                 {
366                 if (shifts2[i])
367                         { c=((c>>2)|(c<<26)); d=((d>>2)|(d<<26)); }
368                 else
369                         { c=((c>>1)|(c<<27)); d=((d>>1)|(d<<27)); }
370                 c&=0x0fffffff;
371                 d&=0x0fffffff;
372                 /* could be a few less shifts but I am to lazy at this
373                  * point in time to investigate */
374                 s=      skb[0][ (c    )&0x3f                ]|
375                         skb[1][((c>> 6)&0x03)|((c>> 7)&0x3c)]|
376                         skb[2][((c>>13)&0x0f)|((c>>14)&0x30)]|
377                         skb[3][((c>>20)&0x01)|((c>>21)&0x06) |
378                                               ((c>>22)&0x38)];
379                 t=      skb[4][ (d    )&0x3f                ]|
380                         skb[5][((d>> 7)&0x03)|((d>> 8)&0x3c)]|
381                         skb[6][ (d>>15)&0x3f                ]|
382                         skb[7][((d>>21)&0x0f)|((d>>22)&0x30)];
383
384                 /* table contained 0213 4657 */
385                 *(k++)=((t<<16)|(s&0x0000ffff))&0xffffffff;
386                 s=     ((s>>16)|(t&0xffff0000));
387                 
388                 s=(s<<4)|(s>>28);
389                 *(k++)=s&0xffffffff;
390                 }
391         return(0);
392         }
393
394 /******************************************************************
395  * modified stuff for crypt.
396  ******************************************************************/
397
398 /* The changes to this macro may help or hinder, depending on the
399  * compiler and the achitecture.  gcc2 always seems to do well :-). 
400  * Inspired by Dana How <how@isl.stanford.edu>
401  * DO NOT use the alternative version on machines with 8 byte longs.
402  */
403 #ifdef ALT_ECB
404 #define D_ENCRYPT(L,R,S) \
405         v=(R^(R>>16)); \
406         u=(v&E0); \
407         v=(v&E1); \
408         u=((u^(u<<16))^R^s[S  ])<<2; \
409         t=(v^(v<<16))^R^s[S+1]; \
410         t=(t>>2)|(t<<30); \
411         L^= \
412         *(unsigned long *)(des_SP+0x0100+((t    )&0xfc))+ \
413         *(unsigned long *)(des_SP+0x0300+((t>> 8)&0xfc))+ \
414         *(unsigned long *)(des_SP+0x0500+((t>>16)&0xfc))+ \
415         *(unsigned long *)(des_SP+0x0700+((t>>24)&0xfc))+ \
416         *(unsigned long *)(des_SP+       ((u    )&0xfc))+ \
417         *(unsigned long *)(des_SP+0x0200+((u>> 8)&0xfc))+ \
418         *(unsigned long *)(des_SP+0x0400+((u>>16)&0xfc))+ \
419         *(unsigned long *)(des_SP+0x0600+((u>>24)&0xfc));
420 #else /* original version */
421 #define D_ENCRYPT(L,R,S)        \
422         v=(R^(R>>16)); \
423         u=(v&E0); \
424         v=(v&E1); \
425         u=(u^(u<<16))^R^s[S  ]; \
426         t=(v^(v<<16))^R^s[S+1]; \
427         t=(t>>4)|(t<<28); \
428         L^=     SPtrans[1][(t    )&0x3f]| \
429                 SPtrans[3][(t>> 8)&0x3f]| \
430                 SPtrans[5][(t>>16)&0x3f]| \
431                 SPtrans[7][(t>>24)&0x3f]| \
432                 SPtrans[0][(u    )&0x3f]| \
433                 SPtrans[2][(u>> 8)&0x3f]| \
434                 SPtrans[4][(u>>16)&0x3f]| \
435                 SPtrans[6][(u>>24)&0x3f];
436 #endif
437
438 unsigned const char con_salt[128]={
439 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
440 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
441 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
442 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
443 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
444 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01,
445 0x02,0x03,0x04,0x05,0x06,0x07,0x08,0x09,
446 0x0A,0x0B,0x05,0x06,0x07,0x08,0x09,0x0A,
447 0x0B,0x0C,0x0D,0x0E,0x0F,0x10,0x11,0x12,
448 0x13,0x14,0x15,0x16,0x17,0x18,0x19,0x1A,
449 0x1B,0x1C,0x1D,0x1E,0x1F,0x20,0x21,0x22,
450 0x23,0x24,0x25,0x20,0x21,0x22,0x23,0x24,
451 0x25,0x26,0x27,0x28,0x29,0x2A,0x2B,0x2C,
452 0x2D,0x2E,0x2F,0x30,0x31,0x32,0x33,0x34,
453 0x35,0x36,0x37,0x38,0x39,0x3A,0x3B,0x3C,
454 0x3D,0x3E,0x3F,0x00,0x00,0x00,0x00,0x00,
455 };
456
457 unsigned const char cov_2char[64]={
458 0x2E,0x2F,0x30,0x31,0x32,0x33,0x34,0x35,
459 0x36,0x37,0x38,0x39,0x41,0x42,0x43,0x44,
460 0x45,0x46,0x47,0x48,0x49,0x4A,0x4B,0x4C,
461 0x4D,0x4E,0x4F,0x50,0x51,0x52,0x53,0x54,
462 0x55,0x56,0x57,0x58,0x59,0x5A,0x61,0x62,
463 0x63,0x64,0x65,0x66,0x67,0x68,0x69,0x6A,
464 0x6B,0x6C,0x6D,0x6E,0x6F,0x70,0x71,0x72,
465 0x73,0x74,0x75,0x76,0x77,0x78,0x79,0x7A
466 };
467
468 /* the salt for classic DES crypt (which is all we implement here)
469    permits [./0-9A-Za-z], since '.' and '/' immediately preceed
470    '0' we don't need individual checks for '.' and '/' 
471 */
472 #define good_for_salt(c) \
473     ((c) >= '.' && (c) <= '9' || (c) >= 'A' && (c) <= 'Z' ||  \
474      (c) >= 'a' && (c) <= 'z')
475
476 char *
477 des_fcrypt(const char *buf, const char *salt, char *buff)
478         {
479         unsigned int i,j,x,y;
480         unsigned long Eswap0,Eswap1;
481         unsigned long out[2],ll;
482         des_cblock key;
483         des_key_schedule ks;
484         unsigned char bb[9];
485         unsigned char *b=bb;
486         unsigned char c,u;
487
488         if (!good_for_salt(salt[0]) || !good_for_salt(salt[1])) {
489             errno = EINVAL;
490             return NULL;
491         }
492
493         /* eay 25/08/92
494          * If you call crypt("pwd","*") as often happens when you
495          * have * as the pwd field in /etc/passwd, the function
496          * returns *\0XXXXXXXXX
497          * The \0 makes the string look like * so the pwd "*" would
498          * crypt to "*".  This was found when replacing the crypt in
499          * our shared libraries.  People found that the disbled
500          * accounts effectivly had no passwd :-(. */
501         x=buff[0]=((salt[0] == '\0')?(char)'A':salt[0]);
502         Eswap0=con_salt[x];
503         x=buff[1]=((salt[1] == '\0')?(char)'A':salt[1]);
504         Eswap1=con_salt[x]<<4;
505
506         for (i=0; i<8; i++)
507                 {
508                 c= *(buf++);
509                 if (!c) break;
510                 key[i]=(char)(c<<1);
511                 }
512         for (; i<8; i++)
513                 key[i]=0;
514
515         des_set_key((des_cblock *)(key),ks);
516         body(&out[0],&out[1],ks,Eswap0,Eswap1);
517
518         ll=out[0]; l2c(ll,b);
519         ll=out[1]; l2c(ll,b);
520         y=0;
521         u=0x80;
522         bb[8]=0;
523         for (i=2; i<13; i++)
524                 {
525                 c=0;
526                 for (j=0; j<6; j++)
527                         {
528                         c<<=1;
529                         if (bb[y] & u) c|=1;
530                         u>>=1;
531                         if (!u)
532                                 {
533                                 y++;
534                                 u=0x80;
535                                 }
536                         }
537                 buff[i]=cov_2char[c];
538                 }
539         buff[13]='\0';
540         return buff;
541         }
542
543 static int 
544 body(   unsigned long *out0,
545         unsigned long *out1,
546         des_key_schedule ks,
547         unsigned long Eswap0,
548         unsigned long Eswap1)
549         {
550         unsigned long l,r,t,u,v;
551 #ifdef ALT_ECB
552         unsigned char *des_SP=(unsigned char *)SPtrans;
553 #endif
554         unsigned long *s;
555         int i,j;
556         unsigned long E0,E1;
557
558         l=0;
559         r=0;
560
561         s=(unsigned long *)ks;
562         E0=Eswap0;
563         E1=Eswap1;
564
565         for (j=0; j<25; j++)
566                 {
567                 for (i=0; i<(ITERATIONS*2); i+=4)
568                         {
569                         D_ENCRYPT(l,r,  i);     /*  1 */
570                         D_ENCRYPT(r,l,  i+2);   /*  2 */
571                         }
572                 t=l;
573                 l=r;
574                 r=t;
575                 }
576         t=r;
577         r=(l>>1)|(l<<31);
578         l=(t>>1)|(t<<31);
579         /* clear the top bits on machines with 8byte longs */
580         l&=0xffffffff;
581         r&=0xffffffff;
582
583         PERM_OP(r,l,t, 1,0x55555555);
584         PERM_OP(l,r,t, 8,0x00ff00ff);
585         PERM_OP(r,l,t, 2,0x33333333);
586         PERM_OP(l,r,t,16,0x0000ffff);
587         PERM_OP(r,l,t, 4,0x0f0f0f0f);
588
589         *out0=l;
590         *out1=r;
591         return(0);
592         }
593