shithub: ft²

ref: c0d8c8d6057aaf54f72a72c632fd0b94f33710e4
dir: /src/mixer/ft2_cubicspline.c/

#include <stdint.h>
#include <stdbool.h>
#include <stdlib.h>
#include "ft2_cubicspline.h"

float *fCubicSplineLUT8 = NULL;
float *fCubicSplineLUT16 = NULL;

/* 4-tap cubic spline interpolation (Horner's method)
**
** This simplified LUT formula has been used in several trackers, but
** I think it originates from ModPlug Tracker (not sure about this).
*/

static void calcCubicLUT(float *fLUTPtr, const double dSale)
{
	for (int32_t i = 0; i < CUBIC_PHASES; i++)
	{
		const double x  = i * (1.0 / CUBIC_PHASES);
		const double x2 = x * x;  // x^2
		const double x3 = x2 * x; // x^3

		*fLUTPtr++ = (float)((-0.5 * x3 + 1.0 * x2 - 0.5 * x) * dSale);
		*fLUTPtr++ = (float)(( 1.5 * x3 - 2.5 * x2 + 1.0)     * dSale);
		*fLUTPtr++ = (float)((-1.5 * x3 + 2.0 * x2 + 0.5 * x) * dSale);
		*fLUTPtr++ = (float)(( 0.5 * x3 - 0.5 * x2)           * dSale);
	}
}

bool calcCubicTable(void)
{
	fCubicSplineLUT8 = (float *)malloc(CUBIC_LUT_LEN * sizeof (float));
	if (fCubicSplineLUT8 == NULL)
		return false;

	fCubicSplineLUT16 = (float *)malloc(CUBIC_LUT_LEN * sizeof (float));
	if (fCubicSplineLUT16 == NULL)
		return false;

	calcCubicLUT(fCubicSplineLUT8, 1.0 / 128.0); // for 8-bit samples
	calcCubicLUT(fCubicSplineLUT16, 1.0 / 32768.0); // for 16-bit samples

	return true;
}

void freeCubicTable(void)
{
	if (fCubicSplineLUT8 != NULL)
	{
		free(fCubicSplineLUT8);
		fCubicSplineLUT8 = NULL;
	}

	if (fCubicSplineLUT16 != NULL)
	{
		free(fCubicSplineLUT16);
		fCubicSplineLUT16 = NULL;
	}
}