Home | History | Annotate | Line # | Download | only in opencrypto
gmac.c revision 1.1
      1  1.1  drochner /* $NetBSD: gmac.c,v 1.1 2011/05/26 21:50:03 drochner Exp $ */
      2  1.1  drochner /* OpenBSD: gmac.c,v 1.3 2011/01/11 15:44:23 deraadt Exp */
      3  1.1  drochner 
      4  1.1  drochner /*
      5  1.1  drochner  * Copyright (c) 2010 Mike Belopuhov <mike (at) vantronix.net>
      6  1.1  drochner  *
      7  1.1  drochner  * Permission to use, copy, modify, and distribute this software for any
      8  1.1  drochner  * purpose with or without fee is hereby granted, provided that the above
      9  1.1  drochner  * copyright notice and this permission notice appear in all copies.
     10  1.1  drochner  *
     11  1.1  drochner  * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
     12  1.1  drochner  * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
     13  1.1  drochner  * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
     14  1.1  drochner  * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
     15  1.1  drochner  * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
     16  1.1  drochner  * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
     17  1.1  drochner  * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
     18  1.1  drochner  */
     19  1.1  drochner 
     20  1.1  drochner /*
     21  1.1  drochner  * This code implements the Message Authentication part of the
     22  1.1  drochner  * Galois/Counter Mode (as being described in the RFC 4543) using
     23  1.1  drochner  * the AES cipher.  FIPS SP 800-38D describes the algorithm details.
     24  1.1  drochner  */
     25  1.1  drochner 
     26  1.1  drochner #include <sys/param.h>
     27  1.1  drochner #include <sys/systm.h>
     28  1.1  drochner 
     29  1.1  drochner #include <crypto/rijndael/rijndael.h>
     30  1.1  drochner #include <opencrypto/gmac.h>
     31  1.1  drochner 
     32  1.1  drochner void	ghash_gfmul(const uint32_t *, const uint32_t *, uint32_t *);
     33  1.1  drochner void	ghash_update(GHASH_CTX *, const uint8_t *, size_t);
     34  1.1  drochner 
     35  1.1  drochner /* Computes a block multiplication in the GF(2^128) */
     36  1.1  drochner void
     37  1.1  drochner ghash_gfmul(const uint32_t *X, const uint32_t *Y, uint32_t *product)
     38  1.1  drochner {
     39  1.1  drochner 	uint32_t	v[4];
     40  1.1  drochner 	uint32_t	z[4] = { 0, 0, 0, 0};
     41  1.1  drochner 	const uint8_t	*x = (const uint8_t *)X;
     42  1.1  drochner 	uint32_t	mul;
     43  1.1  drochner 	int		i;
     44  1.1  drochner 
     45  1.1  drochner 	v[0] = be32toh(Y[0]);
     46  1.1  drochner 	v[1] = be32toh(Y[1]);
     47  1.1  drochner 	v[2] = be32toh(Y[2]);
     48  1.1  drochner 	v[3] = be32toh(Y[3]);
     49  1.1  drochner 
     50  1.1  drochner 	for (i = 0; i < GMAC_BLOCK_LEN * 8; i++) {
     51  1.1  drochner 		/* update Z */
     52  1.1  drochner 		if (x[i >> 3] & (1 << (~i & 7))) {
     53  1.1  drochner 			z[0] ^= v[0];
     54  1.1  drochner 			z[1] ^= v[1];
     55  1.1  drochner 			z[2] ^= v[2];
     56  1.1  drochner 			z[3] ^= v[3];
     57  1.1  drochner 		} /* else: we preserve old values */
     58  1.1  drochner 
     59  1.1  drochner 		/* update V */
     60  1.1  drochner 		mul = v[3] & 1;
     61  1.1  drochner 		v[3] = (v[2] << 31) | (v[3] >> 1);
     62  1.1  drochner 		v[2] = (v[1] << 31) | (v[2] >> 1);
     63  1.1  drochner 		v[1] = (v[0] << 31) | (v[1] >> 1);
     64  1.1  drochner 		v[0] = (v[0] >> 1) ^ (0xe1000000 * mul);
     65  1.1  drochner 	}
     66  1.1  drochner 
     67  1.1  drochner 	product[0] = htobe32(z[0]);
     68  1.1  drochner 	product[1] = htobe32(z[1]);
     69  1.1  drochner 	product[2] = htobe32(z[2]);
     70  1.1  drochner 	product[3] = htobe32(z[3]);
     71  1.1  drochner }
     72  1.1  drochner 
     73  1.1  drochner void
     74  1.1  drochner ghash_update(GHASH_CTX *ctx, const uint8_t *X, size_t len)
     75  1.1  drochner {
     76  1.1  drochner 	uint8_t *s = (uint8_t *)ctx->S;
     77  1.1  drochner 	uint8_t *y = (uint8_t *)ctx->Z;
     78  1.1  drochner 	int i, j;
     79  1.1  drochner 
     80  1.1  drochner 	for (i = 0; i < len / GMAC_BLOCK_LEN; i++) {
     81  1.1  drochner 		for (j = 0; j < GMAC_BLOCK_LEN; j++)
     82  1.1  drochner 			s[j] = y[j] ^ X[j];
     83  1.1  drochner 
     84  1.1  drochner 		ghash_gfmul(ctx->S, ctx->H, ctx->S);
     85  1.1  drochner 
     86  1.1  drochner 		y = s;
     87  1.1  drochner 		X += GMAC_BLOCK_LEN;
     88  1.1  drochner 	}
     89  1.1  drochner 
     90  1.1  drochner 	memcpy(ctx->Z, ctx->S, GMAC_BLOCK_LEN);
     91  1.1  drochner }
     92  1.1  drochner 
     93  1.1  drochner #define AESCTR_NONCESIZE	4
     94  1.1  drochner 
     95  1.1  drochner void
     96  1.1  drochner AES_GMAC_Init(AES_GMAC_CTX *ctx)
     97  1.1  drochner {
     98  1.1  drochner 
     99  1.1  drochner 	memset(ctx, 0, sizeof(AES_GMAC_CTX));
    100  1.1  drochner }
    101  1.1  drochner 
    102  1.1  drochner void
    103  1.1  drochner AES_GMAC_Setkey(AES_GMAC_CTX *ctx, const uint8_t *key, uint16_t klen)
    104  1.1  drochner {
    105  1.1  drochner 	ctx->rounds = rijndaelKeySetupEnc(ctx->K, (const u_char *)key,
    106  1.1  drochner 	    (klen - AESCTR_NONCESIZE) * 8);
    107  1.1  drochner 	/* copy out salt to the counter block */
    108  1.1  drochner 	memcpy(ctx->J, key + klen - AESCTR_NONCESIZE, AESCTR_NONCESIZE);
    109  1.1  drochner 	/* prepare a hash subkey */
    110  1.1  drochner 	rijndaelEncrypt(ctx->K, ctx->rounds, (void *)ctx->ghash.H,
    111  1.1  drochner 			(void *)ctx->ghash.H);
    112  1.1  drochner }
    113  1.1  drochner 
    114  1.1  drochner void
    115  1.1  drochner AES_GMAC_Reinit(AES_GMAC_CTX *ctx, const uint8_t *iv, uint16_t ivlen)
    116  1.1  drochner {
    117  1.1  drochner 	/* copy out IV to the counter block */
    118  1.1  drochner 	memcpy(ctx->J + AESCTR_NONCESIZE, iv, ivlen);
    119  1.1  drochner }
    120  1.1  drochner 
    121  1.1  drochner int
    122  1.1  drochner AES_GMAC_Update(AES_GMAC_CTX *ctx, const uint8_t *data, uint16_t len)
    123  1.1  drochner {
    124  1.1  drochner 	uint32_t	blk[4] = { 0, 0, 0, 0 };
    125  1.1  drochner 	int		plen;
    126  1.1  drochner 
    127  1.1  drochner 	if (len > 0) {
    128  1.1  drochner 		plen = len % GMAC_BLOCK_LEN;
    129  1.1  drochner 		if (len >= GMAC_BLOCK_LEN)
    130  1.1  drochner 			ghash_update(&ctx->ghash, (const uint8_t *)data,
    131  1.1  drochner 				     len - plen);
    132  1.1  drochner 		if (plen) {
    133  1.1  drochner 			memcpy(blk, data + (len - plen), plen);
    134  1.1  drochner 			ghash_update(&ctx->ghash, (uint8_t *)blk,
    135  1.1  drochner 			    GMAC_BLOCK_LEN);
    136  1.1  drochner 		}
    137  1.1  drochner 	}
    138  1.1  drochner 	return (0);
    139  1.1  drochner }
    140  1.1  drochner 
    141  1.1  drochner void
    142  1.1  drochner AES_GMAC_Final(uint8_t digest[GMAC_DIGEST_LEN], AES_GMAC_CTX *ctx)
    143  1.1  drochner {
    144  1.1  drochner 	uint8_t		keystream[GMAC_BLOCK_LEN];
    145  1.1  drochner 	int		i;
    146  1.1  drochner 
    147  1.1  drochner 	/* do one round of GCTR */
    148  1.1  drochner 	ctx->J[GMAC_BLOCK_LEN - 1] = 1;
    149  1.1  drochner 	rijndaelEncrypt(ctx->K, ctx->rounds, ctx->J, keystream);
    150  1.1  drochner 	for (i = 0; i < GMAC_DIGEST_LEN; i++)
    151  1.1  drochner 		digest[i] = ((uint8_t *)ctx->ghash.S)[i] ^ keystream[i];
    152  1.1  drochner 	memset(keystream, 0, sizeof(keystream));
    153  1.1  drochner }
    154