nxdumptool/source/crc32_fast.c

/* Standard CRC32 checksum: fast public domain implementation for
 * little-endian architectures.  Written for compilation with an
 * optimizer set to perform loop unwinding.  Outputs the checksum for
 * each file given as a command line argument.  Invalid file names and
 * files that cause errors are silently skipped.  The program reads
 * from stdin if it is called with no arguments. */

#include <stdio.h>
#include <stdint.h>
#include <stdlib.h>

#include "crc32_fast.h"

u32 crc32_for_byte(u32 r)
{
	for(int j = 0; j < 8; ++j) r = (r & 1? 0: (u32)0xEDB88320L) ^ r >> 1;
	return r ^ (u32)0xFF000000L;
}

/* Any unsigned integer type with at least 32 bits may be used as
 * accumulator type for fast crc32-calulation, but unsigned long is
 * probably the optimal choice for most systems. */
typedef unsigned long accum_t;

void init_tables(u32* table, u32* wtable)
{
	for(u64 i = 0; i < 0x100; ++i) table[i] = crc32_for_byte(i);
	
	for(u64 k = 0; k < sizeof(accum_t); ++k)
	{
		for(u64 w, i = 0; i < 0x100; ++i)
		{
			for(u64 j = w = 0; j < sizeof(accum_t); ++j) w = table[(u8)(j == k? w ^ i: w)] ^ w >> 8;
			wtable[(k << 8) + i] = w ^ (k? wtable[0]: 0);
		}
	}
}

void crc32(const void* data, u64 n_bytes, u32* crc)
{
	static u32 table[0x100], wtable[0x100*sizeof(accum_t)];
	u64 n_accum = n_bytes / sizeof(accum_t);
	
	if (!*table) init_tables(table, wtable);
	for(u64 i = 0; i < n_accum; ++i)
	{
		accum_t a = *crc ^ ((accum_t*)data)[i];
		for(u64 j = *crc = 0; j < sizeof(accum_t); ++j) *crc ^= wtable[(j << 8) + (u8)(a >> 8*j)];
	}
	
	for(u64 i = n_accum*sizeof(accum_t); i < n_bytes; ++i) *crc = table[(u8)*crc ^ ((u8*)data)[i]] ^ *crc >> 8;
}