xref: /openbsd-src/lib/libc/crypt/blowfish.c (revision 5a38ef86d0b61900239c7913d24a05e7b88a58f0)
1 /* $OpenBSD: blowfish.c,v 1.20 2021/11/29 01:04:45 djm Exp $ */
2 /*
3  * Blowfish block cipher for OpenBSD
4  * Copyright 1997 Niels Provos <provos@physnet.uni-hamburg.de>
5  * All rights reserved.
6  *
7  * Implementation advice by David Mazieres <dm@lcs.mit.edu>.
8  *
9  * Redistribution and use in source and binary forms, with or without
10  * modification, are permitted provided that the following conditions
11  * are met:
12  * 1. Redistributions of source code must retain the above copyright
13  *    notice, this list of conditions and the following disclaimer.
14  * 2. Redistributions in binary form must reproduce the above copyright
15  *    notice, this list of conditions and the following disclaimer in the
16  *    documentation and/or other materials provided with the distribution.
17  * 3. The name of the author may not be used to endorse or promote products
18  *    derived from this software without specific prior written permission.
19  *
20  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
21  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
22  * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
23  * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
24  * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
25  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
26  * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
27  * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
28  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
29  * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
30  */
31 
32 /*
33  * This code is derived from section 14.3 and the given source
34  * in section V of Applied Cryptography, second edition.
35  * Blowfish is an unpatented fast block cipher designed by
36  * Bruce Schneier.
37  */
38 
39 #if 0
40 #include <stdio.h>		/* used for debugging */
41 #include <string.h>
42 #endif
43 
44 #include <sys/types.h>
45 #include <blf.h>
46 
47 #undef inline
48 #ifdef __GNUC__
49 #define inline __inline
50 #else				/* !__GNUC__ */
51 #define inline
52 #endif				/* !__GNUC__ */
53 
54 /* Function for Feistel Networks */
55 
56 #define F(s, x) ((((s)[        (((x)>>24)&0xFF)]  \
57 		 + (s)[0x100 + (((x)>>16)&0xFF)]) \
58 		 ^ (s)[0x200 + (((x)>> 8)&0xFF)]) \
59 		 + (s)[0x300 + ( (x)     &0xFF)])
60 
61 #define BLFRND(s,p,i,j,n) (i ^= F(s,j) ^ (p)[n])
62 
63 void
64 Blowfish_encipher(blf_ctx *c, u_int32_t *xl, u_int32_t *xr)
65 {
66 	u_int32_t Xl;
67 	u_int32_t Xr;
68 	u_int32_t *s = c->S[0];
69 	u_int32_t *p = c->P;
70 
71 	Xl = *xl;
72 	Xr = *xr;
73 
74 	Xl ^= p[0];
75 	BLFRND(s, p, Xr, Xl, 1); BLFRND(s, p, Xl, Xr, 2);
76 	BLFRND(s, p, Xr, Xl, 3); BLFRND(s, p, Xl, Xr, 4);
77 	BLFRND(s, p, Xr, Xl, 5); BLFRND(s, p, Xl, Xr, 6);
78 	BLFRND(s, p, Xr, Xl, 7); BLFRND(s, p, Xl, Xr, 8);
79 	BLFRND(s, p, Xr, Xl, 9); BLFRND(s, p, Xl, Xr, 10);
80 	BLFRND(s, p, Xr, Xl, 11); BLFRND(s, p, Xl, Xr, 12);
81 	BLFRND(s, p, Xr, Xl, 13); BLFRND(s, p, Xl, Xr, 14);
82 	BLFRND(s, p, Xr, Xl, 15); BLFRND(s, p, Xl, Xr, 16);
83 
84 	*xl = Xr ^ p[17];
85 	*xr = Xl;
86 }
87 DEF_WEAK(Blowfish_encipher);
88 
89 void
90 Blowfish_decipher(blf_ctx *c, u_int32_t *xl, u_int32_t *xr)
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 = *xl;
98 	Xr = *xr;
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 	*xl = Xr ^ p[0];
111 	*xr = Xl;
112 }
113 DEF_WEAK(Blowfish_decipher);
114 
115 void
116 Blowfish_initstate(blf_ctx *c)
117 {
118 	/* P-box and S-box tables initialized with digits of Pi */
119 
120 	static const blf_ctx initstate =
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 DEF_WEAK(Blowfish_initstate);
394 
395 u_int32_t
396 Blowfish_stream2word(const u_int8_t *data, u_int16_t databytes,
397     u_int16_t *current)
398 {
399 	u_int8_t i;
400 	u_int16_t j;
401 	u_int32_t temp;
402 
403 	temp = 0x00000000;
404 	j = *current;
405 
406 	for (i = 0; i < 4; i++, j++) {
407 		if (j >= databytes)
408 			j = 0;
409 		temp = (temp << 8) | data[j];
410 	}
411 
412 	*current = j;
413 	return temp;
414 }
415 DEF_WEAK(Blowfish_stream2word);
416 
417 void
418 Blowfish_expand0state(blf_ctx *c, const u_int8_t *key, u_int16_t keybytes)
419 {
420 	u_int16_t i;
421 	u_int16_t j;
422 	u_int16_t k;
423 	u_int32_t temp;
424 	u_int32_t datal;
425 	u_int32_t datar;
426 
427 	j = 0;
428 	for (i = 0; i < BLF_N + 2; i++) {
429 		/* Extract 4 int8 to 1 int32 from keystream */
430 		temp = Blowfish_stream2word(key, keybytes, &j);
431 		c->P[i] = c->P[i] ^ temp;
432 	}
433 
434 	j = 0;
435 	datal = 0x00000000;
436 	datar = 0x00000000;
437 	for (i = 0; i < BLF_N + 2; i += 2) {
438 		Blowfish_encipher(c, &datal, &datar);
439 
440 		c->P[i] = datal;
441 		c->P[i + 1] = datar;
442 	}
443 
444 	for (i = 0; i < 4; i++) {
445 		for (k = 0; k < 256; k += 2) {
446 			Blowfish_encipher(c, &datal, &datar);
447 
448 			c->S[i][k] = datal;
449 			c->S[i][k + 1] = datar;
450 		}
451 	}
452 }
453 DEF_WEAK(Blowfish_expand0state);
454 
455 
456 void
457 Blowfish_expandstate(blf_ctx *c, const u_int8_t *data, u_int16_t databytes,
458     const u_int8_t *key, u_int16_t keybytes)
459 {
460 	u_int16_t i;
461 	u_int16_t j;
462 	u_int16_t k;
463 	u_int32_t temp;
464 	u_int32_t datal;
465 	u_int32_t datar;
466 
467 	j = 0;
468 	for (i = 0; i < BLF_N + 2; i++) {
469 		/* Extract 4 int8 to 1 int32 from keystream */
470 		temp = Blowfish_stream2word(key, keybytes, &j);
471 		c->P[i] = c->P[i] ^ temp;
472 	}
473 
474 	j = 0;
475 	datal = 0x00000000;
476 	datar = 0x00000000;
477 	for (i = 0; i < BLF_N + 2; i += 2) {
478 		datal ^= Blowfish_stream2word(data, databytes, &j);
479 		datar ^= Blowfish_stream2word(data, databytes, &j);
480 		Blowfish_encipher(c, &datal, &datar);
481 
482 		c->P[i] = datal;
483 		c->P[i + 1] = datar;
484 	}
485 
486 	for (i = 0; i < 4; i++) {
487 		for (k = 0; k < 256; k += 2) {
488 			datal ^= Blowfish_stream2word(data, databytes, &j);
489 			datar ^= Blowfish_stream2word(data, databytes, &j);
490 			Blowfish_encipher(c, &datal, &datar);
491 
492 			c->S[i][k] = datal;
493 			c->S[i][k + 1] = datar;
494 		}
495 	}
496 
497 }
498 DEF_WEAK(Blowfish_expandstate);
499 
500 void
501 blf_key(blf_ctx *c, const u_int8_t *k, u_int16_t len)
502 {
503 	/* Initialize S-boxes and subkeys with Pi */
504 	Blowfish_initstate(c);
505 
506 	/* Transform S-boxes and subkeys with key */
507 	Blowfish_expand0state(c, k, len);
508 }
509 DEF_WEAK(blf_key);
510 
511 void
512 blf_enc(blf_ctx *c, u_int32_t *data, u_int16_t blocks)
513 {
514 	u_int32_t *d;
515 	u_int16_t i;
516 
517 	d = data;
518 	for (i = 0; i < blocks; i++) {
519 		Blowfish_encipher(c, d, d + 1);
520 		d += 2;
521 	}
522 }
523 DEF_WEAK(blf_enc);
524 
525 void
526 blf_dec(blf_ctx *c, u_int32_t *data, u_int16_t blocks)
527 {
528 	u_int32_t *d;
529 	u_int16_t i;
530 
531 	d = data;
532 	for (i = 0; i < blocks; i++) {
533 		Blowfish_decipher(c, d, d + 1);
534 		d += 2;
535 	}
536 }
537 DEF_WEAK(blf_dec);
538 
539 void
540 blf_ecb_encrypt(blf_ctx *c, u_int8_t *data, u_int32_t len)
541 {
542 	u_int32_t l, r;
543 	u_int32_t i;
544 
545 	for (i = 0; i < len; i += 8) {
546 		l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
547 		r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
548 		Blowfish_encipher(c, &l, &r);
549 		data[0] = l >> 24 & 0xff;
550 		data[1] = l >> 16 & 0xff;
551 		data[2] = l >> 8 & 0xff;
552 		data[3] = l & 0xff;
553 		data[4] = r >> 24 & 0xff;
554 		data[5] = r >> 16 & 0xff;
555 		data[6] = r >> 8 & 0xff;
556 		data[7] = r & 0xff;
557 		data += 8;
558 	}
559 }
560 DEF_WEAK(blf_ecb_encrypt);
561 
562 void
563 blf_ecb_decrypt(blf_ctx *c, u_int8_t *data, u_int32_t len)
564 {
565 	u_int32_t l, r;
566 	u_int32_t i;
567 
568 	for (i = 0; i < len; i += 8) {
569 		l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
570 		r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
571 		Blowfish_decipher(c, &l, &r);
572 		data[0] = l >> 24 & 0xff;
573 		data[1] = l >> 16 & 0xff;
574 		data[2] = l >> 8 & 0xff;
575 		data[3] = l & 0xff;
576 		data[4] = r >> 24 & 0xff;
577 		data[5] = r >> 16 & 0xff;
578 		data[6] = r >> 8 & 0xff;
579 		data[7] = r & 0xff;
580 		data += 8;
581 	}
582 }
583 DEF_WEAK(blf_ecb_decrypt);
584 
585 void
586 blf_cbc_encrypt(blf_ctx *c, u_int8_t *iv, u_int8_t *data, u_int32_t len)
587 {
588 	u_int32_t l, r;
589 	u_int32_t i, j;
590 
591 	for (i = 0; i < len; i += 8) {
592 		for (j = 0; j < 8; j++)
593 			data[j] ^= iv[j];
594 		l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
595 		r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
596 		Blowfish_encipher(c, &l, &r);
597 		data[0] = l >> 24 & 0xff;
598 		data[1] = l >> 16 & 0xff;
599 		data[2] = l >> 8 & 0xff;
600 		data[3] = l & 0xff;
601 		data[4] = r >> 24 & 0xff;
602 		data[5] = r >> 16 & 0xff;
603 		data[6] = r >> 8 & 0xff;
604 		data[7] = r & 0xff;
605 		iv = data;
606 		data += 8;
607 	}
608 }
609 DEF_WEAK(blf_cbc_encrypt);
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;
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 		Blowfish_decipher(c, &l, &r);
624 		data[0] = l >> 24 & 0xff;
625 		data[1] = l >> 16 & 0xff;
626 		data[2] = l >> 8 & 0xff;
627 		data[3] = l & 0xff;
628 		data[4] = r >> 24 & 0xff;
629 		data[5] = r >> 16 & 0xff;
630 		data[6] = r >> 8 & 0xff;
631 		data[7] = r & 0xff;
632 		for (j = 0; j < 8; j++)
633 			data[j] ^= iv[j];
634 		iv -= 8;
635 		data -= 8;
636 	}
637 	l = data[0] << 24 | data[1] << 16 | data[2] << 8 | data[3];
638 	r = data[4] << 24 | data[5] << 16 | data[6] << 8 | data[7];
639 	Blowfish_decipher(c, &l, &r);
640 	data[0] = l >> 24 & 0xff;
641 	data[1] = l >> 16 & 0xff;
642 	data[2] = l >> 8 & 0xff;
643 	data[3] = l & 0xff;
644 	data[4] = r >> 24 & 0xff;
645 	data[5] = r >> 16 & 0xff;
646 	data[6] = r >> 8 & 0xff;
647 	data[7] = r & 0xff;
648 	for (j = 0; j < 8; j++)
649 		data[j] ^= iva[j];
650 }
651 DEF_WEAK(blf_cbc_decrypt);
652 
653 #if 0
654 void
655 report(u_int32_t data[], u_int16_t len)
656 {
657 	u_int16_t i;
658 	for (i = 0; i < len; i += 2)
659 		printf("Block %0hd: %08lx %08lx.\n",
660 		    i / 2, data[i], data[i + 1]);
661 }
662 void
663 main(void)
664 {
665 
666 	blf_ctx c;
667 	char    key[] = "AAAAA";
668 	char    key2[] = "abcdefghijklmnopqrstuvwxyz";
669 
670 	u_int32_t data[10];
671 	u_int32_t data2[] =
672 	{0x424c4f57l, 0x46495348l};
673 
674 	u_int16_t i;
675 
676 	/* First test */
677 	for (i = 0; i < 10; i++)
678 		data[i] = i;
679 
680 	blf_key(&c, (u_int8_t *) key, 5);
681 	blf_enc(&c, data, 5);
682 	blf_dec(&c, data, 1);
683 	blf_dec(&c, data + 2, 4);
684 	printf("Should read as 0 - 9.\n");
685 	report(data, 10);
686 
687 	/* Second test */
688 	blf_key(&c, (u_int8_t *) key2, strlen(key2));
689 	blf_enc(&c, data2, 1);
690 	printf("\nShould read as: 0x324ed0fe 0xf413a203.\n");
691 	report(data2, 2);
692 	blf_dec(&c, data2, 1);
693 	report(data2, 2);
694 }
695 #endif
696