lnsocket

A minimal C library for connecting to the lightning network
git clone git://jb55.com/lnsocket
Log | Files | Refs | Submodules | README | LICENSE

sha512.c (11787B)


      1 
      2 /* MIT (BSD) license - see LICENSE file for details */
      3 /* SHA512 core code translated from the Bitcoin project's C++:
      4  *
      5  * src/crypto/sha512.cpp commit f914f1a746d7f91951c1da262a4a749dd3ebfa71
      6  * Copyright (c) 2014 The Bitcoin Core developers
      7  * Distributed under the MIT software license, see the accompanying
      8  * file COPYING or http://www.opensource.org/licenses/mit-license.php.
      9  */
     10 #include "sha512.h"
     11 #include "endian.h"
     12 #include "compiler.h"
     13 #include <stdbool.h>
     14 #include <assert.h>
     15 #include <string.h>
     16 
     17 static void invalidate_sha512(struct sha512_ctx *ctx)
     18 {
     19 #ifdef CCAN_CRYPTO_SHA512_USE_OPENSSL
     20 	ctx->c.md_len = 0;
     21 #else
     22 	ctx->bytes = (size_t)-1;
     23 #endif
     24 }
     25 
     26 static void check_sha512(struct sha512_ctx *ctx)
     27 {
     28 #ifdef CCAN_CRYPTO_SHA512_USE_OPENSSL
     29 	assert(ctx->c.md_len != 0);
     30 #else
     31 	assert(ctx->bytes != (size_t)-1);
     32 #endif
     33 }
     34 
     35 #ifdef CCAN_CRYPTO_SHA512_USE_OPENSSL
     36 void sha512_init(struct sha512_ctx *ctx)
     37 {
     38 	SHA512_Init(&ctx->c);
     39 }
     40 
     41 void sha512_update(struct sha512_ctx *ctx, const void *p, size_t size)
     42 {
     43 	check_sha512(ctx);
     44 	SHA512_Update(&ctx->c, p, size);
     45 }
     46 
     47 void sha512_done(struct sha512_ctx *ctx, struct sha512 *res)
     48 {
     49 	SHA512_Final(res->u.u8, &ctx->c);
     50 	invalidate_sha512(ctx);
     51 }
     52 #else
     53 static uint64_t Ch(uint64_t x, uint64_t y, uint64_t z)
     54 {
     55 	return z ^ (x & (y ^ z));
     56 }
     57 static uint64_t Maj(uint64_t x, uint64_t y, uint64_t z)
     58 {
     59 	return (x & y) | (z & (x | y));
     60 }
     61 static uint64_t Sigma0(uint64_t x)
     62 {
     63 	return (x >> 28 | x << 36) ^ (x >> 34 | x << 30) ^ (x >> 39 | x << 25);
     64 }
     65 static uint64_t Sigma1(uint64_t x)
     66 {
     67 	return (x >> 14 | x << 50) ^ (x >> 18 | x << 46) ^ (x >> 41 | x << 23);
     68 }
     69 static uint64_t sigma0(uint64_t x)
     70 {
     71 	return (x >> 1 | x << 63) ^ (x >> 8 | x << 56) ^ (x >> 7);
     72 }
     73 static uint64_t sigma1(uint64_t x)
     74 {
     75 	return (x >> 19 | x << 45) ^ (x >> 61 | x << 3) ^ (x >> 6);
     76 }
     77 
     78 /** One round of SHA-512. */
     79 static void Round(uint64_t a, uint64_t b, uint64_t c, uint64_t *d, uint64_t e, uint64_t f, uint64_t g, uint64_t *h, uint64_t k, uint64_t w)
     80 {
     81 	uint64_t t1 = *h + Sigma1(e) + Ch(e, f, g) + k + w;
     82 	uint64_t t2 = Sigma0(a) + Maj(a, b, c);
     83 	*d += t1;
     84 	*h = t1 + t2;
     85 }
     86 
     87 /** Perform one SHA-512 transformation, processing a 128-byte chunk. */
     88 static void Transform(uint64_t *s, const uint64_t *chunk)
     89 {
     90 	uint64_t a = s[0], b = s[1], c = s[2], d = s[3], e = s[4], f = s[5], g = s[6], h = s[7];
     91 	uint64_t w0, w1, w2, w3, w4, w5, w6, w7, w8, w9, w10, w11, w12, w13, w14, w15;
     92 
     93 	Round(a, b, c, &d, e, f, g, &h, 0x428a2f98d728ae22ull, w0 = be64_to_cpu(chunk[0]));
     94 	Round(h, a, b, &c, d, e, f, &g, 0x7137449123ef65cdull, w1 = be64_to_cpu(chunk[1]));
     95 	Round(g, h, a, &b, c, d, e, &f, 0xb5c0fbcfec4d3b2full, w2 = be64_to_cpu(chunk[2]));
     96 	Round(f, g, h, &a, b, c, d, &e, 0xe9b5dba58189dbbcull, w3 = be64_to_cpu(chunk[3]));
     97 	Round(e, f, g, &h, a, b, c, &d, 0x3956c25bf348b538ull, w4 = be64_to_cpu(chunk[4]));
     98 	Round(d, e, f, &g, h, a, b, &c, 0x59f111f1b605d019ull, w5 = be64_to_cpu(chunk[5]));
     99 	Round(c, d, e, &f, g, h, a, &b, 0x923f82a4af194f9bull, w6 = be64_to_cpu(chunk[6]));
    100 	Round(b, c, d, &e, f, g, h, &a, 0xab1c5ed5da6d8118ull, w7 = be64_to_cpu(chunk[7]));
    101 	Round(a, b, c, &d, e, f, g, &h, 0xd807aa98a3030242ull, w8 = be64_to_cpu(chunk[8]));
    102 	Round(h, a, b, &c, d, e, f, &g, 0x12835b0145706fbeull, w9 = be64_to_cpu(chunk[9]));
    103 	Round(g, h, a, &b, c, d, e, &f, 0x243185be4ee4b28cull, w10 = be64_to_cpu(chunk[10]));
    104 	Round(f, g, h, &a, b, c, d, &e, 0x550c7dc3d5ffb4e2ull, w11 = be64_to_cpu(chunk[11]));
    105 	Round(e, f, g, &h, a, b, c, &d, 0x72be5d74f27b896full, w12 = be64_to_cpu(chunk[12]));
    106 	Round(d, e, f, &g, h, a, b, &c, 0x80deb1fe3b1696b1ull, w13 = be64_to_cpu(chunk[13]));
    107 	Round(c, d, e, &f, g, h, a, &b, 0x9bdc06a725c71235ull, w14 = be64_to_cpu(chunk[14]));
    108 	Round(b, c, d, &e, f, g, h, &a, 0xc19bf174cf692694ull, w15 = be64_to_cpu(chunk[15]));
    109 
    110 	Round(a, b, c, &d, e, f, g, &h, 0xe49b69c19ef14ad2ull, w0 += sigma1(w14) + w9 + sigma0(w1));
    111 	Round(h, a, b, &c, d, e, f, &g, 0xefbe4786384f25e3ull, w1 += sigma1(w15) + w10 + sigma0(w2));
    112 	Round(g, h, a, &b, c, d, e, &f, 0x0fc19dc68b8cd5b5ull, w2 += sigma1(w0) + w11 + sigma0(w3));
    113 	Round(f, g, h, &a, b, c, d, &e, 0x240ca1cc77ac9c65ull, w3 += sigma1(w1) + w12 + sigma0(w4));
    114 	Round(e, f, g, &h, a, b, c, &d, 0x2de92c6f592b0275ull, w4 += sigma1(w2) + w13 + sigma0(w5));
    115 	Round(d, e, f, &g, h, a, b, &c, 0x4a7484aa6ea6e483ull, w5 += sigma1(w3) + w14 + sigma0(w6));
    116 	Round(c, d, e, &f, g, h, a, &b, 0x5cb0a9dcbd41fbd4ull, w6 += sigma1(w4) + w15 + sigma0(w7));
    117 	Round(b, c, d, &e, f, g, h, &a, 0x76f988da831153b5ull, w7 += sigma1(w5) + w0 + sigma0(w8));
    118 	Round(a, b, c, &d, e, f, g, &h, 0x983e5152ee66dfabull, w8 += sigma1(w6) + w1 + sigma0(w9));
    119 	Round(h, a, b, &c, d, e, f, &g, 0xa831c66d2db43210ull, w9 += sigma1(w7) + w2 + sigma0(w10));
    120 	Round(g, h, a, &b, c, d, e, &f, 0xb00327c898fb213full, w10 += sigma1(w8) + w3 + sigma0(w11));
    121 	Round(f, g, h, &a, b, c, d, &e, 0xbf597fc7beef0ee4ull, w11 += sigma1(w9) + w4 + sigma0(w12));
    122 	Round(e, f, g, &h, a, b, c, &d, 0xc6e00bf33da88fc2ull, w12 += sigma1(w10) + w5 + sigma0(w13));
    123 	Round(d, e, f, &g, h, a, b, &c, 0xd5a79147930aa725ull, w13 += sigma1(w11) + w6 + sigma0(w14));
    124 	Round(c, d, e, &f, g, h, a, &b, 0x06ca6351e003826full, w14 += sigma1(w12) + w7 + sigma0(w15));
    125 	Round(b, c, d, &e, f, g, h, &a, 0x142929670a0e6e70ull, w15 += sigma1(w13) + w8 + sigma0(w0));
    126 
    127 	Round(a, b, c, &d, e, f, g, &h, 0x27b70a8546d22ffcull, w0 += sigma1(w14) + w9 + sigma0(w1));
    128 	Round(h, a, b, &c, d, e, f, &g, 0x2e1b21385c26c926ull, w1 += sigma1(w15) + w10 + sigma0(w2));
    129 	Round(g, h, a, &b, c, d, e, &f, 0x4d2c6dfc5ac42aedull, w2 += sigma1(w0) + w11 + sigma0(w3));
    130 	Round(f, g, h, &a, b, c, d, &e, 0x53380d139d95b3dfull, w3 += sigma1(w1) + w12 + sigma0(w4));
    131 	Round(e, f, g, &h, a, b, c, &d, 0x650a73548baf63deull, w4 += sigma1(w2) + w13 + sigma0(w5));
    132 	Round(d, e, f, &g, h, a, b, &c, 0x766a0abb3c77b2a8ull, w5 += sigma1(w3) + w14 + sigma0(w6));
    133 	Round(c, d, e, &f, g, h, a, &b, 0x81c2c92e47edaee6ull, w6 += sigma1(w4) + w15 + sigma0(w7));
    134 	Round(b, c, d, &e, f, g, h, &a, 0x92722c851482353bull, w7 += sigma1(w5) + w0 + sigma0(w8));
    135 	Round(a, b, c, &d, e, f, g, &h, 0xa2bfe8a14cf10364ull, w8 += sigma1(w6) + w1 + sigma0(w9));
    136 	Round(h, a, b, &c, d, e, f, &g, 0xa81a664bbc423001ull, w9 += sigma1(w7) + w2 + sigma0(w10));
    137 	Round(g, h, a, &b, c, d, e, &f, 0xc24b8b70d0f89791ull, w10 += sigma1(w8) + w3 + sigma0(w11));
    138 	Round(f, g, h, &a, b, c, d, &e, 0xc76c51a30654be30ull, w11 += sigma1(w9) + w4 + sigma0(w12));
    139 	Round(e, f, g, &h, a, b, c, &d, 0xd192e819d6ef5218ull, w12 += sigma1(w10) + w5 + sigma0(w13));
    140 	Round(d, e, f, &g, h, a, b, &c, 0xd69906245565a910ull, w13 += sigma1(w11) + w6 + sigma0(w14));
    141 	Round(c, d, e, &f, g, h, a, &b, 0xf40e35855771202aull, w14 += sigma1(w12) + w7 + sigma0(w15));
    142 	Round(b, c, d, &e, f, g, h, &a, 0x106aa07032bbd1b8ull, w15 += sigma1(w13) + w8 + sigma0(w0));
    143 
    144 	Round(a, b, c, &d, e, f, g, &h, 0x19a4c116b8d2d0c8ull, w0 += sigma1(w14) + w9 + sigma0(w1));
    145 	Round(h, a, b, &c, d, e, f, &g, 0x1e376c085141ab53ull, w1 += sigma1(w15) + w10 + sigma0(w2));
    146 	Round(g, h, a, &b, c, d, e, &f, 0x2748774cdf8eeb99ull, w2 += sigma1(w0) + w11 + sigma0(w3));
    147 	Round(f, g, h, &a, b, c, d, &e, 0x34b0bcb5e19b48a8ull, w3 += sigma1(w1) + w12 + sigma0(w4));
    148 	Round(e, f, g, &h, a, b, c, &d, 0x391c0cb3c5c95a63ull, w4 += sigma1(w2) + w13 + sigma0(w5));
    149 	Round(d, e, f, &g, h, a, b, &c, 0x4ed8aa4ae3418acbull, w5 += sigma1(w3) + w14 + sigma0(w6));
    150 	Round(c, d, e, &f, g, h, a, &b, 0x5b9cca4f7763e373ull, w6 += sigma1(w4) + w15 + sigma0(w7));
    151 	Round(b, c, d, &e, f, g, h, &a, 0x682e6ff3d6b2b8a3ull, w7 += sigma1(w5) + w0 + sigma0(w8));
    152 	Round(a, b, c, &d, e, f, g, &h, 0x748f82ee5defb2fcull, w8 += sigma1(w6) + w1 + sigma0(w9));
    153 	Round(h, a, b, &c, d, e, f, &g, 0x78a5636f43172f60ull, w9 += sigma1(w7) + w2 + sigma0(w10));
    154 	Round(g, h, a, &b, c, d, e, &f, 0x84c87814a1f0ab72ull, w10 += sigma1(w8) + w3 + sigma0(w11));
    155 	Round(f, g, h, &a, b, c, d, &e, 0x8cc702081a6439ecull, w11 += sigma1(w9) + w4 + sigma0(w12));
    156 	Round(e, f, g, &h, a, b, c, &d, 0x90befffa23631e28ull, w12 += sigma1(w10) + w5 + sigma0(w13));
    157 	Round(d, e, f, &g, h, a, b, &c, 0xa4506cebde82bde9ull, w13 += sigma1(w11) + w6 + sigma0(w14));
    158 	Round(c, d, e, &f, g, h, a, &b, 0xbef9a3f7b2c67915ull, w14 += sigma1(w12) + w7 + sigma0(w15));
    159 	Round(b, c, d, &e, f, g, h, &a, 0xc67178f2e372532bull, w15 += sigma1(w13) + w8 + sigma0(w0));
    160 
    161 	Round(a, b, c, &d, e, f, g, &h, 0xca273eceea26619cull, w0 += sigma1(w14) + w9 + sigma0(w1));
    162 	Round(h, a, b, &c, d, e, f, &g, 0xd186b8c721c0c207ull, w1 += sigma1(w15) + w10 + sigma0(w2));
    163 	Round(g, h, a, &b, c, d, e, &f, 0xeada7dd6cde0eb1eull, w2 += sigma1(w0) + w11 + sigma0(w3));
    164 	Round(f, g, h, &a, b, c, d, &e, 0xf57d4f7fee6ed178ull, w3 += sigma1(w1) + w12 + sigma0(w4));
    165 	Round(e, f, g, &h, a, b, c, &d, 0x06f067aa72176fbaull, w4 += sigma1(w2) + w13 + sigma0(w5));
    166 	Round(d, e, f, &g, h, a, b, &c, 0x0a637dc5a2c898a6ull, w5 += sigma1(w3) + w14 + sigma0(w6));
    167 	Round(c, d, e, &f, g, h, a, &b, 0x113f9804bef90daeull, w6 += sigma1(w4) + w15 + sigma0(w7));
    168 	Round(b, c, d, &e, f, g, h, &a, 0x1b710b35131c471bull, w7 += sigma1(w5) + w0 + sigma0(w8));
    169 	Round(a, b, c, &d, e, f, g, &h, 0x28db77f523047d84ull, w8 += sigma1(w6) + w1 + sigma0(w9));
    170 	Round(h, a, b, &c, d, e, f, &g, 0x32caab7b40c72493ull, w9 += sigma1(w7) + w2 + sigma0(w10));
    171 	Round(g, h, a, &b, c, d, e, &f, 0x3c9ebe0a15c9bebcull, w10 += sigma1(w8) + w3 + sigma0(w11));
    172 	Round(f, g, h, &a, b, c, d, &e, 0x431d67c49c100d4cull, w11 += sigma1(w9) + w4 + sigma0(w12));
    173 	Round(e, f, g, &h, a, b, c, &d, 0x4cc5d4becb3e42b6ull, w12 += sigma1(w10) + w5 + sigma0(w13));
    174 	Round(d, e, f, &g, h, a, b, &c, 0x597f299cfc657e2aull, w13 += sigma1(w11) + w6 + sigma0(w14));
    175 	Round(c, d, e, &f, g, h, a, &b, 0x5fcb6fab3ad6faecull, w14 + sigma1(w12) + w7 + sigma0(w15));
    176 	Round(b, c, d, &e, f, g, h, &a, 0x6c44198c4a475817ull, w15 + sigma1(w13) + w8 + sigma0(w0));
    177 
    178 	s[0] += a;
    179 	s[1] += b;
    180 	s[2] += c;
    181 	s[3] += d;
    182 	s[4] += e;
    183 	s[5] += f;
    184 	s[6] += g;
    185 	s[7] += h;
    186 }
    187 
    188 static void add(struct sha512_ctx *ctx, const void *p, size_t len)
    189 {
    190 	const unsigned char *data = p;
    191 	size_t bufsize = ctx->bytes % 128;
    192 
    193 	if (bufsize + len >= 128) {
    194 		/* Fill the buffer, and process it. */
    195 		memcpy(ctx->buf.u8 + bufsize, data, 128 - bufsize);
    196 		ctx->bytes += 128 - bufsize;
    197 		data += 128 - bufsize;
    198 		len -= 128 - bufsize;
    199 		Transform(ctx->s, ctx->buf.u64);
    200 		bufsize = 0;
    201 	}
    202 
    203 	while (len >= 128) {
    204 		/* Process full chunks directly from the source. */
    205 		if (alignment_ok(data, sizeof(uint64_t)))
    206 			Transform(ctx->s, (const uint64_t *)data);
    207 		else {
    208 			memcpy(ctx->buf.u8, data, sizeof(ctx->buf));
    209 			Transform(ctx->s, ctx->buf.u64);
    210 		}
    211 		ctx->bytes += 128;
    212 		data += 128;
    213 		len -= 128;
    214 	}
    215 
    216 	if (len) {
    217 		/* Fill the buffer with what remains. */
    218 		memcpy(ctx->buf.u8 + bufsize, data, len);
    219 		ctx->bytes += len;
    220 	}
    221 }
    222 
    223 void sha512_init(struct sha512_ctx *ctx)
    224 {
    225 	struct sha512_ctx init = SHA512_INIT;
    226 	*ctx = init;
    227 }
    228 
    229 void sha512_update(struct sha512_ctx *ctx, const void *p, size_t size)
    230 {
    231 	check_sha512(ctx);
    232 	add(ctx, p, size);
    233 }
    234 
    235 void sha512_done(struct sha512_ctx *ctx, struct sha512 *res)
    236 {
    237 	static const unsigned char pad[128] = { 0x80 };
    238 	uint64_t sizedesc[2] = { 0, 0 };
    239 	size_t i;
    240 
    241 	sizedesc[1] = cpu_to_be64((uint64_t)ctx->bytes << 3);
    242 
    243 	/* Add '1' bit to terminate, then all 0 bits, up to next block - 16. */
    244 	add(ctx, pad, 1 + ((256 - 16 - (ctx->bytes % 128) - 1) % 128));
    245 	/* Add number of bits of data (big endian) */
    246 	add(ctx, sizedesc, sizeof(sizedesc));
    247 	for (i = 0; i < sizeof(ctx->s) / sizeof(ctx->s[0]); i++)
    248 		res->u.u64[i] = cpu_to_be64(ctx->s[i]);
    249 	invalidate_sha512(ctx);
    250 }
    251 #endif /* CCAN_CRYPTO_SHA512_USE_OPENSSL */
    252 
    253 void sha512(struct sha512 *sha, const void *p, size_t size)
    254 {
    255 	struct sha512_ctx ctx;
    256 
    257 	sha512_init(&ctx);
    258 	sha512_update(&ctx, p, size);
    259 	sha512_done(&ctx, sha);
    260 	memset(&ctx, 0, sizeof(ctx));
    261 }