rmd160c.c revision 314184
1185029Spjd/* crypto/ripemd/rmd_dgst.c */
2185029Spjd/* Copyright (C) 1995-1998 Eric Young (eay@cryptsoft.com)
3185029Spjd * All rights reserved.
4185029Spjd *
5185029Spjd * This package is an SSL implementation written
6185029Spjd * by Eric Young (eay@cryptsoft.com).
7185029Spjd * The implementation was written so as to conform with Netscapes SSL.
8185029Spjd *
9185029Spjd * This library is free for commercial and non-commercial use as long as
10185029Spjd * the following conditions are aheared to.  The following conditions
11185029Spjd * apply to all code found in this distribution, be it the RC4, RSA,
12185029Spjd * lhash, DES, etc., code; not just the SSL code.  The SSL documentation
13185029Spjd * included with this distribution is covered by the same copyright terms
14185029Spjd * except that the holder is Tim Hudson (tjh@cryptsoft.com).
15185029Spjd *
16185029Spjd * Copyright remains Eric Young's, and as such any Copyright notices in
17185029Spjd * the code are not to be removed.
18185029Spjd * If this package is used in a product, Eric Young should be given attribution
19185029Spjd * as the author of the parts of the library used.
20185029Spjd * This can be in the form of a textual message at program startup or
21185029Spjd * in documentation (online or textual) provided with the package.
22219089Spjd *
23219089Spjd * Redistribution and use in source and binary forms, with or without
24263407Sdelphij * modification, are permitted provided that the following conditions
25185029Spjd * are met:
26185029Spjd * 1. Redistributions of source code must retain the copyright
27263407Sdelphij *    notice, this list of conditions and the following disclaimer.
28263407Sdelphij * 2. Redistributions in binary form must reproduce the above copyright
29185029Spjd *    notice, this list of conditions and the following disclaimer in the
30185029Spjd *    documentation and/or other materials provided with the distribution.
31185029Spjd * 3. All advertising materials mentioning features or use of this software
32185029Spjd *    must display the following acknowledgement:
33185029Spjd *    "This product includes cryptographic software written by
34185029Spjd *     Eric Young (eay@cryptsoft.com)"
35185029Spjd *    The word 'cryptographic' can be left out if the rouines from the library
36185029Spjd *    being used are not cryptographic related :-).
37185029Spjd * 4. If you include any Windows specific code (or a derivative thereof) from
38185029Spjd *    the apps directory (application code) you must include an acknowledgement:
39185029Spjd *    "This product includes software written by Tim Hudson (tjh@cryptsoft.com)"
40185029Spjd *
41185029Spjd * THIS SOFTWARE IS PROVIDED BY ERIC YOUNG ``AS IS'' AND
42185029Spjd * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
43185029Spjd * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
44185029Spjd * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
45185029Spjd * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
46263407Sdelphij * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
47263407Sdelphij * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
48263407Sdelphij * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
49263407Sdelphij * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
50263407Sdelphij * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
51263407Sdelphij * SUCH DAMAGE.
52263407Sdelphij *
53263407Sdelphij * The licence and distribution terms for any publically available version or
54263407Sdelphij * derivative of this code cannot be changed.  i.e. this code cannot simply be
55263407Sdelphij * copied and put under another distribution licence
56263407Sdelphij * [including the GNU Public Licence.]
57263407Sdelphij */
58263407Sdelphij
59263407Sdelphij#include <sys/cdefs.h>
60263407Sdelphij__FBSDID("$FreeBSD: stable/10/lib/libmd/rmd160c.c 314184 2017-02-23 22:10:37Z avg $");
61263407Sdelphij
62263407Sdelphij#include <sys/types.h>
63263407Sdelphij
64263407Sdelphij#include <stdio.h>
65263407Sdelphij#include <string.h>
66263407Sdelphij
67263407Sdelphij#if 0
68185029Spjd#include <machine/ansi.h>	/* we use the __ variants of bit-sized types */
69185029Spjd#endif
70185029Spjd#include <machine/endian.h>
71185029Spjd
72185029Spjd#include "rmd_locl.h"
73185029Spjd
74185029Spjd/*
75185029Spjd * The assembly-language code is not position-independent, so don't
76185029Spjd * try to use it in a shared library.
77185029Spjd */
78185029Spjd#ifdef PIC
79185029Spjd#undef RMD160_ASM
80185029Spjd#endif
81185029Spjd
82185029Spjdchar *RMD160_version="RIPEMD160 part of SSLeay 0.9.0b 11-Oct-1998";
83185029Spjd
84185029Spjd#ifdef RMD160_ASM
85185029Spjdvoid ripemd160_block_x86(RIPEMD160_CTX *c, const u_int32_t *p,int num);
86185029Spjd#define ripemd160_block ripemd160_block_x86
87185029Spjd#else
88185029Spjdvoid ripemd160_block(RIPEMD160_CTX *c, const u_int32_t *p,int num);
89185029Spjd#endif
90185029Spjd
91185029Spjdvoid RIPEMD160_Init(c)
92185029SpjdRIPEMD160_CTX *c;
93185029Spjd	{
94185029Spjd	c->A=RIPEMD160_A;
95185029Spjd	c->B=RIPEMD160_B;
96185029Spjd	c->C=RIPEMD160_C;
97185029Spjd	c->D=RIPEMD160_D;
98185029Spjd	c->E=RIPEMD160_E;
99185029Spjd	c->Nl=0;
100185029Spjd	c->Nh=0;
101185029Spjd	c->num=0;
102185029Spjd	}
103185029Spjd
104185029Spjdvoid RIPEMD160_Update(c, in, len)
105185029SpjdRIPEMD160_CTX *c;
106185029Spjdconst void *in;
107185029Spjdsize_t len;
108185029Spjd	{
109185029Spjd	u_int32_t *p;
110185029Spjd	int sw,sc;
111185029Spjd	u_int32_t l;
112185029Spjd	const unsigned char *data = in;
113185029Spjd
114185029Spjd	if (len == 0) return;
115185029Spjd
116185029Spjd	l=(c->Nl+(len<<3))&0xffffffffL;
117185029Spjd	if (l < c->Nl) /* overflow */
118185029Spjd		c->Nh++;
119185029Spjd	c->Nh+=(len>>29);
120185029Spjd	c->Nl=l;
121185029Spjd
122185029Spjd	if (c->num != 0)
123185029Spjd		{
124185029Spjd		p=c->data;
125185029Spjd		sw=c->num>>2;
126185029Spjd		sc=c->num&0x03;
127185029Spjd
128185029Spjd		if ((c->num+len) >= RIPEMD160_CBLOCK)
129185029Spjd			{
130185029Spjd			l= p[sw];
131185029Spjd			p_c2l(data,l,sc);
132185029Spjd			p[sw++]=l;
133185029Spjd			for (; sw<RIPEMD160_LBLOCK; sw++)
134185029Spjd				{
135185029Spjd				c2l(data,l);
136185029Spjd				p[sw]=l;
137185029Spjd				}
138185029Spjd			len-=(RIPEMD160_CBLOCK-c->num);
139185029Spjd
140185029Spjd			ripemd160_block(c,p,64);
141185029Spjd			c->num=0;
142185029Spjd			/* drop through and do the rest */
143185029Spjd			}
144185029Spjd		else
145185029Spjd			{
146185029Spjd			int ew,ec;
147185029Spjd
148185029Spjd			c->num+=(int)len;
149185029Spjd			if ((sc+len) < 4) /* ugly, add char's to a word */
150185029Spjd				{
151185029Spjd				l= p[sw];
152185029Spjd				p_c2l_p(data,l,sc,len);
153185029Spjd				p[sw]=l;
154185029Spjd				}
155185029Spjd			else
156185029Spjd				{
157185029Spjd				ew=(c->num>>2);
158185029Spjd				ec=(c->num&0x03);
159185029Spjd				l= p[sw];
160185029Spjd				p_c2l(data,l,sc);
161185029Spjd				p[sw++]=l;
162185029Spjd				for (; sw < ew; sw++)
163185029Spjd					{ c2l(data,l); p[sw]=l; }
164185029Spjd				if (ec)
165185029Spjd					{
166185029Spjd					c2l_p(data,l,ec);
167185029Spjd					p[sw]=l;
168185029Spjd					}
169185029Spjd				}
170185029Spjd			return;
171185029Spjd			}
172185029Spjd		}
173185029Spjd	/* we now can process the input data in blocks of RIPEMD160_CBLOCK
174185029Spjd	 * chars and save the leftovers to c->data. */
175185029Spjd#if BYTE_ORDER == LITTLE_ENDIAN
176185029Spjd	if ((((unsigned long)data)%sizeof(u_int32_t)) == 0)
177185029Spjd		{
178185029Spjd		sw=(int)len/RIPEMD160_CBLOCK;
179185029Spjd		if (sw > 0)
180185029Spjd			{
181185029Spjd			sw*=RIPEMD160_CBLOCK;
182185029Spjd			ripemd160_block(c,(u_int32_t *)data,sw);
183185029Spjd			data+=sw;
184185029Spjd			len-=sw;
185185029Spjd			}
186185029Spjd		}
187185029Spjd#endif
188185029Spjd	p=c->data;
189185029Spjd	while (len >= RIPEMD160_CBLOCK)
190185029Spjd		{
191185029Spjd#if BYTE_ORDER == LITTLE_ENDIAN || BYTE_ORDER == BIG_ENDIAN
192185029Spjd		if (p != (u_int32_t *)data)
193185029Spjd			memcpy(p,data,RIPEMD160_CBLOCK);
194185029Spjd		data+=RIPEMD160_CBLOCK;
195185029Spjd#if BYTE_ORDER == BIG_ENDIAN
196185029Spjd		for (sw=(RIPEMD160_LBLOCK/4); sw; sw--)
197185029Spjd			{
198185029Spjd			Endian_Reverse32(p[0]);
199185029Spjd			Endian_Reverse32(p[1]);
200185029Spjd			Endian_Reverse32(p[2]);
201185029Spjd			Endian_Reverse32(p[3]);
202185029Spjd			p+=4;
203185029Spjd			}
204185029Spjd#endif
205185029Spjd#else
206185029Spjd		for (sw=(RIPEMD160_LBLOCK/4); sw; sw--)
207185029Spjd			{
208185029Spjd			c2l(data,l); *(p++)=l;
209185029Spjd			c2l(data,l); *(p++)=l;
210185029Spjd			c2l(data,l); *(p++)=l;
211185029Spjd			c2l(data,l); *(p++)=l;
212185029Spjd			}
213185029Spjd#endif
214185029Spjd		p=c->data;
215185029Spjd		ripemd160_block(c,p,64);
216185029Spjd		len-=RIPEMD160_CBLOCK;
217185029Spjd		}
218185029Spjd	sc=(int)len;
219185029Spjd	c->num=sc;
220185029Spjd	if (sc)
221185029Spjd		{
222185029Spjd		sw=sc>>2;	/* words to copy */
223185029Spjd#if BYTE_ORDER == LITTLE_ENDIAN
224185029Spjd		p[sw]=0;
225185029Spjd		memcpy(p,data,sc);
226185029Spjd#else
227185029Spjd		sc&=0x03;
228185029Spjd		for ( ; sw; sw--)
229185029Spjd			{ c2l(data,l); *(p++)=l; }
230185029Spjd		c2l_p(data,l,sc);
231185029Spjd		*p=l;
232185029Spjd#endif
233		}
234	}
235
236void RIPEMD160_Transform(c,b)
237RIPEMD160_CTX *c;
238unsigned char *b;
239	{
240	u_int32_t p[16];
241#if BYTE_ORDER != LITTLE_ENDIAN
242	u_int32_t *q;
243	int i;
244#endif
245
246#if BYTE_ORDER == BIG_ENDIAN || BYTE_ORDER == LITTLE_ENDIAN
247	memcpy(p,b,64);
248#if BYTE_ORDER == BIG_ENDIAN
249	q=p;
250	for (i=(RIPEMD160_LBLOCK/4); i; i--)
251		{
252		Endian_Reverse32(q[0]);
253		Endian_Reverse32(q[1]);
254		Endian_Reverse32(q[2]);
255		Endian_Reverse32(q[3]);
256		q+=4;
257		}
258#endif
259#else
260	q=p;
261	for (i=(RIPEMD160_LBLOCK/4); i; i--)
262		{
263		u_int32_t l;
264		c2l(b,l); *(q++)=l;
265		c2l(b,l); *(q++)=l;
266		c2l(b,l); *(q++)=l;
267		c2l(b,l); *(q++)=l;
268		}
269#endif
270	ripemd160_block(c,p,64);
271	}
272
273#ifndef RMD160_ASM
274
275void ripemd160_block(ctx, X, num)
276RIPEMD160_CTX *ctx;
277const u_int32_t *X;
278int num;
279	{
280	u_int32_t A,B,C,D,E;
281	u_int32_t a,b,c,d,e;
282
283	for (;;)
284		{
285		A=ctx->A; B=ctx->B; C=ctx->C; D=ctx->D; E=ctx->E;
286
287	RIP1(A,B,C,D,E,WL00,SL00);
288	RIP1(E,A,B,C,D,WL01,SL01);
289	RIP1(D,E,A,B,C,WL02,SL02);
290	RIP1(C,D,E,A,B,WL03,SL03);
291	RIP1(B,C,D,E,A,WL04,SL04);
292	RIP1(A,B,C,D,E,WL05,SL05);
293	RIP1(E,A,B,C,D,WL06,SL06);
294	RIP1(D,E,A,B,C,WL07,SL07);
295	RIP1(C,D,E,A,B,WL08,SL08);
296	RIP1(B,C,D,E,A,WL09,SL09);
297	RIP1(A,B,C,D,E,WL10,SL10);
298	RIP1(E,A,B,C,D,WL11,SL11);
299	RIP1(D,E,A,B,C,WL12,SL12);
300	RIP1(C,D,E,A,B,WL13,SL13);
301	RIP1(B,C,D,E,A,WL14,SL14);
302	RIP1(A,B,C,D,E,WL15,SL15);
303
304	RIP2(E,A,B,C,D,WL16,SL16,KL1);
305	RIP2(D,E,A,B,C,WL17,SL17,KL1);
306	RIP2(C,D,E,A,B,WL18,SL18,KL1);
307	RIP2(B,C,D,E,A,WL19,SL19,KL1);
308	RIP2(A,B,C,D,E,WL20,SL20,KL1);
309	RIP2(E,A,B,C,D,WL21,SL21,KL1);
310	RIP2(D,E,A,B,C,WL22,SL22,KL1);
311	RIP2(C,D,E,A,B,WL23,SL23,KL1);
312	RIP2(B,C,D,E,A,WL24,SL24,KL1);
313	RIP2(A,B,C,D,E,WL25,SL25,KL1);
314	RIP2(E,A,B,C,D,WL26,SL26,KL1);
315	RIP2(D,E,A,B,C,WL27,SL27,KL1);
316	RIP2(C,D,E,A,B,WL28,SL28,KL1);
317	RIP2(B,C,D,E,A,WL29,SL29,KL1);
318	RIP2(A,B,C,D,E,WL30,SL30,KL1);
319	RIP2(E,A,B,C,D,WL31,SL31,KL1);
320
321	RIP3(D,E,A,B,C,WL32,SL32,KL2);
322	RIP3(C,D,E,A,B,WL33,SL33,KL2);
323	RIP3(B,C,D,E,A,WL34,SL34,KL2);
324	RIP3(A,B,C,D,E,WL35,SL35,KL2);
325	RIP3(E,A,B,C,D,WL36,SL36,KL2);
326	RIP3(D,E,A,B,C,WL37,SL37,KL2);
327	RIP3(C,D,E,A,B,WL38,SL38,KL2);
328	RIP3(B,C,D,E,A,WL39,SL39,KL2);
329	RIP3(A,B,C,D,E,WL40,SL40,KL2);
330	RIP3(E,A,B,C,D,WL41,SL41,KL2);
331	RIP3(D,E,A,B,C,WL42,SL42,KL2);
332	RIP3(C,D,E,A,B,WL43,SL43,KL2);
333	RIP3(B,C,D,E,A,WL44,SL44,KL2);
334	RIP3(A,B,C,D,E,WL45,SL45,KL2);
335	RIP3(E,A,B,C,D,WL46,SL46,KL2);
336	RIP3(D,E,A,B,C,WL47,SL47,KL2);
337
338	RIP4(C,D,E,A,B,WL48,SL48,KL3);
339	RIP4(B,C,D,E,A,WL49,SL49,KL3);
340	RIP4(A,B,C,D,E,WL50,SL50,KL3);
341	RIP4(E,A,B,C,D,WL51,SL51,KL3);
342	RIP4(D,E,A,B,C,WL52,SL52,KL3);
343	RIP4(C,D,E,A,B,WL53,SL53,KL3);
344	RIP4(B,C,D,E,A,WL54,SL54,KL3);
345	RIP4(A,B,C,D,E,WL55,SL55,KL3);
346	RIP4(E,A,B,C,D,WL56,SL56,KL3);
347	RIP4(D,E,A,B,C,WL57,SL57,KL3);
348	RIP4(C,D,E,A,B,WL58,SL58,KL3);
349	RIP4(B,C,D,E,A,WL59,SL59,KL3);
350	RIP4(A,B,C,D,E,WL60,SL60,KL3);
351	RIP4(E,A,B,C,D,WL61,SL61,KL3);
352	RIP4(D,E,A,B,C,WL62,SL62,KL3);
353	RIP4(C,D,E,A,B,WL63,SL63,KL3);
354
355	RIP5(B,C,D,E,A,WL64,SL64,KL4);
356	RIP5(A,B,C,D,E,WL65,SL65,KL4);
357	RIP5(E,A,B,C,D,WL66,SL66,KL4);
358	RIP5(D,E,A,B,C,WL67,SL67,KL4);
359	RIP5(C,D,E,A,B,WL68,SL68,KL4);
360	RIP5(B,C,D,E,A,WL69,SL69,KL4);
361	RIP5(A,B,C,D,E,WL70,SL70,KL4);
362	RIP5(E,A,B,C,D,WL71,SL71,KL4);
363	RIP5(D,E,A,B,C,WL72,SL72,KL4);
364	RIP5(C,D,E,A,B,WL73,SL73,KL4);
365	RIP5(B,C,D,E,A,WL74,SL74,KL4);
366	RIP5(A,B,C,D,E,WL75,SL75,KL4);
367	RIP5(E,A,B,C,D,WL76,SL76,KL4);
368	RIP5(D,E,A,B,C,WL77,SL77,KL4);
369	RIP5(C,D,E,A,B,WL78,SL78,KL4);
370	RIP5(B,C,D,E,A,WL79,SL79,KL4);
371
372	a=A; b=B; c=C; d=D; e=E;
373	/* Do other half */
374	A=ctx->A; B=ctx->B; C=ctx->C; D=ctx->D; E=ctx->E;
375
376	RIP5(A,B,C,D,E,WR00,SR00,KR0);
377	RIP5(E,A,B,C,D,WR01,SR01,KR0);
378	RIP5(D,E,A,B,C,WR02,SR02,KR0);
379	RIP5(C,D,E,A,B,WR03,SR03,KR0);
380	RIP5(B,C,D,E,A,WR04,SR04,KR0);
381	RIP5(A,B,C,D,E,WR05,SR05,KR0);
382	RIP5(E,A,B,C,D,WR06,SR06,KR0);
383	RIP5(D,E,A,B,C,WR07,SR07,KR0);
384	RIP5(C,D,E,A,B,WR08,SR08,KR0);
385	RIP5(B,C,D,E,A,WR09,SR09,KR0);
386	RIP5(A,B,C,D,E,WR10,SR10,KR0);
387	RIP5(E,A,B,C,D,WR11,SR11,KR0);
388	RIP5(D,E,A,B,C,WR12,SR12,KR0);
389	RIP5(C,D,E,A,B,WR13,SR13,KR0);
390	RIP5(B,C,D,E,A,WR14,SR14,KR0);
391	RIP5(A,B,C,D,E,WR15,SR15,KR0);
392
393	RIP4(E,A,B,C,D,WR16,SR16,KR1);
394	RIP4(D,E,A,B,C,WR17,SR17,KR1);
395	RIP4(C,D,E,A,B,WR18,SR18,KR1);
396	RIP4(B,C,D,E,A,WR19,SR19,KR1);
397	RIP4(A,B,C,D,E,WR20,SR20,KR1);
398	RIP4(E,A,B,C,D,WR21,SR21,KR1);
399	RIP4(D,E,A,B,C,WR22,SR22,KR1);
400	RIP4(C,D,E,A,B,WR23,SR23,KR1);
401	RIP4(B,C,D,E,A,WR24,SR24,KR1);
402	RIP4(A,B,C,D,E,WR25,SR25,KR1);
403	RIP4(E,A,B,C,D,WR26,SR26,KR1);
404	RIP4(D,E,A,B,C,WR27,SR27,KR1);
405	RIP4(C,D,E,A,B,WR28,SR28,KR1);
406	RIP4(B,C,D,E,A,WR29,SR29,KR1);
407	RIP4(A,B,C,D,E,WR30,SR30,KR1);
408	RIP4(E,A,B,C,D,WR31,SR31,KR1);
409
410	RIP3(D,E,A,B,C,WR32,SR32,KR2);
411	RIP3(C,D,E,A,B,WR33,SR33,KR2);
412	RIP3(B,C,D,E,A,WR34,SR34,KR2);
413	RIP3(A,B,C,D,E,WR35,SR35,KR2);
414	RIP3(E,A,B,C,D,WR36,SR36,KR2);
415	RIP3(D,E,A,B,C,WR37,SR37,KR2);
416	RIP3(C,D,E,A,B,WR38,SR38,KR2);
417	RIP3(B,C,D,E,A,WR39,SR39,KR2);
418	RIP3(A,B,C,D,E,WR40,SR40,KR2);
419	RIP3(E,A,B,C,D,WR41,SR41,KR2);
420	RIP3(D,E,A,B,C,WR42,SR42,KR2);
421	RIP3(C,D,E,A,B,WR43,SR43,KR2);
422	RIP3(B,C,D,E,A,WR44,SR44,KR2);
423	RIP3(A,B,C,D,E,WR45,SR45,KR2);
424	RIP3(E,A,B,C,D,WR46,SR46,KR2);
425	RIP3(D,E,A,B,C,WR47,SR47,KR2);
426
427	RIP2(C,D,E,A,B,WR48,SR48,KR3);
428	RIP2(B,C,D,E,A,WR49,SR49,KR3);
429	RIP2(A,B,C,D,E,WR50,SR50,KR3);
430	RIP2(E,A,B,C,D,WR51,SR51,KR3);
431	RIP2(D,E,A,B,C,WR52,SR52,KR3);
432	RIP2(C,D,E,A,B,WR53,SR53,KR3);
433	RIP2(B,C,D,E,A,WR54,SR54,KR3);
434	RIP2(A,B,C,D,E,WR55,SR55,KR3);
435	RIP2(E,A,B,C,D,WR56,SR56,KR3);
436	RIP2(D,E,A,B,C,WR57,SR57,KR3);
437	RIP2(C,D,E,A,B,WR58,SR58,KR3);
438	RIP2(B,C,D,E,A,WR59,SR59,KR3);
439	RIP2(A,B,C,D,E,WR60,SR60,KR3);
440	RIP2(E,A,B,C,D,WR61,SR61,KR3);
441	RIP2(D,E,A,B,C,WR62,SR62,KR3);
442	RIP2(C,D,E,A,B,WR63,SR63,KR3);
443
444	RIP1(B,C,D,E,A,WR64,SR64);
445	RIP1(A,B,C,D,E,WR65,SR65);
446	RIP1(E,A,B,C,D,WR66,SR66);
447	RIP1(D,E,A,B,C,WR67,SR67);
448	RIP1(C,D,E,A,B,WR68,SR68);
449	RIP1(B,C,D,E,A,WR69,SR69);
450	RIP1(A,B,C,D,E,WR70,SR70);
451	RIP1(E,A,B,C,D,WR71,SR71);
452	RIP1(D,E,A,B,C,WR72,SR72);
453	RIP1(C,D,E,A,B,WR73,SR73);
454	RIP1(B,C,D,E,A,WR74,SR74);
455	RIP1(A,B,C,D,E,WR75,SR75);
456	RIP1(E,A,B,C,D,WR76,SR76);
457	RIP1(D,E,A,B,C,WR77,SR77);
458	RIP1(C,D,E,A,B,WR78,SR78);
459	RIP1(B,C,D,E,A,WR79,SR79);
460
461	D     =ctx->B+c+D;
462	ctx->B=ctx->C+d+E;
463	ctx->C=ctx->D+e+A;
464	ctx->D=ctx->E+a+B;
465	ctx->E=ctx->A+b+C;
466	ctx->A=D;
467
468	X+=16;
469	num-=64;
470	if (num <= 0) break;
471		}
472	}
473#endif
474
475void RIPEMD160_Final(md, c)
476unsigned char *md;
477RIPEMD160_CTX *c;
478	{
479	int i,j;
480	u_int32_t l;
481	u_int32_t *p;
482	static unsigned char end[4]={0x80,0x00,0x00,0x00};
483	unsigned char *cp=end;
484
485	/* c->num should definitly have room for at least one more byte. */
486	p=c->data;
487	j=c->num;
488	i=j>>2;
489
490	/* purify often complains about the following line as an
491	 * Uninitialized Memory Read.  While this can be true, the
492	 * following p_c2l macro will reset l when that case is true.
493	 * This is because j&0x03 contains the number of 'valid' bytes
494	 * already in p[i].  If and only if j&0x03 == 0, the UMR will
495	 * occur but this is also the only time p_c2l will do
496	 * l= *(cp++) instead of l|= *(cp++)
497	 * Many thanks to Alex Tang <altitude@cic.net> for pickup this
498	 * 'potential bug' */
499#ifdef PURIFY
500	if ((j&0x03) == 0) p[i]=0;
501#endif
502	l=p[i];
503	p_c2l(cp,l,j&0x03);
504	p[i]=l;
505	i++;
506	/* i is the next 'undefined word' */
507	if (c->num >= RIPEMD160_LAST_BLOCK)
508		{
509		for (; i<RIPEMD160_LBLOCK; i++)
510			p[i]=0;
511		ripemd160_block(c,p,64);
512		i=0;
513		}
514	for (; i<(RIPEMD160_LBLOCK-2); i++)
515		p[i]=0;
516	p[RIPEMD160_LBLOCK-2]=c->Nl;
517	p[RIPEMD160_LBLOCK-1]=c->Nh;
518	ripemd160_block(c,p,64);
519	cp=md;
520	l=c->A; l2c(l,cp);
521	l=c->B; l2c(l,cp);
522	l=c->C; l2c(l,cp);
523	l=c->D; l2c(l,cp);
524	l=c->E; l2c(l,cp);
525
526	/* clear stuff, ripemd160_block may be leaving some stuff on the stack
527	 * but I'm not worried :-) */
528	c->num=0;
529/*	memset((char *)&c,0,sizeof(c));*/
530	}
531
532#ifdef undef
533int printit(l)
534unsigned long *l;
535	{
536	int i,ii;
537
538	for (i=0; i<2; i++)
539		{
540		for (ii=0; ii<8; ii++)
541			{
542			fprintf(stderr,"%08lx ",l[i*8+ii]);
543			}
544		fprintf(stderr,"\n");
545		}
546	}
547#endif
548
549#ifdef WEAK_REFS
550/* When building libmd, provide weak references. Note: this is not
551   activated in the context of compiling these sources for internal
552   use in libcrypt.
553 */
554#undef RIPEMD160_Init
555__weak_reference(_libmd_RIPEMD160_Init, RIPEMD160_Init);
556#undef RIPEMD160_Update
557__weak_reference(_libmd_RIPEMD160_Update, RIPEMD160_Update);
558#undef RIPEMD160_Final
559__weak_reference(_libmd_RIPEMD160_Final, RIPEMD160_Final);
560#undef RIPEMD160_Transform
561__weak_reference(_libmd_RIPEMD160_Transform, RIPEMD160_Transform);
562#undef RMD160_version
563__weak_reference(_libmd_RMD160_version, RMD160_version);
564#undef ripemd160_block
565__weak_reference(_libmd_ripemd160_block, ripemd160_block);
566#endif
567