/*	$OpenBSD: md_test.c,v 1.4 2025/05/22 03:24:47 joshua Exp $ */
/*
 * Copyright (c) 2022, 2025 Joshua Sing <joshua@joshuasing.dev>
 *
 * Permission to use, copy, modify, and distribute this software for any
 * purpose with or without fee is hereby granted, provided that the above
 * copyright notice and this permission notice appear in all copies.
 *
 * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
 * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
 * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
 * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
 * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
 * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
 * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
 */

#include <openssl/evp.h>
#include <openssl/md4.h>
#include <openssl/md5.h>

#include <stdint.h>
#include <string.h>

#include "test.h"

struct md_test {
	const int algorithm;
	const uint8_t in[128];
	const size_t in_len;
	const uint8_t out[EVP_MAX_MD_SIZE];
};

static const struct md_test md_tests[] = {
	/* MD4 (RFC 1320 test vectors) */
	{
		.algorithm = NID_md4,
		.in = "",
		.in_len = 0,
		.out = {
			0x31, 0xd6, 0xcf, 0xe0, 0xd1, 0x6a, 0xe9, 0x31,
			0xb7, 0x3c, 0x59, 0xd7, 0xe0, 0xc0, 0x89, 0xc0,
		}
	},
	{
		.algorithm = NID_md4,
		.in = "a",
		.in_len = 1,
		.out = {
			0xbd, 0xe5, 0x2c, 0xb3, 0x1d, 0xe3, 0x3e, 0x46,
			0x24, 0x5e, 0x05, 0xfb, 0xdb, 0xd6, 0xfb, 0x24,
		}
	},
	{
		.algorithm = NID_md4,
		.in = "abc",
		.in_len = 3,
		.out = {
			0xa4, 0x48, 0x01, 0x7a, 0xaf, 0x21, 0xd8, 0x52,
			0x5f, 0xc1, 0x0a, 0xe8, 0x7a, 0xa6, 0x72, 0x9d,
		}
	},
	{
		.algorithm = NID_md4,
		.in = "message digest",
		.in_len = 14,
		.out = {
			0xd9, 0x13, 0x0a, 0x81, 0x64, 0x54, 0x9f, 0xe8,
			0x18, 0x87, 0x48, 0x06, 0xe1, 0xc7, 0x01, 0x4b,
		}
	},
	{
		.algorithm = NID_md4,
		.in = "abcdefghijklmnopqrstuvwxyz",
		.in_len = 26,
		.out = {
			0xd7, 0x9e, 0x1c, 0x30, 0x8a, 0xa5, 0xbb, 0xcd,
			0xee, 0xa8, 0xed, 0x63, 0xdf, 0x41, 0x2d, 0xa9,
		}
	},
	{
		.algorithm = NID_md4,
		.in =
		    "ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuv"
		    "wxyz0123456789",
		.in_len = 62,
		.out = {
			0x04, 0x3f, 0x85, 0x82, 0xf2, 0x41, 0xdb, 0x35,
			0x1c, 0xe6, 0x27, 0xe1, 0x53, 0xe7, 0xf0, 0xe4,
		}
	},
	{
		.algorithm = NID_md4,
		.in =
		    "123456789012345678901234567890123456789012345678"
		    "90123456789012345678901234567890",
		.in_len = 80,
		.out = {
			0xe3, 0x3b, 0x4d, 0xdc, 0x9c, 0x38, 0xf2, 0x19,
			0x9c, 0x3e, 0x7b, 0x16, 0x4f, 0xcc, 0x05, 0x36,
		}
	},

	/* MD5 (RFC 1321 test vectors) */
	{
		.algorithm = NID_md5,
		.in = "",
		.in_len = 0,
		.out = {
			0xd4, 0x1d, 0x8c, 0xd9, 0x8f, 0x00, 0xb2, 0x04,
			0xe9, 0x80, 0x09, 0x98, 0xec, 0xf8, 0x42, 0x7e,
		}
	},
	{
		.algorithm = NID_md5,
		.in = "a",
		.in_len = 1,
		.out = {
			0x0c, 0xc1, 0x75, 0xb9, 0xc0, 0xf1, 0xb6, 0xa8,
			0x31, 0xc3, 0x99, 0xe2, 0x69, 0x77, 0x26, 0x61,
		}
	},
	{
		.algorithm = NID_md5,
		.in = "abc",
		.in_len = 3,
		.out = {
			0x90, 0x01, 0x50, 0x98, 0x3c, 0xd2, 0x4f, 0xb0,
			0xd6, 0x96, 0x3f, 0x7d, 0x28, 0xe1, 0x7f, 0x72,
		}
	},
	{
		.algorithm = NID_md5,
		.in = "message digest",
		.in_len = 14,
		.out = {
			0xf9, 0x6b, 0x69, 0x7d, 0x7c, 0xb7, 0x93, 0x8d,
			0x52, 0x5a, 0x2f, 0x31, 0xaa, 0xf1, 0x61, 0xd0,
		}
	},
	{
		.algorithm = NID_md5,
		.in = "abcdefghijklmnopqrstuvwxyz",
		.in_len = 26,
		.out = {
			0xc3, 0xfc, 0xd3, 0xd7, 0x61, 0x92, 0xe4, 0x00,
			0x7d, 0xfb, 0x49, 0x6c, 0xca, 0x67, 0xe1, 0x3b,
		}
	},
	{
		.algorithm = NID_md5,
		.in =
		    "ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuv"
		    "wxyz0123456789",
		.in_len = 62,
		.out = {
			0xd1, 0x74, 0xab, 0x98, 0xd2, 0x77, 0xd9, 0xf5,
			0xa5, 0x61, 0x1c, 0x2c, 0x9f, 0x41, 0x9d, 0x9f,
		}
	},
	{
		.algorithm = NID_md5,
		.in =
		    "123456789012345678901234567890123456789012345678"
		    "90123456789012345678901234567890",
		.in_len = 80,
		.out = {
			0x57, 0xed, 0xf4, 0xa2, 0x2b, 0xe3, 0xc9, 0x55,
			0xac, 0x49, 0xda, 0x2e, 0x21, 0x07, 0xb6, 0x7a,
		}
	},
};

#define N_MD_TESTS (sizeof(md_tests) / sizeof(md_tests[0]))

typedef unsigned char *(*md_hash_func)(const unsigned char *, size_t,
    unsigned char *);

static int
md_hash_from_algorithm(int algorithm, md_hash_func *out_func,
    const EVP_MD **out_md, size_t *out_len)
{
	switch (algorithm) {
	case NID_md4:
		*out_func = MD4;
		*out_md = EVP_md4();
		*out_len = MD4_DIGEST_LENGTH;
		break;
	case NID_md5:
		*out_func = MD5;
		*out_md = EVP_md5();
		*out_len = MD5_DIGEST_LENGTH;
		break;
	default:
		return 0;
	}

	return 1;
}

static void
test_md_tv(struct test *t, const void *arg)
{
	const struct md_test *st = arg;
	md_hash_func md_func;
	const EVP_MD *md;
	EVP_MD_CTX *hash = NULL;
	uint8_t out[EVP_MAX_MD_SIZE];
	size_t in_len, out_len;

	if (!md_hash_from_algorithm(st->algorithm, &md_func, &md, &out_len)) {
		test_errorf(t, "md_hash_from_algorithm: unknown algorithm: %d",
		    st->algorithm);
		goto fail;
	}

	if ((hash = EVP_MD_CTX_new()) == NULL) {
		test_errorf(t, "EVP_MD_CTX_new()");
		goto fail;
	}

	/* Digest */
	memset(out, 0, sizeof(out));
	md_func(st->in, st->in_len, out);
	if (memcmp(st->out, out, out_len) != 0) {
		test_errorf(t, "MD: digest output mismatch");
		test_hexdiff(t, out, out_len, st->out);
	}

	/* EVP single-shot digest */
	memset(out, 0, sizeof(out));
	if (!EVP_Digest(st->in, st->in_len, out, NULL, md, NULL)) {
		test_errorf(t, "EVP_Digest()");
		goto fail;
	}
	if (memcmp(st->out, out, out_len) != 0) {
		test_errorf(t, "EVP_Digest: digest output mismatch");
		test_hexdiff(t, out, out_len, st->out);
	}

	/* EVP digest */
	memset(out, 0, sizeof(out));
	if (!EVP_DigestInit_ex(hash, md, NULL)) {
		test_errorf(t, "EVP_DigestInit_ex()");
		goto fail;
	}

	in_len = st->in_len / 2;
	if (!EVP_DigestUpdate(hash, st->in, in_len)) {
		test_errorf(t, "EVP_DigestUpdate: first half failed");
		goto fail;
	}

	if (!EVP_DigestUpdate(hash, st->in + in_len,
		st->in_len - in_len)) {
		test_errorf(t, "EVP_DigestUpdate: second half failed");
		goto fail;
	}

	if (!EVP_DigestFinal_ex(hash, out, NULL)) {
		test_errorf(t, "EVP_DigestFinal_ex()");
		goto fail;
	}

	if (memcmp(st->out, out, out_len) != 0) {
		test_errorf(t, "EVP: digest output mismatch");
		test_hexdiff(t, out, out_len, st->out);
	}


 fail:
	EVP_MD_CTX_free(hash);
}

static void
test_md(struct test *t, const void *arg)
{
	const struct md_test *st;
	size_t i;
	char *name;

	for (i = 0; i < N_MD_TESTS; i++) {
		st = &md_tests[i];
		if (asprintf(&name, "%s: '%s'", OBJ_nid2sn(st->algorithm), st->in) == -1) {
			test_errorf(t, "create test name");
			return;
		}

		test_run(t, name, test_md_tv, st);
		free(name);
	}
}

static void
test_md5_large(struct test *t, const void *arg)
{
	MD5_CTX ctx;
	uint8_t in[1024];
	uint8_t out[EVP_MAX_MD_SIZE];
	unsigned int out_len;
	size_t in_len;
	size_t i;
	uint8_t want[] = {
		0xd8, 0xbc, 0xae, 0x13, 0xb5, 0x5a, 0xb0, 0xfc,
		0x7f, 0x8a, 0xe1, 0x78, 0x27, 0x8d, 0x44, 0x1b,
	};

	memset(in, 'A', sizeof(in));
	in_len = sizeof(in);

	memset(out, 0, sizeof(out));
	out_len = 16;

	MD5_Init(&ctx);

	for (i = 0; i < (1<<29) + 1; i += in_len) {
		if (!MD5_Update(&ctx, in, in_len)) {
			test_errorf(t, "MD5_Update()");
			return;
		}
	}
	if (!MD5_Final(out, &ctx)) {
		test_errorf(t, "MD5_Final()");
		return;
	}

	if (memcmp(out, want, out_len) != 0) {
		test_errorf(t, "MD5 digest output mismatch");
		test_hexdump(t, out, out_len);
	}
	if (ctx.Nh != 0x1 || ctx.Nl != 0x2000)
		test_errorf(t, "MD5 incorrect bit length");
}

int
main(int argc, char **argv)
{
	struct test *t = test_init();

	test_run(t, "md", test_md, NULL);
	test_run(t, "md5 large", test_md5_large, NULL);

	return test_result(t);
}