libsyncml  0.5.4
sml_md5.c
1 /* parts of this file are :
2  * Written March 1993 by Branko Lankester
3  * Modified June 1993 by Colin Plumb for altered md5.c.
4  * Modified October 1995 by Erik Troan for RPM
5  * Modified again 2005 by Armin Bauer
6  * This code is in the public domain
7  *
8  */
9 
10 #include "syncml.h"
11 #include "syncml_internals.h"
12 #include "sml_error_internals.h"
13 
14 #include <string.h>
15 #include <stdint.h>
16 #include "sml_md5.h"
17 #include <glib.h>
18 
19 typedef struct SmlMD5Context {
20  uint32_t buf[4];
21  uint32_t bits[2];
22  unsigned char in[64];
23  int doByteReverse;
25 
26 void smlMD5Init(SmlMD5Context *ctx);
27 void smlMD5Update(SmlMD5Context *ctx, const char *buf, uint32_t len);
28 void smlMD5Final (SmlMD5Context *ctx, unsigned char digest[16]);
29 static void smlMD5Transform(uint32_t buf[4], const uint32_t in[16]);
30 
31 static gint _ie = 0x44332211;
32 static union _endian { gint i; gchar b[4]; } *_endian = (union _endian *)&_ie;
33 #define IS_BIG_ENDIAN() (_endian->b[0] == '\x44')
34 #define IS_LITTLE_ENDIAN() (_endian->b[0] == '\x11')
35 
36 /* The four core functions - F1 is optimized somewhat */
37 
38 /* #define F1(x, y, z) (x & y | ~x & z) */
39 #define F1(x, y, z) (z ^ (x & (y ^ z)))
40 #define F2(x, y, z) F1(z, x, y)
41 #define F3(x, y, z) (x ^ y ^ z)
42 #define F4(x, y, z) (y ^ (x | ~z))
43 
44 /* This is the central step in the MD5 algorithm. */
45 #define MD5STEP(f, w, x, y, z, data, s) ( w += f(x, y, z) + data, w = w<<s | w>>(32-s), w += x )
46 
47 /*
48  * Note: this code is harmless on little-endian machines.
49  */
50 static void
51 _byte_reverse (guchar *buf, guint32 longs)
52 {
53  guint32 t;
54  do {
55  t = (guint32) ((guint32) buf[3] << 8 | buf[2]) << 16 |
56  ((guint32) buf[1] << 8 | buf[0]);
57  *(guint32 *) buf = t;
58  buf += 4;
59  } while (--longs);
60 }
61 
62 void smlMD5Init(SmlMD5Context *ctx)
63 {
64  ctx->buf[0] = 0x67452301;
65  ctx->buf[1] = 0xefcdab89;
66  ctx->buf[2] = 0x98badcfe;
67  ctx->buf[3] = 0x10325476;
68 
69  ctx->bits[0] = 0;
70  ctx->bits[1] = 0;
71 
72  if (IS_BIG_ENDIAN())
73  ctx->doByteReverse = 1;
74  else
75  ctx->doByteReverse = 0;
76 }
77 
78 void smlMD5Update(SmlMD5Context *ctx, const char *buf, uint32_t len)
79 {
80  guint32 t;
81 
82  /* Update bitcount */
83 
84  t = ctx->bits[0];
85  if ((ctx->bits[0] = t + ((guint32) len << 3)) < t)
86  ctx->bits[1]++; /* Carry from low to high */
87  ctx->bits[1] += len >> 29;
88 
89  t = (t >> 3) & 0x3f; /* Bytes already in shsInfo->data */
90 
91  /* Handle any leading odd-sized chunks */
92 
93  if (t) {
94  guchar *p = (guchar *) ctx->in + t;
95 
96  t = 64 - t;
97  if (len < t) {
98  memcpy (p, buf, len);
99  return;
100  }
101  memcpy (p, buf, t);
102  if (ctx->doByteReverse)
103  _byte_reverse (ctx->in, 16);
104  smlMD5Transform (ctx->buf, (guint32 *) ctx->in);
105  buf += t;
106  len -= t;
107  }
108  /* Process data in 64-byte chunks */
109 
110  while (len >= 64) {
111  memcpy (ctx->in, buf, 64);
112  if (ctx->doByteReverse)
113  _byte_reverse (ctx->in, 16);
114  smlMD5Transform (ctx->buf, (guint32 *) ctx->in);
115  buf += 64;
116  len -= 64;
117  }
118 
119  /* Handle any remaining bytes of data. */
120 
121  memcpy (ctx->in, buf, len);
122 }
123 
124 void smlMD5Final (SmlMD5Context *ctx, unsigned char digest[16])
125 {
126  guint32 count;
127  guchar *p;
128 
129  /* Compute number of bytes mod 64 */
130  count = (ctx->bits[0] >> 3) & 0x3F;
131 
132  /* Set the first char of padding to 0x80. This is safe since there is
133  always at least one byte free */
134  p = ctx->in + count;
135  *p++ = 0x80;
136 
137  /* Bytes of padding needed to make 64 bytes */
138  count = 64 - 1 - count;
139 
140  /* Pad out to 56 mod 64 */
141  if (count < 8) {
142  /* Two lots of padding: Pad the first block to 64 bytes */
143  memset (p, 0, count);
144  if (ctx->doByteReverse)
145  _byte_reverse (ctx->in, 16);
146  smlMD5Transform (ctx->buf, (guint32 *) ctx->in);
147 
148  /* Now fill the next block with 56 bytes */
149  memset (ctx->in, 0, 56);
150  } else {
151  /* Pad block to 56 bytes */
152  memset (p, 0, count - 8);
153  }
154  if (ctx->doByteReverse)
155  _byte_reverse (ctx->in, 14);
156 
157  /* Append length in bits and transform */
158  ((guint32 *) ctx->in)[14] = ctx->bits[0];
159  ((guint32 *) ctx->in)[15] = ctx->bits[1];
160 
161  smlMD5Transform (ctx->buf, (guint32 *) ctx->in);
162  if (ctx->doByteReverse)
163  _byte_reverse ((guchar *) ctx->buf, 4);
164  memcpy (digest, ctx->buf, 16);
165 }
166 
167 static void smlMD5Transform(uint32_t buf[4], const uint32_t in[16])
168 {
169  register guint32 a, b, c, d;
170 
171  a = buf[0];
172  b = buf[1];
173  c = buf[2];
174  d = buf[3];
175 
176  MD5STEP (F1, a, b, c, d, in[0] + 0xd76aa478, 7);
177  MD5STEP (F1, d, a, b, c, in[1] + 0xe8c7b756, 12);
178  MD5STEP (F1, c, d, a, b, in[2] + 0x242070db, 17);
179  MD5STEP (F1, b, c, d, a, in[3] + 0xc1bdceee, 22);
180  MD5STEP (F1, a, b, c, d, in[4] + 0xf57c0faf, 7);
181  MD5STEP (F1, d, a, b, c, in[5] + 0x4787c62a, 12);
182  MD5STEP (F1, c, d, a, b, in[6] + 0xa8304613, 17);
183  MD5STEP (F1, b, c, d, a, in[7] + 0xfd469501, 22);
184  MD5STEP (F1, a, b, c, d, in[8] + 0x698098d8, 7);
185  MD5STEP (F1, d, a, b, c, in[9] + 0x8b44f7af, 12);
186  MD5STEP (F1, c, d, a, b, in[10] + 0xffff5bb1, 17);
187  MD5STEP (F1, b, c, d, a, in[11] + 0x895cd7be, 22);
188  MD5STEP (F1, a, b, c, d, in[12] + 0x6b901122, 7);
189  MD5STEP (F1, d, a, b, c, in[13] + 0xfd987193, 12);
190  MD5STEP (F1, c, d, a, b, in[14] + 0xa679438e, 17);
191  MD5STEP (F1, b, c, d, a, in[15] + 0x49b40821, 22);
192 
193  MD5STEP (F2, a, b, c, d, in[1] + 0xf61e2562, 5);
194  MD5STEP (F2, d, a, b, c, in[6] + 0xc040b340, 9);
195  MD5STEP (F2, c, d, a, b, in[11] + 0x265e5a51, 14);
196  MD5STEP (F2, b, c, d, a, in[0] + 0xe9b6c7aa, 20);
197  MD5STEP (F2, a, b, c, d, in[5] + 0xd62f105d, 5);
198  MD5STEP (F2, d, a, b, c, in[10] + 0x02441453, 9);
199  MD5STEP (F2, c, d, a, b, in[15] + 0xd8a1e681, 14);
200  MD5STEP (F2, b, c, d, a, in[4] + 0xe7d3fbc8, 20);
201  MD5STEP (F2, a, b, c, d, in[9] + 0x21e1cde6, 5);
202  MD5STEP (F2, d, a, b, c, in[14] + 0xc33707d6, 9);
203  MD5STEP (F2, c, d, a, b, in[3] + 0xf4d50d87, 14);
204  MD5STEP (F2, b, c, d, a, in[8] + 0x455a14ed, 20);
205  MD5STEP (F2, a, b, c, d, in[13] + 0xa9e3e905, 5);
206  MD5STEP (F2, d, a, b, c, in[2] + 0xfcefa3f8, 9);
207  MD5STEP (F2, c, d, a, b, in[7] + 0x676f02d9, 14);
208  MD5STEP (F2, b, c, d, a, in[12] + 0x8d2a4c8a, 20);
209 
210  MD5STEP (F3, a, b, c, d, in[5] + 0xfffa3942, 4);
211  MD5STEP (F3, d, a, b, c, in[8] + 0x8771f681, 11);
212  MD5STEP (F3, c, d, a, b, in[11] + 0x6d9d6122, 16);
213  MD5STEP (F3, b, c, d, a, in[14] + 0xfde5380c, 23);
214  MD5STEP (F3, a, b, c, d, in[1] + 0xa4beea44, 4);
215  MD5STEP (F3, d, a, b, c, in[4] + 0x4bdecfa9, 11);
216  MD5STEP (F3, c, d, a, b, in[7] + 0xf6bb4b60, 16);
217  MD5STEP (F3, b, c, d, a, in[10] + 0xbebfbc70, 23);
218  MD5STEP (F3, a, b, c, d, in[13] + 0x289b7ec6, 4);
219  MD5STEP (F3, d, a, b, c, in[0] + 0xeaa127fa, 11);
220  MD5STEP (F3, c, d, a, b, in[3] + 0xd4ef3085, 16);
221  MD5STEP (F3, b, c, d, a, in[6] + 0x04881d05, 23);
222  MD5STEP (F3, a, b, c, d, in[9] + 0xd9d4d039, 4);
223  MD5STEP (F3, d, a, b, c, in[12] + 0xe6db99e5, 11);
224  MD5STEP (F3, c, d, a, b, in[15] + 0x1fa27cf8, 16);
225  MD5STEP (F3, b, c, d, a, in[2] + 0xc4ac5665, 23);
226 
227  MD5STEP (F4, a, b, c, d, in[0] + 0xf4292244, 6);
228  MD5STEP (F4, d, a, b, c, in[7] + 0x432aff97, 10);
229  MD5STEP (F4, c, d, a, b, in[14] + 0xab9423a7, 15);
230  MD5STEP (F4, b, c, d, a, in[5] + 0xfc93a039, 21);
231  MD5STEP (F4, a, b, c, d, in[12] + 0x655b59c3, 6);
232  MD5STEP (F4, d, a, b, c, in[3] + 0x8f0ccc92, 10);
233  MD5STEP (F4, c, d, a, b, in[10] + 0xffeff47d, 15);
234  MD5STEP (F4, b, c, d, a, in[1] + 0x85845dd1, 21);
235  MD5STEP (F4, a, b, c, d, in[8] + 0x6fa87e4f, 6);
236  MD5STEP (F4, d, a, b, c, in[15] + 0xfe2ce6e0, 10);
237  MD5STEP (F4, c, d, a, b, in[6] + 0xa3014314, 15);
238  MD5STEP (F4, b, c, d, a, in[13] + 0x4e0811a1, 21);
239  MD5STEP (F4, a, b, c, d, in[4] + 0xf7537e82, 6);
240  MD5STEP (F4, d, a, b, c, in[11] + 0xbd3af235, 10);
241  MD5STEP (F4, c, d, a, b, in[2] + 0x2ad7d2bb, 15);
242  MD5STEP (F4, b, c, d, a, in[9] + 0xeb86d391, 21);
243 
244  buf[0] += a;
245  buf[1] += b;
246  buf[2] += c;
247  buf[3] += d;
248 }
249 
250 void smlMD5GetDigest (const char *buffer, int buffer_size, unsigned char digest[16])
251 {
252  SmlMD5Context ctx;
253 
254  smlMD5Init (&ctx);
255  smlMD5Update (&ctx, buffer, buffer_size);
256  smlMD5Final (&ctx, digest);
257 
258 }
259 
260 char *smlMD5ToString(unsigned char digest[16], SmlError **error)
261 {
262  CHECK_ERROR_REF
263  char *string = smlTryMalloc0(33, error);
264  if (!string)
265  return NULL;
266 
267  snprintf(string, 33, "%02x%02x%02x%02x%02x%02x%02x%02x%02x%02x%02x%02x%02x%02x%02x%02x", digest[0], digest[1], digest[2], digest[3], digest[4], digest[5], digest[6], digest[7], digest[8], digest[9], digest[10], digest[11], digest[12], digest[13], digest[14], digest[15]);
268  return string;
269 }