[fio.git] / sha256.c

/*
 * Cryptographic API.
 *
 * SHA-256, as specified in
 * http://csrc.nist.gov/cryptval/shs/sha256-384-512.pdf
 *
 * SHA-256 code by Jean-Luc Cooke <jlcooke@certainkey.com>.
 *
 * Copyright (c) Jean-Luc Cooke <jlcooke@certainkey.com>
 * Copyright (c) Andrew McDonald <andrew@mcdonald.org.uk>
 * Copyright (c) 2002 James Morris <jmorris@intercode.com.au>
 *
 * This program is free software; you can redistribute it and/or modify it
 * under the terms of the GNU General Public License as published by the Free
 * Software Foundation; either version 2 of the License, or (at your option) 
 * any later version.
 *
 */
#include <string.h>
#include <inttypes.h>
#include <byteswap.h>
#include <endian.h>

#include "sha256.h"

#if __BYTE_ORDER == __LITTLE_ENDIAN
#define	__be32_to_cpu(x)	__bswap_32(x)
#else
#define __be32_to_cpu(x)	(x)
#endif

#define SHA256_DIGEST_SIZE	32
#define SHA256_HMAC_BLOCK_SIZE	64

static inline uint32_t Ch(uint32_t x, uint32_t y, uint32_t z)
{
	return z ^ (x & (y ^ z));
}

static inline uint32_t Maj(uint32_t x, uint32_t y, uint32_t z)
{
	return (x & y) | (z & (x | y));
}

#define e0(x)       (ror32(x, 2) ^ ror32(x,13) ^ ror32(x,22))
#define e1(x)       (ror32(x, 6) ^ ror32(x,11) ^ ror32(x,25))
#define s0(x)       (ror32(x, 7) ^ ror32(x,18) ^ (x >> 3))
#define s1(x)       (ror32(x,17) ^ ror32(x,19) ^ (x >> 10))

#define H0         0x6a09e667
#define H1         0xbb67ae85
#define H2         0x3c6ef372
#define H3         0xa54ff53a
#define H4         0x510e527f
#define H5         0x9b05688c
#define H6         0x1f83d9ab
#define H7         0x5be0cd19

static inline uint32_t ror32(uint32_t word, unsigned int shift)
{
	 return (word >> shift) | (word << (32 - shift));
}

static inline void LOAD_OP(int I, uint32_t *W, const uint8_t *input)
{
	W[I] = __be32_to_cpu(((uint32_t *)(input))[I]);
}

static inline void BLEND_OP(int I, uint32_t *W)
{
	W[I] = s1(W[I-2]) + W[I-7] + s0(W[I-15]) + W[I-16];
}

static void sha256_transform(uint32_t *state, const uint8_t *input)
{
	uint32_t a, b, c, d, e, f, g, h, t1, t2;
	uint32_t W[64];
	int i;

	/* load the input */
	for (i = 0; i < 16; i++)
		LOAD_OP(i, W, input);

	/* now blend */
	for (i = 16; i < 64; i++)
		BLEND_OP(i, W);

	/* load the state into our registers */
	a=state[0];  b=state[1];  c=state[2];  d=state[3];
	e=state[4];  f=state[5];  g=state[6];  h=state[7];

	/* now iterate */
	t1 = h + e1(e) + Ch(e,f,g) + 0x428a2f98 + W[ 0];
	t2 = e0(a) + Maj(a,b,c);    d+=t1;    h=t1+t2;
	t1 = g + e1(d) + Ch(d,e,f) + 0x71374491 + W[ 1];
	t2 = e0(h) + Maj(h,a,b);    c+=t1;    g=t1+t2;
	t1 = f + e1(c) + Ch(c,d,e) + 0xb5c0fbcf + W[ 2];
	t2 = e0(g) + Maj(g,h,a);    b+=t1;    f=t1+t2;
	t1 = e + e1(b) + Ch(b,c,d) + 0xe9b5dba5 + W[ 3];
	t2 = e0(f) + Maj(f,g,h);    a+=t1;    e=t1+t2;
	t1 = d + e1(a) + Ch(a,b,c) + 0x3956c25b + W[ 4];
	t2 = e0(e) + Maj(e,f,g);    h+=t1;    d=t1+t2;
	t1 = c + e1(h) + Ch(h,a,b) + 0x59f111f1 + W[ 5];
	t2 = e0(d) + Maj(d,e,f);    g+=t1;    c=t1+t2;
	t1 = b + e1(g) + Ch(g,h,a) + 0x923f82a4 + W[ 6];
	t2 = e0(c) + Maj(c,d,e);    f+=t1;    b=t1+t2;
	t1 = a + e1(f) + Ch(f,g,h) + 0xab1c5ed5 + W[ 7];
	t2 = e0(b) + Maj(b,c,d);    e+=t1;    a=t1+t2;

	t1 = h + e1(e) + Ch(e,f,g) + 0xd807aa98 + W[ 8];
	t2 = e0(a) + Maj(a,b,c);    d+=t1;    h=t1+t2;
	t1 = g + e1(d) + Ch(d,e,f) + 0x12835b01 + W[ 9];
	t2 = e0(h) + Maj(h,a,b);    c+=t1;    g=t1+t2;
	t1 = f + e1(c) + Ch(c,d,e) + 0x243185be + W[10];
	t2 = e0(g) + Maj(g,h,a);    b+=t1;    f=t1+t2;
	t1 = e + e1(b) + Ch(b,c,d) + 0x550c7dc3 + W[11];
	t2 = e0(f) + Maj(f,g,h);    a+=t1;    e=t1+t2;
	t1 = d + e1(a) + Ch(a,b,c) + 0x72be5d74 + W[12];
	t2 = e0(e) + Maj(e,f,g);    h+=t1;    d=t1+t2;
	t1 = c + e1(h) + Ch(h,a,b) + 0x80deb1fe + W[13];
	t2 = e0(d) + Maj(d,e,f);    g+=t1;    c=t1+t2;
	t1 = b + e1(g) + Ch(g,h,a) + 0x9bdc06a7 + W[14];
	t2 = e0(c) + Maj(c,d,e);    f+=t1;    b=t1+t2;
	t1 = a + e1(f) + Ch(f,g,h) + 0xc19bf174 + W[15];
	t2 = e0(b) + Maj(b,c,d);    e+=t1;    a=t1+t2;

	t1 = h + e1(e) + Ch(e,f,g) + 0xe49b69c1 + W[16];
	t2 = e0(a) + Maj(a,b,c);    d+=t1;    h=t1+t2;
	t1 = g + e1(d) + Ch(d,e,f) + 0xefbe4786 + W[17];
	t2 = e0(h) + Maj(h,a,b);    c+=t1;    g=t1+t2;
	t1 = f + e1(c) + Ch(c,d,e) + 0x0fc19dc6 + W[18];
	t2 = e0(g) + Maj(g,h,a);    b+=t1;    f=t1+t2;
	t1 = e + e1(b) + Ch(b,c,d) + 0x240ca1cc + W[19];
	t2 = e0(f) + Maj(f,g,h);    a+=t1;    e=t1+t2;
	t1 = d + e1(a) + Ch(a,b,c) + 0x2de92c6f + W[20];
	t2 = e0(e) + Maj(e,f,g);    h+=t1;    d=t1+t2;
	t1 = c + e1(h) + Ch(h,a,b) + 0x4a7484aa + W[21];
	t2 = e0(d) + Maj(d,e,f);    g+=t1;    c=t1+t2;
	t1 = b + e1(g) + Ch(g,h,a) + 0x5cb0a9dc + W[22];
	t2 = e0(c) + Maj(c,d,e);    f+=t1;    b=t1+t2;
	t1 = a + e1(f) + Ch(f,g,h) + 0x76f988da + W[23];
	t2 = e0(b) + Maj(b,c,d);    e+=t1;    a=t1+t2;

	t1 = h + e1(e) + Ch(e,f,g) + 0x983e5152 + W[24];
	t2 = e0(a) + Maj(a,b,c);    d+=t1;    h=t1+t2;
	t1 = g + e1(d) + Ch(d,e,f) + 0xa831c66d + W[25];
	t2 = e0(h) + Maj(h,a,b);    c+=t1;    g=t1+t2;
	t1 = f + e1(c) + Ch(c,d,e) + 0xb00327c8 + W[26];
	t2 = e0(g) + Maj(g,h,a);    b+=t1;    f=t1+t2;
	t1 = e + e1(b) + Ch(b,c,d) + 0xbf597fc7 + W[27];
	t2 = e0(f) + Maj(f,g,h);    a+=t1;    e=t1+t2;
	t1 = d + e1(a) + Ch(a,b,c) + 0xc6e00bf3 + W[28];
	t2 = e0(e) + Maj(e,f,g);    h+=t1;    d=t1+t2;
	t1 = c + e1(h) + Ch(h,a,b) + 0xd5a79147 + W[29];
	t2 = e0(d) + Maj(d,e,f);    g+=t1;    c=t1+t2;
	t1 = b + e1(g) + Ch(g,h,a) + 0x06ca6351 + W[30];
	t2 = e0(c) + Maj(c,d,e);    f+=t1;    b=t1+t2;
	t1 = a + e1(f) + Ch(f,g,h) + 0x14292967 + W[31];
	t2 = e0(b) + Maj(b,c,d);    e+=t1;    a=t1+t2;

	t1 = h + e1(e) + Ch(e,f,g) + 0x27b70a85 + W[32];
	t2 = e0(a) + Maj(a,b,c);    d+=t1;    h=t1+t2;
	t1 = g + e1(d) + Ch(d,e,f) + 0x2e1b2138 + W[33];
	t2 = e0(h) + Maj(h,a,b);    c+=t1;    g=t1+t2;
	t1 = f + e1(c) + Ch(c,d,e) + 0x4d2c6dfc + W[34];
	t2 = e0(g) + Maj(g,h,a);    b+=t1;    f=t1+t2;
	t1 = e + e1(b) + Ch(b,c,d) + 0x53380d13 + W[35];
	t2 = e0(f) + Maj(f,g,h);    a+=t1;    e=t1+t2;
	t1 = d + e1(a) + Ch(a,b,c) + 0x650a7354 + W[36];
	t2 = e0(e) + Maj(e,f,g);    h+=t1;    d=t1+t2;
	t1 = c + e1(h) + Ch(h,a,b) + 0x766a0abb + W[37];
	t2 = e0(d) + Maj(d,e,f);    g+=t1;    c=t1+t2;
	t1 = b + e1(g) + Ch(g,h,a) + 0x81c2c92e + W[38];
	t2 = e0(c) + Maj(c,d,e);    f+=t1;    b=t1+t2;
	t1 = a + e1(f) + Ch(f,g,h) + 0x92722c85 + W[39];
	t2 = e0(b) + Maj(b,c,d);    e+=t1;    a=t1+t2;

	t1 = h + e1(e) + Ch(e,f,g) + 0xa2bfe8a1 + W[40];
	t2 = e0(a) + Maj(a,b,c);    d+=t1;    h=t1+t2;
	t1 = g + e1(d) + Ch(d,e,f) + 0xa81a664b + W[41];
	t2 = e0(h) + Maj(h,a,b);    c+=t1;    g=t1+t2;
	t1 = f + e1(c) + Ch(c,d,e) + 0xc24b8b70 + W[42];
	t2 = e0(g) + Maj(g,h,a);    b+=t1;    f=t1+t2;
	t1 = e + e1(b) + Ch(b,c,d) + 0xc76c51a3 + W[43];
	t2 = e0(f) + Maj(f,g,h);    a+=t1;    e=t1+t2;
	t1 = d + e1(a) + Ch(a,b,c) + 0xd192e819 + W[44];
	t2 = e0(e) + Maj(e,f,g);    h+=t1;    d=t1+t2;
	t1 = c + e1(h) + Ch(h,a,b) + 0xd6990624 + W[45];
	t2 = e0(d) + Maj(d,e,f);    g+=t1;    c=t1+t2;
	t1 = b + e1(g) + Ch(g,h,a) + 0xf40e3585 + W[46];
	t2 = e0(c) + Maj(c,d,e);    f+=t1;    b=t1+t2;
	t1 = a + e1(f) + Ch(f,g,h) + 0x106aa070 + W[47];
	t2 = e0(b) + Maj(b,c,d);    e+=t1;    a=t1+t2;

	t1 = h + e1(e) + Ch(e,f,g) + 0x19a4c116 + W[48];
	t2 = e0(a) + Maj(a,b,c);    d+=t1;    h=t1+t2;
	t1 = g + e1(d) + Ch(d,e,f) + 0x1e376c08 + W[49];
	t2 = e0(h) + Maj(h,a,b);    c+=t1;    g=t1+t2;
	t1 = f + e1(c) + Ch(c,d,e) + 0x2748774c + W[50];
	t2 = e0(g) + Maj(g,h,a);    b+=t1;    f=t1+t2;
	t1 = e + e1(b) + Ch(b,c,d) + 0x34b0bcb5 + W[51];
	t2 = e0(f) + Maj(f,g,h);    a+=t1;    e=t1+t2;
	t1 = d + e1(a) + Ch(a,b,c) + 0x391c0cb3 + W[52];
	t2 = e0(e) + Maj(e,f,g);    h+=t1;    d=t1+t2;
	t1 = c + e1(h) + Ch(h,a,b) + 0x4ed8aa4a + W[53];
	t2 = e0(d) + Maj(d,e,f);    g+=t1;    c=t1+t2;
	t1 = b + e1(g) + Ch(g,h,a) + 0x5b9cca4f + W[54];
	t2 = e0(c) + Maj(c,d,e);    f+=t1;    b=t1+t2;
	t1 = a + e1(f) + Ch(f,g,h) + 0x682e6ff3 + W[55];
	t2 = e0(b) + Maj(b,c,d);    e+=t1;    a=t1+t2;

	t1 = h + e1(e) + Ch(e,f,g) + 0x748f82ee + W[56];
	t2 = e0(a) + Maj(a,b,c);    d+=t1;    h=t1+t2;
	t1 = g + e1(d) + Ch(d,e,f) + 0x78a5636f + W[57];
	t2 = e0(h) + Maj(h,a,b);    c+=t1;    g=t1+t2;
	t1 = f + e1(c) + Ch(c,d,e) + 0x84c87814 + W[58];
	t2 = e0(g) + Maj(g,h,a);    b+=t1;    f=t1+t2;
	t1 = e + e1(b) + Ch(b,c,d) + 0x8cc70208 + W[59];
	t2 = e0(f) + Maj(f,g,h);    a+=t1;    e=t1+t2;
	t1 = d + e1(a) + Ch(a,b,c) + 0x90befffa + W[60];
	t2 = e0(e) + Maj(e,f,g);    h+=t1;    d=t1+t2;
	t1 = c + e1(h) + Ch(h,a,b) + 0xa4506ceb + W[61];
	t2 = e0(d) + Maj(d,e,f);    g+=t1;    c=t1+t2;
	t1 = b + e1(g) + Ch(g,h,a) + 0xbef9a3f7 + W[62];
	t2 = e0(c) + Maj(c,d,e);    f+=t1;    b=t1+t2;
	t1 = a + e1(f) + Ch(f,g,h) + 0xc67178f2 + W[63];
	t2 = e0(b) + Maj(b,c,d);    e+=t1;    a=t1+t2;

	state[0] += a; state[1] += b; state[2] += c; state[3] += d;
	state[4] += e; state[5] += f; state[6] += g; state[7] += h;

	/* clear any sensitive info... */
	a = b = c = d = e = f = g = h = t1 = t2 = 0;
	memset(W, 0, 64 * sizeof(uint32_t));
}

void sha256_init(struct sha256_ctx *sctx)
{
	sctx->state[0] = H0;
	sctx->state[1] = H1;
	sctx->state[2] = H2;
	sctx->state[3] = H3;
	sctx->state[4] = H4;
	sctx->state[5] = H5;
	sctx->state[6] = H6;
	sctx->state[7] = H7;
	sctx->count[0] = sctx->count[1] = 0;
}

void sha256_update(struct sha256_ctx *sctx, const uint8_t *data,
		   unsigned int len)
{
	unsigned int i, index, part_len;

	/* Compute number of bytes mod 128 */
	index = (unsigned int)((sctx->count[0] >> 3) & 0x3f);

	/* Update number of bits */
	if ((sctx->count[0] += (len << 3)) < (len << 3)) {
		sctx->count[1]++;
		sctx->count[1] += (len >> 29);
	}

	part_len = 64 - index;

	/* Transform as many times as possible. */
	if (len >= part_len) {
		memcpy(&sctx->buf[index], data, part_len);
		sha256_transform(sctx->state, sctx->buf);

		for (i = part_len; i + 63 < len; i += 64)
			sha256_transform(sctx->state, &data[i]);
		index = 0;
	} else {
		i = 0;
	}
	
	/* Buffer remaining input */
	memcpy(&sctx->buf[index], &data[i], len-i);
}
Commit	Line	Data
	1	/*
	2	* Cryptographic API.
	3	*
	4	* SHA-256, as specified in
	5	* http://csrc.nist.gov/cryptval/shs/sha256-384-512.pdf
	6	*
	7	* SHA-256 code by Jean-Luc Cooke <jlcooke@certainkey.com>.
	8	*
	9	* Copyright (c) Jean-Luc Cooke <jlcooke@certainkey.com>
	10	* Copyright (c) Andrew McDonald <andrew@mcdonald.org.uk>
	11	* Copyright (c) 2002 James Morris <jmorris@intercode.com.au>
	12	*
	13	* This program is free software; you can redistribute it and/or modify it
	14	* under the terms of the GNU General Public License as published by the Free
	15	* Software Foundation; either version 2 of the License, or (at your option)
	16	* any later version.
	17	*
	18	*/
	19	#include <string.h>
	20	#include <inttypes.h>
	21	#include <byteswap.h>
	22	#include <endian.h>
	23
	24	#include "sha256.h"
	25
	26	#if __BYTE_ORDER == __LITTLE_ENDIAN
	27	#define __be32_to_cpu(x) __bswap_32(x)
	28	#else
	29	#define __be32_to_cpu(x) (x)
	30	#endif
	31
	32	#define SHA256_DIGEST_SIZE 32
	33	#define SHA256_HMAC_BLOCK_SIZE 64
	34
	35	static inline uint32_t Ch(uint32_t x, uint32_t y, uint32_t z)
	36	{
	37	return z ^ (x & (y ^ z));
	38	}
	39
	40	static inline uint32_t Maj(uint32_t x, uint32_t y, uint32_t z)
	41	{
	42	return (x & y) \| (z & (x \| y));
	43	}
	44
	45	#define e0(x) (ror32(x, 2) ^ ror32(x,13) ^ ror32(x,22))
	46	#define e1(x) (ror32(x, 6) ^ ror32(x,11) ^ ror32(x,25))
	47	#define s0(x) (ror32(x, 7) ^ ror32(x,18) ^ (x >> 3))
	48	#define s1(x) (ror32(x,17) ^ ror32(x,19) ^ (x >> 10))
	49
	50	#define H0 0x6a09e667
	51	#define H1 0xbb67ae85
	52	#define H2 0x3c6ef372
	53	#define H3 0xa54ff53a
	54	#define H4 0x510e527f
	55	#define H5 0x9b05688c
	56	#define H6 0x1f83d9ab
	57	#define H7 0x5be0cd19
	58
	59	static inline uint32_t ror32(uint32_t word, unsigned int shift)
	60	{
	61	return (word >> shift) \| (word << (32 - shift));
	62	}
	63
	64	static inline void LOAD_OP(int I, uint32_t W, const uint8_t input)
	65	{
	66	W[I] = __be32_to_cpu(((uint32_t *)(input))[I]);
	67	}
	68
	69	static inline void BLEND_OP(int I, uint32_t *W)
	70	{
	71	W[I] = s1(W[I-2]) + W[I-7] + s0(W[I-15]) + W[I-16];
	72	}
	73
	74	static void sha256_transform(uint32_t state, const uint8_t input)
	75	{
	76	uint32_t a, b, c, d, e, f, g, h, t1, t2;
	77	uint32_t W[64];
	78	int i;
	79
	80	/* load the input */
	81	for (i = 0; i < 16; i++)
	82	LOAD_OP(i, W, input);
	83
	84	/* now blend */
	85	for (i = 16; i < 64; i++)
	86	BLEND_OP(i, W);
	87
	88	/* load the state into our registers */
	89	a=state[0]; b=state[1]; c=state[2]; d=state[3];
	90	e=state[4]; f=state[5]; g=state[6]; h=state[7];
	91
	92	/* now iterate */
	93	t1 = h + e1(e) + Ch(e,f,g) + 0x428a2f98 + W[ 0];
	94	t2 = e0(a) + Maj(a,b,c); d+=t1; h=t1+t2;
	95	t1 = g + e1(d) + Ch(d,e,f) + 0x71374491 + W[ 1];
	96	t2 = e0(h) + Maj(h,a,b); c+=t1; g=t1+t2;
	97	t1 = f + e1(c) + Ch(c,d,e) + 0xb5c0fbcf + W[ 2];
	98	t2 = e0(g) + Maj(g,h,a); b+=t1; f=t1+t2;
	99	t1 = e + e1(b) + Ch(b,c,d) + 0xe9b5dba5 + W[ 3];
	100	t2 = e0(f) + Maj(f,g,h); a+=t1; e=t1+t2;
	101	t1 = d + e1(a) + Ch(a,b,c) + 0x3956c25b + W[ 4];
	102	t2 = e0(e) + Maj(e,f,g); h+=t1; d=t1+t2;
	103	t1 = c + e1(h) + Ch(h,a,b) + 0x59f111f1 + W[ 5];
	104	t2 = e0(d) + Maj(d,e,f); g+=t1; c=t1+t2;
	105	t1 = b + e1(g) + Ch(g,h,a) + 0x923f82a4 + W[ 6];
	106	t2 = e0(c) + Maj(c,d,e); f+=t1; b=t1+t2;
	107	t1 = a + e1(f) + Ch(f,g,h) + 0xab1c5ed5 + W[ 7];
	108	t2 = e0(b) + Maj(b,c,d); e+=t1; a=t1+t2;
	109
	110	t1 = h + e1(e) + Ch(e,f,g) + 0xd807aa98 + W[ 8];
	111	t2 = e0(a) + Maj(a,b,c); d+=t1; h=t1+t2;
	112	t1 = g + e1(d) + Ch(d,e,f) + 0x12835b01 + W[ 9];
	113	t2 = e0(h) + Maj(h,a,b); c+=t1; g=t1+t2;
	114	t1 = f + e1(c) + Ch(c,d,e) + 0x243185be + W[10];
	115	t2 = e0(g) + Maj(g,h,a); b+=t1; f=t1+t2;
	116	t1 = e + e1(b) + Ch(b,c,d) + 0x550c7dc3 + W[11];
	117	t2 = e0(f) + Maj(f,g,h); a+=t1; e=t1+t2;
	118	t1 = d + e1(a) + Ch(a,b,c) + 0x72be5d74 + W[12];
	119	t2 = e0(e) + Maj(e,f,g); h+=t1; d=t1+t2;
	120	t1 = c + e1(h) + Ch(h,a,b) + 0x80deb1fe + W[13];
	121	t2 = e0(d) + Maj(d,e,f); g+=t1; c=t1+t2;
	122	t1 = b + e1(g) + Ch(g,h,a) + 0x9bdc06a7 + W[14];
	123	t2 = e0(c) + Maj(c,d,e); f+=t1; b=t1+t2;
	124	t1 = a + e1(f) + Ch(f,g,h) + 0xc19bf174 + W[15];
	125	t2 = e0(b) + Maj(b,c,d); e+=t1; a=t1+t2;
	126
	127	t1 = h + e1(e) + Ch(e,f,g) + 0xe49b69c1 + W[16];
	128	t2 = e0(a) + Maj(a,b,c); d+=t1; h=t1+t2;
	129	t1 = g + e1(d) + Ch(d,e,f) + 0xefbe4786 + W[17];
	130	t2 = e0(h) + Maj(h,a,b); c+=t1; g=t1+t2;
	131	t1 = f + e1(c) + Ch(c,d,e) + 0x0fc19dc6 + W[18];
	132	t2 = e0(g) + Maj(g,h,a); b+=t1; f=t1+t2;
	133	t1 = e + e1(b) + Ch(b,c,d) + 0x240ca1cc + W[19];
	134	t2 = e0(f) + Maj(f,g,h); a+=t1; e=t1+t2;
	135	t1 = d + e1(a) + Ch(a,b,c) + 0x2de92c6f + W[20];
	136	t2 = e0(e) + Maj(e,f,g); h+=t1; d=t1+t2;
	137	t1 = c + e1(h) + Ch(h,a,b) + 0x4a7484aa + W[21];
	138	t2 = e0(d) + Maj(d,e,f); g+=t1; c=t1+t2;
	139	t1 = b + e1(g) + Ch(g,h,a) + 0x5cb0a9dc + W[22];
	140	t2 = e0(c) + Maj(c,d,e); f+=t1; b=t1+t2;
	141	t1 = a + e1(f) + Ch(f,g,h) + 0x76f988da + W[23];
	142	t2 = e0(b) + Maj(b,c,d); e+=t1; a=t1+t2;
	143
	144	t1 = h + e1(e) + Ch(e,f,g) + 0x983e5152 + W[24];
	145	t2 = e0(a) + Maj(a,b,c); d+=t1; h=t1+t2;
	146	t1 = g + e1(d) + Ch(d,e,f) + 0xa831c66d + W[25];
	147	t2 = e0(h) + Maj(h,a,b); c+=t1; g=t1+t2;
	148	t1 = f + e1(c) + Ch(c,d,e) + 0xb00327c8 + W[26];
	149	t2 = e0(g) + Maj(g,h,a); b+=t1; f=t1+t2;
	150	t1 = e + e1(b) + Ch(b,c,d) + 0xbf597fc7 + W[27];
	151	t2 = e0(f) + Maj(f,g,h); a+=t1; e=t1+t2;
	152	t1 = d + e1(a) + Ch(a,b,c) + 0xc6e00bf3 + W[28];
	153	t2 = e0(e) + Maj(e,f,g); h+=t1; d=t1+t2;
	154	t1 = c + e1(h) + Ch(h,a,b) + 0xd5a79147 + W[29];
	155	t2 = e0(d) + Maj(d,e,f); g+=t1; c=t1+t2;
	156	t1 = b + e1(g) + Ch(g,h,a) + 0x06ca6351 + W[30];
	157	t2 = e0(c) + Maj(c,d,e); f+=t1; b=t1+t2;
	158	t1 = a + e1(f) + Ch(f,g,h) + 0x14292967 + W[31];
	159	t2 = e0(b) + Maj(b,c,d); e+=t1; a=t1+t2;
	160
	161	t1 = h + e1(e) + Ch(e,f,g) + 0x27b70a85 + W[32];
	162	t2 = e0(a) + Maj(a,b,c); d+=t1; h=t1+t2;
	163	t1 = g + e1(d) + Ch(d,e,f) + 0x2e1b2138 + W[33];
	164	t2 = e0(h) + Maj(h,a,b); c+=t1; g=t1+t2;
	165	t1 = f + e1(c) + Ch(c,d,e) + 0x4d2c6dfc + W[34];
	166	t2 = e0(g) + Maj(g,h,a); b+=t1; f=t1+t2;
	167	t1 = e + e1(b) + Ch(b,c,d) + 0x53380d13 + W[35];
	168	t2 = e0(f) + Maj(f,g,h); a+=t1; e=t1+t2;
	169	t1 = d + e1(a) + Ch(a,b,c) + 0x650a7354 + W[36];
	170	t2 = e0(e) + Maj(e,f,g); h+=t1; d=t1+t2;
	171	t1 = c + e1(h) + Ch(h,a,b) + 0x766a0abb + W[37];
	172	t2 = e0(d) + Maj(d,e,f); g+=t1; c=t1+t2;
	173	t1 = b + e1(g) + Ch(g,h,a) + 0x81c2c92e + W[38];
	174	t2 = e0(c) + Maj(c,d,e); f+=t1; b=t1+t2;
	175	t1 = a + e1(f) + Ch(f,g,h) + 0x92722c85 + W[39];
	176	t2 = e0(b) + Maj(b,c,d); e+=t1; a=t1+t2;
	177
	178	t1 = h + e1(e) + Ch(e,f,g) + 0xa2bfe8a1 + W[40];
	179	t2 = e0(a) + Maj(a,b,c); d+=t1; h=t1+t2;
	180	t1 = g + e1(d) + Ch(d,e,f) + 0xa81a664b + W[41];
	181	t2 = e0(h) + Maj(h,a,b); c+=t1; g=t1+t2;
	182	t1 = f + e1(c) + Ch(c,d,e) + 0xc24b8b70 + W[42];
	183	t2 = e0(g) + Maj(g,h,a); b+=t1; f=t1+t2;
	184	t1 = e + e1(b) + Ch(b,c,d) + 0xc76c51a3 + W[43];
	185	t2 = e0(f) + Maj(f,g,h); a+=t1; e=t1+t2;
	186	t1 = d + e1(a) + Ch(a,b,c) + 0xd192e819 + W[44];
	187	t2 = e0(e) + Maj(e,f,g); h+=t1; d=t1+t2;
	188	t1 = c + e1(h) + Ch(h,a,b) + 0xd6990624 + W[45];
	189	t2 = e0(d) + Maj(d,e,f); g+=t1; c=t1+t2;
	190	t1 = b + e1(g) + Ch(g,h,a) + 0xf40e3585 + W[46];
	191	t2 = e0(c) + Maj(c,d,e); f+=t1; b=t1+t2;
	192	t1 = a + e1(f) + Ch(f,g,h) + 0x106aa070 + W[47];
	193	t2 = e0(b) + Maj(b,c,d); e+=t1; a=t1+t2;
	194
	195	t1 = h + e1(e) + Ch(e,f,g) + 0x19a4c116 + W[48];
	196	t2 = e0(a) + Maj(a,b,c); d+=t1; h=t1+t2;
	197	t1 = g + e1(d) + Ch(d,e,f) + 0x1e376c08 + W[49];
	198	t2 = e0(h) + Maj(h,a,b); c+=t1; g=t1+t2;
	199	t1 = f + e1(c) + Ch(c,d,e) + 0x2748774c + W[50];
	200	t2 = e0(g) + Maj(g,h,a); b+=t1; f=t1+t2;
	201	t1 = e + e1(b) + Ch(b,c,d) + 0x34b0bcb5 + W[51];
	202	t2 = e0(f) + Maj(f,g,h); a+=t1; e=t1+t2;
	203	t1 = d + e1(a) + Ch(a,b,c) + 0x391c0cb3 + W[52];
	204	t2 = e0(e) + Maj(e,f,g); h+=t1; d=t1+t2;
	205	t1 = c + e1(h) + Ch(h,a,b) + 0x4ed8aa4a + W[53];
	206	t2 = e0(d) + Maj(d,e,f); g+=t1; c=t1+t2;
	207	t1 = b + e1(g) + Ch(g,h,a) + 0x5b9cca4f + W[54];
	208	t2 = e0(c) + Maj(c,d,e); f+=t1; b=t1+t2;
	209	t1 = a + e1(f) + Ch(f,g,h) + 0x682e6ff3 + W[55];
	210	t2 = e0(b) + Maj(b,c,d); e+=t1; a=t1+t2;
	211
	212	t1 = h + e1(e) + Ch(e,f,g) + 0x748f82ee + W[56];
	213	t2 = e0(a) + Maj(a,b,c); d+=t1; h=t1+t2;
	214	t1 = g + e1(d) + Ch(d,e,f) + 0x78a5636f + W[57];
	215	t2 = e0(h) + Maj(h,a,b); c+=t1; g=t1+t2;
	216	t1 = f + e1(c) + Ch(c,d,e) + 0x84c87814 + W[58];
	217	t2 = e0(g) + Maj(g,h,a); b+=t1; f=t1+t2;
	218	t1 = e + e1(b) + Ch(b,c,d) + 0x8cc70208 + W[59];
	219	t2 = e0(f) + Maj(f,g,h); a+=t1; e=t1+t2;
	220	t1 = d + e1(a) + Ch(a,b,c) + 0x90befffa + W[60];
	221	t2 = e0(e) + Maj(e,f,g); h+=t1; d=t1+t2;
	222	t1 = c + e1(h) + Ch(h,a,b) + 0xa4506ceb + W[61];
	223	t2 = e0(d) + Maj(d,e,f); g+=t1; c=t1+t2;
	224	t1 = b + e1(g) + Ch(g,h,a) + 0xbef9a3f7 + W[62];
	225	t2 = e0(c) + Maj(c,d,e); f+=t1; b=t1+t2;
	226	t1 = a + e1(f) + Ch(f,g,h) + 0xc67178f2 + W[63];
	227	t2 = e0(b) + Maj(b,c,d); e+=t1; a=t1+t2;
	228
	229	state[0] += a; state[1] += b; state[2] += c; state[3] += d;
	230	state[4] += e; state[5] += f; state[6] += g; state[7] += h;
	231
	232	/* clear any sensitive info... */
	233	a = b = c = d = e = f = g = h = t1 = t2 = 0;
	234	memset(W, 0, 64 * sizeof(uint32_t));
	235	}
	236
	237	void sha256_init(struct sha256_ctx *sctx)
	238	{
	239	sctx->state[0] = H0;
	240	sctx->state[1] = H1;
	241	sctx->state[2] = H2;
	242	sctx->state[3] = H3;
	243	sctx->state[4] = H4;
	244	sctx->state[5] = H5;
	245	sctx->state[6] = H6;
	246	sctx->state[7] = H7;
	247	sctx->count[0] = sctx->count[1] = 0;
	248	}
	249
	250	void sha256_update(struct sha256_ctx sctx, const uint8_t data,
	251	unsigned int len)
	252	{
	253	unsigned int i, index, part_len;
	254
	255	/* Compute number of bytes mod 128 */
	256	index = (unsigned int)((sctx->count[0] >> 3) & 0x3f);
	257
	258	/* Update number of bits */
	259	if ((sctx->count[0] += (len << 3)) < (len << 3)) {
	260	sctx->count[1]++;
	261	sctx->count[1] += (len >> 29);
	262	}
	263
	264	part_len = 64 - index;
	265
	266	/* Transform as many times as possible. */
	267	if (len >= part_len) {
	268	memcpy(&sctx->buf[index], data, part_len);
	269	sha256_transform(sctx->state, sctx->buf);
	270
	271	for (i = part_len; i + 63 < len; i += 64)
	272	sha256_transform(sctx->state, &data[i]);
	273	index = 0;
	274	} else {
	275	i = 0;
	276	}
	277
	278	/* Buffer remaining input */
	279	memcpy(&sctx->buf[index], &data[i], len-i);
	280	}