contrib/bearssl/src/hash/md5.c

   1 /*
   2  * Copyright (c) 2016 Thomas Pornin <pornin@bolet.org>
   3  *
   4  * Permission is hereby granted, free of charge, to any person obtaining
   5  * a copy of this software and associated documentation files (the
   6  * "Software"), to deal in the Software without restriction, including
   7  * without limitation the rights to use, copy, modify, merge, publish,
   8  * distribute, sublicense, and/or sell copies of the Software, and to
   9  * permit persons to whom the Software is furnished to do so, subject to
  10  * the following conditions:
  11  *
  12  * The above copyright notice and this permission notice shall be
  13  * included in all copies or substantial portions of the Software.
  14  *
  15  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
  16  * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
  17  * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
  18  * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS
  19  * BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN
  20  * ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
  21  * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
  22  * SOFTWARE.
  23  */
  24
  25 #include "inner.h"
  26
  27 #define F(B, C, D)     ((((C) ^ (D)) & (B)) ^ (D))
  28 #define G(B, C, D)     ((((C) ^ (B)) & (D)) ^ (C))
  29 #define H(B, C, D)     ((B) ^ (C) ^ (D))
  30 #define I(B, C, D)     ((C) ^ ((B) | ~(D)))
  31
  32 #define ROTL(x, n)    (((x) << (n)) | ((x) >> (32 - (n))))
  33
  34 /* see inner.h */
  35 const uint32_t br_md5_IV[4] = {
  36         0x67452301, 0xEFCDAB89, 0x98BADCFE, 0x10325476
  37 };
  38
  39 static const uint32_t K[64] = {
  40         0xD76AA478, 0xE8C7B756, 0x242070DB, 0xC1BDCEEE,
  41         0xF57C0FAF, 0x4787C62A, 0xA8304613, 0xFD469501,
  42         0x698098D8, 0x8B44F7AF, 0xFFFF5BB1, 0x895CD7BE,
  43         0x6B901122, 0xFD987193, 0xA679438E, 0x49B40821,
  44
  45         0xF61E2562, 0xC040B340, 0x265E5A51, 0xE9B6C7AA,
  46         0xD62F105D, 0x02441453, 0xD8A1E681, 0xE7D3FBC8,
  47         0x21E1CDE6, 0xC33707D6, 0xF4D50D87, 0x455A14ED,
  48         0xA9E3E905, 0xFCEFA3F8, 0x676F02D9, 0x8D2A4C8A,
  49
  50         0xFFFA3942, 0x8771F681, 0x6D9D6122, 0xFDE5380C,
  51         0xA4BEEA44, 0x4BDECFA9, 0xF6BB4B60, 0xBEBFBC70,
  52         0x289B7EC6, 0xEAA127FA, 0xD4EF3085, 0x04881D05,
  53         0xD9D4D039, 0xE6DB99E5, 0x1FA27CF8, 0xC4AC5665,
  54
  55         0xF4292244, 0x432AFF97, 0xAB9423A7, 0xFC93A039,
  56         0x655B59C3, 0x8F0CCC92, 0xFFEFF47D, 0x85845DD1,
  57         0x6FA87E4F, 0xFE2CE6E0, 0xA3014314, 0x4E0811A1,
  58         0xF7537E82, 0xBD3AF235, 0x2AD7D2BB, 0xEB86D391
  59 };
  60
  61 static const unsigned char MP[48] = {
  62         1, 6, 11, 0, 5, 10, 15, 4, 9, 14, 3, 8, 13, 2, 7, 12,
  63         5, 8, 11, 14, 1, 4, 7, 10, 13, 0, 3, 6, 9, 12, 15, 2,
  64         0, 7, 14, 5, 12, 3, 10, 1, 8, 15, 6, 13, 4, 11, 2, 9
  65 };
  66
  67 /* see inner.h */
  68 void
  69 br_md5_round(const unsigned char *buf, uint32_t *val)
  70 {
  71         uint32_t m[16];
  72         uint32_t a, b, c, d;
  73         int i;
  74
  75         a = val[0];
  76         b = val[1];
  77         c = val[2];
  78         d = val[3];
  79         /* obsolete
  80         for (i = 0; i < 16; i ++) {
  81                 m[i] = br_dec32le(buf + (i << 2));
  82         }
  83         */
  84         br_range_dec32le(m, 16, buf);
  85
  86         for (i = 0; i < 16; i += 4) {
  87                 a = b + ROTL(a + F(b, c, d) + m[i + 0] + K[i + 0],  7);
  88                 d = a + ROTL(d + F(a, b, c) + m[i + 1] + K[i + 1], 12);
  89                 c = d + ROTL(c + F(d, a, b) + m[i + 2] + K[i + 2], 17);
  90                 b = c + ROTL(b + F(c, d, a) + m[i + 3] + K[i + 3], 22);
  91         }
  92         for (i = 16; i < 32; i += 4) {
  93                 a = b + ROTL(a + G(b, c, d) + m[MP[i - 16]] + K[i + 0],  5);
  94                 d = a + ROTL(d + G(a, b, c) + m[MP[i - 15]] + K[i + 1],  9);
  95                 c = d + ROTL(c + G(d, a, b) + m[MP[i - 14]] + K[i + 2], 14);
  96                 b = c + ROTL(b + G(c, d, a) + m[MP[i - 13]] + K[i + 3], 20);
  97         }
  98         for (i = 32; i < 48; i += 4) {
  99                 a = b + ROTL(a + H(b, c, d) + m[MP[i - 16]] + K[i + 0],  4);
 100                 d = a + ROTL(d + H(a, b, c) + m[MP[i - 15]] + K[i + 1], 11);
 101                 c = d + ROTL(c + H(d, a, b) + m[MP[i - 14]] + K[i + 2], 16);
 102                 b = c + ROTL(b + H(c, d, a) + m[MP[i - 13]] + K[i + 3], 23);
 103         }
 104         for (i = 48; i < 64; i += 4) {
 105                 a = b + ROTL(a + I(b, c, d) + m[MP[i - 16]] + K[i + 0],  6);
 106                 d = a + ROTL(d + I(a, b, c) + m[MP[i - 15]] + K[i + 1], 10);
 107                 c = d + ROTL(c + I(d, a, b) + m[MP[i - 14]] + K[i + 2], 15);
 108                 b = c + ROTL(b + I(c, d, a) + m[MP[i - 13]] + K[i + 3], 21);
 109         }
 110
 111         val[0] += a;
 112         val[1] += b;
 113         val[2] += c;
 114         val[3] += d;
 115 }
 116
 117 /* see bearssl.h */
 118 void
 119 br_md5_init(br_md5_context *cc)
 120 {
 121         cc->vtable = &br_md5_vtable;
 122         memcpy(cc->val, br_md5_IV, sizeof cc->val);
 123         cc->count = 0;
 124 }
 125
 126 /* see bearssl.h */
 127 void
 128 br_md5_update(br_md5_context *cc, const void *data, size_t len)
 129 {
 130         const unsigned char *buf;
 131         size_t ptr;
 132
 133         buf = data;
 134         ptr = (size_t)cc->count & 63;
 135         while (len > 0) {
 136                 size_t clen;
 137
 138                 clen = 64 - ptr;
 139                 if (clen > len) {
 140                         clen = len;
 141                 }
 142                 memcpy(cc->buf + ptr, buf, clen);
 143                 ptr += clen;
 144                 buf += clen;
 145                 len -= clen;
 146                 cc->count += (uint64_t)clen;
 147                 if (ptr == 64) {
 148                         br_md5_round(cc->buf, cc->val);
 149                         ptr = 0;
 150                 }
 151         }
 152 }
 153
 154 /* see bearssl.h */
 155 void
 156 br_md5_out(const br_md5_context *cc, void *dst)
 157 {
 158         unsigned char buf[64];
 159         uint32_t val[4];
 160         size_t ptr;
 161
 162         ptr = (size_t)cc->count & 63;
 163         memcpy(buf, cc->buf, ptr);
 164         memcpy(val, cc->val, sizeof val);
 165         buf[ptr ++] = 0x80;
 166         if (ptr > 56) {
 167                 memset(buf + ptr, 0, 64 - ptr);
 168                 br_md5_round(buf, val);
 169                 memset(buf, 0, 56);
 170         } else {
 171                 memset(buf + ptr, 0, 56 - ptr);
 172         }
 173         br_enc64le(buf + 56, cc->count << 3);
 174         br_md5_round(buf, val);
 175         br_range_enc32le(dst, val, 4);
 176 }
 177
 178 /* see bearssl.h */
 179 uint64_t
 180 br_md5_state(const br_md5_context *cc, void *dst)
 181 {
 182         br_range_enc32le(dst, cc->val, 4);
 183         return cc->count;
 184 }
 185
 186 /* see bearssl.h */
 187 void
 188 br_md5_set_state(br_md5_context *cc, const void *stb, uint64_t count)
 189 {
 190         br_range_dec32le(cc->val, 4, stb);
 191         cc->count = count;
 192 }
 193
 194 /* see bearssl.h */
 195 const br_hash_class br_md5_vtable = {
 196         sizeof(br_md5_context),
 197         BR_HASHDESC_ID(br_md5_ID)
 198                 | BR_HASHDESC_OUT(16)
 199                 | BR_HASHDESC_STATE(16)
 200                 | BR_HASHDESC_LBLEN(6)
 201                 | BR_HASHDESC_MD_PADDING,
 202         (void (*)(const br_hash_class **))&br_md5_init,
 203         (void (*)(const br_hash_class **, const void *, size_t))&br_md5_update,
 204         (void (*)(const br_hash_class *const *, void *))&br_md5_out,
 205         (uint64_t (*)(const br_hash_class *const *, void *))&br_md5_state,
 206         (void (*)(const br_hash_class **, const void *, uint64_t))
 207                 &br_md5_set_state
 208 };