Genesis-Plus-GX/source/sound/eq.c

133 lines
3.2 KiB
C
Raw Normal View History

/*----------------------------------------------------------------------------
2011-07-14 00:49:52 +02:00
//
// 3 Band EQ :)
//
// EQ.C - Main Source file for 3 band EQ
//
// (c) Neil C / Etanza Systems / 2K6
//
// Shouts / Loves / Moans = etanza at lycos dot co dot uk
//
// This work is hereby placed in the public domain for all purposes, including
// use in commercial applications.
//
// The author assumes NO RESPONSIBILITY for any problems caused by the use of
// this software.
//
//----------------------------------------------------------------------------*/
2011-07-14 00:49:52 +02:00
/* NOTES :
2011-07-14 00:49:52 +02:00
//
// - Original filter code by Paul Kellet (musicdsp.pdf)
//
// - Uses 4 first order filters in series, should give 24dB per octave
//
// - Now with P4 Denormal fix :)
//----------------------------------------------------------------------------*/
2011-07-14 00:49:52 +02:00
/* ----------
2011-07-14 00:49:52 +02:00
//| Includes |
// ----------*/
2011-07-14 00:49:52 +02:00
#include <stdio.h>
#include <string.h>
#include <math.h>
#include "eq.h"
#include "macros.h"
2011-07-14 00:49:52 +02:00
/* -----------
2011-07-14 00:49:52 +02:00
//| Constants |
// -----------*/
2011-07-14 00:49:52 +02:00
static double vsa = (1.0 / 4294967295.0); /* Very small amount (Denormal Fix) */
2011-07-14 00:49:52 +02:00
/* ---------------
2011-07-14 00:49:52 +02:00
//| Initialise EQ |
// ---------------*/
2011-07-14 00:49:52 +02:00
/* Recommended frequencies are ...
2011-07-14 00:49:52 +02:00
//
// lowfreq = 880 Hz
// highfreq = 5000 Hz
//
// Set mixfreq to whatever rate your system is using (eg 48Khz)*/
2011-07-14 00:49:52 +02:00
void init_3band_state(EQSTATE * es, int lowfreq, int highfreq, int mixfreq)
{
/* Clear state */
2011-07-14 00:49:52 +02:00
memset(es, 0, sizeof(EQSTATE));
/* Set Low/Mid/High gains to unity */
2011-07-14 00:49:52 +02:00
es->lg = 1.0;
es->mg = 1.0;
es->hg = 1.0;
/* Calculate filter cutoff frequencies */
2011-07-14 00:49:52 +02:00
es->lf = 2 * sin(M_PI * ((double) lowfreq / (double) mixfreq));
es->hf = 2 * sin(M_PI * ((double) highfreq / (double) mixfreq));
}
/* ---------------
2011-07-14 00:49:52 +02:00
//| EQ one sample |
// ---------------*/
2011-07-14 00:49:52 +02:00
/* - sample can be any range you like :)
2011-07-14 00:49:52 +02:00
//
// Note that the output will depend on the gain settings for each band
// (especially the bass) so may require clipping before output, but you
// knew that anyway :)*/
2011-07-14 00:49:52 +02:00
double do_3band(EQSTATE * es, int sample)
{
/* Locals */
2011-07-14 00:49:52 +02:00
double l, m, h; /* Low / Mid / High - Sample Values */
2011-07-14 00:49:52 +02:00
/* Filter #1 (lowpass) */
2011-07-14 00:49:52 +02:00
es->f1p0 += (es->lf * ((double) sample - es->f1p0)) + vsa;
es->f1p1 += (es->lf * (es->f1p0 - es->f1p1));
es->f1p2 += (es->lf * (es->f1p1 - es->f1p2));
es->f1p3 += (es->lf * (es->f1p2 - es->f1p3));
l = es->f1p3;
/* Filter #2 (highpass) */
2011-07-14 00:49:52 +02:00
es->f2p0 += (es->hf * ((double) sample - es->f2p0)) + vsa;
es->f2p1 += (es->hf * (es->f2p0 - es->f2p1));
es->f2p2 += (es->hf * (es->f2p1 - es->f2p2));
es->f2p3 += (es->hf * (es->f2p2 - es->f2p3));
h = es->sdm3 - es->f2p3;
/* Calculate midrange (signal - (low + high)) */
2011-07-14 00:49:52 +02:00
/* m = es->sdm3 - (h + l); */
/* fix from http://www.musicdsp.org/showArchiveComment.php?ArchiveID=236 ? */
2011-07-14 00:49:52 +02:00
m = sample - (h + l);
/* Scale, Combine and store */
2011-07-14 00:49:52 +02:00
l *= es->lg;
m *= es->mg;
h *= es->hg;
/* Shuffle history buffer */
2011-07-14 00:49:52 +02:00
es->sdm3 = es->sdm2;
es->sdm2 = es->sdm1;
es->sdm1 = sample;
/* Return result */
2011-07-14 00:49:52 +02:00
return (int) (l + m + h);
}