xref: /plan9/sys/src/cmd/unix/drawterm/posix-amd64/md5block.c (revision ec59a3ddbfceee0efe34584c2c9981a5e5ff1ec4)
1*ec59a3ddSDavid du Colombier #include <u.h>
2*ec59a3ddSDavid du Colombier #include <libc.h>
3*ec59a3ddSDavid du Colombier #include <libsec.h>
4*ec59a3ddSDavid du Colombier 
5*ec59a3ddSDavid du Colombier /*
6*ec59a3ddSDavid du Colombier  *  rfc1321 requires that I include this.  The code is new.  The constants
7*ec59a3ddSDavid du Colombier  *  all come from the rfc (hence the copyright).  We trade a table for the
8*ec59a3ddSDavid du Colombier  *  macros in rfc.  The total size is a lot less. -- presotto
9*ec59a3ddSDavid du Colombier  *
10*ec59a3ddSDavid du Colombier  *	Copyright (C) 1991-2, RSA Data Security, Inc. Created 1991. All
11*ec59a3ddSDavid du Colombier  *	rights reserved.
12*ec59a3ddSDavid du Colombier  *
13*ec59a3ddSDavid du Colombier  *	License to copy and use this software is granted provided that it
14*ec59a3ddSDavid du Colombier  *	is identified as the "RSA Data Security, Inc. MD5 Message-Digest
15*ec59a3ddSDavid du Colombier  *	Algorithm" in all material mentioning or referencing this software
16*ec59a3ddSDavid du Colombier  *	or this function.
17*ec59a3ddSDavid du Colombier  *
18*ec59a3ddSDavid du Colombier  *	License is also granted to make and use derivative works provided
19*ec59a3ddSDavid du Colombier  *	that such works are identified as "derived from the RSA Data
20*ec59a3ddSDavid du Colombier  *	Security, Inc. MD5 Message-Digest Algorithm" in all material
21*ec59a3ddSDavid du Colombier  *	mentioning or referencing the derived work.
22*ec59a3ddSDavid du Colombier  *
23*ec59a3ddSDavid du Colombier  *	RSA Data Security, Inc. makes no representations concerning either
24*ec59a3ddSDavid du Colombier  *	the merchantability of this software or the suitability of this
25*ec59a3ddSDavid du Colombier  *	software forany particular purpose. It is provided "as is"
26*ec59a3ddSDavid du Colombier  *	without express or implied warranty of any kind.
27*ec59a3ddSDavid du Colombier  *	These notices must be retained in any copies of any part of this
28*ec59a3ddSDavid du Colombier  *	documentation and/or software.
29*ec59a3ddSDavid du Colombier  */
30*ec59a3ddSDavid du Colombier 
31*ec59a3ddSDavid du Colombier /*
32*ec59a3ddSDavid du Colombier  *	Rotate ammounts used in the algorithm
33*ec59a3ddSDavid du Colombier  */
34*ec59a3ddSDavid du Colombier enum
35*ec59a3ddSDavid du Colombier {
36*ec59a3ddSDavid du Colombier 	S11=	7,
37*ec59a3ddSDavid du Colombier 	S12=	12,
38*ec59a3ddSDavid du Colombier 	S13=	17,
39*ec59a3ddSDavid du Colombier 	S14=	22,
40*ec59a3ddSDavid du Colombier 
41*ec59a3ddSDavid du Colombier 	S21=	5,
42*ec59a3ddSDavid du Colombier 	S22=	9,
43*ec59a3ddSDavid du Colombier 	S23=	14,
44*ec59a3ddSDavid du Colombier 	S24=	20,
45*ec59a3ddSDavid du Colombier 
46*ec59a3ddSDavid du Colombier 	S31=	4,
47*ec59a3ddSDavid du Colombier 	S32=	11,
48*ec59a3ddSDavid du Colombier 	S33=	16,
49*ec59a3ddSDavid du Colombier 	S34=	23,
50*ec59a3ddSDavid du Colombier 
51*ec59a3ddSDavid du Colombier 	S41=	6,
52*ec59a3ddSDavid du Colombier 	S42=	10,
53*ec59a3ddSDavid du Colombier 	S43=	15,
54*ec59a3ddSDavid du Colombier 	S44=	21,
55*ec59a3ddSDavid du Colombier };
56*ec59a3ddSDavid du Colombier 
57*ec59a3ddSDavid du Colombier static u32int md5tab[] =
58*ec59a3ddSDavid du Colombier {
59*ec59a3ddSDavid du Colombier 	/* round 1 */
60*ec59a3ddSDavid du Colombier /*[0]*/	0xd76aa478,
61*ec59a3ddSDavid du Colombier 	0xe8c7b756,
62*ec59a3ddSDavid du Colombier 	0x242070db,
63*ec59a3ddSDavid du Colombier 	0xc1bdceee,
64*ec59a3ddSDavid du Colombier 	0xf57c0faf,
65*ec59a3ddSDavid du Colombier 	0x4787c62a,
66*ec59a3ddSDavid du Colombier 	0xa8304613,
67*ec59a3ddSDavid du Colombier 	0xfd469501,
68*ec59a3ddSDavid du Colombier 	0x698098d8,
69*ec59a3ddSDavid du Colombier 	0x8b44f7af,
70*ec59a3ddSDavid du Colombier 	0xffff5bb1,
71*ec59a3ddSDavid du Colombier 	0x895cd7be,
72*ec59a3ddSDavid du Colombier 	0x6b901122,
73*ec59a3ddSDavid du Colombier 	0xfd987193,
74*ec59a3ddSDavid du Colombier 	0xa679438e,
75*ec59a3ddSDavid du Colombier 	0x49b40821,
76*ec59a3ddSDavid du Colombier 
77*ec59a3ddSDavid du Colombier 	/* round 2 */
78*ec59a3ddSDavid du Colombier /*[16]*/0xf61e2562,
79*ec59a3ddSDavid du Colombier 	0xc040b340,
80*ec59a3ddSDavid du Colombier 	0x265e5a51,
81*ec59a3ddSDavid du Colombier 	0xe9b6c7aa,
82*ec59a3ddSDavid du Colombier 	0xd62f105d,
83*ec59a3ddSDavid du Colombier 	 0x2441453,
84*ec59a3ddSDavid du Colombier 	0xd8a1e681,
85*ec59a3ddSDavid du Colombier 	0xe7d3fbc8,
86*ec59a3ddSDavid du Colombier 	0x21e1cde6,
87*ec59a3ddSDavid du Colombier 	0xc33707d6,
88*ec59a3ddSDavid du Colombier 	0xf4d50d87,
89*ec59a3ddSDavid du Colombier 	0x455a14ed,
90*ec59a3ddSDavid du Colombier 	0xa9e3e905,
91*ec59a3ddSDavid du Colombier 	0xfcefa3f8,
92*ec59a3ddSDavid du Colombier 	0x676f02d9,
93*ec59a3ddSDavid du Colombier 	0x8d2a4c8a,
94*ec59a3ddSDavid du Colombier 
95*ec59a3ddSDavid du Colombier 	/* round 3 */
96*ec59a3ddSDavid du Colombier /*[32]*/0xfffa3942,
97*ec59a3ddSDavid du Colombier 	0x8771f681,
98*ec59a3ddSDavid du Colombier 	0x6d9d6122,
99*ec59a3ddSDavid du Colombier 	0xfde5380c,
100*ec59a3ddSDavid du Colombier 	0xa4beea44,
101*ec59a3ddSDavid du Colombier 	0x4bdecfa9,
102*ec59a3ddSDavid du Colombier 	0xf6bb4b60,
103*ec59a3ddSDavid du Colombier 	0xbebfbc70,
104*ec59a3ddSDavid du Colombier 	0x289b7ec6,
105*ec59a3ddSDavid du Colombier 	0xeaa127fa,
106*ec59a3ddSDavid du Colombier 	0xd4ef3085,
107*ec59a3ddSDavid du Colombier 	 0x4881d05,
108*ec59a3ddSDavid du Colombier 	0xd9d4d039,
109*ec59a3ddSDavid du Colombier 	0xe6db99e5,
110*ec59a3ddSDavid du Colombier 	0x1fa27cf8,
111*ec59a3ddSDavid du Colombier 	0xc4ac5665,
112*ec59a3ddSDavid du Colombier 
113*ec59a3ddSDavid du Colombier 	/* round 4 */
114*ec59a3ddSDavid du Colombier /*[48]*/0xf4292244,
115*ec59a3ddSDavid du Colombier 	0x432aff97,
116*ec59a3ddSDavid du Colombier 	0xab9423a7,
117*ec59a3ddSDavid du Colombier 	0xfc93a039,
118*ec59a3ddSDavid du Colombier 	0x655b59c3,
119*ec59a3ddSDavid du Colombier 	0x8f0ccc92,
120*ec59a3ddSDavid du Colombier 	0xffeff47d,
121*ec59a3ddSDavid du Colombier 	0x85845dd1,
122*ec59a3ddSDavid du Colombier 	0x6fa87e4f,
123*ec59a3ddSDavid du Colombier 	0xfe2ce6e0,
124*ec59a3ddSDavid du Colombier 	0xa3014314,
125*ec59a3ddSDavid du Colombier 	0x4e0811a1,
126*ec59a3ddSDavid du Colombier 	0xf7537e82,
127*ec59a3ddSDavid du Colombier 	0xbd3af235,
128*ec59a3ddSDavid du Colombier 	0x2ad7d2bb,
129*ec59a3ddSDavid du Colombier 	0xeb86d391,
130*ec59a3ddSDavid du Colombier };
131*ec59a3ddSDavid du Colombier 
132*ec59a3ddSDavid du Colombier static void decode(u32int*, uchar*, ulong);
133*ec59a3ddSDavid du Colombier extern void _md5block(uchar *p, ulong len, u32int *s);
134*ec59a3ddSDavid du Colombier 
135*ec59a3ddSDavid du Colombier void
_md5block(uchar * p,ulong len,u32int * s)136*ec59a3ddSDavid du Colombier _md5block(uchar *p, ulong len, u32int *s)
137*ec59a3ddSDavid du Colombier {
138*ec59a3ddSDavid du Colombier 	u32int a, b, c, d, sh;
139*ec59a3ddSDavid du Colombier 	u32int *t;
140*ec59a3ddSDavid du Colombier 	uchar *end;
141*ec59a3ddSDavid du Colombier 	u32int x[16];
142*ec59a3ddSDavid du Colombier 
143*ec59a3ddSDavid du Colombier 	for(end = p+len; p < end; p += 64){
144*ec59a3ddSDavid du Colombier 		a = s[0];
145*ec59a3ddSDavid du Colombier 		b = s[1];
146*ec59a3ddSDavid du Colombier 		c = s[2];
147*ec59a3ddSDavid du Colombier 		d = s[3];
148*ec59a3ddSDavid du Colombier 
149*ec59a3ddSDavid du Colombier 		decode(x, p, 64);
150*ec59a3ddSDavid du Colombier 
151*ec59a3ddSDavid du Colombier 		t = md5tab;
152*ec59a3ddSDavid du Colombier 		sh = 0;
153*ec59a3ddSDavid du Colombier 		for(; sh != 16; t += 4){
154*ec59a3ddSDavid du Colombier 			a += ((c ^ d) & b) ^ d;
155*ec59a3ddSDavid du Colombier 			a += x[sh] + t[0];
156*ec59a3ddSDavid du Colombier 			a = (a << S11) | (a >> (32 - S11));
157*ec59a3ddSDavid du Colombier 			a += b;
158*ec59a3ddSDavid du Colombier 
159*ec59a3ddSDavid du Colombier 			d += ((b ^ c) & a) ^ c;
160*ec59a3ddSDavid du Colombier 			d += x[sh + 1] + t[1];
161*ec59a3ddSDavid du Colombier 			d = (d << S12) | (d >> (32 - S12));
162*ec59a3ddSDavid du Colombier 			d += a;
163*ec59a3ddSDavid du Colombier 
164*ec59a3ddSDavid du Colombier 			c += ((a ^ b) & d) ^ b;
165*ec59a3ddSDavid du Colombier 			c += x[sh + 2] + t[2];
166*ec59a3ddSDavid du Colombier 			c = (c << S13) | (c >> (32 - S13));
167*ec59a3ddSDavid du Colombier 			c += d;
168*ec59a3ddSDavid du Colombier 
169*ec59a3ddSDavid du Colombier 			b += ((d ^ a) & c) ^ a;
170*ec59a3ddSDavid du Colombier 			b += x[sh + 3] + t[3];
171*ec59a3ddSDavid du Colombier 			b = (b << S14) | (b >> (32 - S14));
172*ec59a3ddSDavid du Colombier 			b += c;
173*ec59a3ddSDavid du Colombier 
174*ec59a3ddSDavid du Colombier 			sh += 4;
175*ec59a3ddSDavid du Colombier 		}
176*ec59a3ddSDavid du Colombier 		sh = 1;
177*ec59a3ddSDavid du Colombier 		for(; sh != 1+20*4; t += 4){
178*ec59a3ddSDavid du Colombier 			a += ((b ^ c) & d) ^ c;
179*ec59a3ddSDavid du Colombier 			a += x[sh & 0xf] + t[0];
180*ec59a3ddSDavid du Colombier 			a = (a << S21) | (a >> (32 - S21));
181*ec59a3ddSDavid du Colombier 			a += b;
182*ec59a3ddSDavid du Colombier 
183*ec59a3ddSDavid du Colombier 			d += ((a ^ b) & c) ^ b;
184*ec59a3ddSDavid du Colombier 			d += x[(sh + 5) & 0xf] + t[1];
185*ec59a3ddSDavid du Colombier 			d = (d << S22) | (d >> (32 - S22));
186*ec59a3ddSDavid du Colombier 			d += a;
187*ec59a3ddSDavid du Colombier 
188*ec59a3ddSDavid du Colombier 			c += ((d ^ a) & b) ^ a;
189*ec59a3ddSDavid du Colombier 			c += x[(sh + 10) & 0xf] + t[2];
190*ec59a3ddSDavid du Colombier 			c = (c << S23) | (c >> (32 - S23));
191*ec59a3ddSDavid du Colombier 			c += d;
192*ec59a3ddSDavid du Colombier 
193*ec59a3ddSDavid du Colombier 			b += ((c ^ d) & a) ^ d;
194*ec59a3ddSDavid du Colombier 			b += x[(sh + 15) & 0xf] + t[3];
195*ec59a3ddSDavid du Colombier 			b = (b << S24) | (b >> (32 - S24));
196*ec59a3ddSDavid du Colombier 			b += c;
197*ec59a3ddSDavid du Colombier 
198*ec59a3ddSDavid du Colombier 			sh += 20;
199*ec59a3ddSDavid du Colombier 		}
200*ec59a3ddSDavid du Colombier 		sh = 5;
201*ec59a3ddSDavid du Colombier 		for(; sh != 5+12*4; t += 4){
202*ec59a3ddSDavid du Colombier 			a += b ^ c ^ d;
203*ec59a3ddSDavid du Colombier 			a += x[sh & 0xf] + t[0];
204*ec59a3ddSDavid du Colombier 			a = (a << S31) | (a >> (32 - S31));
205*ec59a3ddSDavid du Colombier 			a += b;
206*ec59a3ddSDavid du Colombier 
207*ec59a3ddSDavid du Colombier 			d += a ^ b ^ c;
208*ec59a3ddSDavid du Colombier 			d += x[(sh + 3) & 0xf] + t[1];
209*ec59a3ddSDavid du Colombier 			d = (d << S32) | (d >> (32 - S32));
210*ec59a3ddSDavid du Colombier 			d += a;
211*ec59a3ddSDavid du Colombier 
212*ec59a3ddSDavid du Colombier 			c += d ^ a ^ b;
213*ec59a3ddSDavid du Colombier 			c += x[(sh + 6) & 0xf] + t[2];
214*ec59a3ddSDavid du Colombier 			c = (c << S33) | (c >> (32 - S33));
215*ec59a3ddSDavid du Colombier 			c += d;
216*ec59a3ddSDavid du Colombier 
217*ec59a3ddSDavid du Colombier 			b += c ^ d ^ a;
218*ec59a3ddSDavid du Colombier 			b += x[(sh + 9) & 0xf] + t[3];
219*ec59a3ddSDavid du Colombier 			b = (b << S34) | (b >> (32 - S34));
220*ec59a3ddSDavid du Colombier 			b += c;
221*ec59a3ddSDavid du Colombier 
222*ec59a3ddSDavid du Colombier 			sh += 12;
223*ec59a3ddSDavid du Colombier 		}
224*ec59a3ddSDavid du Colombier 		sh = 0;
225*ec59a3ddSDavid du Colombier 		for(; sh != 28*4; t += 4){
226*ec59a3ddSDavid du Colombier 			a += c ^ (b | ~d);
227*ec59a3ddSDavid du Colombier 			a += x[sh & 0xf] + t[0];
228*ec59a3ddSDavid du Colombier 			a = (a << S41) | (a >> (32 - S41));
229*ec59a3ddSDavid du Colombier 			a += b;
230*ec59a3ddSDavid du Colombier 
231*ec59a3ddSDavid du Colombier 			d += b ^ (a | ~c);
232*ec59a3ddSDavid du Colombier 			d += x[(sh + 7) & 0xf] + t[1];
233*ec59a3ddSDavid du Colombier 			d = (d << S42) | (d >> (32 - S42));
234*ec59a3ddSDavid du Colombier 			d += a;
235*ec59a3ddSDavid du Colombier 
236*ec59a3ddSDavid du Colombier 			c += a ^ (d | ~b);
237*ec59a3ddSDavid du Colombier 			c += x[(sh + 14) & 0xf] + t[2];
238*ec59a3ddSDavid du Colombier 			c = (c << S43) | (c >> (32 - S43));
239*ec59a3ddSDavid du Colombier 			c += d;
240*ec59a3ddSDavid du Colombier 
241*ec59a3ddSDavid du Colombier 			b += d ^ (c | ~a);
242*ec59a3ddSDavid du Colombier 			b += x[(sh + 21) & 0xf] + t[3];
243*ec59a3ddSDavid du Colombier 			b = (b << S44) | (b >> (32 - S44));
244*ec59a3ddSDavid du Colombier 			b += c;
245*ec59a3ddSDavid du Colombier 
246*ec59a3ddSDavid du Colombier 			sh += 28;
247*ec59a3ddSDavid du Colombier 		}
248*ec59a3ddSDavid du Colombier 
249*ec59a3ddSDavid du Colombier 		s[0] += a;
250*ec59a3ddSDavid du Colombier 		s[1] += b;
251*ec59a3ddSDavid du Colombier 		s[2] += c;
252*ec59a3ddSDavid du Colombier 		s[3] += d;
253*ec59a3ddSDavid du Colombier 	}
254*ec59a3ddSDavid du Colombier }
255*ec59a3ddSDavid du Colombier 
256*ec59a3ddSDavid du Colombier /*
257*ec59a3ddSDavid du Colombier  *	decodes input (uchar) into output (u32int). Assumes len is
258*ec59a3ddSDavid du Colombier  *	a multiple of 4.
259*ec59a3ddSDavid du Colombier  */
260*ec59a3ddSDavid du Colombier static void
decode(u32int * output,uchar * input,ulong len)261*ec59a3ddSDavid du Colombier decode(u32int *output, uchar *input, ulong len)
262*ec59a3ddSDavid du Colombier {
263*ec59a3ddSDavid du Colombier 	uchar *e;
264*ec59a3ddSDavid du Colombier 
265*ec59a3ddSDavid du Colombier 	for(e = input+len; input < e; input += 4)
266*ec59a3ddSDavid du Colombier 		*output++ = input[0] | (input[1] << 8) |
267*ec59a3ddSDavid du Colombier 			(input[2] << 16) | (input[3] << 24);
268*ec59a3ddSDavid du Colombier }
269