des_locl.h 16 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444
  1. /* crypto/des/des_locl.h */
  2. /* Copyright (C) 1995-1997 Eric Young (eay@cryptsoft.com)
  3. * All rights reserved.
  4. *
  5. * This package is an SSL implementation written
  6. * by Eric Young (eay@cryptsoft.com).
  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 (tjh@cryptsoft.com).
  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 (eay@cryptsoft.com)"
  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 (tjh@cryptsoft.com)"
  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))
  148. # define ROTATE(a,n) (_lrotr(a,n))
  149. # elif defined(__ICC)
  150. # define ROTATE(a,n) (_rotr(a,n))
  151. # elif defined(__GNUC__) && __GNUC__>=2 && !defined(__STRICT_ANSI__) && !defined(OPENSSL_NO_ASM) && !defined(OPENSSL_NO_INLINE_ASM) && !defined(PEDANTIC)
  152. # if defined(__i386) || defined(__i386__) || defined(__x86_64) || defined(__x86_64__)
  153. # define ROTATE(a,n) ({ register unsigned int ret; \
  154. asm ("rorl %1,%0" \
  155. : "=r"(ret) \
  156. : "I"(n),"0"(a) \
  157. : "cc"); \
  158. ret; \
  159. })
  160. # endif
  161. # endif
  162. # ifndef ROTATE
  163. # define ROTATE(a,n) (((a)>>(n))+((a)<<(32-(n))))
  164. # endif
  165. /*
  166. * Don't worry about the LOAD_DATA() stuff, that is used by fcrypt() to add
  167. * it's little bit to the front
  168. */
  169. # ifdef DES_FCRYPT
  170. # define LOAD_DATA_tmp(R,S,u,t,E0,E1) \
  171. { DES_LONG tmp; LOAD_DATA(R,S,u,t,E0,E1,tmp); }
  172. # define LOAD_DATA(R,S,u,t,E0,E1,tmp) \
  173. t=R^(R>>16L); \
  174. u=t&E0; t&=E1; \
  175. tmp=(u<<16); u^=R^s[S ]; u^=tmp; \
  176. tmp=(t<<16); t^=R^s[S+1]; t^=tmp
  177. # else
  178. # define LOAD_DATA_tmp(a,b,c,d,e,f) LOAD_DATA(a,b,c,d,e,f,g)
  179. # define LOAD_DATA(R,S,u,t,E0,E1,tmp) \
  180. u=R^s[S ]; \
  181. t=R^s[S+1]
  182. # endif
  183. /*
  184. * The changes to this macro may help or hinder, depending on the compiler
  185. * and the architecture. gcc2 always seems to do well :-). Inspired by Dana
  186. * How <how@isl.stanford.edu> DO NOT use the alternative version on machines
  187. * with 8 byte longs. It does not seem to work on the Alpha, even when
  188. * DES_LONG is 4 bytes, probably an issue of accessing non-word aligned
  189. * objects :-(
  190. */
  191. # ifdef DES_PTR
  192. /*
  193. * It recently occurred to me that 0^0^0^0^0^0^0 == 0, so there is no reason
  194. * to not xor all the sub items together. This potentially saves a register
  195. * since things can be xored directly into L
  196. */
  197. # if defined(DES_RISC1) || defined(DES_RISC2)
  198. # ifdef DES_RISC1
  199. # define D_ENCRYPT(LL,R,S) { \
  200. unsigned int u1,u2,u3; \
  201. LOAD_DATA(R,S,u,t,E0,E1,u1); \
  202. u2=(int)u>>8L; \
  203. u1=(int)u&0xfc; \
  204. u2&=0xfc; \
  205. t=ROTATE(t,4); \
  206. u>>=16L; \
  207. LL^= *(const DES_LONG *)(des_SP +u1); \
  208. LL^= *(const DES_LONG *)(des_SP+0x200+u2); \
  209. u3=(int)(u>>8L); \
  210. u1=(int)u&0xfc; \
  211. u3&=0xfc; \
  212. LL^= *(const DES_LONG *)(des_SP+0x400+u1); \
  213. LL^= *(const DES_LONG *)(des_SP+0x600+u3); \
  214. u2=(int)t>>8L; \
  215. u1=(int)t&0xfc; \
  216. u2&=0xfc; \
  217. t>>=16L; \
  218. LL^= *(const DES_LONG *)(des_SP+0x100+u1); \
  219. LL^= *(const DES_LONG *)(des_SP+0x300+u2); \
  220. u3=(int)t>>8L; \
  221. u1=(int)t&0xfc; \
  222. u3&=0xfc; \
  223. LL^= *(const DES_LONG *)(des_SP+0x500+u1); \
  224. LL^= *(const DES_LONG *)(des_SP+0x700+u3); }
  225. # endif
  226. # ifdef DES_RISC2
  227. # define D_ENCRYPT(LL,R,S) { \
  228. unsigned int u1,u2,s1,s2; \
  229. LOAD_DATA(R,S,u,t,E0,E1,u1); \
  230. u2=(int)u>>8L; \
  231. u1=(int)u&0xfc; \
  232. u2&=0xfc; \
  233. t=ROTATE(t,4); \
  234. LL^= *(const DES_LONG *)(des_SP +u1); \
  235. LL^= *(const DES_LONG *)(des_SP+0x200+u2); \
  236. s1=(int)(u>>16L); \
  237. s2=(int)(u>>24L); \
  238. s1&=0xfc; \
  239. s2&=0xfc; \
  240. LL^= *(const DES_LONG *)(des_SP+0x400+s1); \
  241. LL^= *(const DES_LONG *)(des_SP+0x600+s2); \
  242. u2=(int)t>>8L; \
  243. u1=(int)t&0xfc; \
  244. u2&=0xfc; \
  245. LL^= *(const DES_LONG *)(des_SP+0x100+u1); \
  246. LL^= *(const DES_LONG *)(des_SP+0x300+u2); \
  247. s1=(int)(t>>16L); \
  248. s2=(int)(t>>24L); \
  249. s1&=0xfc; \
  250. s2&=0xfc; \
  251. LL^= *(const DES_LONG *)(des_SP+0x500+s1); \
  252. LL^= *(const DES_LONG *)(des_SP+0x700+s2); }
  253. # endif
  254. # else
  255. # define D_ENCRYPT(LL,R,S) { \
  256. LOAD_DATA_tmp(R,S,u,t,E0,E1); \
  257. t=ROTATE(t,4); \
  258. LL^= \
  259. *(const DES_LONG *)(des_SP +((u )&0xfc))^ \
  260. *(const DES_LONG *)(des_SP+0x200+((u>> 8L)&0xfc))^ \
  261. *(const DES_LONG *)(des_SP+0x400+((u>>16L)&0xfc))^ \
  262. *(const DES_LONG *)(des_SP+0x600+((u>>24L)&0xfc))^ \
  263. *(const DES_LONG *)(des_SP+0x100+((t )&0xfc))^ \
  264. *(const DES_LONG *)(des_SP+0x300+((t>> 8L)&0xfc))^ \
  265. *(const DES_LONG *)(des_SP+0x500+((t>>16L)&0xfc))^ \
  266. *(const DES_LONG *)(des_SP+0x700+((t>>24L)&0xfc)); }
  267. # endif
  268. # else /* original version */
  269. # if defined(DES_RISC1) || defined(DES_RISC2)
  270. # ifdef DES_RISC1
  271. # define D_ENCRYPT(LL,R,S) {\
  272. unsigned int u1,u2,u3; \
  273. LOAD_DATA(R,S,u,t,E0,E1,u1); \
  274. u>>=2L; \
  275. t=ROTATE(t,6); \
  276. u2=(int)u>>8L; \
  277. u1=(int)u&0x3f; \
  278. u2&=0x3f; \
  279. u>>=16L; \
  280. LL^=DES_SPtrans[0][u1]; \
  281. LL^=DES_SPtrans[2][u2]; \
  282. u3=(int)u>>8L; \
  283. u1=(int)u&0x3f; \
  284. u3&=0x3f; \
  285. LL^=DES_SPtrans[4][u1]; \
  286. LL^=DES_SPtrans[6][u3]; \
  287. u2=(int)t>>8L; \
  288. u1=(int)t&0x3f; \
  289. u2&=0x3f; \
  290. t>>=16L; \
  291. LL^=DES_SPtrans[1][u1]; \
  292. LL^=DES_SPtrans[3][u2]; \
  293. u3=(int)t>>8L; \
  294. u1=(int)t&0x3f; \
  295. u3&=0x3f; \
  296. LL^=DES_SPtrans[5][u1]; \
  297. LL^=DES_SPtrans[7][u3]; }
  298. # endif
  299. # ifdef DES_RISC2
  300. # define D_ENCRYPT(LL,R,S) {\
  301. unsigned int u1,u2,s1,s2; \
  302. LOAD_DATA(R,S,u,t,E0,E1,u1); \
  303. u>>=2L; \
  304. t=ROTATE(t,6); \
  305. u2=(int)u>>8L; \
  306. u1=(int)u&0x3f; \
  307. u2&=0x3f; \
  308. LL^=DES_SPtrans[0][u1]; \
  309. LL^=DES_SPtrans[2][u2]; \
  310. s1=(int)u>>16L; \
  311. s2=(int)u>>24L; \
  312. s1&=0x3f; \
  313. s2&=0x3f; \
  314. LL^=DES_SPtrans[4][s1]; \
  315. LL^=DES_SPtrans[6][s2]; \
  316. u2=(int)t>>8L; \
  317. u1=(int)t&0x3f; \
  318. u2&=0x3f; \
  319. LL^=DES_SPtrans[1][u1]; \
  320. LL^=DES_SPtrans[3][u2]; \
  321. s1=(int)t>>16; \
  322. s2=(int)t>>24L; \
  323. s1&=0x3f; \
  324. s2&=0x3f; \
  325. LL^=DES_SPtrans[5][s1]; \
  326. LL^=DES_SPtrans[7][s2]; }
  327. # endif
  328. # else
  329. # define D_ENCRYPT(LL,R,S) {\
  330. LOAD_DATA_tmp(R,S,u,t,E0,E1); \
  331. t=ROTATE(t,4); \
  332. LL^=\
  333. DES_SPtrans[0][(u>> 2L)&0x3f]^ \
  334. DES_SPtrans[2][(u>>10L)&0x3f]^ \
  335. DES_SPtrans[4][(u>>18L)&0x3f]^ \
  336. DES_SPtrans[6][(u>>26L)&0x3f]^ \
  337. DES_SPtrans[1][(t>> 2L)&0x3f]^ \
  338. DES_SPtrans[3][(t>>10L)&0x3f]^ \
  339. DES_SPtrans[5][(t>>18L)&0x3f]^ \
  340. DES_SPtrans[7][(t>>26L)&0x3f]; }
  341. # endif
  342. # endif
  343. /*-
  344. * IP and FP
  345. * The problem is more of a geometric problem that random bit fiddling.
  346. 0 1 2 3 4 5 6 7 62 54 46 38 30 22 14 6
  347. 8 9 10 11 12 13 14 15 60 52 44 36 28 20 12 4
  348. 16 17 18 19 20 21 22 23 58 50 42 34 26 18 10 2
  349. 24 25 26 27 28 29 30 31 to 56 48 40 32 24 16 8 0
  350. 32 33 34 35 36 37 38 39 63 55 47 39 31 23 15 7
  351. 40 41 42 43 44 45 46 47 61 53 45 37 29 21 13 5
  352. 48 49 50 51 52 53 54 55 59 51 43 35 27 19 11 3
  353. 56 57 58 59 60 61 62 63 57 49 41 33 25 17 9 1
  354. The output has been subject to swaps of the form
  355. 0 1 -> 3 1 but the odd and even bits have been put into
  356. 2 3 2 0
  357. different words. The main trick is to remember that
  358. t=((l>>size)^r)&(mask);
  359. r^=t;
  360. l^=(t<<size);
  361. can be used to swap and move bits between words.
  362. So l = 0 1 2 3 r = 16 17 18 19
  363. 4 5 6 7 20 21 22 23
  364. 8 9 10 11 24 25 26 27
  365. 12 13 14 15 28 29 30 31
  366. becomes (for size == 2 and mask == 0x3333)
  367. t = 2^16 3^17 -- -- l = 0 1 16 17 r = 2 3 18 19
  368. 6^20 7^21 -- -- 4 5 20 21 6 7 22 23
  369. 10^24 11^25 -- -- 8 9 24 25 10 11 24 25
  370. 14^28 15^29 -- -- 12 13 28 29 14 15 28 29
  371. Thanks for hints from Richard Outerbridge - he told me IP&FP
  372. could be done in 15 xor, 10 shifts and 5 ands.
  373. When I finally started to think of the problem in 2D
  374. I first got ~42 operations without xors. When I remembered
  375. how to use xors :-) I got it to its final state.
  376. */
  377. # define PERM_OP(a,b,t,n,m) ((t)=((((a)>>(n))^(b))&(m)),\
  378. (b)^=(t),\
  379. (a)^=((t)<<(n)))
  380. # define IP(l,r) \
  381. { \
  382. register DES_LONG tt; \
  383. PERM_OP(r,l,tt, 4,0x0f0f0f0fL); \
  384. PERM_OP(l,r,tt,16,0x0000ffffL); \
  385. PERM_OP(r,l,tt, 2,0x33333333L); \
  386. PERM_OP(l,r,tt, 8,0x00ff00ffL); \
  387. PERM_OP(r,l,tt, 1,0x55555555L); \
  388. }
  389. # define FP(l,r) \
  390. { \
  391. register DES_LONG tt; \
  392. PERM_OP(l,r,tt, 1,0x55555555L); \
  393. PERM_OP(r,l,tt, 8,0x00ff00ffL); \
  394. PERM_OP(l,r,tt, 2,0x33333333L); \
  395. PERM_OP(r,l,tt,16,0x0000ffffL); \
  396. PERM_OP(l,r,tt, 4,0x0f0f0f0fL); \
  397. }
  398. extern const DES_LONG DES_SPtrans[8][64];
  399. void fcrypt_body(DES_LONG *out, DES_key_schedule *ks,
  400. DES_LONG Eswap0, DES_LONG Eswap1);
  401. # ifdef OPENSSL_SMALL_FOOTPRINT
  402. # undef DES_UNROLL
  403. # endif
  404. #endif