xref: /openbsd/sys/crypto/blf.c (revision 09467b48)
1 /*	$OpenBSD: blf.c,v 1.7 2007/11/26 09:28:34 martynas Exp $	*/
2 
3 /*
4  * Blowfish block cipher for OpenBSD
5  * Copyright 1997 Niels Provos <provos@physnet.uni-hamburg.de>
6  * All rights reserved.
7  *
8  * Implementation advice by David Mazieres <dm@lcs.mit.edu>.
9  *
10  * Redistribution and use in source and binary forms, with or without
11  * modification, are permitted provided that the following conditions
12  * are met:
13  * 1. Redistributions of source code must retain the above copyright
14  *    notice, this list of conditions and the following disclaimer.
15  * 2. Redistributions in binary form must reproduce the above copyright
16  *    notice, this list of conditions and the following disclaimer in the
17  *    documentation and/or other materials provided with the distribution.
18  * 3. All advertising materials mentioning features or use of this software
19  *    must display the following acknowledgement:
20  *      This product includes software developed by Niels Provos.
21  * 4. The name of the author may not be used to endorse or promote products
22  *    derived from this software without specific prior written permission.
23  *
24  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
25  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
26  * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
27  * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
28  * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
29  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
30  * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
31  * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
32  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
33  * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
34  */
35 
36 /*
37  * This code is derived from section 14.3 and the given source
38  * in section V of Applied Cryptography, second edition.
39  * Blowfish is an unpatented fast block cipher designed by
40  * Bruce Schneier.
41  */
42 
43 #include <sys/param.h>
44 #include <sys/systm.h>
45 
46 #include <crypto/blf.h>
47 
48 #undef inline
49 #ifdef __GNUC__
50 #define inline __inline
51 #else				/* !__GNUC__ */
52 #define inline
53 #endif				/* !__GNUC__ */
54 
55 /* Function for Feistel Networks */
56 
57 #define F(s, x) ((((s)[        (((x)>>24)&0xFF)]  \
58 		 + (s)[0x100 + (((x)>>16)&0xFF)]) \
59 		 ^ (s)[0x200 + (((x)>> 8)&0xFF)]) \
60 		 + (s)[0x300 + ( (x)     &0xFF)])
61 
62 #define BLFRND(s,p,i,j,n) (i ^= F(s,j) ^ (p)[n])
63 
64 void
65 Blowfish_encipher(blf_ctx *c, u_int32_t *x)
66 {
67 	u_int32_t Xl;
68 	u_int32_t Xr;
69 	u_int32_t *s = c->S[0];
70 	u_int32_t *p = c->P;
71 
72 	Xl = x[0];
73 	Xr = x[1];
74 
75 	Xl ^= p[0];
76 	BLFRND(s, p, Xr, Xl, 1); BLFRND(s, p, Xl, Xr, 2);
77 	BLFRND(s, p, Xr, Xl, 3); BLFRND(s, p, Xl, Xr, 4);
78 	BLFRND(s, p, Xr, Xl, 5); BLFRND(s, p, Xl, Xr, 6);
79 	BLFRND(s, p, Xr, Xl, 7); BLFRND(s, p, Xl, Xr, 8);
80 	BLFRND(s, p, Xr, Xl, 9); BLFRND(s, p, Xl, Xr, 10);
81 	BLFRND(s, p, Xr, Xl, 11); BLFRND(s, p, Xl, Xr, 12);
82 	BLFRND(s, p, Xr, Xl, 13); BLFRND(s, p, Xl, Xr, 14);
83 	BLFRND(s, p, Xr, Xl, 15); BLFRND(s, p, Xl, Xr, 16);
84 
85 	x[0] = Xr ^ p[17];
86 	x[1] = Xl;
87 }
88 
89 void
90 Blowfish_decipher(blf_ctx *c, u_int32_t *x)
91 {
92 	u_int32_t Xl;
93 	u_int32_t Xr;
94 	u_int32_t *s = c->S[0];
95 	u_int32_t *p = c->P;
96 
97 	Xl = x[0];
98 	Xr = x[1];
99 
100 	Xl ^= p[17];
101 	BLFRND(s, p, Xr, Xl, 16); BLFRND(s, p, Xl, Xr, 15);
102 	BLFRND(s, p, Xr, Xl, 14); BLFRND(s, p, Xl, Xr, 13);
103 	BLFRND(s, p, Xr, Xl, 12); BLFRND(s, p, Xl, Xr, 11);
104 	BLFRND(s, p, Xr, Xl, 10); BLFRND(s, p, Xl, Xr, 9);
105 	BLFRND(s, p, Xr, Xl, 8); BLFRND(s, p, Xl, Xr, 7);
106 	BLFRND(s, p, Xr, Xl, 6); BLFRND(s, p, Xl, Xr, 5);
107 	BLFRND(s, p, Xr, Xl, 4); BLFRND(s, p, Xl, Xr, 3);
108 	BLFRND(s, p, Xr, Xl, 2); BLFRND(s, p, Xl, Xr, 1);
109 
110 	x[0] = Xr ^ p[0];
111 	x[1] = Xl;
112 }
113 
114 void
115 Blowfish_initstate(blf_ctx *c)
116 {
117 	/* P-box and S-box tables initialized with digits of Pi */
118 
119 	static const blf_ctx initstate =
120 
121 	{ {
122 		{
123 			0xd1310ba6, 0x98dfb5ac, 0x2ffd72db, 0xd01adfb7,
124 			0xb8e1afed, 0x6a267e96, 0xba7c9045, 0xf12c7f99,
125 			0x24a19947, 0xb3916cf7, 0x0801f2e2, 0x858efc16,
126 			0x636920d8, 0x71574e69, 0xa458fea3, 0xf4933d7e,
127 			0x0d95748f, 0x728eb658, 0x718bcd58, 0x82154aee,
128 			0x7b54a41d, 0xc25a59b5, 0x9c30d539, 0x2af26013,
129 			0xc5d1b023, 0x286085f0, 0xca417918, 0xb8db38ef,
130 			0x8e79dcb0, 0x603a180e, 0x6c9e0e8b, 0xb01e8a3e,
131 			0xd71577c1, 0xbd314b27, 0x78af2fda, 0x55605c60,
132 			0xe65525f3, 0xaa55ab94, 0x57489862, 0x63e81440,
133 			0x55ca396a, 0x2aab10b6, 0xb4cc5c34, 0x1141e8ce,
134 			0xa15486af, 0x7c72e993, 0xb3ee1411, 0x636fbc2a,
135 			0x2ba9c55d, 0x741831f6, 0xce5c3e16, 0x9b87931e,
136 			0xafd6ba33, 0x6c24cf5c, 0x7a325381, 0x28958677,
137 			0x3b8f4898, 0x6b4bb9af, 0xc4bfe81b, 0x66282193,
138 			0x61d809cc, 0xfb21a991, 0x487cac60, 0x5dec8032,
139 			0xef845d5d, 0xe98575b1, 0xdc262302, 0xeb651b88,
140 			0x23893e81, 0xd396acc5, 0x0f6d6ff3, 0x83f44239,
141 			0x2e0b4482, 0xa4842004, 0x69c8f04a, 0x9e1f9b5e,
142 			0x21c66842, 0xf6e96c9a, 0x670c9c61, 0xabd388f0,
143 			0x6a51a0d2, 0xd8542f68, 0x960fa728, 0xab5133a3,
144 			0x6eef0b6c, 0x137a3be4, 0xba3bf050, 0x7efb2a98,
145 			0xa1f1651d, 0x39af0176, 0x66ca593e, 0x82430e88,
146 			0x8cee8619, 0x456f9fb4, 0x7d84a5c3, 0x3b8b5ebe,
147 			0xe06f75d8, 0x85c12073, 0x401a449f, 0x56c16aa6,
148 			0x4ed3aa62, 0x363f7706, 0x1bfedf72, 0x429b023d,
149 			0x37d0d724, 0xd00a1248, 0xdb0fead3, 0x49f1c09b,
150 			0x075372c9, 0x80991b7b, 0x25d479d8, 0xf6e8def7,
151 			0xe3fe501a, 0xb6794c3b, 0x976ce0bd, 0x04c006ba,
152 			0xc1a94fb6, 0x409f60c4, 0x5e5c9ec2, 0x196a2463,
153 			0x68fb6faf, 0x3e6c53b5, 0x1339b2eb, 0x3b52ec6f,
154 			0x6dfc511f, 0x9b30952c, 0xcc814544, 0xaf5ebd09,
155 			0xbee3d004, 0xde334afd, 0x660f2807, 0x192e4bb3,
156 			0xc0cba857, 0x45c8740f, 0xd20b5f39, 0xb9d3fbdb,
157 			0x5579c0bd, 0x1a60320a, 0xd6a100c6, 0x402c7279,
158 			0x679f25fe, 0xfb1fa3cc, 0x8ea5e9f8, 0xdb3222f8,
159 			0x3c7516df, 0xfd616b15, 0x2f501ec8, 0xad0552ab,
160 			0x323db5fa, 0xfd238760, 0x53317b48, 0x3e00df82,
161 			0x9e5c57bb, 0xca6f8ca0, 0x1a87562e, 0xdf1769db,
162 			0xd542a8f6, 0x287effc3, 0xac6732c6, 0x8c4f5573,
163 			0x695b27b0, 0xbbca58c8, 0xe1ffa35d, 0xb8f011a0,
164 			0x10fa3d98, 0xfd2183b8, 0x4afcb56c, 0x2dd1d35b,
165 			0x9a53e479, 0xb6f84565, 0xd28e49bc, 0x4bfb9790,
166 			0xe1ddf2da, 0xa4cb7e33, 0x62fb1341, 0xcee4c6e8,
167 			0xef20cada, 0x36774c01, 0xd07e9efe, 0x2bf11fb4,
168 			0x95dbda4d, 0xae909198, 0xeaad8e71, 0x6b93d5a0,
169 			0xd08ed1d0, 0xafc725e0, 0x8e3c5b2f, 0x8e7594b7,
170 			0x8ff6e2fb, 0xf2122b64, 0x8888b812, 0x900df01c,
171 			0x4fad5ea0, 0x688fc31c, 0xd1cff191, 0xb3a8c1ad,
172 			0x2f2f2218, 0xbe0e1777, 0xea752dfe, 0x8b021fa1,
173 			0xe5a0cc0f, 0xb56f74e8, 0x18acf3d6, 0xce89e299,
174 			0xb4a84fe0, 0xfd13e0b7, 0x7cc43b81, 0xd2ada8d9,
175 			0x165fa266, 0x80957705, 0x93cc7314, 0x211a1477,
176 			0xe6ad2065, 0x77b5fa86, 0xc75442f5, 0xfb9d35cf,
177 			0xebcdaf0c, 0x7b3e89a0, 0xd6411bd3, 0xae1e7e49,
178 			0x00250e2d, 0x2071b35e, 0x226800bb, 0x57b8e0af,
179 			0x2464369b, 0xf009b91e, 0x5563911d, 0x59dfa6aa,
180 			0x78c14389, 0xd95a537f, 0x207d5ba2, 0x02e5b9c5,
181 			0x83260376, 0x6295cfa9, 0x11c81968, 0x4e734a41,
182 			0xb3472dca, 0x7b14a94a, 0x1b510052, 0x9a532915,
183 			0xd60f573f, 0xbc9bc6e4, 0x2b60a476, 0x81e67400,
184 			0x08ba6fb5, 0x571be91f, 0xf296ec6b, 0x2a0dd915,
185 			0xb6636521, 0xe7b9f9b6, 0xff34052e, 0xc5855664,
186 		0x53b02d5d, 0xa99f8fa1, 0x08ba4799, 0x6e85076a},
187 		{
188 			0x4b7a70e9, 0xb5b32944, 0xdb75092e, 0xc4192623,
189 			0xad6ea6b0, 0x49a7df7d, 0x9cee60b8, 0x8fedb266,
190 			0xecaa8c71, 0x699a17ff, 0x5664526c, 0xc2b19ee1,
191 			0x193602a5, 0x75094c29, 0xa0591340, 0xe4183a3e,
192 			0x3f54989a, 0x5b429d65, 0x6b8fe4d6, 0x99f73fd6,
193 			0xa1d29c07, 0xefe830f5, 0x4d2d38e6, 0xf0255dc1,
194 			0x4cdd2086, 0x8470eb26, 0x6382e9c6, 0x021ecc5e,
195 			0x09686b3f, 0x3ebaefc9, 0x3c971814, 0x6b6a70a1,
196 			0x687f3584, 0x52a0e286, 0xb79c5305, 0xaa500737,
197 			0x3e07841c, 0x7fdeae5c, 0x8e7d44ec, 0x5716f2b8,
198 			0xb03ada37, 0xf0500c0d, 0xf01c1f04, 0x0200b3ff,
199 			0xae0cf51a, 0x3cb574b2, 0x25837a58, 0xdc0921bd,
200 			0xd19113f9, 0x7ca92ff6, 0x94324773, 0x22f54701,
201 			0x3ae5e581, 0x37c2dadc, 0xc8b57634, 0x9af3dda7,
202 			0xa9446146, 0x0fd0030e, 0xecc8c73e, 0xa4751e41,
203 			0xe238cd99, 0x3bea0e2f, 0x3280bba1, 0x183eb331,
204 			0x4e548b38, 0x4f6db908, 0x6f420d03, 0xf60a04bf,
205 			0x2cb81290, 0x24977c79, 0x5679b072, 0xbcaf89af,
206 			0xde9a771f, 0xd9930810, 0xb38bae12, 0xdccf3f2e,
207 			0x5512721f, 0x2e6b7124, 0x501adde6, 0x9f84cd87,
208 			0x7a584718, 0x7408da17, 0xbc9f9abc, 0xe94b7d8c,
209 			0xec7aec3a, 0xdb851dfa, 0x63094366, 0xc464c3d2,
210 			0xef1c1847, 0x3215d908, 0xdd433b37, 0x24c2ba16,
211 			0x12a14d43, 0x2a65c451, 0x50940002, 0x133ae4dd,
212 			0x71dff89e, 0x10314e55, 0x81ac77d6, 0x5f11199b,
213 			0x043556f1, 0xd7a3c76b, 0x3c11183b, 0x5924a509,
214 			0xf28fe6ed, 0x97f1fbfa, 0x9ebabf2c, 0x1e153c6e,
215 			0x86e34570, 0xeae96fb1, 0x860e5e0a, 0x5a3e2ab3,
216 			0x771fe71c, 0x4e3d06fa, 0x2965dcb9, 0x99e71d0f,
217 			0x803e89d6, 0x5266c825, 0x2e4cc978, 0x9c10b36a,
218 			0xc6150eba, 0x94e2ea78, 0xa5fc3c53, 0x1e0a2df4,
219 			0xf2f74ea7, 0x361d2b3d, 0x1939260f, 0x19c27960,
220 			0x5223a708, 0xf71312b6, 0xebadfe6e, 0xeac31f66,
221 			0xe3bc4595, 0xa67bc883, 0xb17f37d1, 0x018cff28,
222 			0xc332ddef, 0xbe6c5aa5, 0x65582185, 0x68ab9802,
223 			0xeecea50f, 0xdb2f953b, 0x2aef7dad, 0x5b6e2f84,
224 			0x1521b628, 0x29076170, 0xecdd4775, 0x619f1510,
225 			0x13cca830, 0xeb61bd96, 0x0334fe1e, 0xaa0363cf,
226 			0xb5735c90, 0x4c70a239, 0xd59e9e0b, 0xcbaade14,
227 			0xeecc86bc, 0x60622ca7, 0x9cab5cab, 0xb2f3846e,
228 			0x648b1eaf, 0x19bdf0ca, 0xa02369b9, 0x655abb50,
229 			0x40685a32, 0x3c2ab4b3, 0x319ee9d5, 0xc021b8f7,
230 			0x9b540b19, 0x875fa099, 0x95f7997e, 0x623d7da8,
231 			0xf837889a, 0x97e32d77, 0x11ed935f, 0x16681281,
232 			0x0e358829, 0xc7e61fd6, 0x96dedfa1, 0x7858ba99,
233 			0x57f584a5, 0x1b227263, 0x9b83c3ff, 0x1ac24696,
234 			0xcdb30aeb, 0x532e3054, 0x8fd948e4, 0x6dbc3128,
235 			0x58ebf2ef, 0x34c6ffea, 0xfe28ed61, 0xee7c3c73,
236 			0x5d4a14d9, 0xe864b7e3, 0x42105d14, 0x203e13e0,
237 			0x45eee2b6, 0xa3aaabea, 0xdb6c4f15, 0xfacb4fd0,
238 			0xc742f442, 0xef6abbb5, 0x654f3b1d, 0x41cd2105,
239 			0xd81e799e, 0x86854dc7, 0xe44b476a, 0x3d816250,
240 			0xcf62a1f2, 0x5b8d2646, 0xfc8883a0, 0xc1c7b6a3,
241 			0x7f1524c3, 0x69cb7492, 0x47848a0b, 0x5692b285,
242 			0x095bbf00, 0xad19489d, 0x1462b174, 0x23820e00,
243 			0x58428d2a, 0x0c55f5ea, 0x1dadf43e, 0x233f7061,
244 			0x3372f092, 0x8d937e41, 0xd65fecf1, 0x6c223bdb,
245 			0x7cde3759, 0xcbee7460, 0x4085f2a7, 0xce77326e,
246 			0xa6078084, 0x19f8509e, 0xe8efd855, 0x61d99735,
247 			0xa969a7aa, 0xc50c06c2, 0x5a04abfc, 0x800bcadc,
248 			0x9e447a2e, 0xc3453484, 0xfdd56705, 0x0e1e9ec9,
249 			0xdb73dbd3, 0x105588cd, 0x675fda79, 0xe3674340,
250 			0xc5c43465, 0x713e38d8, 0x3d28f89e, 0xf16dff20,
251 		0x153e21e7, 0x8fb03d4a, 0xe6e39f2b, 0xdb83adf7},
252 		{
253 			0xe93d5a68, 0x948140f7, 0xf64c261c, 0x94692934,
254 			0x411520f7, 0x7602d4f7, 0xbcf46b2e, 0xd4a20068,
255 			0xd4082471, 0x3320f46a, 0x43b7d4b7, 0x500061af,
256 			0x1e39f62e, 0x97244546, 0x14214f74, 0xbf8b8840,
257 			0x4d95fc1d, 0x96b591af, 0x70f4ddd3, 0x66a02f45,
258 			0xbfbc09ec, 0x03bd9785, 0x7fac6dd0, 0x31cb8504,
259 			0x96eb27b3, 0x55fd3941, 0xda2547e6, 0xabca0a9a,
260 			0x28507825, 0x530429f4, 0x0a2c86da, 0xe9b66dfb,
261 			0x68dc1462, 0xd7486900, 0x680ec0a4, 0x27a18dee,
262 			0x4f3ffea2, 0xe887ad8c, 0xb58ce006, 0x7af4d6b6,
263 			0xaace1e7c, 0xd3375fec, 0xce78a399, 0x406b2a42,
264 			0x20fe9e35, 0xd9f385b9, 0xee39d7ab, 0x3b124e8b,
265 			0x1dc9faf7, 0x4b6d1856, 0x26a36631, 0xeae397b2,
266 			0x3a6efa74, 0xdd5b4332, 0x6841e7f7, 0xca7820fb,
267 			0xfb0af54e, 0xd8feb397, 0x454056ac, 0xba489527,
268 			0x55533a3a, 0x20838d87, 0xfe6ba9b7, 0xd096954b,
269 			0x55a867bc, 0xa1159a58, 0xcca92963, 0x99e1db33,
270 			0xa62a4a56, 0x3f3125f9, 0x5ef47e1c, 0x9029317c,
271 			0xfdf8e802, 0x04272f70, 0x80bb155c, 0x05282ce3,
272 			0x95c11548, 0xe4c66d22, 0x48c1133f, 0xc70f86dc,
273 			0x07f9c9ee, 0x41041f0f, 0x404779a4, 0x5d886e17,
274 			0x325f51eb, 0xd59bc0d1, 0xf2bcc18f, 0x41113564,
275 			0x257b7834, 0x602a9c60, 0xdff8e8a3, 0x1f636c1b,
276 			0x0e12b4c2, 0x02e1329e, 0xaf664fd1, 0xcad18115,
277 			0x6b2395e0, 0x333e92e1, 0x3b240b62, 0xeebeb922,
278 			0x85b2a20e, 0xe6ba0d99, 0xde720c8c, 0x2da2f728,
279 			0xd0127845, 0x95b794fd, 0x647d0862, 0xe7ccf5f0,
280 			0x5449a36f, 0x877d48fa, 0xc39dfd27, 0xf33e8d1e,
281 			0x0a476341, 0x992eff74, 0x3a6f6eab, 0xf4f8fd37,
282 			0xa812dc60, 0xa1ebddf8, 0x991be14c, 0xdb6e6b0d,
283 			0xc67b5510, 0x6d672c37, 0x2765d43b, 0xdcd0e804,
284 			0xf1290dc7, 0xcc00ffa3, 0xb5390f92, 0x690fed0b,
285 			0x667b9ffb, 0xcedb7d9c, 0xa091cf0b, 0xd9155ea3,
286 			0xbb132f88, 0x515bad24, 0x7b9479bf, 0x763bd6eb,
287 			0x37392eb3, 0xcc115979, 0x8026e297, 0xf42e312d,
288 			0x6842ada7, 0xc66a2b3b, 0x12754ccc, 0x782ef11c,
289 			0x6a124237, 0xb79251e7, 0x06a1bbe6, 0x4bfb6350,
290 			0x1a6b1018, 0x11caedfa, 0x3d25bdd8, 0xe2e1c3c9,
291 			0x44421659, 0x0a121386, 0xd90cec6e, 0xd5abea2a,
292 			0x64af674e, 0xda86a85f, 0xbebfe988, 0x64e4c3fe,
293 			0x9dbc8057, 0xf0f7c086, 0x60787bf8, 0x6003604d,
294 			0xd1fd8346, 0xf6381fb0, 0x7745ae04, 0xd736fccc,
295 			0x83426b33, 0xf01eab71, 0xb0804187, 0x3c005e5f,
296 			0x77a057be, 0xbde8ae24, 0x55464299, 0xbf582e61,
297 			0x4e58f48f, 0xf2ddfda2, 0xf474ef38, 0x8789bdc2,
298 			0x5366f9c3, 0xc8b38e74, 0xb475f255, 0x46fcd9b9,
299 			0x7aeb2661, 0x8b1ddf84, 0x846a0e79, 0x915f95e2,
300 			0x466e598e, 0x20b45770, 0x8cd55591, 0xc902de4c,
301 			0xb90bace1, 0xbb8205d0, 0x11a86248, 0x7574a99e,
302 			0xb77f19b6, 0xe0a9dc09, 0x662d09a1, 0xc4324633,
303 			0xe85a1f02, 0x09f0be8c, 0x4a99a025, 0x1d6efe10,
304 			0x1ab93d1d, 0x0ba5a4df, 0xa186f20f, 0x2868f169,
305 			0xdcb7da83, 0x573906fe, 0xa1e2ce9b, 0x4fcd7f52,
306 			0x50115e01, 0xa70683fa, 0xa002b5c4, 0x0de6d027,
307 			0x9af88c27, 0x773f8641, 0xc3604c06, 0x61a806b5,
308 			0xf0177a28, 0xc0f586e0, 0x006058aa, 0x30dc7d62,
309 			0x11e69ed7, 0x2338ea63, 0x53c2dd94, 0xc2c21634,
310 			0xbbcbee56, 0x90bcb6de, 0xebfc7da1, 0xce591d76,
311 			0x6f05e409, 0x4b7c0188, 0x39720a3d, 0x7c927c24,
312 			0x86e3725f, 0x724d9db9, 0x1ac15bb4, 0xd39eb8fc,
313 			0xed545578, 0x08fca5b5, 0xd83d7cd3, 0x4dad0fc4,
314 			0x1e50ef5e, 0xb161e6f8, 0xa28514d9, 0x6c51133c,
315 			0x6fd5c7e7, 0x56e14ec4, 0x362abfce, 0xddc6c837,
316 		0xd79a3234, 0x92638212, 0x670efa8e, 0x406000e0},
317 		{
318 			0x3a39ce37, 0xd3faf5cf, 0xabc27737, 0x5ac52d1b,
319 			0x5cb0679e, 0x4fa33742, 0xd3822740, 0x99bc9bbe,
320 			0xd5118e9d, 0xbf0f7315, 0xd62d1c7e, 0xc700c47b,
321 			0xb78c1b6b, 0x21a19045, 0xb26eb1be, 0x6a366eb4,
322 			0x5748ab2f, 0xbc946e79, 0xc6a376d2, 0x6549c2c8,
323 			0x530ff8ee, 0x468dde7d, 0xd5730a1d, 0x4cd04dc6,
324 			0x2939bbdb, 0xa9ba4650, 0xac9526e8, 0xbe5ee304,
325 			0xa1fad5f0, 0x6a2d519a, 0x63ef8ce2, 0x9a86ee22,
326 			0xc089c2b8, 0x43242ef6, 0xa51e03aa, 0x9cf2d0a4,
327 			0x83c061ba, 0x9be96a4d, 0x8fe51550, 0xba645bd6,
328 			0x2826a2f9, 0xa73a3ae1, 0x4ba99586, 0xef5562e9,
329 			0xc72fefd3, 0xf752f7da, 0x3f046f69, 0x77fa0a59,
330 			0x80e4a915, 0x87b08601, 0x9b09e6ad, 0x3b3ee593,
331 			0xe990fd5a, 0x9e34d797, 0x2cf0b7d9, 0x022b8b51,
332 			0x96d5ac3a, 0x017da67d, 0xd1cf3ed6, 0x7c7d2d28,
333 			0x1f9f25cf, 0xadf2b89b, 0x5ad6b472, 0x5a88f54c,
334 			0xe029ac71, 0xe019a5e6, 0x47b0acfd, 0xed93fa9b,
335 			0xe8d3c48d, 0x283b57cc, 0xf8d56629, 0x79132e28,
336 			0x785f0191, 0xed756055, 0xf7960e44, 0xe3d35e8c,
337 			0x15056dd4, 0x88f46dba, 0x03a16125, 0x0564f0bd,
338 			0xc3eb9e15, 0x3c9057a2, 0x97271aec, 0xa93a072a,
339 			0x1b3f6d9b, 0x1e6321f5, 0xf59c66fb, 0x26dcf319,
340 			0x7533d928, 0xb155fdf5, 0x03563482, 0x8aba3cbb,
341 			0x28517711, 0xc20ad9f8, 0xabcc5167, 0xccad925f,
342 			0x4de81751, 0x3830dc8e, 0x379d5862, 0x9320f991,
343 			0xea7a90c2, 0xfb3e7bce, 0x5121ce64, 0x774fbe32,
344 			0xa8b6e37e, 0xc3293d46, 0x48de5369, 0x6413e680,
345 			0xa2ae0810, 0xdd6db224, 0x69852dfd, 0x09072166,
346 			0xb39a460a, 0x6445c0dd, 0x586cdecf, 0x1c20c8ae,
347 			0x5bbef7dd, 0x1b588d40, 0xccd2017f, 0x6bb4e3bb,
348 			0xdda26a7e, 0x3a59ff45, 0x3e350a44, 0xbcb4cdd5,
349 			0x72eacea8, 0xfa6484bb, 0x8d6612ae, 0xbf3c6f47,
350 			0xd29be463, 0x542f5d9e, 0xaec2771b, 0xf64e6370,
351 			0x740e0d8d, 0xe75b1357, 0xf8721671, 0xaf537d5d,
352 			0x4040cb08, 0x4eb4e2cc, 0x34d2466a, 0x0115af84,
353 			0xe1b00428, 0x95983a1d, 0x06b89fb4, 0xce6ea048,
354 			0x6f3f3b82, 0x3520ab82, 0x011a1d4b, 0x277227f8,
355 			0x611560b1, 0xe7933fdc, 0xbb3a792b, 0x344525bd,
356 			0xa08839e1, 0x51ce794b, 0x2f32c9b7, 0xa01fbac9,
357 			0xe01cc87e, 0xbcc7d1f6, 0xcf0111c3, 0xa1e8aac7,
358 			0x1a908749, 0xd44fbd9a, 0xd0dadecb, 0xd50ada38,
359 			0x0339c32a, 0xc6913667, 0x8df9317c, 0xe0b12b4f,
360 			0xf79e59b7, 0x43f5bb3a, 0xf2d519ff, 0x27d9459c,
361 			0xbf97222c, 0x15e6fc2a, 0x0f91fc71, 0x9b941525,
362 			0xfae59361, 0xceb69ceb, 0xc2a86459, 0x12baa8d1,
363 			0xb6c1075e, 0xe3056a0c, 0x10d25065, 0xcb03a442,
364 			0xe0ec6e0e, 0x1698db3b, 0x4c98a0be, 0x3278e964,
365 			0x9f1f9532, 0xe0d392df, 0xd3a0342b, 0x8971f21e,
366 			0x1b0a7441, 0x4ba3348c, 0xc5be7120, 0xc37632d8,
367 			0xdf359f8d, 0x9b992f2e, 0xe60b6f47, 0x0fe3f11d,
368 			0xe54cda54, 0x1edad891, 0xce6279cf, 0xcd3e7e6f,
369 			0x1618b166, 0xfd2c1d05, 0x848fd2c5, 0xf6fb2299,
370 			0xf523f357, 0xa6327623, 0x93a83531, 0x56cccd02,
371 			0xacf08162, 0x5a75ebb5, 0x6e163697, 0x88d273cc,
372 			0xde966292, 0x81b949d0, 0x4c50901b, 0x71c65614,
373 			0xe6c6c7bd, 0x327a140a, 0x45e1d006, 0xc3f27b9a,
374 			0xc9aa53fd, 0x62a80f00, 0xbb25bfe2, 0x35bdd2f6,
375 			0x71126905, 0xb2040222, 0xb6cbcf7c, 0xcd769c2b,
376 			0x53113ec0, 0x1640e3d3, 0x38abbd60, 0x2547adf0,
377 			0xba38209c, 0xf746ce76, 0x77afa1c5, 0x20756060,
378 			0x85cbfe4e, 0x8ae88dd8, 0x7aaaf9b0, 0x4cf9aa7e,
379 			0x1948c25c, 0x02fb8a8c, 0x01c36ae4, 0xd6ebe1f9,
380 			0x90d4f869, 0xa65cdea0, 0x3f09252d, 0xc208e69f,
381 		0xb74e6132, 0xce77e25b, 0x578fdfe3, 0x3ac372e6}
382 	},
383 	{
384 		0x243f6a88, 0x85a308d3, 0x13198a2e, 0x03707344,
385 		0xa4093822, 0x299f31d0, 0x082efa98, 0xec4e6c89,
386 		0x452821e6, 0x38d01377, 0xbe5466cf, 0x34e90c6c,
387 		0xc0ac29b7, 0xc97c50dd, 0x3f84d5b5, 0xb5470917,
388 		0x9216d5d9, 0x8979fb1b
389 	} };
390 
391 	*c = initstate;
392 }
393 
394 u_int32_t
395 Blowfish_stream2word(const u_int8_t *data, u_int16_t databytes,
396     u_int16_t *current)
397 {
398 	u_int8_t i;
399 	u_int16_t j;
400 	u_int32_t temp;
401 
402 	temp = 0x00000000;
403 	j = *current;
404 
405 	for (i = 0; i < 4; i++, j++) {
406 		if (j >= databytes)
407 			j = 0;
408 		temp = (temp << 8) | data[j];
409 	}
410 
411 	*current = j;
412 	return temp;
413 }
414 
415 void
416 Blowfish_expand0state(blf_ctx *c, const u_int8_t *key, u_int16_t keybytes)
417 {
418 	u_int16_t i;
419 	u_int16_t j;
420 	u_int16_t k;
421 	u_int32_t temp;
422 	u_int32_t data[2];
423 
424 	j = 0;
425 	for (i = 0; i < BLF_N + 2; i++) {
426 		/* Extract 4 int8 to 1 int32 from keystream */
427 		temp = Blowfish_stream2word(key, keybytes, &j);
428 		c->P[i] = c->P[i] ^ temp;
429 	}
430 
431 	j = 0;
432 	data[0] = 0x00000000;
433 	data[1] = 0x00000000;
434 	for (i = 0; i < BLF_N + 2; i += 2) {
435 		Blowfish_encipher(c, data);
436 
437 		c->P[i] = data[0];
438 		c->P[i + 1] = data[1];
439 	}
440 
441 	for (i = 0; i < 4; i++) {
442 		for (k = 0; k < 256; k += 2) {
443 			Blowfish_encipher(c, data);
444 
445 			c->S[i][k] = data[0];
446 			c->S[i][k + 1] = data[1];
447 		}
448 	}
449 }
450 
451 
452 void
453 Blowfish_expandstate(blf_ctx *c, const u_int8_t *data, u_int16_t databytes,
454     const u_int8_t *key, u_int16_t keybytes)
455 {
456 	u_int16_t i;
457 	u_int16_t j;
458 	u_int16_t k;
459 	u_int32_t temp;
460 	u_int32_t d[2];
461 
462 	j = 0;
463 	for (i = 0; i < BLF_N + 2; i++) {
464 		/* Extract 4 int8 to 1 int32 from keystream */
465 		temp = Blowfish_stream2word(key, keybytes, &j);
466 		c->P[i] = c->P[i] ^ temp;
467 	}
468 
469 	j = 0;
470 	d[0] = 0x00000000;
471 	d[1] = 0x00000000;
472 	for (i = 0; i < BLF_N + 2; i += 2) {
473 		d[0] ^= Blowfish_stream2word(data, databytes, &j);
474 		d[1] ^= Blowfish_stream2word(data, databytes, &j);
475 		Blowfish_encipher(c, d);
476 
477 		c->P[i] = d[0];
478 		c->P[i + 1] = d[1];
479 	}
480 
481 	for (i = 0; i < 4; i++) {
482 		for (k = 0; k < 256; k += 2) {
483 			d[0]^= Blowfish_stream2word(data, databytes, &j);
484 			d[1] ^= Blowfish_stream2word(data, databytes, &j);
485 			Blowfish_encipher(c, d);
486 
487 			c->S[i][k] = d[0];
488 			c->S[i][k + 1] = d[1];
489 		}
490 	}
491 
492 }
493 
494 void
495 blf_key(blf_ctx *c, const u_int8_t *k, u_int16_t len)
496 {
497 	/* Initialize S-boxes and subkeys with Pi */
498 	Blowfish_initstate(c);
499 
500 	/* Transform S-boxes and subkeys with key */
501 	Blowfish_expand0state(c, k, len);
502 }
503 
504 void
505 blf_enc(blf_ctx *c, u_int32_t *data, u_int16_t blocks)
506 {
507 	u_int32_t *d;
508 	u_int16_t i;
509 
510 	d = data;
511 	for (i = 0; i < blocks; i++) {
512 		Blowfish_encipher(c, d);
513 		d += 2;
514 	}
515 }
516 
517 void
518 blf_dec(blf_ctx *c, u_int32_t *data, u_int16_t blocks)
519 {
520 	u_int32_t *d;
521 	u_int16_t i;
522 
523 	d = data;
524 	for (i = 0; i < blocks; i++) {
525 		Blowfish_decipher(c, d);
526 		d += 2;
527 	}
528 }
529 
530 void
531 blf_ecb_encrypt(blf_ctx *c, u_int8_t *data, u_int32_t len)
532 {
533 	u_int32_t l, r, d[2];
534 	u_int32_t i;
535 
536 	for (i = 0; i < len; i += 8) {
537 		l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
538 		r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
539 		d[0] = l;
540 		d[1] = r;
541 		Blowfish_encipher(c, d);
542 		l = d[0];
543 		r = d[1];
544 		data[0] = l >> 24 & 0xff;
545 		data[1] = l >> 16 & 0xff;
546 		data[2] = l >> 8 & 0xff;
547 		data[3] = l & 0xff;
548 		data[4] = r >> 24 & 0xff;
549 		data[5] = r >> 16 & 0xff;
550 		data[6] = r >> 8 & 0xff;
551 		data[7] = r & 0xff;
552 		data += 8;
553 	}
554 }
555 
556 void
557 blf_ecb_decrypt(blf_ctx *c, u_int8_t *data, u_int32_t len)
558 {
559 	u_int32_t l, r, d[2];
560 	u_int32_t i;
561 
562 	for (i = 0; i < len; i += 8) {
563 		l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
564 		r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
565 		d[0] = l;
566 		d[1] = r;
567 		Blowfish_decipher(c, d);
568 		l = d[0];
569 		r = d[1];
570 		data[0] = l >> 24 & 0xff;
571 		data[1] = l >> 16 & 0xff;
572 		data[2] = l >> 8 & 0xff;
573 		data[3] = l & 0xff;
574 		data[4] = r >> 24 & 0xff;
575 		data[5] = r >> 16 & 0xff;
576 		data[6] = r >> 8 & 0xff;
577 		data[7] = r & 0xff;
578 		data += 8;
579 	}
580 }
581 
582 void
583 blf_cbc_encrypt(blf_ctx *c, u_int8_t *iv, u_int8_t *data, u_int32_t len)
584 {
585 	u_int32_t l, r, d[2];
586 	u_int32_t i, j;
587 
588 	for (i = 0; i < len; i += 8) {
589 		for (j = 0; j < 8; j++)
590 			data[j] ^= iv[j];
591 		l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
592 		r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
593 		d[0] = l;
594 		d[1] = r;
595 		Blowfish_encipher(c, d);
596 		l = d[0];
597 		r = d[1];
598 		data[0] = l >> 24 & 0xff;
599 		data[1] = l >> 16 & 0xff;
600 		data[2] = l >> 8 & 0xff;
601 		data[3] = l & 0xff;
602 		data[4] = r >> 24 & 0xff;
603 		data[5] = r >> 16 & 0xff;
604 		data[6] = r >> 8 & 0xff;
605 		data[7] = r & 0xff;
606 		iv = data;
607 		data += 8;
608 	}
609 }
610 
611 void
612 blf_cbc_decrypt(blf_ctx *c, u_int8_t *iva, u_int8_t *data, u_int32_t len)
613 {
614 	u_int32_t l, r, d[2];
615 	u_int8_t *iv;
616 	u_int32_t i, j;
617 
618 	iv = data + len - 16;
619 	data = data + len - 8;
620 	for (i = len - 8; i >= 8; i -= 8) {
621 		l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
622 		r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
623 		d[0] = l;
624 		d[1] = r;
625 		Blowfish_decipher(c, d);
626 		l = d[0];
627 		r = d[1];
628 		data[0] = l >> 24 & 0xff;
629 		data[1] = l >> 16 & 0xff;
630 		data[2] = l >> 8 & 0xff;
631 		data[3] = l & 0xff;
632 		data[4] = r >> 24 & 0xff;
633 		data[5] = r >> 16 & 0xff;
634 		data[6] = r >> 8 & 0xff;
635 		data[7] = r & 0xff;
636 		for (j = 0; j < 8; j++)
637 			data[j] ^= iv[j];
638 		iv -= 8;
639 		data -= 8;
640 	}
641 	l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
642 	r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
643 	d[0] = l;
644 	d[1] = r;
645 	Blowfish_decipher(c, d);
646 	l = d[0];
647 	r = d[1];
648 	data[0] = l >> 24 & 0xff;
649 	data[1] = l >> 16 & 0xff;
650 	data[2] = l >> 8 & 0xff;
651 	data[3] = l & 0xff;
652 	data[4] = r >> 24 & 0xff;
653 	data[5] = r >> 16 & 0xff;
654 	data[6] = r >> 8 & 0xff;
655 	data[7] = r & 0xff;
656 	for (j = 0; j < 8; j++)
657 		data[j] ^= iva[j];
658 }
659