/*
 * Copyright (C) 2015 Michael Brown <mbrown@fensystems.co.uk>.
 *
 * 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 any later version.
 *
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
 * 02110-1301, USA.
 *
 * You can also choose to distribute this program under the terms of
 * the Unmodified Binary Distribution Licence (as given in the file
 * COPYING.UBDL), provided that you have satisfied its requirements.
 */

FILE_LICENCE ( GPL2_OR_LATER_OR_UBDL );

/** @file
 *
 * SHA-512 algorithm
 *
 */

#include <stdint.h>
#include <string.h>
#include <byteswap.h>
#include <assert.h>
#include <ipxe/rotate.h>
#include <ipxe/crypto.h>
#include <ipxe/sha512.h>

/** SHA-512 variables */
struct sha512_variables {
	/* This layout matches that of struct sha512_digest_data,
	 * allowing for efficient endianness-conversion,
	 */
	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 w[SHA512_ROUNDS];
} __attribute__ (( packed ));

/** SHA-512 constants */
static const uint64_t k[SHA512_ROUNDS] = {
	0x428a2f98d728ae22ULL, 0x7137449123ef65cdULL, 0xb5c0fbcfec4d3b2fULL,
	0xe9b5dba58189dbbcULL, 0x3956c25bf348b538ULL, 0x59f111f1b605d019ULL,
	0x923f82a4af194f9bULL, 0xab1c5ed5da6d8118ULL, 0xd807aa98a3030242ULL,
	0x12835b0145706fbeULL, 0x243185be4ee4b28cULL, 0x550c7dc3d5ffb4e2ULL,
	0x72be5d74f27b896fULL, 0x80deb1fe3b1696b1ULL, 0x9bdc06a725c71235ULL,
	0xc19bf174cf692694ULL, 0xe49b69c19ef14ad2ULL, 0xefbe4786384f25e3ULL,
	0x0fc19dc68b8cd5b5ULL, 0x240ca1cc77ac9c65ULL, 0x2de92c6f592b0275ULL,
	0x4a7484aa6ea6e483ULL, 0x5cb0a9dcbd41fbd4ULL, 0x76f988da831153b5ULL,
	0x983e5152ee66dfabULL, 0xa831c66d2db43210ULL, 0xb00327c898fb213fULL,
	0xbf597fc7beef0ee4ULL, 0xc6e00bf33da88fc2ULL, 0xd5a79147930aa725ULL,
	0x06ca6351e003826fULL, 0x142929670a0e6e70ULL, 0x27b70a8546d22ffcULL,
	0x2e1b21385c26c926ULL, 0x4d2c6dfc5ac42aedULL, 0x53380d139d95b3dfULL,
	0x650a73548baf63deULL, 0x766a0abb3c77b2a8ULL, 0x81c2c92e47edaee6ULL,
	0x92722c851482353bULL, 0xa2bfe8a14cf10364ULL, 0xa81a664bbc423001ULL,
	0xc24b8b70d0f89791ULL, 0xc76c51a30654be30ULL, 0xd192e819d6ef5218ULL,
	0xd69906245565a910ULL, 0xf40e35855771202aULL, 0x106aa07032bbd1b8ULL,
	0x19a4c116b8d2d0c8ULL, 0x1e376c085141ab53ULL, 0x2748774cdf8eeb99ULL,
	0x34b0bcb5e19b48a8ULL, 0x391c0cb3c5c95a63ULL, 0x4ed8aa4ae3418acbULL,
	0x5b9cca4f7763e373ULL, 0x682e6ff3d6b2b8a3ULL, 0x748f82ee5defb2fcULL,
	0x78a5636f43172f60ULL, 0x84c87814a1f0ab72ULL, 0x8cc702081a6439ecULL,
	0x90befffa23631e28ULL, 0xa4506cebde82bde9ULL, 0xbef9a3f7b2c67915ULL,
	0xc67178f2e372532bULL, 0xca273eceea26619cULL, 0xd186b8c721c0c207ULL,
	0xeada7dd6cde0eb1eULL, 0xf57d4f7fee6ed178ULL, 0x06f067aa72176fbaULL,
	0x0a637dc5a2c898a6ULL, 0x113f9804bef90daeULL, 0x1b710b35131c471bULL,
	0x28db77f523047d84ULL, 0x32caab7b40c72493ULL, 0x3c9ebe0a15c9bebcULL,
	0x431d67c49c100d4cULL, 0x4cc5d4becb3e42b6ULL, 0x597f299cfc657e2aULL,
	0x5fcb6fab3ad6faecULL, 0x6c44198c4a475817ULL
};

/** SHA-512 initial digest values */
static const struct sha512_digest sha512_init_digest = {
	.h = {
		cpu_to_be64 ( 0x6a09e667f3bcc908ULL ),
		cpu_to_be64 ( 0xbb67ae8584caa73bULL ),
		cpu_to_be64 ( 0x3c6ef372fe94f82bULL ),
		cpu_to_be64 ( 0xa54ff53a5f1d36f1ULL ),
		cpu_to_be64 ( 0x510e527fade682d1ULL ),
		cpu_to_be64 ( 0x9b05688c2b3e6c1fULL ),
		cpu_to_be64 ( 0x1f83d9abfb41bd6bULL ),
		cpu_to_be64 ( 0x5be0cd19137e2179ULL ),
	},
};

/**
 * Initialise SHA-512 family algorithm
 *
 * @v context		SHA-512 context
 * @v init		Initial digest values
 * @v digestsize	Digest size
 */
void sha512_family_init ( struct sha512_context *context,
			  const struct sha512_digest *init,
			  size_t digestsize ) {

	context->len = 0;
	context->digestsize = digestsize;
	memcpy ( &context->ddq.dd.digest, init,
		 sizeof ( context->ddq.dd.digest ) );
}

/**
 * Initialise SHA-512 algorithm
 *
 * @v ctx		SHA-512 context
 */
static void sha512_init ( void *ctx ) {
	struct sha512_context *context = ctx;

	sha512_family_init ( context, &sha512_init_digest,
			     sizeof ( struct sha512_digest ) );
}

/**
 * Calculate SHA-512 digest of accumulated data
 *
 * @v context		SHA-512 context
 */
static void sha512_digest ( struct sha512_context *context ) {
        union {
		union sha512_digest_data_qwords ddq;
		struct sha512_variables v;
	} u;
	uint64_t *a = &u.v.a;
	uint64_t *b = &u.v.b;
	uint64_t *c = &u.v.c;
	uint64_t *d = &u.v.d;
	uint64_t *e = &u.v.e;
	uint64_t *f = &u.v.f;
	uint64_t *g = &u.v.g;
	uint64_t *h = &u.v.h;
	uint64_t *w = u.v.w;
	uint64_t s0;
	uint64_t s1;
	uint64_t maj;
	uint64_t t1;
	uint64_t t2;
	uint64_t ch;
	unsigned int i;

	/* Sanity checks */
	assert ( ( context->len % sizeof ( context->ddq.dd.data ) ) == 0 );
	build_assert ( &u.ddq.dd.digest.h[0] == a );
	build_assert ( &u.ddq.dd.digest.h[1] == b );
	build_assert ( &u.ddq.dd.digest.h[2] == c );
	build_assert ( &u.ddq.dd.digest.h[3] == d );
	build_assert ( &u.ddq.dd.digest.h[4] == e );
	build_assert ( &u.ddq.dd.digest.h[5] == f );
	build_assert ( &u.ddq.dd.digest.h[6] == g );
	build_assert ( &u.ddq.dd.digest.h[7] == h );
	build_assert ( &u.ddq.dd.data.qword[0] == w );

	DBGC ( context, "SHA512 digesting:\n" );
	DBGC_HDA ( context, 0, &context->ddq.dd.digest,
		   sizeof ( context->ddq.dd.digest ) );
	DBGC_HDA ( context, context->len, &context->ddq.dd.data,
		   sizeof ( context->ddq.dd.data ) );

	/* Convert h[0..7] to host-endian, and initialise a, b, c, d,
	 * e, f, g, h, and w[0..15]
	 */
	for ( i = 0 ; i < ( sizeof ( u.ddq.qword ) /
			    sizeof ( u.ddq.qword[0] ) ) ; i++ ) {
		be64_to_cpus ( &context->ddq.qword[i] );
		u.ddq.qword[i] = context->ddq.qword[i];
	}

	/* Initialise w[16..79] */
	for ( i = 16 ; i < SHA512_ROUNDS ; i++ ) {
		s0 = ( ror64 ( w[i-15], 1 ) ^ ror64 ( w[i-15], 8 ) ^
		       ( w[i-15] >> 7 ) );
		s1 = ( ror64 ( w[i-2], 19 ) ^ ror64 ( w[i-2], 61 ) ^
		       ( w[i-2] >> 6 ) );
		w[i] = ( w[i-16] + s0 + w[i-7] + s1 );
	}

	/* Main loop */
	for ( i = 0 ; i < SHA512_ROUNDS ; i++ ) {
		s0 = ( ror64 ( *a, 28 ) ^ ror64 ( *a, 34 ) ^ ror64 ( *a, 39 ) );
		maj = ( ( *a & *b ) ^ ( *a & *c ) ^ ( *b & *c ) );
		t2 = ( s0 + maj );
		s1 = ( ror64 ( *e, 14 ) ^ ror64 ( *e, 18 ) ^ ror64 ( *e, 41 ) );
		ch = ( ( *e & *f ) ^ ( (~*e) & *g ) );
		t1 = ( *h + s1 + ch + k[i] + w[i] );
		*h = *g;
		*g = *f;
		*f = *e;
		*e = ( *d + t1 );
		*d = *c;
		*c = *b;
		*b = *a;
		*a = ( t1 + t2 );
		DBGC2 ( context, "%2d : %016llx %016llx %016llx %016llx "
			"%016llx %016llx %016llx %016llx\n",
			i, *a, *b, *c, *d, *e, *f, *g, *h );
	}

	/* Add chunk to hash and convert back to big-endian */
	for ( i = 0 ; i < 8 ; i++ ) {
		context->ddq.dd.digest.h[i] =
			cpu_to_be64 ( context->ddq.dd.digest.h[i] +
				      u.ddq.dd.digest.h[i] );
	}

	DBGC ( context, "SHA512 digested:\n" );
	DBGC_HDA ( context, 0, &context->ddq.dd.digest,
		   sizeof ( context->ddq.dd.digest ) );
}

/**
 * Accumulate data with SHA-512 algorithm
 *
 * @v ctx		SHA-512 context
 * @v data		Data
 * @v len		Length of data
 */
void sha512_update ( void *ctx, const void *data, size_t len ) {
	struct sha512_context *context = ctx;
	const uint8_t *byte = data;
	size_t offset;

	/* Accumulate data a byte at a time, performing the digest
	 * whenever we fill the data buffer
	 */
	while ( len-- ) {
		offset = ( context->len % sizeof ( context->ddq.dd.data ) );
		context->ddq.dd.data.byte[offset] = *(byte++);
		context->len++;
		if ( ( context->len % sizeof ( context->ddq.dd.data ) ) == 0 )
			sha512_digest ( context );
	}
}

/**
 * Generate SHA-512 digest
 *
 * @v ctx		SHA-512 context
 * @v out		Output buffer
 */
void sha512_final ( void *ctx, void *out ) {
	struct sha512_context *context = ctx;
	uint64_t len_bits_hi;
	uint64_t len_bits_lo;
	uint8_t pad;

	/* Record length before pre-processing */
	len_bits_hi = 0;
	len_bits_lo = cpu_to_be64 ( ( ( uint64_t ) context->len ) * 8 );

	/* Pad with a single "1" bit followed by as many "0" bits as required */
	pad = 0x80;
	do {
		sha512_update ( ctx, &pad, sizeof ( pad ) );
		pad = 0x00;
	} while ( ( context->len % sizeof ( context->ddq.dd.data ) ) !=
		  offsetof ( typeof ( context->ddq.dd.data ), final.len_hi ) );

	/* Append length (in bits) */
	sha512_update ( ctx, &len_bits_hi, sizeof ( len_bits_hi ) );
	sha512_update ( ctx, &len_bits_lo, sizeof ( len_bits_lo ) );
	assert ( ( context->len % sizeof ( context->ddq.dd.data ) ) == 0 );

	/* Copy out final digest */
	memcpy ( out, &context->ddq.dd.digest, context->digestsize );
}

/** SHA-512 algorithm */
struct digest_algorithm sha512_algorithm = {
	.name		= "sha512",
	.ctxsize	= sizeof ( struct sha512_context ),
	.blocksize	= sizeof ( union sha512_block ),
	.digestsize	= sizeof ( struct sha512_digest ),
	.init		= sha512_init,
	.update		= sha512_update,
	.final		= sha512_final,
};