des_locl.h 16 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441
  1. /* crypto/des/des_locl.h */
  2. /* Copyright (C) 1995-1997 Eric Young ([email protected])
  3. * All rights reserved.
  4. *
  5. * This package is an SSL implementation written
  6. * by Eric Young ([email protected]).
  7. * The implementation was written so as to conform with Netscapes SSL.
  8. *
  9. * This library is free for commercial and non-commercial use as long as
  10. * the following conditions are aheared to. The following conditions
  11. * apply to all code found in this distribution, be it the RC4, RSA,
  12. * lhash, DES, etc., code; not just the SSL code. The SSL documentation
  13. * included with this distribution is covered by the same copyright terms
  14. * except that the holder is Tim Hudson ([email protected]).
  15. *
  16. * Copyright remains Eric Young's, and as such any Copyright notices in
  17. * the code are not to be removed.
  18. * If this package is used in a product, Eric Young should be given attribution
  19. * as the author of the parts of the library used.
  20. * This can be in the form of a textual message at program startup or
  21. * in documentation (online or textual) provided with the package.
  22. *
  23. * Redistribution and use in source and binary forms, with or without
  24. * modification, are permitted provided that the following conditions
  25. * are met:
  26. * 1. Redistributions of source code must retain the copyright
  27. * notice, this list of conditions and the following disclaimer.
  28. * 2. Redistributions in binary form must reproduce the above copyright
  29. * notice, this list of conditions and the following disclaimer in the
  30. * documentation and/or other materials provided with the distribution.
  31. * 3. All advertising materials mentioning features or use of this software
  32. * must display the following acknowledgement:
  33. * "This product includes cryptographic software written by
  34. * Eric Young ([email protected])"
  35. * The word 'cryptographic' can be left out if the rouines from the library
  36. * being used are not cryptographic related :-).
  37. * 4. If you include any Windows specific code (or a derivative thereof) from
  38. * the apps directory (application code) you must include an acknowledgement:
  39. * "This product includes software written by Tim Hudson ([email protected])"
  40. *
  41. * THIS SOFTWARE IS PROVIDED BY ERIC YOUNG ``AS IS'' AND
  42. * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  43. * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  44. * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
  45. * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  46. * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  47. * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  48. * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  49. * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  50. * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  51. * SUCH DAMAGE.
  52. *
  53. * The licence and distribution terms for any publically available version or
  54. * derivative of this code cannot be changed. i.e. this code cannot simply be
  55. * copied and put under another distribution licence
  56. * [including the GNU Public Licence.]
  57. */
  58. #ifndef HEADER_DES_LOCL_H
  59. # define HEADER_DES_LOCL_H
  60. # include <openssl/e_os2.h>
  61. # if defined(OPENSSL_SYS_WIN32)
  62. # ifndef OPENSSL_SYS_MSDOS
  63. # define OPENSSL_SYS_MSDOS
  64. # endif
  65. # endif
  66. # include <stdio.h>
  67. # include <stdlib.h>
  68. # ifndef OPENSSL_SYS_MSDOS
  69. # if !defined(OPENSSL_SYS_VMS) || defined(__DECC)
  70. # ifdef OPENSSL_UNISTD
  71. # include OPENSSL_UNISTD
  72. # else
  73. # include <unistd.h>
  74. # endif
  75. # include <math.h>
  76. # endif
  77. # endif
  78. # include <openssl/des.h>
  79. # ifdef OPENSSL_SYS_MSDOS /* Visual C++ 2.1 (Windows NT/95) */
  80. # include <stdlib.h>
  81. # include <errno.h>
  82. # include <time.h>
  83. # include <io.h>
  84. # endif
  85. # if defined(__STDC__) || defined(OPENSSL_SYS_VMS) || defined(M_XENIX) || defined(OPENSSL_SYS_MSDOS)
  86. # include <string.h>
  87. # endif
  88. # ifdef OPENSSL_BUILD_SHLIBCRYPTO
  89. # undef OPENSSL_EXTERN
  90. # define OPENSSL_EXTERN OPENSSL_EXPORT
  91. # endif
  92. # define ITERATIONS 16
  93. # define HALF_ITERATIONS 8
  94. /* used in des_read and des_write */
  95. # define MAXWRITE (1024*16)
  96. # define BSIZE (MAXWRITE+4)
  97. # define c2l(c,l) (l =((DES_LONG)(*((c)++))) , \
  98. l|=((DES_LONG)(*((c)++)))<< 8L, \
  99. l|=((DES_LONG)(*((c)++)))<<16L, \
  100. l|=((DES_LONG)(*((c)++)))<<24L)
  101. /* NOTE - c is not incremented as per c2l */
  102. # define c2ln(c,l1,l2,n) { \
  103. c+=n; \
  104. l1=l2=0; \
  105. switch (n) { \
  106. case 8: l2 =((DES_LONG)(*(--(c))))<<24L; \
  107. case 7: l2|=((DES_LONG)(*(--(c))))<<16L; \
  108. case 6: l2|=((DES_LONG)(*(--(c))))<< 8L; \
  109. case 5: l2|=((DES_LONG)(*(--(c)))); \
  110. case 4: l1 =((DES_LONG)(*(--(c))))<<24L; \
  111. case 3: l1|=((DES_LONG)(*(--(c))))<<16L; \
  112. case 2: l1|=((DES_LONG)(*(--(c))))<< 8L; \
  113. case 1: l1|=((DES_LONG)(*(--(c)))); \
  114. } \
  115. }
  116. # define l2c(l,c) (*((c)++)=(unsigned char)(((l) )&0xff), \
  117. *((c)++)=(unsigned char)(((l)>> 8L)&0xff), \
  118. *((c)++)=(unsigned char)(((l)>>16L)&0xff), \
  119. *((c)++)=(unsigned char)(((l)>>24L)&0xff))
  120. /*
  121. * replacements for htonl and ntohl since I have no idea what to do when
  122. * faced with machines with 8 byte longs.
  123. */
  124. # define HDRSIZE 4
  125. # define n2l(c,l) (l =((DES_LONG)(*((c)++)))<<24L, \
  126. l|=((DES_LONG)(*((c)++)))<<16L, \
  127. l|=((DES_LONG)(*((c)++)))<< 8L, \
  128. l|=((DES_LONG)(*((c)++))))
  129. # define l2n(l,c) (*((c)++)=(unsigned char)(((l)>>24L)&0xff), \
  130. *((c)++)=(unsigned char)(((l)>>16L)&0xff), \
  131. *((c)++)=(unsigned char)(((l)>> 8L)&0xff), \
  132. *((c)++)=(unsigned char)(((l) )&0xff))
  133. /* NOTE - c is not incremented as per l2c */
  134. # define l2cn(l1,l2,c,n) { \
  135. c+=n; \
  136. switch (n) { \
  137. case 8: *(--(c))=(unsigned char)(((l2)>>24L)&0xff); \
  138. case 7: *(--(c))=(unsigned char)(((l2)>>16L)&0xff); \
  139. case 6: *(--(c))=(unsigned char)(((l2)>> 8L)&0xff); \
  140. case 5: *(--(c))=(unsigned char)(((l2) )&0xff); \
  141. case 4: *(--(c))=(unsigned char)(((l1)>>24L)&0xff); \
  142. case 3: *(--(c))=(unsigned char)(((l1)>>16L)&0xff); \
  143. case 2: *(--(c))=(unsigned char)(((l1)>> 8L)&0xff); \
  144. case 1: *(--(c))=(unsigned char)(((l1) )&0xff); \
  145. } \
  146. }
  147. # if (defined(OPENSSL_SYS_WIN32) && defined(_MSC_VER)) || defined(__ICC)
  148. # define ROTATE(a,n) (_lrotr(a,n))
  149. # elif defined(__GNUC__) && __GNUC__>=2 && !defined(__STRICT_ANSI__) && !defined(OPENSSL_NO_ASM) && !defined(OPENSSL_NO_INLINE_ASM) && !defined(PEDANTIC)
  150. # if defined(__i386) || defined(__i386__) || defined(__x86_64) || defined(__x86_64__)
  151. # define ROTATE(a,n) ({ register unsigned int ret; \
  152. asm ("rorl %1,%0" \
  153. : "=r"(ret) \
  154. : "I"(n),"0"(a) \
  155. : "cc"); \
  156. ret; \
  157. })
  158. # endif
  159. # endif
  160. # ifndef ROTATE
  161. # define ROTATE(a,n) (((a)>>(n))+((a)<<(32-(n))))
  162. # endif
  163. /*
  164. * Don't worry about the LOAD_DATA() stuff, that is used by fcrypt() to add
  165. * it's little bit to the front
  166. */
  167. # ifdef DES_FCRYPT
  168. # define LOAD_DATA_tmp(R,S,u,t,E0,E1) \
  169. { DES_LONG tmp; LOAD_DATA(R,S,u,t,E0,E1,tmp); }
  170. # define LOAD_DATA(R,S,u,t,E0,E1,tmp) \
  171. t=R^(R>>16L); \
  172. u=t&E0; t&=E1; \
  173. tmp=(u<<16); u^=R^s[S ]; u^=tmp; \
  174. tmp=(t<<16); t^=R^s[S+1]; t^=tmp
  175. # else
  176. # define LOAD_DATA_tmp(a,b,c,d,e,f) LOAD_DATA(a,b,c,d,e,f,g)
  177. # define LOAD_DATA(R,S,u,t,E0,E1,tmp) \
  178. u=R^s[S ]; \
  179. t=R^s[S+1]
  180. # endif
  181. /*
  182. * The changes to this macro may help or hinder, depending on the compiler
  183. * and the architecture. gcc2 always seems to do well :-). Inspired by Dana
  184. * How <[email protected]> DO NOT use the alternative version on machines
  185. * with 8 byte longs. It does not seem to work on the Alpha, even when
  186. * DES_LONG is 4 bytes, probably an issue of accessing non-word aligned
  187. * objects :-(
  188. */
  189. # ifdef DES_PTR
  190. /*
  191. * It recently occurred to me that 0^0^0^0^0^0^0 == 0, so there is no reason
  192. * to not xor all the sub items together. This potentially saves a register
  193. * since things can be xored directly into L
  194. */
  195. # if defined(DES_RISC1) || defined(DES_RISC2)
  196. # ifdef DES_RISC1
  197. # define D_ENCRYPT(LL,R,S) { \
  198. unsigned int u1,u2,u3; \
  199. LOAD_DATA(R,S,u,t,E0,E1,u1); \
  200. u2=(int)u>>8L; \
  201. u1=(int)u&0xfc; \
  202. u2&=0xfc; \
  203. t=ROTATE(t,4); \
  204. u>>=16L; \
  205. LL^= *(const DES_LONG *)(des_SP +u1); \
  206. LL^= *(const DES_LONG *)(des_SP+0x200+u2); \
  207. u3=(int)(u>>8L); \
  208. u1=(int)u&0xfc; \
  209. u3&=0xfc; \
  210. LL^= *(const DES_LONG *)(des_SP+0x400+u1); \
  211. LL^= *(const DES_LONG *)(des_SP+0x600+u3); \
  212. u2=(int)t>>8L; \
  213. u1=(int)t&0xfc; \
  214. u2&=0xfc; \
  215. t>>=16L; \
  216. LL^= *(const DES_LONG *)(des_SP+0x100+u1); \
  217. LL^= *(const DES_LONG *)(des_SP+0x300+u2); \
  218. u3=(int)t>>8L; \
  219. u1=(int)t&0xfc; \
  220. u3&=0xfc; \
  221. LL^= *(const DES_LONG *)(des_SP+0x500+u1); \
  222. LL^= *(const DES_LONG *)(des_SP+0x700+u3); }
  223. # endif
  224. # ifdef DES_RISC2
  225. # define D_ENCRYPT(LL,R,S) { \
  226. unsigned int u1,u2,s1,s2; \
  227. LOAD_DATA(R,S,u,t,E0,E1,u1); \
  228. u2=(int)u>>8L; \
  229. u1=(int)u&0xfc; \
  230. u2&=0xfc; \
  231. t=ROTATE(t,4); \
  232. LL^= *(const DES_LONG *)(des_SP +u1); \
  233. LL^= *(const DES_LONG *)(des_SP+0x200+u2); \
  234. s1=(int)(u>>16L); \
  235. s2=(int)(u>>24L); \
  236. s1&=0xfc; \
  237. s2&=0xfc; \
  238. LL^= *(const DES_LONG *)(des_SP+0x400+s1); \
  239. LL^= *(const DES_LONG *)(des_SP+0x600+s2); \
  240. u2=(int)t>>8L; \
  241. u1=(int)t&0xfc; \
  242. u2&=0xfc; \
  243. LL^= *(const DES_LONG *)(des_SP+0x100+u1); \
  244. LL^= *(const DES_LONG *)(des_SP+0x300+u2); \
  245. s1=(int)(t>>16L); \
  246. s2=(int)(t>>24L); \
  247. s1&=0xfc; \
  248. s2&=0xfc; \
  249. LL^= *(const DES_LONG *)(des_SP+0x500+s1); \
  250. LL^= *(const DES_LONG *)(des_SP+0x700+s2); }
  251. # endif
  252. # else
  253. # define D_ENCRYPT(LL,R,S) { \
  254. LOAD_DATA_tmp(R,S,u,t,E0,E1); \
  255. t=ROTATE(t,4); \
  256. LL^= \
  257. *(const DES_LONG *)(des_SP +((u )&0xfc))^ \
  258. *(const DES_LONG *)(des_SP+0x200+((u>> 8L)&0xfc))^ \
  259. *(const DES_LONG *)(des_SP+0x400+((u>>16L)&0xfc))^ \
  260. *(const DES_LONG *)(des_SP+0x600+((u>>24L)&0xfc))^ \
  261. *(const DES_LONG *)(des_SP+0x100+((t )&0xfc))^ \
  262. *(const DES_LONG *)(des_SP+0x300+((t>> 8L)&0xfc))^ \
  263. *(const DES_LONG *)(des_SP+0x500+((t>>16L)&0xfc))^ \
  264. *(const DES_LONG *)(des_SP+0x700+((t>>24L)&0xfc)); }
  265. # endif
  266. # else /* original version */
  267. # if defined(DES_RISC1) || defined(DES_RISC2)
  268. # ifdef DES_RISC1
  269. # define D_ENCRYPT(LL,R,S) {\
  270. unsigned int u1,u2,u3; \
  271. LOAD_DATA(R,S,u,t,E0,E1,u1); \
  272. u>>=2L; \
  273. t=ROTATE(t,6); \
  274. u2=(int)u>>8L; \
  275. u1=(int)u&0x3f; \
  276. u2&=0x3f; \
  277. u>>=16L; \
  278. LL^=DES_SPtrans[0][u1]; \
  279. LL^=DES_SPtrans[2][u2]; \
  280. u3=(int)u>>8L; \
  281. u1=(int)u&0x3f; \
  282. u3&=0x3f; \
  283. LL^=DES_SPtrans[4][u1]; \
  284. LL^=DES_SPtrans[6][u3]; \
  285. u2=(int)t>>8L; \
  286. u1=(int)t&0x3f; \
  287. u2&=0x3f; \
  288. t>>=16L; \
  289. LL^=DES_SPtrans[1][u1]; \
  290. LL^=DES_SPtrans[3][u2]; \
  291. u3=(int)t>>8L; \
  292. u1=(int)t&0x3f; \
  293. u3&=0x3f; \
  294. LL^=DES_SPtrans[5][u1]; \
  295. LL^=DES_SPtrans[7][u3]; }
  296. # endif
  297. # ifdef DES_RISC2
  298. # define D_ENCRYPT(LL,R,S) {\
  299. unsigned int u1,u2,s1,s2; \
  300. LOAD_DATA(R,S,u,t,E0,E1,u1); \
  301. u>>=2L; \
  302. t=ROTATE(t,6); \
  303. u2=(int)u>>8L; \
  304. u1=(int)u&0x3f; \
  305. u2&=0x3f; \
  306. LL^=DES_SPtrans[0][u1]; \
  307. LL^=DES_SPtrans[2][u2]; \
  308. s1=(int)u>>16L; \
  309. s2=(int)u>>24L; \
  310. s1&=0x3f; \
  311. s2&=0x3f; \
  312. LL^=DES_SPtrans[4][s1]; \
  313. LL^=DES_SPtrans[6][s2]; \
  314. u2=(int)t>>8L; \
  315. u1=(int)t&0x3f; \
  316. u2&=0x3f; \
  317. LL^=DES_SPtrans[1][u1]; \
  318. LL^=DES_SPtrans[3][u2]; \
  319. s1=(int)t>>16; \
  320. s2=(int)t>>24L; \
  321. s1&=0x3f; \
  322. s2&=0x3f; \
  323. LL^=DES_SPtrans[5][s1]; \
  324. LL^=DES_SPtrans[7][s2]; }
  325. # endif
  326. # else
  327. # define D_ENCRYPT(LL,R,S) {\
  328. LOAD_DATA_tmp(R,S,u,t,E0,E1); \
  329. t=ROTATE(t,4); \
  330. LL^=\
  331. DES_SPtrans[0][(u>> 2L)&0x3f]^ \
  332. DES_SPtrans[2][(u>>10L)&0x3f]^ \
  333. DES_SPtrans[4][(u>>18L)&0x3f]^ \
  334. DES_SPtrans[6][(u>>26L)&0x3f]^ \
  335. DES_SPtrans[1][(t>> 2L)&0x3f]^ \
  336. DES_SPtrans[3][(t>>10L)&0x3f]^ \
  337. DES_SPtrans[5][(t>>18L)&0x3f]^ \
  338. DES_SPtrans[7][(t>>26L)&0x3f]; }
  339. # endif
  340. # endif
  341. /*-
  342. * IP and FP
  343. * The problem is more of a geometric problem that random bit fiddling.
  344. 0 1 2 3 4 5 6 7 62 54 46 38 30 22 14 6
  345. 8 9 10 11 12 13 14 15 60 52 44 36 28 20 12 4
  346. 16 17 18 19 20 21 22 23 58 50 42 34 26 18 10 2
  347. 24 25 26 27 28 29 30 31 to 56 48 40 32 24 16 8 0
  348. 32 33 34 35 36 37 38 39 63 55 47 39 31 23 15 7
  349. 40 41 42 43 44 45 46 47 61 53 45 37 29 21 13 5
  350. 48 49 50 51 52 53 54 55 59 51 43 35 27 19 11 3
  351. 56 57 58 59 60 61 62 63 57 49 41 33 25 17 9 1
  352. The output has been subject to swaps of the form
  353. 0 1 -> 3 1 but the odd and even bits have been put into
  354. 2 3 2 0
  355. different words. The main trick is to remember that
  356. t=((l>>size)^r)&(mask);
  357. r^=t;
  358. l^=(t<<size);
  359. can be used to swap and move bits between words.
  360. So l = 0 1 2 3 r = 16 17 18 19
  361. 4 5 6 7 20 21 22 23
  362. 8 9 10 11 24 25 26 27
  363. 12 13 14 15 28 29 30 31
  364. becomes (for size == 2 and mask == 0x3333)
  365. t = 2^16 3^17 -- -- l = 0 1 16 17 r = 2 3 18 19
  366. 6^20 7^21 -- -- 4 5 20 21 6 7 22 23
  367. 10^24 11^25 -- -- 8 9 24 25 10 11 24 25
  368. 14^28 15^29 -- -- 12 13 28 29 14 15 28 29
  369. Thanks for hints from Richard Outerbridge - he told me IP&FP
  370. could be done in 15 xor, 10 shifts and 5 ands.
  371. When I finally started to think of the problem in 2D
  372. I first got ~42 operations without xors. When I remembered
  373. how to use xors :-) I got it to its final state.
  374. */
  375. # define PERM_OP(a,b,t,n,m) ((t)=((((a)>>(n))^(b))&(m)),\
  376. (b)^=(t),\
  377. (a)^=((t)<<(n)))
  378. # define IP(l,r) \
  379. { \
  380. register DES_LONG tt; \
  381. PERM_OP(r,l,tt, 4,0x0f0f0f0fL); \
  382. PERM_OP(l,r,tt,16,0x0000ffffL); \
  383. PERM_OP(r,l,tt, 2,0x33333333L); \
  384. PERM_OP(l,r,tt, 8,0x00ff00ffL); \
  385. PERM_OP(r,l,tt, 1,0x55555555L); \
  386. }
  387. # define FP(l,r) \
  388. { \
  389. register DES_LONG tt; \
  390. PERM_OP(l,r,tt, 1,0x55555555L); \
  391. PERM_OP(r,l,tt, 8,0x00ff00ffL); \
  392. PERM_OP(l,r,tt, 2,0x33333333L); \
  393. PERM_OP(r,l,tt,16,0x0000ffffL); \
  394. PERM_OP(l,r,tt, 4,0x0f0f0f0fL); \
  395. }
  396. extern const DES_LONG DES_SPtrans[8][64];
  397. void fcrypt_body(DES_LONG *out, DES_key_schedule *ks,
  398. DES_LONG Eswap0, DES_LONG Eswap1);
  399. # ifdef OPENSSL_SMALL_FOOTPRINT
  400. # undef DES_UNROLL
  401. # endif
  402. #endif