gra-projekt/Implementierung/src/md2_impls/md2_1.c

78 lines
2.4 KiB
C
Raw Normal View History

2022-07-15 11:32:22 +02:00
#include "../../lib/md2_impls/md2_1.h"
2022-07-18 21:22:03 +02:00
#include <immintrin.h>
2022-07-15 11:32:22 +02:00
2022-07-20 17:46:35 +02:00
#include "../../lib/md2_impls/md2_common.h"
2022-07-18 21:22:03 +02:00
static uint8_t PADDING[17][16] = {
{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{3, 3, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{4, 4, 4, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{5, 5, 5, 5, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{6, 6, 6, 6, 6, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{7, 7, 7, 7, 7, 7, 7, 0, 0, 0, 0, 0, 0, 0, 0, 0},
{8, 8, 8, 8, 8, 8, 8, 8, 0, 0, 0, 0, 0, 0, 0, 0},
{9, 9, 9, 9, 9, 9, 9, 9, 9, 0, 0, 0, 0, 0, 0, 0},
{10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 0, 0, 0, 0, 0, 0},
{11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 0, 0, 0, 0, 0},
{12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 0, 0, 0, 0},
{13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 0, 0, 0},
{14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 0, 0},
{15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 0},
{16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16}};
2022-07-20 17:46:35 +02:00
void md2_hash_1(size_t len, const uint8_t buf[len], uint8_t out[16]) {
2022-07-15 11:32:22 +02:00
// === step 1 ===
2022-07-18 21:22:03 +02:00
int paddingNeeded = 16 - (len & 15);
2022-07-15 11:32:22 +02:00
len += paddingNeeded;
2022-07-20 17:46:35 +02:00
uint8_t *newBuf = aligned_alloc(16, sizeof(uint8_t) * (len + 16));
if (newBuf == NULL) {
2022-07-20 12:11:16 +02:00
return;
}
2022-07-18 21:22:03 +02:00
2022-07-20 17:46:35 +02:00
for (size_t i = 0; i < 16; i++) {
2022-07-18 21:22:03 +02:00
newBuf[len + i] = 0;
}
2022-07-15 11:32:22 +02:00
memcpy(newBuf, buf, len - paddingNeeded);
2022-07-18 21:22:03 +02:00
memcpy(newBuf + len - paddingNeeded, PADDING + paddingNeeded, paddingNeeded);
2022-07-15 11:32:22 +02:00
// === step 2 ===
2022-07-20 23:23:18 +02:00
md2_checksum(len, newBuf);
2022-07-15 11:32:22 +02:00
// === step 3 ===
2022-07-18 21:22:03 +02:00
uint8_t *messageDigestBuf = aligned_alloc(16, sizeof(uint8_t) * 48);
if (messageDigestBuf == NULL) {
2022-07-20 12:11:16 +02:00
return;
}
2022-07-18 21:22:03 +02:00
for (size_t i = 0; i < 48; i++) {
messageDigestBuf[i] = 0;
}
2022-07-15 11:32:22 +02:00
// === step 4 ===
2022-07-18 21:22:03 +02:00
__m128i vx;
__m128i vy;
2022-07-20 17:46:35 +02:00
for (size_t i = 0; i <= (len + 16) / 16 - 1; i++) {
FIRST_LOOP_START_MARK
2022-07-24 23:16:32 +02:00
// first loop replaced by 5 SIMD-Operations
2022-07-20 17:46:35 +02:00
vx = _mm_load_si128((__m128i *)(newBuf + i * 16));
_mm_store_si128((__m128i *)(messageDigestBuf + 16), vx);
vy = _mm_load_si128((__m128i *)(messageDigestBuf));
2022-07-18 21:22:03 +02:00
vy = _mm_xor_si128(vy, vx);
2022-07-20 17:46:35 +02:00
_mm_store_si128((__m128i *)(messageDigestBuf + 32), vy);
FIRST_LOOP_END_MARK
2022-07-15 11:32:22 +02:00
2022-07-20 23:23:18 +02:00
md2_second_loop(messageDigestBuf);
2022-07-15 11:32:22 +02:00
}
2022-07-20 17:46:35 +02:00
END_MARK
2022-07-15 11:32:22 +02:00
memcpy(out, messageDigestBuf, 16);
free(messageDigestBuf);
free(newBuf);
}