|
|
//========= Copyright Valve Corporation, All rights reserved. ============//
//
// Purpose:
//
// $NoKeywords: $
//
//=============================================================================//
// snd_dsp.c -- audio processing routines
#include "audio_pch.h"
#include "snd_mix_buf.h"
#include "iprediction.h"
#include "../../common.h" // for parsing routines
#include "vstdlib/random.h"
// memdbgon must be the last include file in a .cpp file!!!
#include "tier0/memdbgon.h"
#define SIGN(d) ((d)<0?-1:1)
#define ABS(a) abs(a)
#define MSEC_TO_SAMPS(a) (((a)*SOUND_DMA_SPEED) / 1000) // convert milliseconds to # samples in equivalent time
#define SEC_TO_SAMPS(a) ((a)*SOUND_DMA_SPEED) // convert seconds to # samples in equivalent time
// Suppress the noisy warnings caused by CLIP_DSP
#if defined(__clang__)
#pragma GCC diagnostic ignored "-Wself-assign"
#endif
#define CLIP_DSP(x) (x)
extern ConVar das_debug;
#define SOUND_MS_PER_FT 1 // sound travels approx 1 foot per millisecond
#define ROOM_MAX_SIZE 1000 // max size in feet of room simulation for dsp
void DSP_ReleaseMemory( void ); bool DSP_LoadPresetFile( void );
extern float Gain_To_dB ( float gain ); extern float dB_To_Gain ( float dB ); extern float Gain_To_Amplitude ( float gain ); extern float Amplitude_To_Gain ( float amplitude );
extern bool g_bdas_room_init; extern bool g_bdas_init_nodes;
//===============================================================================
//
// Digital Signal Processing algorithms for audio FX.
//
// KellyB 2/18/03
//===============================================================================
// Performance notes:
// DSP processing should take no more than 3ms total time per frame to remain on par with hl1
// Assume a min frame rate of 24fps = 42ms per frame
// at 24fps, to maintain 44.1khz output rate, we must process about 1840 mono samples per frame.
// So we must process 1840 samples in 3ms.
// on a 1Ghz CPU (mid-low end CPU) 3ms provides roughly 3,000,000 cycles.
// Thus we have 3e6 / 1840 = 1630 cycles per sample.
#define PBITS 12 // parameter bits
#define PMAX ((1 << PBITS)) // parameter max
// crossfade from y2 to y1 at point r (0 < r < PMAX)
#define XFADE(y1,y2,r) ((y2) + ( ( ((y1) - (y2)) * (r) ) >> PBITS) )
// exponential crossfade from y2 to y1 at point r (0 < r < PMAX)
#define XFADE_EXP(y1, y2, r) ((y2) + ((((((y1) - (y2)) * (r) ) >> PBITS) * (r)) >> PBITS) )
/////////////////////
// dsp helpers
/////////////////////
// reverse delay pointer
inline void DlyPtrReverse (int dlysize, int *psamps, int **ppsamp) { // when *ppsamp = psamps - 1, it wraps around to *ppsamp = psamps + dlysize
if ( *ppsamp < psamps ) *ppsamp += dlysize + 1; }
// advance delay pointer
inline void DlyPtrForward (int dlysize, int *psamps, int **ppsamp) { // when *ppsamp = psamps + dlysize + 1, it wraps around to *ppsamp = psamps
if ( *ppsamp > psamps + dlysize ) *ppsamp -= dlysize + 1; }
// Infinite Impulse Response (feedback) filter, cannonical form
// returns single sample 'out' for current input value 'in'
// in: input sample
// psamp: internal state array, dimension max(cdenom,cnumer) + 1
// cnumer,cdenom: numerator and denominator filter orders
// denom,numer: cdenom+1 dimensional arrays of filter params
//
// for cdenom = 4:
//
// 1 psamp0(n) numer0
// in(n)--->(+)--(*)---.------(*)---->(+)---> out(n)
// ^ | ^
// | [Delay d] |
// | | |
// | -denom1 |psamp1 numer1 |
// ----(*)---.------(*)-------
// ^ | ^
// | [Delay d] |
// | | |
// | -denom2 |psamp2 numer2 |
// ----(*)---.------(*)-------
// ^ | ^
// | [Delay d] |
// | | |
// | -denom3 |psamp3 numer3 |
// ----(*)---.------(*)-------
// ^ | ^
// | [Delay d] |
// | | |
// | -denom4 |psamp4 numer4 |
// ----(*)---.------(*)-------
//
// for each input sample in:
// psamp0 = in - denom1*psamp1 - denom2*psamp2 - ...
// out = numer0*psamp0 + numer1*psamp1 + ...
// psampi = psampi-1, i = cmax, cmax-1, ..., 1
inline int IIRFilter_Update_OrderN ( int cdenom, int *denom, int cnumer, int *numer, int *psamp, int in ) { int cmax, i; int out; int in0; out = 0; in0 = in;
cmax = max ( cdenom, cnumer ); // add input values
// for (i = 1; i <= cdenom; i++)
// psamp[0] -= ( denom[i] * psamp[i] ) >> PBITS;
switch (cdenom) { case 12: in0 -= ( denom[12] * psamp[12] ) >> PBITS; case 11: in0 -= ( denom[11] * psamp[11] ) >> PBITS; case 10: in0 -= ( denom[10] * psamp[10] ) >> PBITS; case 9: in0 -= ( denom[9] * psamp[9] ) >> PBITS; case 8: in0 -= ( denom[8] * psamp[8] ) >> PBITS; case 7: in0 -= ( denom[7] * psamp[7] ) >> PBITS; case 6: in0 -= ( denom[6] * psamp[6] ) >> PBITS; case 5: in0 -= ( denom[5] * psamp[5] ) >> PBITS; case 4: in0 -= ( denom[4] * psamp[4] ) >> PBITS; case 3: in0 -= ( denom[3] * psamp[3] ) >> PBITS; case 2: in0 -= ( denom[2] * psamp[2] ) >> PBITS; default: case 1: in0 -= ( denom[1] * psamp[1] ) >> PBITS; }
psamp[0] = in0; // add output values
//for (i = 0; i <= cnumer; i++)
// out += ( numer[i] * psamp[i] ) >> PBITS;
switch (cnumer) { case 12: out += ( numer[12] * psamp[12] ) >> PBITS; case 11: out += ( numer[11] * psamp[11] ) >> PBITS; case 10: out += ( numer[10] * psamp[10] ) >> PBITS; case 9: out += ( numer[9] * psamp[9] ) >> PBITS; case 8: out += ( numer[8] * psamp[8] ) >> PBITS; case 7: out += ( numer[7] * psamp[7] ) >> PBITS; case 6: out += ( numer[6] * psamp[6] ) >> PBITS; case 5: out += ( numer[5] * psamp[5] ) >> PBITS; case 4: out += ( numer[4] * psamp[4] ) >> PBITS; case 3: out += ( numer[3] * psamp[3] ) >> PBITS; case 2: out += ( numer[2] * psamp[2] ) >> PBITS; default: case 1: out += ( numer[1] * psamp[1] ) >> PBITS; case 0: out += ( numer[0] * psamp[0] ) >> PBITS; } // update internal state (reverse order)
for (i = cmax; i >= 1; i--) psamp[i] = psamp[i-1];
// return current output sample
return out; }
// 1st order filter - faster version
inline int IIRFilter_Update_Order1 ( int *denom, int cnumer, int *numer, int *psamp, int in ) { int out; if (!psamp[0] && !psamp[1] && !in) return 0;
psamp[0] = in - (( denom[1] * psamp[1] ) >> PBITS);
out = ( ( numer[1] * psamp[1] ) + ( numer[0] * psamp[0] ) ) >> PBITS;
psamp[1] = psamp[0];
return out; }
// return 'tdelay' delayed sample from delay buffer
// dlysize: delay samples
// psamps: head of delay buffer psamps[0...dlysize]
// psamp: current data pointer
// sdly: 0...dlysize
inline int GetDly ( int dlysize, int *psamps, int *psamp, int tdelay ) { int *pout;
pout = psamp + tdelay; if ( pout <= (psamps + dlysize)) return *pout; else return *(pout - dlysize - 1); }
// update the delay buffer pointer
// dlysize: delay samples
// psamps: head of delay buffer psamps[0...dlysize]
// ppsamp: data pointer
inline void DlyUpdate ( int dlysize, int *psamps, int **ppsamp ) { // decrement pointer and fix up on buffer boundary
// when *ppsamp = psamps-1, it wraps around to *ppsamp = psamps+dlysize
(*ppsamp)--; DlyPtrReverse ( dlysize, psamps, ppsamp ); }
// simple delay with feedback, no filter in feedback line.
// delaysize: delay line size in samples
// tdelay: tap from this location - <= delaysize
// psamps: delay line buffer pointer of dimension delaysize+1
// ppsamp: circular pointer, must be init to &psamps[0] before first call
// fbgain: feedback value, 0-PMAX (normalized to 0.0-1.0)
// outgain: gain
// in: input sample
// psamps0(n) outgain
// in(n)--->(+)--------.-----(*)-> out(n)
// ^ |
// | [Delay d]
// | |
// | fbgain |Wd(n)
// ----(*)---.
inline int ReverbSimple ( int delaysize, int tdelay, int *psamps, int **ppsamp, int fbgain, int outgain, int in ) { int out, sD;
// get current delay output
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay );
// calculate output + delay * gain
out = in + (( fbgain * sD ) >> PBITS); // write to delay
**ppsamp = out;
// advance internal delay pointers
DlyUpdate ( delaysize, psamps, ppsamp );
return ( (out * outgain) >> PBITS ); }
inline int ReverbSimple_xfade ( int delaysize, int tdelay, int tdelaynew, int xf, int *psamps, int **ppsamp, int fbgain, int outgain, int in ) { int out, sD; int sDnew;
// crossfade from tdelay to tdelaynew samples. xfade is 0..PMAX
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay ); sDnew = GetDly ( delaysize, psamps, *ppsamp, tdelaynew ); sD = sD + (((sDnew - sD) * xf) >> PBITS);
out = in + (( fbgain * sD ) >> PBITS); **ppsamp = out; DlyUpdate ( delaysize, psamps, ppsamp );
return ( (out * outgain) >> PBITS ); }
// multitap simple reverb
// NOTE: tdelay3 > tdelay2 > tdelay1 > t0
// NOTE: fbgain * 4 < 1!
inline int ReverbSimple_multitap ( int delaysize, int tdelay0, int tdelay1, int tdelay2, int tdelay3, int *psamps, int **ppsamp, int fbgain, int outgain, int in ) { int s1, s2, s3, s4, sum; s1 = GetDly ( delaysize, psamps, *ppsamp, tdelay0 ); s2 = GetDly ( delaysize, psamps, *ppsamp, tdelay1 ); s3 = GetDly ( delaysize, psamps, *ppsamp, tdelay2 ); s4 = GetDly ( delaysize, psamps, *ppsamp, tdelay3 );
sum = s1 + s2 + s3 + s4;
// write to delay
**ppsamp = in + ((s4 * fbgain) >> PBITS); // update delay pointers
DlyUpdate ( delaysize, psamps, ppsamp );
return ( ((sum + in) * outgain ) >> PBITS ); }
// modulate smallest tap delay only
inline int ReverbSimple_multitap_xfade ( int delaysize, int tdelay0, int tdelaynew, int xf, int tdelay1, int tdelay2, int tdelay3, int *psamps, int **ppsamp, int fbgain, int outgain, int in ) { int s1, s2, s3, s4, sum; int sD, sDnew; // crossfade from tdelay to tdelaynew tap. xfade is 0..PMAX
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay3 ); sDnew = GetDly ( delaysize, psamps, *ppsamp, tdelaynew );
s4 = sD + (((sDnew - sD) * xf) >> PBITS);
s1 = GetDly ( delaysize, psamps, *ppsamp, tdelay0 ); s2 = GetDly ( delaysize, psamps, *ppsamp, tdelay1 ); s3 = GetDly ( delaysize, psamps, *ppsamp, tdelay2 );
sum = s1 + s2 + s3 + s4;
// write to delay
**ppsamp = in + ((s4 * fbgain) >> PBITS);
// update delay pointers
DlyUpdate ( delaysize, psamps, ppsamp );
return ( ((sum + in) * outgain ) >> PBITS ); }
// straight delay, no feedback
//
// delaysize: delay line size in samples
// tdelay: tap from this location - <= delaysize
// psamps: delay line buffer pointer of dimension delaysize+1
// ppsamp: circular pointer, must be init to &psamps[0] before first call
// in: input sample
//
// in(n)--->[Delay d]---> out(n)
//
inline int DelayLinear ( int delaysize, int tdelay, int *psamps, int **ppsamp, int in ) { int out;
out = GetDly ( delaysize, psamps, *ppsamp, tdelay );
**ppsamp = in;
DlyUpdate ( delaysize, psamps, ppsamp );
return ( out ); }
// crossfade delay values from tdelay to tdelaynew, with xfade1 for tdelay and xfade2 for tdelaynew. xfade = 0...PMAX
inline int DelayLinear_xfade ( int delaysize, int tdelay, int tdelaynew, int xf, int *psamps, int **ppsamp, int in ) { int out; int outnew;
out = GetDly ( delaysize, psamps, *ppsamp, tdelay ); outnew = GetDly ( delaysize, psamps, *ppsamp, tdelaynew );
out = out + (((outnew - out) * xf) >> PBITS);
**ppsamp = in; DlyUpdate ( delaysize, psamps, ppsamp );
return ( out ); }
// lowpass reverberator, replace feedback multiplier 'fbgain' in
// reverberator with a low pass filter
// delaysize: delay line size in samples
// tdelay: tap from this location - <= delaysize
// psamps: delay line buffer pointer of dimension delaysize+1
// ppsamp: circular pointer, must be init to &w[0] before first call
// fbgain: feedback gain (built into filter gain)
// outgain: output gain
// cnumer: filter order
// numer: filter numerator, 0-PMAX (normalized to 0.0-1.0), cnumer+1 dimensional
// denom: filter denominator, 0-PMAX (normalized to 0.0-1.0), cnumer+1 dimensional
// pfsamps: filter state, cnumer+1 dimensional
// in: input sample
// psamps0(n) outgain
// in(n)--->(+)--------------.----(*)--> out(n)
// ^ |
// | [Delay d]
// | |
// | fbgain |Wd(n)
// --(*)--[Filter])-
inline int DelayLowpass ( int delaysize, int tdelay, int *psamps, int **ppsamp, int fbgain, int outgain, int *denom, int Ll, int *numer, int *pfsamps, int in ) { int out, sD;
// delay output is filter input
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay );
// filter output, with feedback 'fbgain' baked into filter params
out = in + IIRFilter_Update_Order1 ( denom, Ll, numer, pfsamps, sD ); // write to delay
**ppsamp = out; // update delay pointers
DlyUpdate ( delaysize, psamps, ppsamp );
// output with gain
return ( (out * outgain) >> PBITS ); }
inline int DelayLowpass_xfade ( int delaysize, int tdelay, int tdelaynew, int xf, int *psamps, int **ppsamp, int fbgain, int outgain, int *denom, int Ll, int *numer, int *pfsamps, int in ) { int out, sD; int sDnew;
// crossfade from tdelay to tdelaynew tap. xfade is 0..PMAX
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay ); sDnew = GetDly ( delaysize, psamps, *ppsamp, tdelaynew ); sD = sD + (((sDnew - sD) * xf) >> PBITS);
// filter output with feedback 'fbgain' baked into filter params
out = in + IIRFilter_Update_Order1 ( denom, Ll, numer, pfsamps, sD );
// write to delay
**ppsamp = out;
// update delay ptrs
DlyUpdate ( delaysize, psamps, ppsamp );
// output with gain
return ( (out * outgain) >> PBITS ); }
// delay is multitap tdelay0,tdelay1,tdelay2,tdelay3
// NOTE: tdelay3 > tdelay2 > tdelay1 > tdelay0
// NOTE: fbgain * 4 < 1!
inline int DelayLowpass_multitap ( int delaysize, int tdelay0, int tdelay1, int tdelay2, int tdelay3, int *psamps, int **ppsamp, int fbgain, int outgain, int *denom, int Ll, int *numer, int *pfsamps, int in ) { int s0, s1, s2, s3, s4, sum; s1 = GetDly ( delaysize, psamps, *ppsamp, tdelay0 ); s2 = GetDly ( delaysize, psamps, *ppsamp, tdelay1 ); s3 = GetDly ( delaysize, psamps, *ppsamp, tdelay2 ); s4 = GetDly ( delaysize, psamps, *ppsamp, tdelay3 );
sum = s1 + s2 + s3 + s4;
s0 = in + IIRFilter_Update_Order1 ( denom, Ll, numer, pfsamps, s4 ); // write to delay
**ppsamp = s0;
// update delay ptrs
DlyUpdate ( delaysize, psamps, ppsamp );
return ( ((sum + in) * outgain ) >> PBITS ); }
inline int DelayLowpass_multitap_xfade ( int delaysize, int tdelay0, int tdelaynew, int xf, int tdelay1, int tdelay2, int tdelay3, int *psamps, int **ppsamp, int fbgain, int outgain, int *denom, int Ll, int *numer, int *pfsamps, int in ) { int s0, s1, s2, s3, s4, sum;
int sD, sDnew;
// crossfade from tdelay to tdelaynew tap. xfade is 0..PMAX
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay3 ); sDnew = GetDly ( delaysize, psamps, *ppsamp, tdelaynew ); s4 = sD + (((sDnew - sD) * xf) >> PBITS);
s1 = GetDly ( delaysize, psamps, *ppsamp, tdelay0 ); s2 = GetDly ( delaysize, psamps, *ppsamp, tdelay1 ); s3 = GetDly ( delaysize, psamps, *ppsamp, tdelay2 ); sum = s1 + s2 + s3 + s4;
s0 = in + IIRFilter_Update_Order1 ( denom, Ll, numer, pfsamps, s4 ); **ppsamp = s0; DlyUpdate ( delaysize, psamps, ppsamp );
return ( ((sum + in) * outgain ) >> PBITS ); }
// linear delay with lowpass filter on delay output and gain stage
// delaysize: delay line size in samples
// tdelay: delay tap from this location - <= delaysize
// psamps: delay line buffer pointer of dimension delaysize+1
// ppsamp: circular pointer, must init &psamps[0] before first call
// fbgain: feedback gain (ignored)
// outgain: output gain
// cnumer: filter order
// numer: filter numerator, 0-PMAX (normalized to 0.0-1.0), cnumer+1 dimensional
// denom: filter denominator, 0-PMAX (normalized to 0.0-1.0), cnumer+1 dimensional
// pfsamps: filter state, cnumer+1 dimensional
// in: input sample
// in(n)--->[Delay d]--->[Filter]-->(*outgain)---> out(n)
inline int DelayLinear_lowpass ( int delaysize, int tdelay, int *psamps, int **ppsamp, int fbgain, int outgain, int *denom, int cnumer, int *numer, int *pfsamps, int in ) { int out, sD; // delay output is filter input
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay );
// calc filter output
out = IIRFilter_Update_Order1 ( denom, cnumer, numer, pfsamps, sD );
// input sample to delay input
**ppsamp = in;
// update delay pointers
DlyUpdate ( delaysize, psamps, ppsamp );
// output with gain
return ( (out * outgain) >> PBITS ); }
inline int DelayLinear_lowpass_xfade ( int delaysize, int tdelay, int tdelaynew, int xf, int *psamps, int **ppsamp, int fbgain, int outgain, int *denom, int cnumer, int *numer, int *pfsamps, int in ) { int out, sD; int sDnew;
// crossfade from tdelay to tdelaynew tap. xfade is 0..PMAX
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay ); sDnew = GetDly ( delaysize, psamps, *ppsamp, tdelaynew ); sD = sD + (((sDnew - sD) * xf) >> PBITS);
out = IIRFilter_Update_Order1 ( denom, cnumer, numer, pfsamps, sD );
**ppsamp = in;
DlyUpdate ( delaysize, psamps, ppsamp );
return ( (out * outgain) >> PBITS ); }
// classic allpass reverb
// delaysize: delay line size in samples
// tdelay: tap from this location - <= D
// psamps: delay line buffer pointer of dimension delaysize+1
// ppsamp: circular pointer, must be init to &psamps[0] before first call
// fbgain: feedback value, 0-PMAX (normalized to 0.0-1.0)
// outgain: gain
// psamps0(n) -fbgain outgain
// in(n)--->(+)--------.-----(*)-->(+)--(*)-> out(n)
// ^ | ^
// | [Delay d] |
// | | |
// | fbgain |psampsd(n) |
// ----(*)---.-------------
//
// for each input sample 'in':
// psamps0 = in + fbgain * psampsd
// y = -fbgain * psamps0 + psampsd
// delay (d, psamps) - psamps is the delay buffer array
//
// or, using circular delay, for each input sample 'in':
//
// Sd = GetDly (delaysize,psamps,ppsamp,delaysize)
// S0 = in + fbgain*Sd
// y = -fbgain*S0 + Sd
// *ppsamp = S0
// DlyUpdate(delaysize, psamps, &ppsamp)
inline int DelayAllpass ( int delaysize, int tdelay, int *psamps, int **ppsamp, int fbgain, int outgain, int in ) { int out, s0, sD; sD = GetDly ( delaysize, psamps, *ppsamp, tdelay ); s0 = in + (( fbgain * sD ) >> PBITS);
out = ( ( -fbgain * s0 ) >> PBITS ) + sD; **ppsamp = s0; DlyUpdate ( delaysize, psamps, ppsamp );
return ( (out * outgain) >> PBITS ); }
inline int DelayAllpass_xfade ( int delaysize, int tdelay, int tdelaynew, int xf, int *psamps, int **ppsamp, int fbgain, int outgain, int in ) { int out, s0, sD; int sDnew;
// crossfade from t to tnew tap. xfade is 0..PMAX
sD = GetDly ( delaysize, psamps, *ppsamp, tdelay ); sDnew = GetDly ( delaysize, psamps, *ppsamp, tdelaynew ); sD = sD + (((sDnew - sD) * xf) >> PBITS);
s0 = in + (( fbgain * sD ) >> PBITS);
out = ( ( -fbgain * s0 ) >> PBITS ) + sD; **ppsamp = s0; DlyUpdate ( delaysize, psamps, ppsamp );
return ( (out * outgain) >> PBITS ); }
///////////////////////////////////////////////////////////////////////////////////
// fixed point math for real-time wave table traversing, pitch shifting, resampling
///////////////////////////////////////////////////////////////////////////////////
#define FIX20_BITS 20 // 20 bits of fractional part
#define FIX20_SCALE (1 << FIX20_BITS)
#define FIX20_INTMAX ((1 << (32 - FIX20_BITS))-1) // maximum step integer
#define FLOAT_TO_FIX20(a) ((int)((a) * (float)FIX20_SCALE)) // convert float to fixed point
#define INT_TO_FIX20(a) (((int)(a)) << FIX20_BITS) // convert int to fixed point
#define FIX20_TO_FLOAT(a) ((float)(a) / (float)FIX20_SCALE) // convert fix20 to float
#define FIX20_INTPART(a) (((int)(a)) >> FIX20_BITS) // get integer part of fixed point
#define FIX20_FRACPART(a) ((a) - (((a) >> FIX20_BITS) << FIX20_BITS)) // get fractional part of fixed point
#define FIX20_FRACTION(a,b) (FIX(a)/(b)) // convert int a to fixed point, divide by b
typedef int fix20int;
/////////////////////////////////
// DSP processor parameter block
/////////////////////////////////
// NOTE: these prototypes must match the XXX_Params ( prc_t *pprc ) and XXX_GetNext ( XXX_t *p, int x ) functions
typedef void * (*prc_Param_t)( void *pprc ); // individual processor allocation functions
typedef int (*prc_GetNext_t) ( void *pdata, int x ); // get next function for processor
typedef int (*prc_GetNextN_t) ( void *pdata, portable_samplepair_t *pbuffer, int SampleCount, int op); // batch version of getnext
typedef void (*prc_Free_t) ( void *pdata ); // free function for processor
typedef void (*prc_Mod_t) (void *pdata, float v); // modulation function for processor
#define OP_LEFT 0 // batch process left channel in place
#define OP_RIGHT 1 // batch process right channel in place
#define OP_LEFT_DUPLICATE 2 // batch process left channel in place, duplicate to right channel
#define PRC_NULL 0 // pass through - must be 0
#define PRC_DLY 1 // simple feedback reverb
#define PRC_RVA 2 // parallel reverbs
#define PRC_FLT 3 // lowpass or highpass filter
#define PRC_CRS 4 // chorus
#define PRC_PTC 5 // pitch shifter
#define PRC_ENV 6 // adsr envelope
#define PRC_LFO 7 // lfo
#define PRC_EFO 8 // envelope follower
#define PRC_MDY 9 // mod delay
#define PRC_DFR 10 // diffusor - n series allpass delays
#define PRC_AMP 11 // amplifier with distortion
#define QUA_LO 0 // quality of filter or reverb. Must be 0,1,2,3.
#define QUA_MED 1
#define QUA_HI 2
#define QUA_VHI 3
#define QUA_MAX QUA_VHI
#define CPRCPARAMS 16 // up to 16 floating point params for each processor type
// processor definition - one for each running instance of a dsp processor
struct prc_t { int type; // PRC type
float prm[CPRCPARAMS]; // dsp processor parameters - array of floats
prc_Param_t pfnParam; // allocation function - takes ptr to prc, returns ptr to specialized data struct for proc type
prc_GetNext_t pfnGetNext; // get next function
prc_GetNextN_t pfnGetNextN; // batch version of get next
prc_Free_t pfnFree; // free function
prc_Mod_t pfnMod; // modulation function
void *pdata; // processor state data - ie: pdly, pflt etc.
};
// processor parameter ranges - for validating parameters during allocation of new processor
typedef struct prm_rng_t { int iprm; // parameter index
float lo; // min value of parameter
float hi; // max value of parameter
} prm_rng_s;
void PRC_CheckParams ( prc_t *pprc, prm_rng_t *prng );
///////////
// Filters
///////////
#define CFLTS 64 // max number of filters simultaneously active
#define FLT_M 12 // max order of any filter
#define FLT_LP 0 // lowpass filter
#define FLT_HP 1 // highpass filter
#define FLT_BP 2 // bandpass filter
#define FTR_MAX FLT_BP
// flt parameters
struct flt_t { bool fused; // true if slot in use
int b[FLT_M+1]; // filter numerator parameters (convert 0.0-1.0 to 0-PMAX representation)
int a[FLT_M+1]; // filter denominator parameters (convert 0.0-1.0 to 0-PMAX representation)
int w[FLT_M+1]; // filter state - samples (dimension of max (M, L))
int L; // filter order numerator (dimension of a[M+1])
int M; // filter order denominator (dimension of b[L+1])
int N; // # of series sections - 1 (0 = 1 section, 1 = 2 sections etc)
flt_t *pf1; // series cascaded versions of filter
flt_t *pf2; flt_t *pf3; };
// flt flts
flt_t flts[CFLTS];
void FLT_Init ( flt_t *pf ) { if ( pf ) Q_memset ( pf, 0, sizeof (flt_t) ); } void FLT_InitAll ( void ) { for ( int i = 0 ; i < CFLTS; i++ ) FLT_Init ( &flts[i] ); }
void FLT_Free ( flt_t *pf ) { if ( pf ) { if (pf->pf1) Q_memset ( pf->pf1, 0, sizeof (flt_t) ); if (pf->pf2) Q_memset ( pf->pf2, 0, sizeof (flt_t) ); if (pf->pf3) Q_memset ( pf->pf3, 0, sizeof (flt_t) ); Q_memset ( pf, 0, sizeof (flt_t) ); } }
void FLT_FreeAll ( void ) { for (int i = 0 ; i < CFLTS; i++) FLT_Free ( &flts[i] ); }
// find a free filter from the filter pool
// initialize filter numerator, denominator b[0..M], a[0..L]
// gain scales filter numerator
// N is # of series sections - 1
flt_t * FLT_Alloc ( int N, int M, int L, int *a, int *b, float gain ) { int i, j; flt_t *pf = NULL; for (i = 0; i < CFLTS; i++) { if ( !flts[i].fused ) { pf = &flts[i];
// transfer filter params into filter struct
pf->M = M; pf->L = L; pf->N = N;
for (j = 0; j <= M; j++) pf->a[j] = a[j];
for (j = 0; j <= L; j++) pf->b[j] = (int)((float)(b[j]) * gain);
pf->pf1 = NULL; pf->pf2 = NULL; pf->pf3 = NULL;
pf->fused = true; break; } }
Assert(pf); // make sure we're not trying to alloc more than CFLTS flts
return pf; }
// convert filter params cutoff and type into
// iir transfer function params M, L, a[], b[]
// iir filter, 1st order, transfer function is H(z) = b0 + b1 Z^-1 / a0 + a1 Z^-1
// or H(z) = b0 - b1 Z^-1 / a0 + a1 Z^-1 for lowpass
// design cutoff filter at 3db (.5 gain) p579
void FLT_Design_3db_IIR ( float cutoff, float ftype, int *pM, int *pL, int *a, int *b ) { // ftype: FLT_LP, FLT_HP, FLT_BP
double Wc = 2.0 * M_PI * cutoff / SOUND_DMA_SPEED; // radians per sample
double Oc; double fa; double fb;
// calculations:
// Wc = 2pi * fc/44100 convert to radians
// Oc = tan (Wc/2) * Gc / sqt ( 1 - Gc^2) get analog version, low pass
// Oc = tan (Wc/2) * (sqt (1 - Gc^2)) / Gc analog version, high pass
// Gc = 10 ^ (-Ac/20) gain at cutoff. Ac = 3db, so Gc^2 = 0.5
// a = ( 1 - Oc ) / ( 1 + Oc )
// b = ( 1 - a ) / 2
Oc = tan ( Wc / 2.0 );
fa = ( 1.0 - Oc ) / ( 1.0 + Oc );
fb = ( 1.0 - fa ) / 2.0;
if ( ftype == FLT_HP ) fb = ( 1.0 + fa ) / 2.0;
a[0] = 0; // a0 always ignored
a[1] = (int)( -fa * PMAX ); // quantize params down to 0-PMAX >> PBITS
b[0] = (int)( fb * PMAX ); b[1] = b[0]; if ( ftype == FLT_HP ) b[1] = -b[1];
*pM = *pL = 1;
return; }
// filter parameter order
typedef enum { flt_iftype, flt_icutoff, flt_iqwidth, flt_iquality, flt_igain, flt_cparam // # of params
} flt_e;
// filter parameter ranges
prm_rng_t flt_rng[] = {
{flt_cparam, 0, 0}, // first entry is # of parameters
{flt_iftype, 0, FTR_MAX}, // filter type FLT_LP, FLT_HP, FLT_BP
{flt_icutoff, 10, 22050}, // cutoff frequency in hz at -3db gain
{flt_iqwidth, 0, 11025}, // width of BP (cut in starts at cutoff)
{flt_iquality, 0, QUA_MAX}, // QUA_LO, _MED, _HI, _VHI = # of series sections
{flt_igain, 0.0, 10.0}, // output gain 0-10.0
};
// convert prc float params to iir filter params, alloc filter and return ptr to it
// filter quality set by prc quality - 0,1,2
flt_t * FLT_Params ( prc_t *pprc ) { float qual = pprc->prm[flt_iquality]; float cutoff = pprc->prm[flt_icutoff]; float ftype = pprc->prm[flt_iftype]; float qwidth = pprc->prm[flt_iqwidth]; float gain = pprc->prm[flt_igain];
int L = 0; // numerator order
int M = 0; // denominator order
int b[FLT_M+1]; // numerator params 0..PMAX
int b_scaled[FLT_M+1]; // gain scaled numerator
int a[FLT_M+1]; // denominator params 0..PMAX
int L_bp = 0; // bandpass numerator order
int M_bp = 0; // bandpass denominator order
int b_bp[FLT_M+1]; // bandpass numerator params 0..PMAX
int b_bp_scaled[FLT_M+1]; // gain scaled numerator
int a_bp[FLT_M+1]; // bandpass denominator params 0..PMAX
int N; // # of series sections
bool bpass = false;
// if qwidth > 0 then alloc bandpass filter (pf is lowpass)
if ( qwidth > 0.0 ) bpass = true; if (bpass) { ftype = FLT_LP; }
// low pass and highpass filter design
// 1st order IIR filter, 3db cutoff at fc
if ( bpass ) { // highpass section
FLT_Design_3db_IIR ( cutoff, FLT_HP, &M_bp, &L_bp, a_bp, b_bp ); M_bp = clamp (M_bp, 1, FLT_M); L_bp = clamp (L_bp, 1, FLT_M); cutoff += qwidth; } // lowpass section
FLT_Design_3db_IIR ( cutoff, (int)ftype, &M, &L, a, b ); M = clamp (M, 1, FLT_M); L = clamp (L, 1, FLT_M);
// quality = # of series sections - 1
N = clamp ((int)qual, 0, 3);
// make sure we alloc at least 2 filters
if (bpass) N = max(N, 1);
flt_t *pf0 = NULL; flt_t *pf1 = NULL; flt_t *pf2 = NULL; flt_t *pf3 = NULL; // scale b numerators with gain - only scale for first filter if series filters
for (int i = 0; i < FLT_M; i++) { b_bp_scaled[i] = (int)((float)(b_bp[i]) * gain ); b_scaled[i] = (int)((float)(b[i]) * gain ); }
if (bpass) { // 1st filter is lowpass
pf0 = FLT_Alloc ( N, M_bp, L_bp, a_bp, b_bp_scaled, 1.0 ); } else { pf0 = FLT_Alloc ( N, M, L, a, b_scaled, 1.0 ); }
// allocate series filters
if (pf0) { switch (N) { case 3: // alloc last filter as lowpass also if FLT_BP
if (bpass) pf3 = FLT_Alloc ( 0, M_bp, L_bp, a_bp, b_bp, 1.0 ); else pf3 = FLT_Alloc ( 0, M, L, a, b, 1.0 ); case 2: pf2 = FLT_Alloc ( 0, M, L, a, b, 1.0 ); case 1: pf1 = FLT_Alloc ( 0, M, L, a, b, 1.0 ); case 0: break; } pf0->pf1 = pf1; pf0->pf2 = pf2; pf0->pf3 = pf3; }
return pf0; }
inline void * FLT_VParams ( void *p ) { PRC_CheckParams( (prc_t *)p, flt_rng); return (void *) FLT_Params ((prc_t *)p); }
inline void FLT_Mod ( void *p, float v ) { return; }
// get next filter value for filter pf and input x
inline int FLT_GetNext ( flt_t *pf, int x ) { flt_t *pf1; flt_t *pf2; flt_t *pf3; int y;
switch( pf->N ) { default: case 0: return IIRFilter_Update_Order1(pf->a, pf->L, pf->b, pf->w, x); case 1: pf1 = pf->pf1;
y = IIRFilter_Update_Order1(pf->a, pf->L, pf->b, pf->w, x); return IIRFilter_Update_Order1(pf1->a, pf1->L, pf1->b, pf1->w, y); case 2: pf1 = pf->pf1; pf2 = pf->pf2; y = IIRFilter_Update_Order1(pf->a, pf->L, pf->b, pf->w, x); y = IIRFilter_Update_Order1(pf1->a, pf1->L, pf1->b, pf1->w, y); return IIRFilter_Update_Order1(pf2->a, pf2->L, pf2->b, pf2->w, y); case 3: pf1 = pf->pf1; pf2 = pf->pf2; pf3 = pf->pf3;
y = IIRFilter_Update_Order1(pf->a, pf->L, pf->b, pf->w, x); y = IIRFilter_Update_Order1(pf1->a, pf1->L, pf1->b, pf1->w, y); y = IIRFilter_Update_Order1(pf2->a, pf2->L, pf2->b, pf2->w, y); return IIRFilter_Update_Order1(pf3->a, pf3->L, pf3->b, pf3->w, y); } }
// batch version for performance
inline void FLT_GetNextN( flt_t *pflt, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = FLT_GetNext( pflt, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = FLT_GetNext( pflt, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = FLT_GetNext( pflt, pb->left ); pb++; } return; } }
///////////////////////////////////////////////////////////////////////////
// Positional updaters for pitch shift etc
///////////////////////////////////////////////////////////////////////////
// looping position within a wav, with integer and fractional parts
// used for pitch shifting, upsampling/downsampling
// 20 bits of fraction, 8+ bits of integer
struct pos_t {
fix20int step; // wave table whole and fractional step value
fix20int cstep; // current cummulative step value
int pos; // current position within wav table
int D; // max dimension of array w[0...D] ie: # of samples = D+1
};
// circular wrap of pointer p, relative to array w
// D max buffer index w[0...D] (count of samples in buffer is D+1)
// i circular index
inline void POS_Wrap ( int D, int *i ) { if ( *i > D ) *i -= D + 1; // when *pi = D + 1, it wraps around to *pi = 0
if ( *i < 0 ) *i += D + 1; // when *pi = - 1, it wraps around to *pi = D
}
// set initial update value - fstep can have no more than 8 bits of integer and 20 bits of fract
// D is array max dimension w[0...D] (ie: size D+1)
// w is ptr to array
// p is ptr to pos_t to initialize
inline void POS_Init( pos_t *p, int D, float fstep ) { float step = fstep;
// make sure int part of step is capped at fix20_intmax
if ((int)step > FIX20_INTMAX) step = (step - (int)step) + FIX20_INTMAX;
p->step = FLOAT_TO_FIX20(step); // convert fstep to fixed point
p->cstep = 0; p->pos = 0; // current update value
p->D = D; // always init to end value, in case we're stepping backwards
}
// change step value - this is an instantaneous change, not smoothed.
inline void POS_ChangeVal( pos_t *p, float fstepnew ) { p->step = FLOAT_TO_FIX20( fstepnew ); // convert fstep to fixed point
}
// return current integer position, then update internal position value
inline int POS_GetNext ( pos_t *p ) { //float f = FIX20_TO_FLOAT(p->cstep);
//int i1 = FIX20_INTPART(p->cstep);
//float f1 = FIX20_TO_FLOAT(FIX20_FRACPART(p->cstep));
//float f2 = FIX20_TO_FLOAT(p->step);
p->cstep += p->step; // update accumulated fraction step value (fixed point)
p->pos += FIX20_INTPART( p->cstep ); // update pos with integer part of accumulated step
p->cstep = FIX20_FRACPART( p->cstep ); // throw away the integer part of accumulated step
// wrap pos around either end of buffer if needed
POS_Wrap(p->D, &(p->pos));
// make sure returned position is within array bounds
Assert (p->pos <= p->D);
return p->pos; }
// oneshot position within wav
struct pos_one_t { pos_t p; // pos_t
bool fhitend; // flag indicating we hit end of oneshot wav
};
// set initial update value - fstep can have no more than 8 bits of integer and 20 bits of fract
// one shot position - play only once, don't wrap, when hit end of buffer, return last position
inline void POS_ONE_Init( pos_one_t *p1, int D, float fstep ) { POS_Init( &p1->p, D, fstep ) ; p1->fhitend = false; }
// return current integer position, then update internal position value
inline int POS_ONE_GetNext ( pos_one_t *p1 ) { int pos; pos_t *p0;
pos = p1->p.pos; // return current position
if (p1->fhitend) return pos;
p0 = &(p1->p); p0->cstep += p0->step; // update accumulated fraction step value (fixed point)
p0->pos += FIX20_INTPART( p0->cstep ); // update pos with integer part of accumulated step
//p0->cstep = SIGN(p0->cstep) * FIX20_FRACPART( p0->cstep );
p0->cstep = FIX20_FRACPART( p0->cstep ); // throw away the integer part of accumulated step
// if we wrapped, stop updating, always return last position
// if step value is 0, return hit end
if (!p0->step || p0->pos < 0 || p0->pos >= p0->D ) p1->fhitend = true; else pos = p0->pos;
// make sure returned value is within array bounds
Assert ( pos <= p0->D );
return pos; }
/////////////////////
// Reverbs and delays
/////////////////////
#define CDLYS 128 // max delay lines active. Also used for lfos.
#define DLY_PLAIN 0 // single feedback loop
#define DLY_ALLPASS 1 // feedback and feedforward loop - flat frequency response (diffusor)
#define DLY_LOWPASS 2 // lowpass filter in feedback loop
#define DLY_LINEAR 3 // linear delay, no feedback, unity gain
#define DLY_FLINEAR 4 // linear delay with lowpass filter and output gain
#define DLY_LOWPASS_4TAP 5 // lowpass filter in feedback loop, 4 delay taps
#define DLY_PLAIN_4TAP 6 // single feedback loop, 4 delay taps
#define DLY_MAX DLY_PLAIN_4TAP
#define DLY_HAS_MULTITAP(a) ((a) == DLY_LOWPASS_4TAP || (a) == DLY_PLAIN_4TAP)
#define DLY_HAS_FILTER(a) ((a) == DLY_FLINEAR || (a) == DLY_LOWPASS || (a) == DLY_LOWPASS_4TAP)
#define DLY_TAP_FEEDBACK_GAIN 0.25 // drop multitap feedback to compensate for sum of taps in dly_*multitap()
#define DLY_NORMALIZING_REDUCTION_MAX 0.25 // don't reduce gain (due to feedback) below N% of original gain
// delay line
struct dly_t { bool fused; // true if dly is in use
int type; // delay type
int D; // delay size, in samples
int t; // current tap, <= D
int tnew; // crossfading to tnew
int xf; // crossfade value of t (0..PMAX)
int t1,t2,t3; // additional taps for multi-tap delays
int a1,a2,a3; // feedback values for taps
int D0; // original delay size (only relevant if calling DLY_ChangeVal)
int *p; // circular buffer pointer
int *w; // array of samples
int a; // feedback value 0..PMAX,normalized to 0-1.0
int b; // gain value 0..PMAX, normalized to 0-1.0
flt_t *pflt; // pointer to filter, if type DLY_LOWPASS
};
dly_t dlys[CDLYS]; // delay lines
void DLY_Init ( dly_t *pdly ) { if ( pdly ) Q_memset( pdly, 0, sizeof (dly_t)); } void DLY_InitAll ( void ) { for (int i = 0 ; i < CDLYS; i++) DLY_Init ( &dlys[i] ); } void DLY_Free ( dly_t *pdly ) { // free memory buffer
if ( pdly ) { FLT_Free ( pdly->pflt );
if ( pdly->w ) { delete[] pdly->w; } // free dly slot
Q_memset ( pdly, 0, sizeof (dly_t) ); } }
void DLY_FreeAll ( void ) { for (int i = 0; i < CDLYS; i++ ) DLY_Free ( &dlys[i] ); }
// return adjusted feedback value for given dly
// such that decay time is same as that for dmin and fbmin
// dmin - minimum delay
// fbmin - minimum feedback
// dly - delay to match decay to dmin, fbmin
float DLY_NormalizeFeedback ( int dmin, float fbmin, int dly ) { // minimum decay time T to -60db for a simple reverb is:
// Tmin = (ln 10^-3 / Ln fbmin) * (Dmin / fs)
// where fs = sample frequency
// similarly,
// Tdly = (ln 10^-3 / Ln fb) * (D / fs)
// setting Tdly = Tmin and solving for fb gives:
// D / Dmin = ln fb / ln fbmin
// since y^x = z gives x = ln z / ln y
// fb = fbmin ^ (D/Dmin)
float fb = powf (fbmin, (float)dly / (float) dmin);
return fb; }
// set up 'b' gain parameter of feedback delay to
// compensate for gain caused by feedback 'fb'.
void DLY_SetNormalizingGain ( dly_t *pdly, int feedback ) { // compute normalized gain, set as output gain
// calculate gain of delay line with feedback, and use it to
// reduce output. ie: force delay line with feedback to unity gain
// for constant input x with feedback fb:
// out = x + x*fb + x * fb^2 + x * fb^3...
// gain = out/x
// so gain = 1 + fb + fb^2 + fb^3...
// which, by the miracle of geometric series, equates to 1/1-fb
// thus, gain = 1/(1-fb)
float fgain = 0; float gain; int b; float fb = (float)feedback;
fb = fb / (float)PMAX; fb = fpmin(fb, 0.999f); // if b is 0, set b to PMAX (1)
b = pdly->b ? pdly->b : PMAX;
fgain = 1.0 / (1.0 - fb);
// compensating gain - multiply rva output by gain then >> PBITS
gain = (int)((1.0 / fgain) * PMAX);
gain = gain * 4; // compensate for fact that gain calculation is for +/- 32767 amplitude wavs
// ie: ok to allow a bit more gain because most wavs are not at theoretical peak amplitude at all times
// limit gain reduction to N% PMAX
gain = clamp (gain, (float)(PMAX * DLY_NORMALIZING_REDUCTION_MAX), (float)PMAX); gain = ((float)b/(float)PMAX) * gain; // scale final gain by pdly->b.
pdly->b = (int)gain; }
void DLY_ChangeTaps ( dly_t *pdly, int t0, int t1, int t2, int t3 );
// allocate a new delay line
// D number of samples to delay
// a feedback value (0-PMAX normalized to 0.0-1.0)
// b gain value (0-PMAX normalized to 0.0-1.0) - this is folded into the filter fb params
// if DLY_LOWPASS or DLY_FLINEAR:
// L - numerator order of filter
// M - denominator order of filter
// fb - numerator params, M+1
// fa - denominator params, L+1
dly_t * DLY_AllocLP ( int D, int a, int b, int type, int M, int L, int *fa, int *fb ) { int *w; int i; dly_t *pdly = NULL; int feedback;
// find open slot
for (i = 0; i < CDLYS; i++) { if (!dlys[i].fused) { pdly = &dlys[i]; DLY_Init( pdly ); break; } } if ( i == CDLYS ) { DevMsg ("DSP: Warning, failed to allocate delay line.\n" ); return NULL; // all delay lines in use
}
// save original feedback value
feedback = a;
// adjust feedback a, gain b if delay is multitap unit
if ( DLY_HAS_MULTITAP(type) ) { // split output gain over 4 taps
b = (int)((float)(b) * DLY_TAP_FEEDBACK_GAIN); }
if ( DLY_HAS_FILTER(type) ) { // alloc lowpass iir_filter
// delay feedback gain is built into filter gain
float gain = (float)a / (float)(PMAX); pdly->pflt = FLT_Alloc( 0, M, L, fa, fb, gain ); if ( !pdly->pflt ) { DevMsg ("DSP: Warning, failed to allocate filter for delay line.\n" ); return NULL; } }
// alloc delay memory
w = new int[D+1]; if ( !w ) { Warning( "Sound DSP: Failed to lock.\n"); FLT_Free ( pdly->pflt ); return NULL; } // clear delay array
Q_memset (w, 0, sizeof(int) * (D+1)); // init values
pdly->type = type; pdly->D = D; pdly->t = D; // set delay tap to full delay
pdly->tnew = D; pdly->xf = 0; pdly->D0 = D; pdly->p = w; // init circular pointer to head of buffer
pdly->w = w; pdly->a = min( a, PMAX - 1 ); // do not allow 100% feedback
pdly->b = b; pdly->fused = true;
if ( type == DLY_LINEAR || type == DLY_FLINEAR ) { // linear delay has no feedback and unity gain
pdly->a = 0; pdly->b = PMAX; } else { // adjust b to compensate for feedback gain of steady state max input
DLY_SetNormalizingGain( pdly, feedback ); }
if ( DLY_HAS_MULTITAP(type) ) { // initially set up all taps to same value - caller uses DLY_ChangeTaps to change values
DLY_ChangeTaps( pdly, D, D, D, D ); }
return (pdly); }
// allocate lowpass or allpass delay
dly_t * DLY_Alloc( int D, int a, int b, int type ) { return DLY_AllocLP( D, a, b, type, 0, 0, 0, 0 ); }
// Allocate new delay, convert from float params in prc preset to internal parameters
// Uses filter params in prc if delay is type lowpass
// delay parameter order
typedef enum {
dly_idtype, // NOTE: first 8 params must match those in mdy_e
dly_idelay, dly_ifeedback, dly_igain,
dly_iftype, dly_icutoff, dly_iqwidth, dly_iquality,
dly_itap1, dly_itap2, dly_itap3,
dly_cparam
} dly_e;
// delay parameter ranges
prm_rng_t dly_rng[] = {
{dly_cparam, 0, 0}, // first entry is # of parameters
// delay params
{dly_idtype, 0, DLY_MAX}, // delay type DLY_PLAIN, DLY_LOWPASS, DLY_ALLPASS etc
{dly_idelay, -1.0, 1000.0}, // delay in milliseconds (-1 forces auto dsp to set delay value from room size)
{dly_ifeedback, 0.0, 0.99}, // feedback 0-1.0
{dly_igain, 0.0, 10.0}, // final gain of output stage, 0-10.0
// filter params if dly type DLY_LOWPASS or DLY_FLINEAR
{dly_iftype, 0, FTR_MAX}, {dly_icutoff, 10.0, 22050.0}, {dly_iqwidth, 100.0, 11025.0}, {dly_iquality, 0, QUA_MAX}, // note: -1 flag tells auto dsp to get value directly from room size
{dly_itap1, -1.0, 1000.0}, // delay in milliseconds NOTE: delay > tap3 > tap2 > tap1
{dly_itap2, -1.0, 1000.0}, // delay in milliseconds
{dly_itap3, -1.0, 1000.0}, // delay in milliseconds
};
dly_t * DLY_Params ( prc_t *pprc ) { dly_t *pdly = NULL; int D, a, b; float delay = fabs(pprc->prm[dly_idelay]); float feedback = pprc->prm[dly_ifeedback]; float gain = pprc->prm[dly_igain]; int type = pprc->prm[dly_idtype];
float ftype = pprc->prm[dly_iftype]; float cutoff = pprc->prm[dly_icutoff]; float qwidth = pprc->prm[dly_iqwidth]; float qual = pprc->prm[dly_iquality];
float t1 = fabs(pprc->prm[dly_itap1]); float t2 = fabs(pprc->prm[dly_itap2]); float t3 = fabs(pprc->prm[dly_itap3]);
D = MSEC_TO_SAMPS(delay); // delay samples
a = feedback * PMAX; // feedback
b = gain * PMAX; // gain
switch ( (int) type ) { case DLY_PLAIN: case DLY_PLAIN_4TAP: case DLY_ALLPASS: case DLY_LINEAR: pdly = DLY_Alloc( D, a, b, type ); break;
case DLY_FLINEAR: case DLY_LOWPASS: case DLY_LOWPASS_4TAP: { // set up dummy lowpass filter to convert params
prc_t prcf;
prcf.prm[flt_iquality] = qual; // 0,1,2 - (0 or 1 low quality implies faster execution time)
prcf.prm[flt_icutoff] = cutoff; prcf.prm[flt_iftype] = ftype; prcf.prm[flt_iqwidth] = qwidth; prcf.prm[flt_igain] = 1.0;
flt_t *pflt = (flt_t *)FLT_Params ( &prcf ); if ( !pflt ) { DevMsg ("DSP: Warning, failed to allocate filter.\n" ); return NULL; }
pdly = DLY_AllocLP ( D, a, b, type, pflt->M, pflt->L, pflt->a, pflt->b ); FLT_Free ( pflt ); break; } }
// set up multi-tap delays
if ( pdly && DLY_HAS_MULTITAP((int)type) ) DLY_ChangeTaps( pdly, D, MSEC_TO_SAMPS(t1), MSEC_TO_SAMPS(t2), MSEC_TO_SAMPS(t3) );
return pdly; }
inline void * DLY_VParams ( void *p ) { PRC_CheckParams( (prc_t *)p, dly_rng ); return (void *) DLY_Params ((prc_t *)p); }
// get next value from delay line, move x into delay line
inline int DLY_GetNext ( dly_t *pdly, int x ) { switch (pdly->type) { default: case DLY_PLAIN: return ReverbSimple( pdly->D, pdly->t, pdly->w, &pdly->p, pdly->a, pdly->b, x ); case DLY_ALLPASS: return DelayAllpass( pdly->D, pdly->t, pdly->w, &pdly->p, pdly->a, pdly->b, x ); case DLY_LOWPASS: return DelayLowpass( pdly->D, pdly->t, pdly->w, &(pdly->p), pdly->a, pdly->b, pdly->pflt->a, pdly->pflt->L, pdly->pflt->b, pdly->pflt->w, x ); case DLY_LINEAR: return DelayLinear( pdly->D, pdly->t, pdly->w, &pdly->p, x ); case DLY_FLINEAR: return DelayLinear_lowpass( pdly->D, pdly->t, pdly->w, &(pdly->p), pdly->a, pdly->b, pdly->pflt->a, pdly->pflt->L, pdly->pflt->b, pdly->pflt->w, x ); case DLY_PLAIN_4TAP: return ReverbSimple_multitap( pdly->D, pdly->t, pdly->t1, pdly->t2, pdly->t3, pdly->w, &pdly->p, pdly->a, pdly->b, x ); case DLY_LOWPASS_4TAP: return DelayLowpass_multitap( pdly->D, pdly->t, pdly->t1, pdly->t2,pdly->t3, pdly->w, &(pdly->p), pdly->a, pdly->b, pdly->pflt->a, pdly->pflt->L, pdly->pflt->b, pdly->pflt->w, x ); } }
inline int DLY_GetNextXfade ( dly_t *pdly, int x ) {
switch (pdly->type) { default: case DLY_PLAIN: return ReverbSimple_xfade( pdly->D, pdly->t, pdly->tnew, pdly->xf, pdly->w, &pdly->p, pdly->a, pdly->b, x ); case DLY_ALLPASS: return DelayAllpass_xfade( pdly->D, pdly->t, pdly->tnew, pdly->xf, pdly->w, &pdly->p, pdly->a, pdly->b, x ); case DLY_LOWPASS: return DelayLowpass_xfade( pdly->D, pdly->t, pdly->tnew, pdly->xf, pdly->w, &(pdly->p), pdly->a, pdly->b, pdly->pflt->a, pdly->pflt->L, pdly->pflt->b, pdly->pflt->w, x ); case DLY_LINEAR: return DelayLinear_xfade( pdly->D, pdly->t, pdly->tnew, pdly->xf, pdly->w, &pdly->p, x ); case DLY_FLINEAR: return DelayLinear_lowpass_xfade( pdly->D, pdly->t, pdly->tnew, pdly->xf, pdly->w, &(pdly->p), pdly->a, pdly->b, pdly->pflt->a, pdly->pflt->L, pdly->pflt->b, pdly->pflt->w, x ); case DLY_PLAIN_4TAP: return ReverbSimple_multitap_xfade( pdly->D, pdly->t, pdly->tnew, pdly->xf, pdly->t1, pdly->t2, pdly->t3, pdly->w, &pdly->p, pdly->a, pdly->b, x ); case DLY_LOWPASS_4TAP: return DelayLowpass_multitap_xfade( pdly->D, pdly->t, pdly->tnew, pdly->xf, pdly->t1, pdly->t2, pdly->t3, pdly->w, &(pdly->p), pdly->a, pdly->b, pdly->pflt->a, pdly->pflt->L, pdly->pflt->b, pdly->pflt->w, x ); } }
// batch version for performance
// UNDONE: a) unwind this more - pb increments by 2 to avoid pb->left or pb->right deref.
// UNDONE: b) all filter and delay params are dereferenced outside of DLY_GetNext and passed as register values
// UNDONE: c) pull case statement in dly_getnext out, so loop directly calls the inline dly_*() routine.
inline void DLY_GetNextN( dly_t *pdly, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = DLY_GetNext( pdly, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = DLY_GetNext( pdly, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = DLY_GetNext( pdly, pb->left ); pb++; } return; } }
// get tap on t'th sample in delay - don't update buffer pointers, this is done via DLY_GetNext
// Only valid for DLY_LINEAR.
inline int DLY_GetTap ( dly_t *pdly, int t ) { return GetDly (pdly->D, pdly->w, pdly->p, t ); }
#define SWAP(a,b,t) {(t) = (a); (a) = (b); (b) = (t);}
// make instantaneous change to tap values t0..t3
// all values of t must be less than original delay D
// only processed for DLY_LOWPASS_4TAP & DLY_PLAIN_4TAP
// NOTE: pdly->a feedback must have been set before this call!
void DLY_ChangeTaps ( dly_t *pdly, int t0, int t1, int t2, int t3 ) { if (!pdly) return; int temp;
// sort taps to make sure t3 > t2 > t1 > t0 !
for (int i = 0; i < 4; i++) { if (t0 > t1) SWAP(t0, t1, temp); if (t1 > t2) SWAP(t1, t2, temp); if (t2 > t3) SWAP(t2, t3, temp); }
pdly->t = min ( t0, pdly->D0 ); pdly->t1 = min ( t1, pdly->D0 ); pdly->t2 = min ( t2, pdly->D0 ); pdly->t3 = min ( t3, pdly->D0 );
}
// make instantaneous change for first delay tap 't' to new delay value.
// t tap value must be <= original D (ie: we don't do any reallocation here)
void DLY_ChangeVal ( dly_t *pdly, int t ) { // never set delay > original delay
pdly->t = min ( t, pdly->D0 ); }
// ignored - use MDY_ for modulatable delay
inline void DLY_Mod ( void *p, float v ) { return; }
/////////////////////////////////////////////////////////////////////////////
// Ramp - used for varying smoothly between int parameters ie: modulation delays
/////////////////////////////////////////////////////////////////////////////
struct rmp_t { int initval; // initial ramp value
int target; // final ramp value
int sign; // increasing (1) or decreasing (-1) ramp
int yprev; // previous output value
bool fhitend; // true if hit end of ramp
bool bEndAtTime; // if true, fhitend is true when ramp time is hit (even if target not hit)
// if false, then fhitend is true only when target is hit
pos_one_t ps; // current ramp output
};
// ramp smoothly between initial value and target value in approx 'ramptime' seconds.
// (initial value may be greater or less than target value)
// never changes output by more than +1 or -1 (which can cause the ramp to take longer to complete than ramptime - see bEndAtTime)
// called once per sample while ramping
// ramptime - duration of ramp in seconds
// initval - initial ramp value
// targetval - target ramp value
// if bEndAtTime is true, then RMP_HitEnd returns true when ramp time is reached, EVEN IF TARGETVAL IS NOT REACHED
// if bEndAtTime is false, then RMP_HitEnd returns true when targetval is reached, EVEN IF DELTA IN RAMP VALUES IS > +/- 1
void RMP_Init( rmp_t *prmp, float ramptime, int initval, int targetval, bool bEndAtTime ) { int rise; int run;
if (prmp) Q_memset( prmp, 0, sizeof (rmp_t) ); else return; run = (int) (ramptime * SOUND_DMA_SPEED); // 'samples' in ramp
rise = (targetval - initval); // height of ramp
// init fixed point iterator to iterate along the height of the ramp 'rise'
// always iterates from 0..'rise', increasing in value
POS_ONE_Init( &prmp->ps, ABS( rise ), ABS((float) rise) / ((float) run) ); prmp->yprev = initval; prmp->initval = initval; prmp->target = targetval; prmp->sign = SIGN( rise ); prmp->bEndAtTime = bEndAtTime;
}
// continues from current position to new target position
void RMP_SetNext( rmp_t *prmp, float ramptime, int targetval ) { RMP_Init ( prmp, ramptime, prmp->yprev, targetval, prmp->bEndAtTime ); }
inline bool RMP_HitEnd ( rmp_t *prmp ) { return prmp->fhitend; }
inline void RMP_SetEnd ( rmp_t *prmp ) { prmp->fhitend = true; }
// get next ramp value & update ramp, if bEndAtTime is true, never varies by more than +1 or -1 between calls
// when ramp hits target value, it thereafter always returns last value
inline int RMP_GetNext( rmp_t *prmp ) { int y; int d; // if we hit ramp end, return last value
if (prmp->fhitend) return prmp->yprev; // get next integer position in ramp height.
d = POS_ONE_GetNext( &prmp->ps ); if ( prmp->ps.fhitend ) prmp->fhitend = true;
// increase or decrease from initval, depending on ramp sign
if ( prmp->sign > 0 ) y = prmp->initval + d; else y = prmp->initval - d;
// if bEndAtTime is true, only update current height by a max of +1 or -1
// this also means that for short ramp times, we may not hit target
if (prmp->bEndAtTime) { if ( ABS( y - prmp->yprev ) >= 1 ) prmp->yprev += prmp->sign; } else { // always hits target - but varies by more than +/- 1
prmp->yprev = y; }
return prmp->yprev; }
// get current ramp value, don't update ramp
inline int RMP_GetCurrent( rmp_t *prmp ) { return prmp->yprev; }
//////////////
// mod delay
//////////////
// modulate delay time anywhere from 0..D using MDY_ChangeVal. no output glitches (uses RMP)
#define CMDYS 64 // max # of mod delays active (steals from delays)
struct mdy_t { bool fused;
bool fchanging; // true if modulating to new delay value
dly_t *pdly; // delay
float ramptime; // ramp 'glide' time - time in seconds to change between values
int mtime; // time in samples between delay changes. 0 implies no self-modulating
int mtimecur; // current time in samples until next delay change
float depth; // modulate delay from D to D - (D*depth) depth 0-1.0
int mix; // PMAX as % processed fx signal mix
rmp_t rmp_interp; // interpolation ramp 0...PMAX
bool bPhaseInvert; // if true, invert phase of output
};
mdy_t mdys[CMDYS];
void MDY_Init( mdy_t *pmdy ) { if (pmdy) Q_memset( pmdy, 0, sizeof (mdy_t) ); }; void MDY_Free( mdy_t *pmdy ) { if (pmdy) { DLY_Free (pmdy->pdly); Q_memset( pmdy, 0, sizeof (mdy_t) ); } }; void MDY_InitAll() { for (int i = 0; i < CMDYS; i++) MDY_Init( &mdys[i] ); }; void MDY_FreeAll() { for (int i = 0; i < CMDYS; i++) MDY_Free( &mdys[i] ); };
// allocate mod delay, given previously allocated dly (NOTE: mod delay only sweeps tap 0, not t1,t2 or t3)
// ramptime is time in seconds for delay to change from dcur to dnew
// modtime is time in seconds between modulations. 0 if no self-modulation
// depth is 0-1.0 multiplier, new delay values when modulating are Dnew = randomlong (D - D*depth, D)
// mix - 0-1.0, default 1.0 for 100% fx mix - pans between input signal and fx signal
mdy_t *MDY_Alloc ( dly_t *pdly, float ramptime, float modtime, float depth, float mix ) { int i; mdy_t *pmdy;
if ( !pdly ) return NULL;
for (i = 0; i < CMDYS; i++) { if ( !mdys[i].fused ) { pmdy = &mdys[i];
MDY_Init ( pmdy );
pmdy->pdly = pdly;
if ( !pmdy->pdly ) { DevMsg ("DSP: Warning, failed to allocate delay for mod delay.\n" ); return NULL; }
pmdy->fused = true; pmdy->ramptime = ramptime; pmdy->mtime = SEC_TO_SAMPS(modtime); pmdy->mtimecur = pmdy->mtime; pmdy->depth = depth; pmdy->mix = int ( PMAX * mix ); pmdy->bPhaseInvert = false;
return pmdy; } }
DevMsg ("DSP: Warning, failed to allocate mod delay.\n" ); return NULL; }
// change to new delay tap value t samples, ramp linearly over ramptime seconds
void MDY_ChangeVal ( mdy_t *pmdy, int t ) { // if D > original delay value, cap at original value
t = min (pmdy->pdly->D0, t); pmdy->fchanging = true;
// init interpolation ramp - always hit target
RMP_Init ( &pmdy->rmp_interp, pmdy->ramptime, 0, PMAX, false );
// init delay xfade values
pmdy->pdly->tnew = t; pmdy->pdly->xf = 0; }
// interpolate between current and target delay values
inline int MDY_GetNext( mdy_t *pmdy, int x ) { int xout; if ( !pmdy->fchanging ) { // not modulating...
xout = DLY_GetNext( pmdy->pdly, x );
if ( !pmdy->mtime ) { // return right away if not modulating (not changing and not self modulating)
goto mdy_return; } } else { // modulating...
xout = DLY_GetNextXfade( pmdy->pdly, x ); // get xfade ramp & set up delay xfade value for next call to DLY_GetNextXfade()
pmdy->pdly->xf = RMP_GetNext( &pmdy->rmp_interp ); // 0...PMAX
if ( RMP_HitEnd( &pmdy->rmp_interp ) ) { // done. set delay tap & value = target
DLY_ChangeVal( pmdy->pdly, pmdy->pdly->tnew );
pmdy->pdly->t = pmdy->pdly->tnew;
pmdy->fchanging = false; } } // if self-modulating and timer has expired, get next change
if ( pmdy->mtime && !pmdy->mtimecur-- ) { pmdy->mtimecur = pmdy->mtime;
int D0 = pmdy->pdly->D0; int Dnew; float D1;
// modulate between 0 and 100% of d0
D1 = (float)D0 * (1.0 - pmdy->depth);
Dnew = RandomInt( (int)D1, D0 );
// set up modulation to new value
MDY_ChangeVal ( pmdy, Dnew ); }
mdy_return:
// reverse phase of output
if ( pmdy->bPhaseInvert ) xout = -xout;
// 100% fx mix
if ( pmdy->mix == PMAX) return xout; // special case 50/50 mix
if ( pmdy->mix == PMAX / 2) return ( (xout + x) >> 1 );
// return mix of input and processed signal
return ( x + (((xout - x) * pmdy->mix) >> PBITS) ); }
// batch version for performance
// UNDONE: unwind MDY_GetNext so that it directly calls DLY_GetNextN:
// UNDONE: a) if not currently modulating and never self-modulating, then just unwind like DLY_GetNext
// UNDONE: b) if not currently modulating, figure out how many samples N until self-modulation timer kicks in again
// and stream out N samples just like DLY_GetNext
inline void MDY_GetNextN( mdy_t *pmdy, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = MDY_GetNext( pmdy, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = MDY_GetNext( pmdy, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = MDY_GetNext( pmdy, pb->left ); pb++; } return; } }
// parameter order
typedef enum {
mdy_idtype, // NOTE: first 8 params must match params in dly_e
mdy_idelay, mdy_ifeedback, mdy_igain,
mdy_iftype, mdy_icutoff, mdy_iqwidth, mdy_iquality, mdy_imodrate, mdy_imoddepth, mdy_imodglide,
mdy_imix, mdy_ibxfade,
mdy_cparam
} mdy_e;
// parameter ranges
prm_rng_t mdy_rng[] = {
{mdy_cparam, 0, 0}, // first entry is # of parameters
// delay params
{mdy_idtype, 0, DLY_MAX}, // delay type DLY_PLAIN, DLY_LOWPASS, DLY_ALLPASS
{mdy_idelay, 0.0, 1000.0}, // delay in milliseconds
{mdy_ifeedback, 0.0, 0.99}, // feedback 0-1.0
{mdy_igain, 0.0, 1.0}, // final gain of output stage, 0-1.0
// filter params if mdy type DLY_LOWPASS
{mdy_iftype, 0, FTR_MAX}, {mdy_icutoff, 10.0, 22050.0}, {mdy_iqwidth, 100.0, 11025.0}, {mdy_iquality, 0, QUA_MAX}, {mdy_imodrate, 0.01, 200.0}, // frequency at which delay values change to new random value. 0 is no self-modulation
{mdy_imoddepth, 0.0, 1.0}, // how much delay changes (decreases) from current value (0-1.0)
{mdy_imodglide, 0.01, 100.0}, // glide time between dcur and dnew in milliseconds
{mdy_imix, 0.0, 1.0} // 1.0 = full fx mix, 0.5 = 50% fx, 50% dry
};
// convert user parameters to internal parameters, allocate and return
mdy_t * MDY_Params ( prc_t *pprc ) { mdy_t *pmdy; dly_t *pdly;
float ramptime = pprc->prm[mdy_imodglide] / 1000.0; // get ramp time in seconds
float modtime = 0.0f; if ( pprc->prm[mdy_imodrate] != 0.0f ) { modtime = 1.0 / pprc->prm[mdy_imodrate]; // time between modulations in seconds
} float depth = pprc->prm[mdy_imoddepth]; // depth of modulations 0-1.0
float mix = pprc->prm[mdy_imix];
// alloc plain, allpass or lowpass delay
pdly = DLY_Params( pprc ); if ( !pdly ) return NULL;
pmdy = MDY_Alloc ( pdly, ramptime, modtime, depth, mix ); return pmdy; }
inline void * MDY_VParams ( void *p ) { PRC_CheckParams ( (prc_t *)p, mdy_rng ); return (void *) MDY_Params ((prc_t *)p); }
// v is +/- 0-1.0
// change current delay value 0..D
void MDY_Mod ( mdy_t *pmdy, float v ) {
int D0 = pmdy->pdly->D0; // base delay value
float v2;
// if v is < -2.0 then delay is v + 10.0
// invert phase of output. hack.
if ( v < -2.0 ) { v = v + 10.0; pmdy->bPhaseInvert = true; } else { pmdy->bPhaseInvert = false; }
v2 = -(v + 1.0)/2.0; // v2 varies -1.0-0.0
// D0 varies 0..D0
D0 = D0 + (int)((float)D0 * v2);
// change delay
MDY_ChangeVal( pmdy, D0 );
return; }
///////////////////
// Parallel reverbs
///////////////////
// Reverb A
// M parallel reverbs, mixed to mono output
#define CRVAS 64 // max number of parallel series reverbs active
#define CRVA_DLYS 12 // max number of delays making up reverb_a
struct rva_t { bool fused; int m; // number of parallel plain or lowpass delays
int fparallel; // true if filters in parallel with delays, otherwise single output filter
flt_t *pflt; // series filters
dly_t *pdlys[CRVA_DLYS]; // array of pointers to delays
mdy_t *pmdlys[CRVA_DLYS]; // array of pointers to mod delays
bool fmoddly; // true if using mod delays
};
rva_t rvas[CRVAS];
void RVA_Init ( rva_t *prva ) { if ( prva ) Q_memset (prva, 0, sizeof (rva_t)); } void RVA_InitAll( void ) { for (int i = 0; i < CRVAS; i++) RVA_Init ( &rvas[i] ); }
// free parallel series reverb
void RVA_Free( rva_t *prva ) { int i;
if ( prva ) { // free all delays
for (i = 0; i < CRVA_DLYS; i++) DLY_Free ( prva->pdlys[i] ); // zero all ptrs to delays in mdy array
for (i = 0; i < CRVA_DLYS; i++) { if ( prva->pmdlys[i] ) prva->pmdlys[i]->pdly = NULL; }
// free all mod delays
for (i = 0; i < CRVA_DLYS; i++) MDY_Free ( prva->pmdlys[i] ); FLT_Free( prva->pflt ); Q_memset( prva, 0, sizeof (rva_t) ); } }
void RVA_FreeAll( void ) { for (int i = 0; i < CRVAS; i++) RVA_Free( &rvas[i] ); }
// create parallel reverb - m parallel reverbs summed
// D array of CRVB_DLYS reverb delay sizes max sample index w[0...D] (ie: D+1 samples)
// a array of reverb feedback parms for parallel reverbs (CRVB_P_DLYS)
// if a[i] < 0 then this is a predelay - use DLY_FLINEAR instead of DLY_LOWPASS
// b array of CRVB_P_DLYS - mix params for parallel reverbs
// m - number of parallel delays
// pflt - filter template, to be used by all parallel delays
// fparallel - true if filter operates in parallel with delays, otherwise filter output only
// fmoddly - > 0 if delays are all mod delays (milliseconds of delay modulation)
// fmodrate - # of delay repetitions between changes to mod delay
// ftaps - if > 0, use 4 taps per reverb delay unit (increases density) tap = D - n*ftaps n = 0,1,2,3
rva_t * RVA_Alloc ( int *D, int *a, int *b, int m, flt_t *pflt, int fparallel, float fmoddly, float fmodrate, float ftaps ) { int i; int dtype; rva_t *prva; flt_t *pflt2 = NULL; bool btaps = ftaps > 0.0;
// find open slot
for ( i = 0; i < CRVAS; i++ ) { if ( !rvas[i].fused ) break; }
// return null if no free slots
if (i == CRVAS) { DevMsg ("DSP: Warning, failed to allocate reverb.\n" ); return NULL; } prva = &rvas[i];
// if series filter specified, alloc two series filters
if ( pflt && !fparallel) { // use filter data as template for a filter on output (2 cascaded filters)
pflt2 = FLT_Alloc (0, pflt->M, pflt->L, pflt->a, pflt->b, 1.0);
if (!pflt2) { DevMsg ("DSP: Warning, failed to allocate flt for reverb.\n" ); return NULL; }
pflt2->pf1 = FLT_Alloc (0, pflt->M, pflt->L, pflt->a, pflt->b, 1.0); pflt2->N = 1; }
// allocate parallel delays
for (i = 0; i < m; i++) { // set delay type
if ( pflt && fparallel ) // if a[i] param is < 0, allocate delay as predelay instead of feedback delay
dtype = a[i] < 0 ? DLY_FLINEAR : DLY_LOWPASS; else // if no filter specified, alloc as plain or multitap plain delay
dtype = btaps ? DLY_PLAIN_4TAP : DLY_PLAIN;
if ( dtype == DLY_LOWPASS && btaps ) dtype = DLY_LOWPASS_4TAP;
// if filter specified and parallel specified, alloc 1 filter per delay
if ( DLY_HAS_FILTER(dtype) ) prva->pdlys[i] = DLY_AllocLP( D[i], abs(a[i]), b[i], dtype, pflt->M, pflt->L, pflt->a, pflt->b ); else prva->pdlys[i] = DLY_Alloc( D[i], abs(a[i]), b[i], dtype );
if ( DLY_HAS_MULTITAP(dtype) ) { // set up delay taps to increase density around delay value.
// value of ftaps is the seed for all tap values
float t1 = max((double)MSEC_TO_SAMPS(5), D[i] * (1.0 - ftaps * 3.141592) ); float t2 = max((double)MSEC_TO_SAMPS(7), D[i] * (1.0 - ftaps * 1.697043) ); float t3 = max((double)MSEC_TO_SAMPS(10), D[i] * (1.0 - ftaps * 0.96325) );
DLY_ChangeTaps( prva->pdlys[i], (int)t1, (int)t2, (int)t3, D[i] ); } }
if ( fmoddly > 0.0 ) { // alloc mod delays, using previously alloc'd delays
// ramptime is time in seconds for delay to change from dcur to dnew
// modtime is time in seconds between modulations. 0 if no self-modulation
// depth is 0-1.0 multiplier, new delay values when modulating are Dnew = randomlong (D - D*depth, D)
float ramptime; float modtime; float depth;
for (i = 0; i < m; i++) { int Do = prva->pdlys[i]->D;
modtime = (float)Do / (float)(SOUND_DMA_SPEED); // seconds per delay
depth = (fmoddly * 0.001f) / modtime; // convert milliseconds to 'depth' %
depth = clamp (depth, 0.01f, 0.99f); modtime = modtime * fmodrate; // modulate every N delay passes
ramptime = fpmin(20.0f/1000.0f, modtime / 2); // ramp between delay values in N ms
prva->pmdlys[i] = MDY_Alloc( prva->pdlys[i], ramptime, modtime, depth, 1.0 ); }
prva->fmoddly = true; } // if we failed to alloc any reverb, free all, return NULL
for (i = 0; i < m; i++) { if ( !prva->pdlys[i] ) { FLT_Free( pflt2 ); RVA_Free( prva ); DevMsg ("DSP: Warning, failed to allocate delay for reverb.\n" ); return NULL; } }
prva->fused = true; prva->m = m; prva->fparallel = fparallel; prva->pflt = pflt2; return prva; }
// parallel reverberator
//
// for each input sample x do:
// x0 = plain(D0,w0,&p0,a0,x)
// x1 = plain(D1,w1,&p1,a1,x)
// x2 = plain(D2,w2,&p2,a2,x)
// x3 = plain(D3,w3,&p3,a3,x)
// y = b0*x0 + b1*x1 + b2*x2 + b3*x3
//
// rgdly - array of M delays:
// D - Delay values (typical - 29, 37, 44, 50, 27, 31)
// w - array of delayed values
// p - array of pointers to circular delay line pointers
// a - array of M feedback values (typical - all equal, like 0.75 * PMAX)
// b - array of M gain values for plain reverb outputs (1, .9, .8, .7)
// xin - input value
// if fparallel, filters are built into delays,
// otherwise, filter is in feedback loop
int g_MapIntoPBITSDivInt[] = { 0, PMAX/1, PMAX/2, PMAX/3, PMAX/4, PMAX/5, PMAX/6, PMAX/7, PMAX/8, PMAX/9, PMAX/10, PMAX/11,PMAX/12,PMAX/13,PMAX/14,PMAX/15,PMAX/16, };
inline int RVA_GetNext( rva_t *prva, int x ) { int m = prva->m; int y = 0;
if ( prva->fmoddly ) { // get output of parallel mod delays
for (int i = 0; i < m; i++ ) y += MDY_GetNext( prva->pmdlys[i], x ); } else { // get output of parallel delays
for (int i = 0; i < m; i++ ) y += DLY_GetNext( prva->pdlys[i], x ); }
// PERFORMANCE: y/m is now baked into the 'b' gain params for each delay ( b = b/m )
// y = (y * g_MapIntoPBITSDivInt[m]) >> PBITS;
if ( prva->fparallel ) return y;
// run series filters if present
if ( prva->pflt ) { y = FLT_GetNext( prva->pflt, y); } return y; }
// batch version for performance
// UNDONE: unwind RVA_GetNextN so that it directly calls DLY_GetNextN or MDY_GetNextN
inline void RVA_GetNextN( rva_t *prva, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = RVA_GetNext( prva, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = RVA_GetNext( prva, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = RVA_GetNext( prva, pb->left ); pb++; } return; } }
// reverb parameter order
typedef enum {
// parameter order
rva_size_max, rva_size_min,
rva_inumdelays, rva_ifeedback, rva_igain, rva_icutoff, rva_ifparallel, rva_imoddly, rva_imodrate,
rva_width, rva_depth, rva_height,
rva_fbwidth, rva_fbdepth, rva_fbheight,
rva_iftaps,
rva_cparam // # of params
} rva_e;
// filter parameter ranges
prm_rng_t rva_rng[] = {
{rva_cparam, 0, 0}, // first entry is # of parameters
// reverb params
{rva_size_max, 0.0, 1000.0}, // max room delay in milliseconds
{rva_size_min, 0.0, 1000.0}, // min room delay in milliseconds
{rva_inumdelays,1.0, 12.0}, // controls # of parallel or series delays
{rva_ifeedback, 0.0, 1.0}, // feedback of delays
{rva_igain, 0.0, 10.0}, // output gain
// filter params for each parallel reverb (quality set to 0 for max execution speed)
{rva_icutoff, 10, 22050}, {rva_ifparallel, 0, 1}, // if 1, then all filters operate in parallel with delays. otherwise filter output only
{rva_imoddly, 0.0, 50.0}, // if > 0 then all delays are modulating delays, mod param controls milliseconds of mod depth
{rva_imodrate, 0.0, 10.0}, // how many delay repetitions pass between mod changes to delayl
// override params - for more detailed description of room
// note: width/depth/height < 0 only for some automatic dsp presets
{rva_width, -1000.0, 1000.0}, // 0-1000.0 millisec (room width in feet) - used instead of size if non-zero
{rva_depth, -1000.0, 1000.0}, // 0-1000.0 room depth in feet - used instead of size if non-zero
{rva_height, -1000.0, 1000.0}, // 0-1000.0 room height in feet - used instead of size if non-zero
{rva_fbwidth, -1.0, 1.0}, // 0-1.0 material reflectivity - used as feedback param instead of decay if non-zero
{rva_fbdepth, -1.0, 1.0}, // 0-1.0 material reflectivity - used as feedback param instead of decay if non-zero
{rva_fbheight, -1.0, 1.0}, // 0-1.0 material reflectivity - used as feedback param instead of decay if non-zero
// if < 0, a predelay is allocated, then feedback is -1*param given
{rva_iftaps, 0.0, 0.333} // if > 0, use 3 extra taps with delay values = d * (1 - faps*n) n = 0,1,2,3
};
#define RVA_BASEM 1 // base number of parallel delays
// nominal delay and feedback values. More delays = more density.
#define RVADLYSMAX 49
float rvadlys[] = {18, 23, 28, 33, 42, 21, 26, 36, 39, 45, 47, 30}; float rvafbs[] = {0.9, 0.9, 0.9, 0.85, 0.8, 0.9, 0.9, 0.85, 0.8, 0.8, 0.8, 0.85};
#define SWAP(a,b,t) {(t) = (a); (a) = (b); (b) = (t);}
#define RVA_MIN_SEPARATION 7 // minimum separation between reverbs, in ms.
// Construct D,a,b delay arrays given array of length,width,height sizes and feedback values
// rgd[] array of delay values in milliseconds (feet)
// rgf[] array of feedback values 0..1
// m # of parallel reverbs to construct
// D[] array of output delay values for parallel reverbs
// a[] array of output feedback values
// b[] array of output gain values = 1/m
// gain - output gain
// feedback - default feedback if rgf members are 0
void RVA_ConstructDelays( float *rgd, float *rgf, int m, int *D, int *a, int *b, float gain, float feedback ) {
int i; float r; int d; float t, d1, d2, dm; bool bpredelay; // sort descending, so rgd[0] is largest delay & rgd[2] is smallest
if (rgd[2] > rgd[1]) { SWAP(rgd[2], rgd[1], t); SWAP(rgf[2], rgf[1], t); } if (rgd[1] > rgd[0]) { SWAP(rgd[0], rgd[1], t); SWAP(rgf[0], rgf[1], t); } if (rgd[2] > rgd[1]) { SWAP(rgd[2], rgd[1], t); SWAP(rgf[2], rgf[1], t); }
// if all feedback values 0, use default feedback
if (rgf[0] == 0.0 && rgf[1] == 0.0 && rgf[2] == 0.0 ) { // use feedback param for all
rgf[0] = rgf[1] = rgf[2] = feedback;
// adjust feedback down for larger delays so that decay is constant for all delays
rgf[0] = DLY_NormalizeFeedback( rgd[2], rgf[2], rgd[0] ); rgf[1] = DLY_NormalizeFeedback( rgd[2], rgf[2], rgd[1] );
} // make sure all reverbs are different by at least RVA_MIN_SEPARATION * m/3 m is 3,6,9 or 12
int dmin = (m/3) * RVA_MIN_SEPARATION;
d1 = rgd[1] - rgd[2]; if (d1 <= dmin) rgd[1] += (dmin-d1); // make difference = dmin
d2 = rgd[0] - rgd[1]; if (d2 <= dmin) rgd[0] += (dmin-d1); // make difference = dmin
for ( i = 0; i < m; i++ ) { // reverberations due to room width, depth, height
// assume sound moves at approx 1ft/ms
int j = (int)(fmod ((float)i, 3.0f)); // j counts 0,1,2 0,1,2 0,1..
d = (int)rgd[j]; r = fabs(rgf[j]);
bpredelay = ((rgf[j] < 0) && i < 3);
// re-use predelay values as reverb values:
if (rgf[j] < 0 && !bpredelay) d = max((int)(rgd[j] / 4.0), RVA_MIN_SEPARATION);
if (i < 3) dm = 0.0; else dm = max( (double)(RVA_MIN_SEPARATION * (i/3)), ((i/3) * ((float)d * 0.18)) );
d += (int)dm; D[i] = MSEC_TO_SAMPS(d);
// D[i] = MSEC_TO_SAMPS(d + ((i/3) * RVA_MIN_SEPARATION)); // (i/3) counts 0,0,0 1,1,1 2,2,2 ... separate all reverbs by 5ms
// feedback - due to wall/floor/ceiling reflectivity
a[i] = (int) min (0.999 * PMAX, (double)PMAX * r);
if (bpredelay) a[i] = -a[i]; // flag delay as predelay
b[i] = (int)((float)(gain * PMAX) / (float)m); } }
void RVA_PerfTest() { double time1, time2;
int i; int k; int j; int m; int a[100];
time1 = Plat_FloatTime();
for (m = 0; m < 1000; m++) { for (i = 0, j = 10000; i < 10000; i++, j--) { // j = j % 6;
// k = (i * j) >> PBITS;
k = i / ((j % 6) + 1); } }
time2 = Plat_FloatTime(); DevMsg("divide = %2.5f \n", (time2-time1));
for (i=1;i<10;i++) a[i] = PMAX / i;
time1 = Plat_FloatTime();
for (m = 0; m < 1000; m++) { for (i = 0, j = 10000; i < 10000; i++, j--) { k = (i * a[(j % 6) + 1] ) >> PBITS; } }
time2 = Plat_FloatTime();
DevMsg("shift & multiply = %2.5f \n", (time2-time1)); }
rva_t * RVA_Params ( prc_t *pprc ) { rva_t *prva;
float size_max = pprc->prm[rva_size_max]; // max delay size
float size_min = pprc->prm[rva_size_min]; // min delay size
float numdelays = pprc->prm[rva_inumdelays]; // controls # of parallel delays
float feedback = pprc->prm[rva_ifeedback]; // 0-1.0 controls feedback parameters
float gain = pprc->prm[rva_igain]; // 0-10.0 controls output gain
float cutoff = pprc->prm[rva_icutoff]; // filter cutoff
float fparallel = pprc->prm[rva_ifparallel]; // if true, all filters are in delay feedback paths - otherwise single flt on output
float fmoddly = pprc->prm[rva_imoddly]; // if > 0, milliseconds of delay mod depth
float fmodrate = pprc->prm[rva_imodrate]; // if fmoddly > 0, # of delay repetitions between modulations
float width = fabs(pprc->prm[rva_width]); // 0-1000 controls size of 1/3 of delays - used instead of size if non-zero
float depth = fabs(pprc->prm[rva_depth]); // 0-1000 controls size of 1/3 of delays - used instead of size if non-zero
float height = fabs(pprc->prm[rva_height]); // 0-1000 controls size of 1/3 of delays - used instead of size if non-zero
float fbwidth = pprc->prm[rva_fbwidth]; // feedback parameter for walls 0..2
float fbdepth = pprc->prm[rva_fbdepth]; // feedback parameter for floor
float fbheight = pprc->prm[rva_fbheight]; // feedback parameter for ceiling
float ftaps = pprc->prm[rva_iftaps]; // if > 0 increase reverb density using 3 extra taps d = (1.0 - ftaps * n) n = 0,1,2,3
// RVA_PerfTest();
// D array of CRVB_DLYS reverb delay sizes max sample index w[0...D] (ie: D+1 samples)
// a array of reverb feedback parms for parallel delays
// b array of CRVB_P_DLYS - mix params for parallel reverbs
// m - number of parallel delays
int D[CRVA_DLYS]; int a[CRVA_DLYS]; int b[CRVA_DLYS]; int m;
// limit # delays 1-12
m = clamp (numdelays, (float)RVA_BASEM, (float)CRVA_DLYS);
// set up D (delay) a (feedback) b (gain) arrays
if ( int(width) || int(height) || int(depth) ) { // if width, height, depth given, use values as simple delays
float rgd[3]; float rgfb[3];
// force m to 3, 6, 9 or 12
if (m < 3) m = 3; if (m > 3 && m < 6) m = 6; if (m > 6 && m < 9) m = 9; if (m > 9) m = 12;
rgd[0] = width; rgfb[0] = fbwidth; rgd[1] = depth; rgfb[1] = fbdepth; rgd[2] = height; rgfb[2] = fbheight;
RVA_ConstructDelays( rgd, rgfb, m, D, a, b, gain, feedback ); } else { // use size parameter instead of width/depth/height
for ( int i = 0; i < m; i++ ) { // delays of parallel reverb. D[0] = size_min.
D[i] = MSEC_TO_SAMPS( size_min + (int)( ((float)(size_max - size_min) / (float)m) * (float)i) ); // feedback and gain of parallel reverb
if (i == 0) { // set feedback for smallest delay
a[i] = (int) min (0.999 * PMAX, (double)PMAX * feedback ); } else { // adjust feedback down for larger delays so that decay time is constant
a[i] = (int) min (0.999 * PMAX, (double)PMAX * DLY_NormalizeFeedback( D[0], feedback, D[i] ) ); } b[i] = (int) ((float)(gain * PMAX) / (float)m); } }
// add filter
flt_t *pflt = NULL;
if ( cutoff ) {
// set up dummy lowpass filter to convert params
prc_t prcf;
prcf.prm[flt_iquality] = QUA_LO; // force filter to low quality for faster execution time
prcf.prm[flt_icutoff] = cutoff; prcf.prm[flt_iftype] = FLT_LP; prcf.prm[flt_iqwidth] = 0; prcf.prm[flt_igain] = 1.0;
pflt = (flt_t *)FLT_Params ( &prcf ); } prva = RVA_Alloc ( D, a, b, m, pflt, fparallel, fmoddly, fmodrate, ftaps );
FLT_Free( pflt );
return prva; }
inline void * RVA_VParams ( void *p ) { PRC_CheckParams ( (prc_t *)p, rva_rng ); return (void *) RVA_Params ((prc_t *)p); }
inline void RVA_Mod ( void *p, float v ) { return; }
////////////
// Diffusor
///////////
// (N series allpass reverbs)
#define CDFRS 64 // max number of series reverbs active
#define CDFR_DLYS 16 // max number of delays making up diffusor
struct dfr_t { bool fused; int n; // series allpass delays
int w[CDFR_DLYS]; // internal state array for series allpass filters
dly_t *pdlys[CDFR_DLYS]; // array of pointers to delays
};
dfr_t dfrs[CDFRS];
void DFR_Init ( dfr_t *pdfr ) { if ( pdfr ) Q_memset (pdfr, 0, sizeof (dfr_t)); } void DFR_InitAll( void ) { for (int i = 0; i < CDFRS; i++) DFR_Init ( &dfrs[i] ); }
// free parallel series reverb
void DFR_Free( dfr_t *pdfr ) { if ( pdfr ) { // free all delays
for (int i = 0; i < CDFR_DLYS; i++) DLY_Free ( pdfr->pdlys[i] ); Q_memset( pdfr, 0, sizeof (dfr_t) ); } }
void DFR_FreeAll( void ) { for (int i = 0; i < CDFRS; i++) DFR_Free( &dfrs[i] ); }
// create n series allpass reverbs
// D array of CRVB_DLYS reverb delay sizes max sample index w[0...D] (ie: D+1 samples)
// a array of reverb feedback parms for series delays
// b array of gain params for parallel reverbs
// n - number of series delays
dfr_t * DFR_Alloc ( int *D, int *a, int *b, int n ) { int i; dfr_t *pdfr;
// find open slot
for (i = 0; i < CDFRS; i++) { if (!dfrs[i].fused) break; } // return null if no free slots
if (i == CDFRS) { DevMsg ("DSP: Warning, failed to allocate diffusor.\n" ); return NULL; } pdfr = &dfrs[i];
DFR_Init( pdfr );
// alloc reverbs
for (i = 0; i < n; i++) pdfr->pdlys[i] = DLY_Alloc( D[i], a[i], b[i], DLY_ALLPASS ); // if we failed to alloc any reverb, free all, return NULL
for (i = 0; i < n; i++) { if ( !pdfr->pdlys[i]) { DFR_Free( pdfr ); DevMsg ("DSP: Warning, failed to allocate delay for diffusor.\n" ); return NULL; } } pdfr->fused = true; pdfr->n = n;
return pdfr; }
// series reverberator
inline int DFR_GetNext( dfr_t *pdfr, int x ) { int i; int y; dly_t *pdly;
y = x; for (i = 0; i < pdfr->n; i++) { pdly = pdfr->pdlys[i]; y = DelayAllpass( pdly->D, pdly->t, pdly->w, &pdly->p, pdly->a, pdly->b, y ); }
return y; }
// batch version for performance
inline void DFR_GetNextN( dfr_t *pdfr, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = DFR_GetNext( pdfr, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = DFR_GetNext( pdfr, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = DFR_GetNext( pdfr, pb->left ); pb++; } return; } }
#define DFR_BASEN 1 // base number of series allpass delays
// nominal diffusor delay and feedback values
float dfrdlys[] = {13, 19, 26, 21, 32, 36, 38, 16, 24, 28, 41, 35, 10, 46, 50, 27}; float dfrfbs[] = {1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0};
// diffusor parameter order
typedef enum {
// parameter order
dfr_isize, dfr_inumdelays, dfr_ifeedback, dfr_igain,
dfr_cparam // # of params
} dfr_e;
// diffusor parameter ranges
prm_rng_t dfr_rng[] = {
{dfr_cparam, 0, 0}, // first entry is # of parameters
{dfr_isize, 0.0, 1.0}, // 0-1.0 scales all delays
{dfr_inumdelays,0.0, 4.0}, // 0-4.0 controls # of series delays
{dfr_ifeedback, 0.0, 1.0}, // 0-1.0 scales all feedback parameters
{dfr_igain, 0.0, 10.0}, // 0-1.0 scales all feedback parameters
};
dfr_t * DFR_Params ( prc_t *pprc ) { dfr_t *pdfr; int i; int s; float size = pprc->prm[dfr_isize]; // 0-1.0 scales all delays
float numdelays = pprc->prm[dfr_inumdelays]; // 0-4.0 controls # of series delays
float feedback = pprc->prm[dfr_ifeedback]; // 0-1.0 scales all feedback parameters
float gain = pprc->prm[dfr_igain]; // 0-10.0 controls output gain
// D array of CRVB_DLYS reverb delay sizes max sample index w[0...D] (ie: D+1 samples)
// a array of reverb feedback parms for series delays (CRVB_S_DLYS)
// b gain of each reverb section
// n - number of series delays
int D[CDFR_DLYS]; int a[CDFR_DLYS]; int b[CDFR_DLYS]; int n;
if (gain == 0.0) gain = 1.0;
// get # series diffusors
// limit m, n to half max number of delays
n = clamp (Float2Int(numdelays), DFR_BASEN, CDFR_DLYS/2);
// compute delays for diffusors
for (i = 0; i < n; i++) { s = (int)( dfrdlys[i] * size );
// delay of diffusor
D[i] = MSEC_TO_SAMPS(s);
// feedback and gain of diffusor
a[i] = min (0.999 * PMAX, (double)(dfrfbs[i] * PMAX * feedback)); b[i] = (int) ( (float)(gain * (float)PMAX) ); }
pdfr = DFR_Alloc ( D, a, b, n );
return pdfr; }
inline void * DFR_VParams ( void *p ) { PRC_CheckParams ((prc_t *)p, dfr_rng); return (void *) DFR_Params ((prc_t *)p); }
inline void DFR_Mod ( void *p, float v ) { return; }
//////////////////////
// LFO wav definitions
//////////////////////
#define CLFOSAMPS 512 // samples per wav table - single cycle only
#define LFOBITS 14 // bits of peak amplitude of lfo wav
#define LFOAMP ((1<<LFOBITS)-1) // peak amplitude of lfo wav
//types of lfo wavs
#define LFO_SIN 0 // sine wav
#define LFO_TRI 1 // triangle wav
#define LFO_SQR 2 // square wave, 50% duty cycle
#define LFO_SAW 3 // forward saw wav
#define LFO_RND 4 // random wav
#define LFO_LOG_IN 5 // logarithmic fade in
#define LFO_LOG_OUT 6 // logarithmic fade out
#define LFO_LIN_IN 7 // linear fade in
#define LFO_LIN_OUT 8 // linear fade out
#define LFO_MAX LFO_LIN_OUT
#define CLFOWAV 9 // number of LFO wav tables
struct lfowav_t // lfo or envelope wave table
{ int type; // lfo type
dly_t *pdly; // delay holds wav values and step pointers
};
lfowav_t lfowavs[CLFOWAV];
// deallocate lfo wave table. Called only when sound engine exits.
void LFOWAV_Free( lfowav_t *plw ) { // free delay
if ( plw ) DLY_Free( plw->pdly );
Q_memset( plw, 0, sizeof (lfowav_t) ); }
// deallocate all lfo wave tables. Called only when sound engine exits.
void LFOWAV_FreeAll( void ) { for ( int i = 0; i < CLFOWAV; i++ ) LFOWAV_Free( &lfowavs[i] ); }
// fill lfo array w with count samples of lfo type 'type'
// all lfo wavs except fade out, rnd, and log_out should start with 0 output
void LFOWAV_Fill( int *w, int count, int type ) { int i,x; switch (type) { default: case LFO_SIN: // sine wav, all values 0 <= x <= LFOAMP, initial value = 0
for (i = 0; i < count; i++ ) { x = ( int )( (float)(LFOAMP) * sinf( (2.0 * M_PI_F * (float)i / (float)count ) + (M_PI_F * 1.5) ) ); w[i] = (x + LFOAMP)/2; } break; case LFO_TRI: // triangle wav, all values 0 <= x <= LFOAMP, initial value = 0
for (i = 0; i < count; i++) { w[i] = ( int ) ( (float)(2 * LFOAMP * i ) / (float)(count) ); if ( i > count / 2 ) w[i] = ( int ) ( (float) (2 * LFOAMP) - (float)( 2 * LFOAMP * i ) / (float)(count) ); } break; case LFO_SQR: // square wave, 50% duty cycle, all values 0 <= x <= LFOAMP, initial value = 0
for (i = 0; i < count; i++) w[i] = i > count / 2 ? 0 : LFOAMP; break; case LFO_SAW: // forward saw wav, aall values 0 <= x <= LFOAMP, initial value = 0
for (i = 0; i < count; i++) w[i] = ( int ) ( (float)(LFOAMP) * (float)i / (float)(count) ); break; case LFO_RND: // random wav, all values 0 <= x <= LFOAMP
for (i = 0; i < count; i++) w[i] = ( int ) ( RandomInt(0, LFOAMP) ); break; case LFO_LOG_IN: // logarithmic fade in, all values 0 <= x <= LFOAMP, initial value = 0
for (i = 0; i < count; i++) w[i] = ( int ) ( (float)(LFOAMP) * powf( (float)i / (float)count, 2)); break; case LFO_LOG_OUT: // logarithmic fade out, all values 0 <= x <= LFOAMP, initial value = LFOAMP
for (i = 0; i < count; i++) w[i] = ( int ) ( (float)(LFOAMP) * powf( 1.0 - ((float)i / (float)count), 2 )); break; case LFO_LIN_IN: // linear fade in, all values 0 <= x <= LFOAMP, initial value = 0
for (i = 0; i < count; i++) w[i] = ( int ) ( (float)(LFOAMP) * (float)i / (float)(count) ); break; case LFO_LIN_OUT: // linear fade out, all values 0 <= x <= LFOAMP, initial value = LFOAMP
for (i = 0; i < count; i++) w[i] = LFOAMP - ( int ) ( (float)(LFOAMP) * (float)i / (float)(count) ); break; } }
// allocate all lfo wave tables. Called only when sound engine loads.
void LFOWAV_InitAll() { int i; dly_t *pdly;
Q_memset( lfowavs, 0, sizeof( lfowavs ) );
// alloc space for each lfo wav type
for (i = 0; i < CLFOWAV; i++) { pdly = DLY_Alloc( CLFOSAMPS, 0, 0 , DLY_PLAIN); lfowavs[i].pdly = pdly; lfowavs[i].type = i;
LFOWAV_Fill( pdly->w, CLFOSAMPS, i ); } // if any dlys fail to alloc, free all
for (i = 0; i < CLFOWAV; i++) { if ( !lfowavs[i].pdly ) LFOWAV_FreeAll(); } }
////////////////////////////////////////
// LFO iterators - one shot and looping
////////////////////////////////////////
#define CLFO 16 // max active lfos (this steals from active delays)
struct lfo_t { bool fused; // true if slot take
dly_t *pdly; // delay points to lfo wav within lfowav_t (don't free this)
int gain;
float f; // playback frequency in hz
pos_t pos; // current position within wav table, looping
pos_one_t pos1; // current position within wav table, one shot
int foneshot; // true - one shot only, don't repeat
};
lfo_t lfos[CLFO];
void LFO_Init( lfo_t *plfo ) { if ( plfo ) Q_memset( plfo, 0, sizeof (lfo_t) ); } void LFO_InitAll( void ) { for (int i = 0; i < CLFO; i++) LFO_Init(&lfos[i]); } void LFO_Free( lfo_t *plfo ) { if ( plfo ) Q_memset( plfo, 0, sizeof (lfo_t) ); } void LFO_FreeAll( void ) { for (int i = 0; i < CLFO; i++) LFO_Free(&lfos[i]); }
// get step value given desired playback frequency
inline float LFO_HzToStep ( float freqHz ) { float lfoHz;
// calculate integer and fractional step values,
// assume an update rate of SOUND_DMA_SPEED samples/sec
// 1 cycle/CLFOSAMPS * SOUND_DMA_SPEED samps/sec = cycles/sec = current lfo rate
//
// lforate * X = freqHz so X = freqHz/lforate = update rate
lfoHz = (float)(SOUND_DMA_SPEED) / (float)(CLFOSAMPS);
return freqHz / lfoHz; }
// return pointer to new lfo
lfo_t * LFO_Alloc( int wtype, float freqHz, bool foneshot, float gain ) { int i; int type = min ( CLFOWAV - 1, wtype ); float lfostep; for (i = 0; i < CLFO; i++) if (!lfos[i].fused) { lfo_t *plfo = &lfos[i];
LFO_Init( plfo );
plfo->fused = true; plfo->pdly = lfowavs[type].pdly; // pdly in lfo points to wav table data in lfowavs
plfo->f = freqHz; plfo->foneshot = foneshot; plfo->gain = gain * PMAX;
lfostep = LFO_HzToStep( freqHz );
// init positional pointer (ie: fixed point updater for controlling pitch of lfo)
if ( !foneshot ) POS_Init(&(plfo->pos), plfo->pdly->D, lfostep ); else POS_ONE_Init(&(plfo->pos1), plfo->pdly->D,lfostep );
return plfo; } DevMsg ("DSP: Warning, failed to allocate LFO.\n" ); return NULL; }
// get next lfo value
// Value returned is 0..LFOAMP. can be normalized by shifting right by LFOBITS
// To play back at correct passed in frequency, routien should be
// called once for every output sample (ie: at SOUND_DMA_SPEED)
// x is dummy param
inline int LFO_GetNext( lfo_t *plfo, int x ) { int i;
// get current position
if ( !plfo->foneshot ) i = POS_GetNext( &plfo->pos ); else i = POS_ONE_GetNext( &plfo->pos1 );
// return current sample
if (plfo->gain == PMAX) return plfo->pdly->w[i]; else return (plfo->pdly->w[i] * plfo->gain ) >> PBITS; }
// batch version for performance
inline void LFO_GetNextN( lfo_t *plfo, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = LFO_GetNext( plfo, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = LFO_GetNext( plfo, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = LFO_GetNext( plfo, pb->left ); pb++; } return; } }
// uses lfowav, rate, foneshot
typedef enum {
// parameter order
lfo_iwav, lfo_irate, lfo_ifoneshot, lfo_igain, lfo_cparam // # of params
} lfo_e;
// parameter ranges
prm_rng_t lfo_rng[] = {
{lfo_cparam, 0, 0}, // first entry is # of parameters
{lfo_iwav, 0.0, LFO_MAX}, // lfo type to use (LFO_SIN, LFO_RND...)
{lfo_irate, 0.0, 16000.0}, // modulation rate in hz. for MDY, 1/rate = 'glide' time in seconds
{lfo_ifoneshot, 0.0, 1.0}, // 1.0 if lfo is oneshot
{lfo_igain, 0.0, 10.0}, // output gain
};
lfo_t * LFO_Params ( prc_t *pprc ) { lfo_t *plfo; bool foneshot = pprc->prm[lfo_ifoneshot] > 0 ? true : false; float gain = pprc->prm[lfo_igain];
plfo = LFO_Alloc ( pprc->prm[lfo_iwav], pprc->prm[lfo_irate], foneshot, gain );
return plfo; }
void LFO_ChangeVal ( lfo_t *plfo, float fhz ) { float fstep = LFO_HzToStep( fhz );
// change lfo playback rate to new frequency fhz
if ( plfo->foneshot ) POS_ChangeVal( &plfo->pos, fstep ); else POS_ChangeVal( &plfo->pos1.p, fstep ); }
inline void * LFO_VParams ( void *p ) { PRC_CheckParams ( (prc_t *)p, lfo_rng ); return (void *) LFO_Params ((prc_t *)p); }
// v is +/- 0-1.0
// v changes current lfo frequency up/down by +/- v%
inline void LFO_Mod ( lfo_t *plfo, float v ) { float fhz; float fhznew;
fhz = plfo->f; fhznew = fhz * (1.0 + v);
LFO_ChangeVal ( plfo, fhznew );
return; }
////////////////////////////////////////
// Time Compress/expand with pitch shift
////////////////////////////////////////
// realtime pitch shift - ie: pitch shift without change to playback rate
#define CPTCS 64
struct ptc_t { bool fused; dly_t *pdly_in; // input buffer space
dly_t *pdly_out; // output buffer space
int *pin; // input buffer (pdly_in->w)
int *pout; // output buffer (pdly_out->w)
int cin; // # samples in input buffer
int cout; // # samples in output buffer
int cxfade; // # samples in crossfade segment
int ccut; // # samples to cut
int cduplicate; // # samples to duplicate (redundant - same as ccut)
int iin; // current index into input buffer (reading)
pos_one_t psn; // stepping index through output buffer
bool fdup; // true if duplicating, false if cutting
float fstep; // pitch shift & time compress/expand
};
ptc_t ptcs[CPTCS];
void PTC_Init( ptc_t *pptc ) { if (pptc) Q_memset( pptc, 0, sizeof (ptc_t) ); }; void PTC_Free( ptc_t *pptc ) { if (pptc) { DLY_Free (pptc->pdly_in); DLY_Free (pptc->pdly_out);
Q_memset( pptc, 0, sizeof (ptc_t) ); } }; void PTC_InitAll() { for (int i = 0; i < CPTCS; i++) PTC_Init( &ptcs[i] ); }; void PTC_FreeAll() { for (int i = 0; i < CPTCS; i++) PTC_Free( &ptcs[i] ); };
// Time compressor/expander with pitch shift (ie: pitch changes, playback rate does not)
//
// Algorithm:
// 1) Duplicate or discard chunks of sound to provide tslice * fstep seconds of sound.
// (The user-selectable size of the buffer to process is tslice milliseconds in length)
// 2) Resample this compressed/expanded buffer at fstep to produce a pitch shifted
// output with the same duration as the input (ie: #samples out = # samples in, an
// obvious requirement for realtime inline processing).
// timeslice is size in milliseconds of full buffer to process.
// timeslice * fstep is the size of the expanded/compressed buffer
// timexfade is length in milliseconds of crossfade region between duplicated or cut sections
// fstep is % expanded/compressed sound normalized to 0.01-2.0 (1% - 200%)
// input buffer:
// iin-->
// [0... tslice ...D] input samples 0...D (D is NEWEST sample)
// [0... ...n][m... tseg ...D] region to be cut or duplicated m...D
// [0... [p..txf1..n][m... tseg ...D] fade in region 1 txf1 p...n
// [0... ...n][m..[q..txf2..D] fade out region 2 txf2 q...D
// pitch up: duplicate into output buffer: tdup = tseg
// [0... ...n][m... tdup ...D][m... tdup ...D] output buffer size with duplicate region
// [0... ...n][m..[p...xf1..n][m... tdup ...D] fade in p...n while fading out q...D
// [0... ...n][m..[q...xf2..D][m... tdup ...D]
// [0... ...n][m..[.XFADE...n][m... tdup ...D] final duplicated output buffer - resample at fstep
// pitch down: cut into output buffer: tcut = tseg
// [0... ...n][m... tcut ...D] input samples with cut region delineated m...D
// [0... ...n] output buffer size after cut
// [0... [q..txf2...D] fade in txf1 q...D while fade out txf2 p...n
// [0... [.XFADE ...D] final cut output buffer - resample at fstep
ptc_t * PTC_Alloc( float timeslice, float timexfade, float fstep ) { int i; ptc_t *pptc; float tout; int cin, cout; float tslice = timeslice; float txfade = timexfade; float tcutdup;
// find time compressor slot
for ( i = 0; i < CPTCS; i++ ) { if ( !ptcs[i].fused ) break; } if ( i == CPTCS ) { DevMsg ("DSP: Warning, failed to allocate pitch shifter.\n" ); return NULL; }
pptc = &ptcs[i]; PTC_Init ( pptc );
// get size of region to cut or duplicate
tcutdup = abs((fstep - 1.0) * timeslice);
// to prevent buffer overruns:
// make sure timeslice is greater than cut/dup time
tslice = max ( (double)tslice, 1.1 * tcutdup);
// make sure xfade time smaller than cut/dup time, and smaller than (timeslice-cutdup) time
txfade = min ( (double)txfade, 0.9 * tcutdup ); txfade = min ( (double)txfade, 0.9 * (tslice - tcutdup));
pptc->cxfade = MSEC_TO_SAMPS( txfade ); pptc->ccut = MSEC_TO_SAMPS( tcutdup ); pptc->cduplicate = MSEC_TO_SAMPS( tcutdup ); // alloc delay lines (buffers)
tout = tslice * fstep;
cin = MSEC_TO_SAMPS( tslice ); cout = MSEC_TO_SAMPS( tout ); pptc->pdly_in = DLY_Alloc( cin, 0, 1, DLY_LINEAR ); // alloc input buffer
pptc->pdly_out = DLY_Alloc( cout, 0, 1, DLY_LINEAR); // alloc output buffer
if ( !pptc->pdly_in || !pptc->pdly_out ) { PTC_Free( pptc ); DevMsg ("DSP: Warning, failed to allocate delay for pitch shifter.\n" ); return NULL; }
// buffer pointers
pptc->pin = pptc->pdly_in->w; pptc->pout = pptc->pdly_out->w;
// input buffer index
pptc->iin = 0;
// output buffer index
POS_ONE_Init ( &pptc->psn, cout, fstep );
// if fstep > 1.0 we're pitching shifting up, so fdup = true
pptc->fdup = fstep > 1.0 ? true : false; pptc->cin = cin; pptc->cout = cout;
pptc->fstep = fstep; pptc->fused = true;
return pptc; }
// linear crossfader
// yfadein - instantaneous value fading in
// ydafeout -instantaneous value fading out
// nsamples - duration in #samples of fade
// isample - index in to fade 0...nsamples-1
inline int xfade ( int yfadein, int yfadeout, int nsamples, int isample ) { int yout; int m = (isample << PBITS ) / nsamples;
// yout = ((yfadein * m) >> PBITS) + ((yfadeout * (PMAX - m)) >> PBITS);
yout = (yfadeout + (yfadein - yfadeout) * m ) >> PBITS;
return yout; }
// w - pointer to start of input buffer samples
// v - pointer to start of output buffer samples
// cin - # of input buffer samples
// cout = # of output buffer samples
// cxfade = # of crossfade samples
// cduplicate = # of samples in duplicate/cut segment
void TimeExpand( int *w, int *v, int cin, int cout, int cxfade, int cduplicate ) { int i,j; int m; int p; int q; int D;
// input buffer
// xfade source duplicate
// [0...........][p.......n][m...........D]
// output buffer
// xfade region duplicate
// [0.....................n][m..[q.......D][m...........D]
// D - index of last sample in input buffer
// m - index of 1st sample in duplication region
// p - index of 1st sample of crossfade source
// q - index of 1st sample in crossfade region
D = cin - 1; m = cin - cduplicate; p = m - cxfade; q = cin - cxfade;
// copy up to crossfade region
for (i = 0; i < q; i++) v[i] = w[i]; // crossfade region
j = p;
for (i = q; i <= D; i++) v[i] = xfade (w[j++], w[i], cxfade, i-q); // fade out p..n, fade in q..D
// duplicate region
j = D+1;
for (i = m; i <= D; i++) v[j++] = w[i];
}
// cut ccut samples from end of input buffer, crossfade end of cut section
// with end of remaining section
// w - pointer to start of input buffer samples
// v - pointer to start of output buffer samples
// cin - # of input buffer samples
// cout = # of output buffer samples
// cxfade = # of crossfade samples
// ccut = # of samples in cut segment
void TimeCompress( int *w, int *v, int cin, int cout, int cxfade, int ccut ) { int i,j; int m; int p; int q; int D;
// input buffer
// xfade source
// [0.....................n][m..[p.......D]
// xfade region cut
// [0...........][q.......n][m...........D]
// output buffer
// xfade to source
// [0...........][p.......D]
// D - index of last sample in input buffer
// m - index of 1st sample in cut region
// p - index of 1st sample of crossfade source
// q - index of 1st sample in crossfade region
D = cin - 1; m = cin - ccut; p = cin - cxfade; q = m - cxfade;
// copy up to crossfade region
for (i = 0; i < q; i++) v[i] = w[i]; // crossfade region
j = p;
for (i = q; i < m; i++) v[i] = xfade (w[j++], w[i], cxfade, i-q); // fade out p..n, fade in q..D
// skip rest of input buffer
}
// get next sample
// put input sample into input (delay) buffer
// get output sample from output buffer, step by fstep %
// output buffer is time expanded or compressed version of previous input buffer
inline int PTC_GetNext( ptc_t *pptc, int x ) { int iout, xout; bool fhitend = false;
// write x into input buffer
Assert (pptc->iin < pptc->cin);
pptc->pin[pptc->iin] = x;
pptc->iin++; // check for end of input buffer
if ( pptc->iin >= pptc->cin ) fhitend = true;
// read sample from output buffer, resampling at fstep
iout = POS_ONE_GetNext( &pptc->psn ); Assert (iout < pptc->cout); xout = pptc->pout[iout]; if ( fhitend ) { // if hit end of input buffer (ie: input buffer is full)
// reset input buffer pointer
// reset output buffer pointer
// rebuild entire output buffer (TimeCompress/TimeExpand)
pptc->iin = 0;
POS_ONE_Init( &pptc->psn, pptc->cout, pptc->fstep );
if ( pptc->fdup ) TimeExpand ( pptc->pin, pptc->pout, pptc->cin, pptc->cout, pptc->cxfade, pptc->cduplicate ); else TimeCompress ( pptc->pin, pptc->pout, pptc->cin, pptc->cout, pptc->cxfade, pptc->ccut ); }
return xout; }
// batch version for performance
inline void PTC_GetNextN( ptc_t *pptc, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = PTC_GetNext( pptc, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = PTC_GetNext( pptc, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = PTC_GetNext( pptc, pb->left ); pb++; } return; } }
// change time compression to new value
// fstep is new value
// ramptime is how long change takes in seconds (ramps smoothly), 0 for no ramp
void PTC_ChangeVal( ptc_t *pptc, float fstep, float ramptime ) { // UNDONE: ignored
// UNDONE: just realloc time compressor with new fstep
}
// uses pitch:
// 1.0 = playback normal rate
// 0.5 = cut 50% of sound (2x playback)
// 1.5 = add 50% sound (0.5x playback)
typedef enum {
// parameter order
ptc_ipitch, ptc_itimeslice, ptc_ixfade, ptc_cparam // # of params
} ptc_e;
// diffusor parameter ranges
prm_rng_t ptc_rng[] = {
{ptc_cparam, 0, 0}, // first entry is # of parameters
{ptc_ipitch, 0.1, 4.0}, // 0-n.0 where 1.0 = 1 octave up and 0.5 is one octave down
{ptc_itimeslice, 20.0, 300.0}, // in milliseconds - size of sound chunk to analyze and cut/duplicate - 100ms nominal
{ptc_ixfade, 1.0, 200.0}, // in milliseconds - size of crossfade region between spliced chunks - 20ms nominal
};
ptc_t * PTC_Params ( prc_t *pprc ) { ptc_t *pptc;
float pitch = pprc->prm[ptc_ipitch]; float timeslice = pprc->prm[ptc_itimeslice]; float txfade = pprc->prm[ptc_ixfade];
pptc = PTC_Alloc( timeslice, txfade, pitch ); return pptc; }
inline void * PTC_VParams ( void *p ) { PRC_CheckParams ( (prc_t *)p, ptc_rng ); return (void *) PTC_Params ((prc_t *)p); }
// change to new pitch value
// v is +/- 0-1.0
// v changes current pitch up/down by +/- v%
void PTC_Mod ( ptc_t *pptc, float v ) { float fstep; float fstepnew;
fstep = pptc->fstep; fstepnew = fstep * (1.0 + v);
PTC_ChangeVal( pptc, fstepnew, 0.01 ); }
////////////////////
// ADSR envelope
////////////////////
#define CENVS 64 // max # of envelopes active
#define CENVRMPS 4 // A, D, S, R
#define ENV_LIN 0 // linear a,d,s,r
#define ENV_EXP 1 // exponential a,d,s,r
#define ENV_MAX ENV_EXP
#define ENV_BITS 14 // bits of resolution of ramp
struct env_t { bool fused;
bool fhitend; // true if done
bool fexp; // true if exponential ramps
int ienv; // current ramp
rmp_t rmps[CENVRMPS]; // ramps
};
env_t envs[CENVS];
void ENV_Init( env_t *penv ) { if (penv) Q_memset( penv, 0, sizeof (env_t) ); }; void ENV_Free( env_t *penv ) { if (penv) Q_memset( penv, 0, sizeof (env_t) ); }; void ENV_InitAll() { for (int i = 0; i < CENVS; i++) ENV_Init( &envs[i] ); }; void ENV_FreeAll() { for (int i = 0; i < CENVS; i++) ENV_Free( &envs[i] ); };
// allocate ADSR envelope
// all times are in seconds
// amp1 - attack amplitude multiplier 0-1.0
// amp2 - sustain amplitude multiplier 0-1.0
// amp3 - end of sustain amplitude multiplier 0-1.0
env_t *ENV_Alloc ( int type, float famp1, float famp2, float famp3, float attack, float decay, float sustain, float release, bool fexp) { int i; env_t *penv;
for (i = 0; i < CENVS; i++) { if ( !envs[i].fused ) {
int amp1 = famp1 * (1 << ENV_BITS); // ramp resolution
int amp2 = famp2 * (1 << ENV_BITS); int amp3 = famp3 * (1 << ENV_BITS);
penv = &envs[i]; ENV_Init (penv);
// UNDONE: ignoring type = ENV_EXP - use oneshot LFOS instead with sawtooth/exponential
// set up ramps
RMP_Init( &penv->rmps[0], attack, 0, amp1, true ); RMP_Init( &penv->rmps[1], decay, amp1, amp2, true ); RMP_Init( &penv->rmps[2], sustain, amp2, amp3, true ); RMP_Init( &penv->rmps[3], release, amp3, 0, true );
penv->ienv = 0; penv->fused = true; penv->fhitend = false; penv->fexp = fexp; return penv; } } DevMsg ("DSP: Warning, failed to allocate envelope.\n" ); return NULL; }
inline int ENV_GetNext( env_t *penv, int x ) { if ( !penv->fhitend ) { int i; int y;
i = penv->ienv; y = RMP_GetNext ( &penv->rmps[i] ); // check for next ramp
if ( penv->rmps[i].fhitend ) i++;
penv->ienv = i;
// check for end of all ramps
if ( i > 3) penv->fhitend = true;
// multiply input signal by ramp
if (penv->fexp) return (((x * y) >> ENV_BITS) * y) >> ENV_BITS; else return (x * y) >> ENV_BITS; }
return 0; }
// batch version for performance
inline void ENV_GetNextN( env_t *penv, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = ENV_GetNext( penv, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = ENV_GetNext( penv, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = ENV_GetNext( penv, pb->left ); pb++; } return; } }
// uses lfowav, amp1, amp2, amp3, attack, decay, sustain, release
// lfowav is type, currently ignored - ie: LFO_LIN_IN, LFO_LOG_IN
// parameter order
typedef enum { env_itype, env_iamp1, env_iamp2, env_iamp3, env_iattack, env_idecay, env_isustain, env_irelease, env_ifexp, env_cparam // # of params
} env_e;
// parameter ranges
prm_rng_t env_rng[] = {
{env_cparam, 0, 0}, // first entry is # of parameters
{env_itype, 0.0,ENV_MAX}, // ENV_LINEAR, ENV_LOG - currently ignored
{env_iamp1, 0.0, 1.0}, // attack peak amplitude 0-1.0
{env_iamp2, 0.0, 1.0}, // decay target amplitued 0-1.0
{env_iamp3, 0.0, 1.0}, // sustain target amplitude 0-1.0
{env_iattack, 0.0, 20000.0}, // attack time in milliseconds
{env_idecay, 0.0, 20000.0}, // envelope decay time in milliseconds
{env_isustain, 0.0, 20000.0}, // sustain time in milliseconds
{env_irelease, 0.0, 20000.0}, // release time in milliseconds
{env_ifexp, 0.0, 1.0}, // 1.0 if exponential ramps
};
env_t * ENV_Params ( prc_t *pprc ) { env_t *penv;
float type = pprc->prm[env_itype]; float amp1 = pprc->prm[env_iamp1]; float amp2 = pprc->prm[env_iamp2]; float amp3 = pprc->prm[env_iamp3]; float attack = pprc->prm[env_iattack]/1000.0; float decay = pprc->prm[env_idecay]/1000.0; float sustain = pprc->prm[env_isustain]/1000.0; float release = pprc->prm[env_irelease]/1000.0; float fexp = pprc->prm[env_ifexp]; bool bexp;
bexp = fexp > 0.0 ? 1 : 0; penv = ENV_Alloc ( type, amp1, amp2, amp3, attack, decay, sustain, release, bexp ); return penv; }
inline void * ENV_VParams ( void *p ) { PRC_CheckParams( (prc_t *)p, env_rng ); return (void *) ENV_Params ((prc_t *)p); }
inline void ENV_Mod ( void *p, float v ) { return; }
//////////////////////////
// Gate & envelope follower
//////////////////////////
#define CEFOS 64 // max # of envelope followers active
struct efo_t { bool fused;
int xout; // current output value
// gate params
bool bgate; // if true, gate function is on
bool bgateon; // if true, gate is on
bool bexp; // if true, use exponential fade out
int thresh; // amplitude threshold for gate on
int thresh_off; // amplitidue threshold for gate off
float attack_time; // gate attack time in seconds
float decay_time; // gate decay time in seconds
rmp_t rmp_attack; // gate on ramp - attack
rmp_t rmp_decay; // gate off ramp - decay
};
efo_t efos[CEFOS];
void EFO_Init( efo_t *pefo ) { if (pefo) Q_memset( pefo, 0, sizeof (efo_t) ); }; void EFO_Free( efo_t *pefo ) { if (pefo) Q_memset( pefo, 0, sizeof (efo_t) ); }; void EFO_InitAll() { for (int i = 0; i < CEFOS; i++) EFO_Init( &efos[i] ); }; void EFO_FreeAll() { for (int i = 0; i < CEFOS; i++) EFO_Free( &efos[i] ); };
// return true when gate is off AND decay ramp has hit end
inline bool EFO_GateOff( efo_t *pefo ) { return ( !pefo->bgateon && RMP_HitEnd( &pefo->rmp_decay ) ); }
// allocate enveloper follower
#define EFO_HYST_AMP 1000 // hysteresis amplitude
efo_t *EFO_Alloc ( float threshold, float attack_sec, float decay_sec, bool bexp ) { int i; efo_t *pefo;
for (i = 0; i < CEFOS; i++) { if ( !efos[i].fused ) { pefo = &efos[i];
EFO_Init ( pefo );
pefo->xout = 0; pefo->fused = true; // init gate params
pefo->bgate = threshold > 0.0;
if (pefo->bgate) { pefo->attack_time = attack_sec; pefo->decay_time = decay_sec;
RMP_Init( &pefo->rmp_attack, attack_sec, 0, PMAX, false); RMP_Init( &pefo->rmp_decay, decay_sec, PMAX, 0, false); RMP_SetEnd( &pefo->rmp_attack ); RMP_SetEnd( &pefo->rmp_decay );
pefo->thresh = threshold; pefo->thresh_off = max(1.f, threshold - EFO_HYST_AMP); pefo->bgateon = false; pefo->bexp = bexp; } return pefo; } }
DevMsg ("DSP: Warning, failed to allocate envelope follower.\n" ); return NULL; }
// values of L for CEFO_BITS_DIVIDE: L = (1 - 1/(1 << CEFO_BITS_DIVIDE))
// 1 L = 0.5
// 2 L = 0.75
// 3 L = 0.875
// 4 L = 0.9375
// 5 L = 0.96875
// 6 L = 0.984375
// 7 L = 0.9921875
// 8 L = 0.99609375
// 9 L = 0.998046875
// 10 L = 0.9990234375
// 11 L = 0.99951171875
// 12 L = 0.999755859375
// decay time constant for values of L, for E = 10^-3 = 60dB of attenuation
//
// Neff = Ln E / Ln L = -6.9077552 / Ln L
//
// 1 L = 0.5 Neff = 10 samples
// 2 L = 0.75 Neff = 24
// 3 L = 0.875 Neff = 51
// 4 L = 0.9375 Neff = 107
// 5 L = 0.96875 Neff = 217
// 6 L = 0.984375 Neff = 438
// 7 L = 0.9921875 Neff = 880
// 8 L = 0.99609375 Neff = 1764
// 9 L = 0.998046875 Neff = 3533
// 10 L = 0.9990234375 Neff = 7070
// 11 L = 0.99951171875 Neff = 14143
// 12 L = 0.999755859375 Neff = 28290
#define CEFO_BITS 11 // 14143 samples in gate window (3hz)
inline int EFO_GetNext( efo_t *pefo, int x ) { int r; int xa = abs(x); int xdif;
// get envelope:
// Cn = L * Cn-1 + ( 1 - L ) * |x|
// which simplifies to:
// Cn = |x| + (Cn-1 - |x|) * L
// for 0 < L < 1
// increasing L increases time to rise or fall to a new input level
// so: increasing CEFO_BITS_DIVIDE increases rise/fall time
// where: L = (1 - 1/(1 << CEFO_BITS))
// xdif = Cn-1 - |x|
// so: xdif * L = xdif - xdif / (1 << CEFO_BITS) = ((xdif << CEFO_BITS) - xdif ) >> CEFO_BITS
xdif = pefo->xout - xa;
pefo->xout = xa + (((xdif << CEFO_BITS) - xdif) >> CEFO_BITS);
if ( pefo->bgate ) { // gate
bool bgateon_prev = pefo->bgateon;
// gate hysteresis
if (bgateon_prev) // gate was on - it's off only if amp drops below thresh_off
pefo->bgateon = ( pefo->xout >= pefo->thresh_off ); else // gate was off - it's on only if amp > thresh
pefo->bgateon = ( pefo->xout >= pefo->thresh ); if ( pefo->bgateon ) { // gate is on
if ( bgateon_prev && RMP_HitEnd( &pefo->rmp_attack )) return x; // gate is fully on
if ( !bgateon_prev ) { // gate just turned on, start ramp attack
// start attack from previous decay ramp if active
r = RMP_HitEnd( &pefo->rmp_decay ) ? 0 : RMP_GetNext( &pefo->rmp_decay ); RMP_SetEnd( &pefo->rmp_decay);
// DevMsg ("GATE ON \n");
RMP_Init( &pefo->rmp_attack, pefo->attack_time, r, PMAX, false); return (x * r) >> PBITS; } if ( !RMP_HitEnd( &pefo->rmp_attack ) ) { r = RMP_GetNext( &pefo->rmp_attack );
// gate is on and ramping up
return (x * r) >> PBITS; } } else { // gate is fully off
if ( !bgateon_prev && RMP_HitEnd( &pefo->rmp_decay)) return 0;
if ( bgateon_prev ) { // gate just turned off, start ramp decay
// start decay from previous attack ramp if active
r = RMP_HitEnd( &pefo->rmp_attack ) ? PMAX : RMP_GetNext( &pefo->rmp_attack ); RMP_SetEnd( &pefo->rmp_attack);
RMP_Init( &pefo->rmp_decay, pefo->decay_time, r, 0, false); // DevMsg ("GATE OFF \n");
// if exponential set, gate has exponential ramp down. Otherwise linear ramp down.
if ( pefo->bexp ) return ( (((x * r) >> PBITS) * r ) >> PBITS ); else return (x * r) >> PBITS;
} else if ( !RMP_HitEnd( &pefo->rmp_decay ) ) { // gate is off and ramping down
r = RMP_GetNext( &pefo->rmp_decay );
// if exponential set, gate has exponential ramp down. Otherwise linear ramp down.
if ( pefo->bexp ) return ( (((x * r) >> PBITS) * r ) >> PBITS ); else return (x * r) >> PBITS; } }
return x; }
return pefo->xout; }
// batch version for performance
inline void EFO_GetNextN( efo_t *pefo, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = EFO_GetNext( pefo, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = EFO_GetNext( pefo, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = EFO_GetNext( pefo, pb->left ); pb++; } return; } } // parameter order
typedef enum { efo_ithreshold, efo_iattack, efo_idecay, efo_iexp,
efo_cparam // # of params
} efo_e;
// parameter ranges
prm_rng_t efo_rng[] = {
{efo_cparam, 0, 0}, // first entry is # of parameters
{efo_ithreshold, -140.0, 0.0}, // gate threshold in db. if 0.0 then no gate.
{efo_iattack, 0.0, 20000.0}, // attack time in milliseconds
{efo_idecay, 0.0, 20000.0}, // envelope decay time in milliseconds
{efo_iexp, 0.0, 1.0}, // if 1, use exponential decay ramp (for more realistic reverb tail)
};
efo_t * EFO_Params ( prc_t *pprc ) { efo_t *penv;
float threshold = Gain_To_Amplitude( dB_To_Gain(pprc->prm[efo_ithreshold]) ); float attack = pprc->prm[efo_iattack]/1000.0; float decay = pprc->prm[efo_idecay]/1000.0; float fexp = pprc->prm[efo_iexp]; bool bexp;
// check for no gate
if ( pprc->prm[efo_ithreshold] == 0.0 ) threshold = 0.0;
bexp = fexp > 0.0 ? 1 : 0;
penv = EFO_Alloc ( threshold, attack, decay, bexp ); return penv; }
inline void * EFO_VParams ( void *p ) { PRC_CheckParams( (prc_t *)p, efo_rng ); return (void *) EFO_Params ((prc_t *)p); }
inline void EFO_Mod ( void *p, float v ) { return; }
///////////////////////////////////////////
// Chorus - lfo modulated delay
///////////////////////////////////////////
#define CCRSS 64 // max number chorus' active
struct crs_t { bool fused;
mdy_t *pmdy; // modulatable delay
lfo_t *plfo; // modulating lfo
int lfoprev; // previous modulator value from lfo
};
crs_t crss[CCRSS];
void CRS_Init( crs_t *pcrs ) { if (pcrs) Q_memset( pcrs, 0, sizeof (crs_t) ); }; void CRS_Free( crs_t *pcrs ) { if (pcrs) { MDY_Free ( pcrs->pmdy ); LFO_Free ( pcrs->plfo ); Q_memset( pcrs, 0, sizeof (crs_t) ); } }
void CRS_InitAll() { for (int i = 0; i < CCRSS; i++) CRS_Init( &crss[i] ); } void CRS_FreeAll() { for (int i = 0; i < CCRSS; i++) CRS_Free( &crss[i] ); }
// fstep is base pitch shift, ie: floating point step value, where 1.0 = +1 octave, 0.5 = -1 octave
// lfotype is LFO_SIN, LFO_RND, LFO_TRI etc (LFO_RND for chorus, LFO_SIN for flange)
// fHz is modulation frequency in Hz
// depth is modulation depth, 0-1.0
// mix is mix of chorus and clean signal
#define CRS_DELAYMAX 100 // max milliseconds of sweepable delay
#define CRS_RAMPTIME 5 // milliseconds to ramp between new delay values
crs_t * CRS_Alloc( int lfotype, float fHz, float fdepth, float mix ) { int i; crs_t *pcrs; dly_t *pdly; mdy_t *pmdy; lfo_t *plfo; float ramptime; int D;
// find free chorus slot
for ( i = 0; i < CCRSS; i++ ) { if ( !crss[i].fused ) break; }
if ( i == CCRSS ) { DevMsg ("DSP: Warning, failed to allocate chorus.\n" ); return NULL; }
pcrs = &crss[i];
CRS_Init ( pcrs );
D = fdepth * MSEC_TO_SAMPS(CRS_DELAYMAX); // sweep from 0 - n milliseconds
ramptime = (float) CRS_RAMPTIME / 1000.0; // # milliseconds to ramp between new values
pdly = DLY_Alloc ( D, 0, 1, DLY_LINEAR );
pmdy = MDY_Alloc ( pdly, ramptime, 0.0, 0.0, mix );
plfo = LFO_Alloc ( lfotype, fHz, false, 1.0 ); if ( !plfo || !pmdy ) { LFO_Free ( plfo ); MDY_Free ( pmdy ); DevMsg ("DSP: Warning, failed to allocate lfo or mdy for chorus.\n" ); return NULL; }
pcrs->pmdy = pmdy; pcrs->plfo = plfo; pcrs->fused = true;
return pcrs; }
// return next chorused sample (modulated delay) mixed with input sample
inline int CRS_GetNext( crs_t *pcrs, int x ) { int l; int y; // get current mod delay value
y = MDY_GetNext ( pcrs->pmdy, x );
// get next lfo value for modulation
// note: lfo must return 0 as first value
l = LFO_GetNext ( pcrs->plfo, x );
// if modulator has changed, change mdy
if ( l != pcrs->lfoprev ) { // calculate new tap starts at D)
int D = pcrs->pmdy->pdly->D0; int tap;
// lfo should always output values 0 <= l <= LFOMAX
if (l < 0) l = 0;
tap = D - ((l * D) >> LFOBITS);
MDY_ChangeVal ( pcrs->pmdy, tap );
pcrs->lfoprev = l; }
return y; }
// batch version for performance
inline void CRS_GetNextN( crs_t *pcrs, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = CRS_GetNext( pcrs, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = CRS_GetNext( pcrs, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = CRS_GetNext( pcrs, pb->left ); pb++; } return; } }
// parameter order
typedef enum {
crs_ilfotype, crs_irate, crs_idepth, crs_imix, crs_cparam
} crs_e;
// parameter ranges
prm_rng_t crs_rng[] = {
{crs_cparam, 0, 0}, // first entry is # of parameters
{crs_ilfotype, 0, LFO_MAX}, // lfotype is LFO_SIN, LFO_RND, LFO_TRI etc (LFO_RND for chorus, LFO_SIN for flange)
{crs_irate, 0.0, 1000.0}, // rate is modulation frequency in Hz
{crs_idepth, 0.0, 1.0}, // depth is modulation depth, 0-1.0
{crs_imix, 0.0, 1.0}, // mix is mix of chorus and clean signal
};
// uses pitch, lfowav, rate, depth
crs_t * CRS_Params ( prc_t *pprc ) { crs_t *pcrs; pcrs = CRS_Alloc ( pprc->prm[crs_ilfotype], pprc->prm[crs_irate], pprc->prm[crs_idepth], pprc->prm[crs_imix] );
return pcrs; }
inline void * CRS_VParams ( void *p ) { PRC_CheckParams ( (prc_t *)p, crs_rng ); return (void *) CRS_Params ((prc_t *)p); }
inline void CRS_Mod ( void *p, float v ) { return; }
////////////////////////////////////////////////////
// amplifier - modulatable gain, distortion
////////////////////////////////////////////////////
#define CAMPS 64 // max number amps active
#define AMPSLEW 10 // milliseconds of slew time between gain changes
struct amp_t { bool fused;
int gain; // amplification 0-6.0 * PMAX
int gain_max; // original gain setting
int distmix; // 0-1.0 mix of distortion with clean * PMAX
int vfeed; // 0-1.0 feedback with distortion * PMAX
int vthresh; // amplitude of clipping threshold 0..32768
bool fchanging; // true if modulating to new amp value
float ramptime; // ramp 'glide' time - time in seconds to change between values
int mtime; // time in samples between amp changes. 0 implies no self-modulating
int mtimecur; // current time in samples until next amp change
int depth; // modulate amp from A to A - (A*depth) depth 0-1.0
bool brand; // if true, use random modulation otherwise alternate btwn max/min
rmp_t rmp_interp; // interpolation ramp 0...PMAX
};
amp_t amps[CAMPS];
void AMP_Init( amp_t *pamp ) { if (pamp) Q_memset( pamp, 0, sizeof (amp_t) ); }; void AMP_Free( amp_t *pamp ) { if (pamp) { Q_memset( pamp, 0, sizeof (amp_t) ); } }
void AMP_InitAll() { for (int i = 0; i < CAMPS; i++) AMP_Init( &s[i] ); } void AMP_FreeAll() { for (int i = 0; i < CAMPS; i++) AMP_Free( &s[i] ); }
amp_t * AMP_Alloc( float gain, float vthresh, float distmix, float vfeed, float ramptime, float modtime, float depth, bool brand ) { int i; amp_t *pamp;
// find free amp slot
for ( i = 0; i < CAMPS; i++ ) { if ( !amps[i].fused ) break; }
if ( i == CAMPS ) { DevMsg ("DSP: Warning, failed to allocate amp.\n" ); return NULL; }
pamp = &s[i];
AMP_Init ( pamp );
pamp->fused = true;
pamp->gain = gain * PMAX; pamp->gain_max = gain * PMAX; pamp->distmix = distmix * PMAX; pamp->vfeed = vfeed * PMAX; pamp->vthresh = vthresh * 32767.0;
// modrate, 0.01, 200.0}, // frequency at which amplitude values change to new random value. 0 is no self-modulation
// moddepth, 0.0, 1.0}, // how much amplitude changes (decreases) from current value (0-1.0)
// modglide, 0.01, 100.0}, // glide time between mapcur and ampnew in milliseconds
pamp->ramptime = ramptime; pamp->mtime = SEC_TO_SAMPS(modtime); pamp->mtimecur = pamp->mtime; pamp->depth = depth * PMAX; pamp->brand = brand;
return pamp; }
// return next amplified sample
inline int AMP_GetNext( amp_t *pamp, int x ) { int y = x; int d;
// if distortion is on, add distortion, feedback
if ( pamp->vthresh < PMAX && pamp->distmix ) { int vthresh = pamp->vthresh;
/* if ( pamp->vfeed > 0.0 )
{ // UNDONE: feedback
} */ // clip distort
d = ( y > vthresh ? vthresh : ( y < -vthresh ? -vthresh : y));
// mix distorted with clean (1.0 = full distortion)
if ( pamp->distmix < PMAX ) y = y + (((d - y) * pamp->distmix ) >> PBITS); else y = d; }
// get output for current gain value
int xout = (y * pamp->gain) >> PBITS; if ( !pamp->fchanging && !pamp->mtime ) { // if not modulating and not self modulating, return right away
return xout; }
if (pamp->fchanging) { // modulating...
// get next gain value
pamp->gain = RMP_GetNext( &pamp->rmp_interp ); // 0...next gain
if ( RMP_HitEnd( &pamp->rmp_interp ) ) { // done.
pamp->fchanging = false; } } // if self-modulating and timer has expired, get next change
if ( pamp->mtime && !pamp->mtimecur-- ) { pamp->mtimecur = pamp->mtime;
int gain_new; int G1; int G2 = pamp->gain_max;
// modulate between 0 and 100% of gain_max
G1 = pamp->gain_max - ((pamp->gain_max * pamp->depth) >> PBITS);
if (pamp->brand) { gain_new = RandomInt( min(G1,G2), max(G1,G2) ); } else { // alternate between min & max
gain_new = (pamp->gain == G1 ? G2 : G1); }
// set up modulation to new value
pamp->fchanging = true;
// init gain ramp - always hit target
RMP_Init ( &pamp->rmp_interp, pamp->ramptime, pamp->gain, gain_new, false ); }
return xout;
}
// batch version for performance
inline void AMP_GetNextN( amp_t *pamp, portable_samplepair_t *pbuffer, int SampleCount, int op ) { int count = SampleCount; portable_samplepair_t *pb = pbuffer; switch (op) { default: case OP_LEFT: while (count--) { pb->left = AMP_GetNext( pamp, pb->left ); pb++; } return; case OP_RIGHT: while (count--) { pb->right = AMP_GetNext( pamp, pb->right ); pb++; } return; case OP_LEFT_DUPLICATE: while (count--) { pb->left = pb->right = AMP_GetNext( pamp, pb->left ); pb++; } return; } }
inline void AMP_Mod( amp_t *pamp, float v ) { }
// parameter order
typedef enum {
amp_gain, amp_vthresh, amp_distmix, amp_vfeed, amp_imodrate, amp_imoddepth, amp_imodglide, amp_irand, amp_cparam
} amp_e;
// parameter ranges
prm_rng_t amp_rng[] = {
{amp_cparam, 0, 0}, // first entry is # of parameters
{amp_gain, 0.0, 1000.0}, // amplification
{amp_vthresh, 0.0, 1.0}, // threshold for distortion (1.0 = no distortion)
{amp_distmix, 0.0, 1.0}, // mix of clean and distortion (1.0 = full distortion, 0.0 = full clean)
{amp_vfeed, 0.0, 1.0}, // distortion feedback
{amp_imodrate, 0.0, 200.0}, // frequency at which amplitude values change to new random value. 0 is no self-modulation
{amp_imoddepth, 0.0, 1.0}, // how much amplitude changes (decreases) from current value (0-1.0)
{amp_imodglide, 0.01, 100.0}, // glide time between mapcur and ampnew in milliseconds
{amp_irand, 0.0, 1.0}, // if 1, use random modulation otherwise alternate from max-min-max
};
amp_t * AMP_Params ( prc_t *pprc ) { amp_t *pamp;
float ramptime = 0.0; float modtime = 0.0; float depth = 0.0; float rand = pprc->prm[amp_irand]; bool brand;
if (pprc->prm[amp_imodrate] > 0.0) { ramptime = pprc->prm[amp_imodglide] / 1000.0; // get ramp time in seconds
modtime = 1.0 / max((double)pprc->prm[amp_imodrate], 0.01); // time between modulations in seconds
depth = pprc->prm[amp_imoddepth]; // depth of modulations 0-1.0
}
brand = rand > 0.0 ? 1 : 0;
pamp = AMP_Alloc ( pprc->prm[amp_gain], pprc->prm[amp_vthresh], pprc->prm[amp_distmix], pprc->prm[amp_vfeed], ramptime, modtime, depth, brand );
return pamp; }
inline void * AMP_VParams ( void *p ) { PRC_CheckParams ( (prc_t *)p, amp_rng ); return (void *) AMP_Params ((prc_t *)p); }
/////////////////
// NULL processor
/////////////////
struct nul_t { int type; };
nul_t nuls[] = {{0}};
void NULL_Init ( nul_t *pnul ) { } void NULL_InitAll( ) { } void NULL_Free ( nul_t *pnul ) { } void NULL_FreeAll ( ) { } nul_t *NULL_Alloc ( ) { return &nuls[0]; }
inline int NULL_GetNext ( void *p, int x) { return x; }
inline void NULL_GetNextN( nul_t *pnul, portable_samplepair_t *pbuffer, int SampleCount, int op ) { return; }
inline void NULL_Mod ( void *p, float v ) { return; }
inline void * NULL_VParams ( void *p ) { return (void *) (&nuls[0]); }
//////////////////////////
// DSP processors presets - see dsp_presets.txt
//////////////////////////
// init array of processors - first store pfnParam, pfnGetNext and pfnFree functions for type,
// then call the pfnParam function to initialize each processor
// prcs - an array of prc structures, all with initialized params
// count - number of elements in the array
// returns false if failed to init one or more processors
bool PRC_InitAll( prc_t *prcs, int count ) { int i; prc_Param_t pfnParam; // allocation function - takes ptr to prc, returns ptr to specialized data struct for proc type
prc_GetNext_t pfnGetNext; // get next function
prc_GetNextN_t pfnGetNextN; // get next function, batch version
prc_Free_t pfnFree; prc_Mod_t pfnMod;
bool fok = true;;
if ( count == 0 ) count = 1;
// set up pointers to XXX_Free, XXX_GetNext and XXX_Params functions
for (i = 0; i < count; i++) { switch (prcs[i].type) { default: case PRC_NULL: pfnFree = (prc_Free_t)NULL_Free; pfnGetNext = (prc_GetNext_t)NULL_GetNext; pfnGetNextN = (prc_GetNextN_t)NULL_GetNextN; pfnParam = NULL_VParams; pfnMod = (prc_Mod_t)NULL_Mod; break; case PRC_DLY: pfnFree = (prc_Free_t)DLY_Free; pfnGetNext = (prc_GetNext_t)DLY_GetNext; pfnGetNextN = (prc_GetNextN_t)DLY_GetNextN; pfnParam = DLY_VParams; pfnMod = (prc_Mod_t)DLY_Mod; break; case PRC_RVA: pfnFree = (prc_Free_t)RVA_Free; pfnGetNext = (prc_GetNext_t)RVA_GetNext; pfnGetNextN = (prc_GetNextN_t)RVA_GetNextN; pfnParam = RVA_VParams; pfnMod = (prc_Mod_t)RVA_Mod; break; case PRC_FLT: pfnFree = (prc_Free_t)FLT_Free; pfnGetNext = (prc_GetNext_t)FLT_GetNext; pfnGetNextN = (prc_GetNextN_t)FLT_GetNextN; pfnParam = FLT_VParams; pfnMod = (prc_Mod_t)FLT_Mod; break; case PRC_CRS: pfnFree = (prc_Free_t)CRS_Free; pfnGetNext = (prc_GetNext_t)CRS_GetNext; pfnGetNextN = (prc_GetNextN_t)CRS_GetNextN; pfnParam = CRS_VParams; pfnMod = (prc_Mod_t)CRS_Mod; break; case PRC_PTC: pfnFree = (prc_Free_t)PTC_Free; pfnGetNext = (prc_GetNext_t)PTC_GetNext; pfnGetNextN = (prc_GetNextN_t)PTC_GetNextN; pfnParam = PTC_VParams; pfnMod = (prc_Mod_t)PTC_Mod; break; case PRC_ENV: pfnFree = (prc_Free_t)ENV_Free; pfnGetNext = (prc_GetNext_t)ENV_GetNext; pfnGetNextN = (prc_GetNextN_t)ENV_GetNextN; pfnParam = ENV_VParams; pfnMod = (prc_Mod_t)ENV_Mod; break; case PRC_LFO: pfnFree = (prc_Free_t)LFO_Free; pfnGetNext = (prc_GetNext_t)LFO_GetNext; pfnGetNextN = (prc_GetNextN_t)LFO_GetNextN; pfnParam = LFO_VParams; pfnMod = (prc_Mod_t)LFO_Mod; break; case PRC_EFO: pfnFree = (prc_Free_t)EFO_Free; pfnGetNext = (prc_GetNext_t)EFO_GetNext; pfnGetNextN = (prc_GetNextN_t)EFO_GetNextN; pfnParam = EFO_VParams; pfnMod = (prc_Mod_t)EFO_Mod; break; case PRC_MDY: pfnFree = (prc_Free_t)MDY_Free; pfnGetNext = (prc_GetNext_t)MDY_GetNext; pfnGetNextN = (prc_GetNextN_t)MDY_GetNextN; pfnParam = MDY_VParams; pfnMod = (prc_Mod_t)MDY_Mod; break; case PRC_DFR: pfnFree = (prc_Free_t)DFR_Free; pfnGetNext = (prc_GetNext_t)DFR_GetNext; pfnGetNextN = (prc_GetNextN_t)DFR_GetNextN; pfnParam = DFR_VParams; pfnMod = (prc_Mod_t)DFR_Mod; break; case PRC_AMP: pfnFree = (prc_Free_t)AMP_Free; pfnGetNext = (prc_GetNext_t)AMP_GetNext; pfnGetNextN = (prc_GetNextN_t)AMP_GetNextN; pfnParam = AMP_VParams; pfnMod = (prc_Mod_t)AMP_Mod; break; }
// set up function pointers
prcs[i].pfnParam = pfnParam; prcs[i].pfnGetNext = pfnGetNext; prcs[i].pfnGetNextN = pfnGetNextN; prcs[i].pfnFree = pfnFree; prcs[i].pfnMod = pfnMod;
// call param function, store pdata for the processor type
prcs[i].pdata = pfnParam ( (void *) (&prcs[i]) );
if ( !prcs[i].pdata ) fok = false; }
return fok; }
// free individual processor's data
void PRC_Free ( prc_t *pprc ) { if ( pprc->pfnFree && pprc->pdata ) pprc->pfnFree ( pprc->pdata ); }
// free all processors for supplied array
// prcs - array of processors
// count - elements in array
void PRC_FreeAll ( prc_t *prcs, int count ) { for (int i = 0; i < count; i++) PRC_Free( &prcs[i] ); }
// get next value for processor - (usually called directly by PSET_GetNext)
inline int PRC_GetNext ( prc_t *pprc, int x ) { return pprc->pfnGetNext ( pprc->pdata, x ); }
// automatic parameter range limiting
// force parameters between specified min/max in param_rng
void PRC_CheckParams ( prc_t *pprc, prm_rng_t *prng ) { // first entry in param_rng is # of parameters
int cprm = prng[0].iprm;
for (int i = 0; i < cprm; i++) { // if parameter is 0.0, always allow it (this is 'off' for most params)
if ( pprc->prm[i] != 0.0 && (pprc->prm[i] > prng[i+1].hi || pprc->prm[i] < prng[i+1].lo) ) { DevMsg ("DSP: Warning, clamping out of range parameter.\n" ); pprc->prm[i] = clamp (pprc->prm[i], prng[i+1].lo, prng[i+1].hi); } } }
// DSP presets
// A dsp preset comprises one or more dsp processors in linear, parallel or feedback configuration
// preset configurations
//
#define PSET_SIMPLE 0
// x(n)--->P(0)--->y(n)
#define PSET_LINEAR 1
// x(n)--->P(0)-->P(1)-->...P(m)--->y(n)
#define PSET_PARALLEL2 5
// x(n)--->P(0)-->(+)-->y(n)
// ^
// |
// x(n)--->P(1)-----
#define PSET_PARALLEL4 6
// x(n)--->P(0)-->P(1)-->(+)-->y(n)
// ^
// |
// x(n)--->P(2)-->P(3)-----
#define PSET_PARALLEL5 7
// x(n)--->P(0)-->P(1)-->(+)-->P(4)-->y(n)
// ^
// |
// x(n)--->P(2)-->P(3)-----
#define PSET_FEEDBACK 8
// x(n)-P(0)--(+)-->P(1)-->P(2)---->y(n)
// ^ |
// | v
// -----P(4)<--P(3)--
#define PSET_FEEDBACK3 9
// x(n)---(+)-->P(0)--------->y(n)
// ^ |
// | v
// -----P(2)<--P(1)--
#define PSET_FEEDBACK4 10
// x(n)---(+)-->P(0)-------->P(3)--->y(n)
// ^ |
// | v
// ---P(2)<--P(1)--
#define PSET_MOD 11
//
// x(n)------>P(1)--P(2)--P(3)--->y(n)
// ^
// x(n)------>P(0)....:
#define PSET_MOD2 12
//
// x(n)-------P(1)-->y(n)
// ^
// x(n)-->P(0)..:
#define PSET_MOD3 13
//
// x(n)-------P(1)-->P(2)-->y(n)
// ^
// x(n)-->P(0)..:
#define CPSETS 64 // max number of presets simultaneously active
#define CPSET_PRCS 5 // max # of processors per dsp preset
#define CPSET_STATES (CPSET_PRCS+3) // # of internal states
// NOTE: do not reorder members of pset_t - g_psettemplates relies on it!!!
struct pset_t { int type; // preset configuration type
int cprcs; // number of processors for this preset
prc_t prcs[CPSET_PRCS]; // processor preset data
float mix_min; // min dsp mix at close range
float mix_max; // max dsp mix at long range
float db_min; // if sndlvl of a new sound is < db_min, reduce mix_min/max by db_mixdrop
float db_mixdrop; // reduce mix_min/max by n% if sndlvl of new sound less than db_min
float duration; // if > 0, duration of preset in seconds (duration 0 = infinite)
float fade; // fade out time, exponential fade
int csamp_duration; // duration counter # samples
int w[CPSET_STATES]; // internal states
int fused; };
pset_t psets[CPSETS];
pset_t *g_psettemplates = NULL; int g_cpsettemplates = 0;
// returns true if preset will expire after duration
bool PSET_IsOneShot( pset_t *ppset ) { return ppset->duration > 0.0; }
// return true if preset is no longer active - duration has expired
bool PSET_HasExpired( pset_t *ppset ) { if (!PSET_IsOneShot( ppset )) return false;
return ppset->csamp_duration <= 0; }
// if preset is oneshot, update duration counter by SampleCount samples
void PSET_UpdateDuration( pset_t *ppset, int SampleCount ) { if ( PSET_IsOneShot( ppset ) ) { // if oneshot preset and not expired, decrement sample count
if (ppset->csamp_duration > 0) ppset->csamp_duration -= SampleCount; } }
// A dsp processor (prc) performs a single-sample function, such as pitch shift, delay, reverb, filter
// init a preset - just clear state array
void PSET_Init( pset_t *ppset ) { // clear state array
if (ppset) Q_memset( ppset->w, 0, sizeof (int) * (CPSET_STATES) ); }
// clear runtime slots
void PSET_InitAll( void ) { for (int i = 0; i < CPSETS; i++) Q_memset( &psets[i], 0, sizeof(pset_t)); }
// free the preset - free all processors
void PSET_Free( pset_t *ppset ) { if (ppset) { // free processors
PRC_FreeAll ( ppset->prcs, ppset->cprcs );
// clear
Q_memset( ppset, 0, sizeof (pset_t)); } }
void PSET_FreeAll() { for (int i = 0; i < CPSETS; i++) PSET_Free( &psets[i] ); };
// return preset struct, given index into preset template array
// NOTE: should not ever be more than 2 or 3 of these active simultaneously
pset_t * PSET_Alloc ( int ipsettemplate ) { pset_t *ppset; bool fok;
// don't excede array bounds
if ( ipsettemplate >= g_cpsettemplates) ipsettemplate = 0;
// find free slot
int i = 0; for (i = 0; i < CPSETS; i++) { if ( !psets[i].fused ) break; }
if ( i == CPSETS ) return NULL;
if (das_debug.GetInt()) { int nSlots = 0; for ( int j = 0; j < CPSETS; j++) { if ( psets[j].fused ) nSlots++; } DevMsg("total preset slots used: %d \n", nSlots ); }
ppset = &psets[i]; // clear preset
Q_memset(ppset, 0, sizeof(pset_t));
// copy template into preset
*ppset = g_psettemplates[ipsettemplate];
ppset->fused = true;
// clear state array
PSET_Init ( ppset ); // init all processors, set up processor function pointers
fok = PRC_InitAll( ppset->prcs, ppset->cprcs );
if ( !fok ) { // failed to init one or more processors
Warning( "Sound DSP: preset failed to init.\n"); PRC_FreeAll ( ppset->prcs, ppset->cprcs ); return NULL; }
// if preset has duration, setup duration sample counter
if ( PSET_IsOneShot( ppset ) ) { ppset->csamp_duration = SEC_TO_SAMPS( ppset->duration ); }
return ppset; }
// batch version of PSET_GetNext for linear array of processors. For performance.
// ppset - preset array
// pbuffer - input sample data
// SampleCount - size of input buffer
// OP: OP_LEFT - process left channel in place
// OP_RIGHT - process right channel in place
// OP_LEFT_DUPLICATe - process left channel, duplicate into right
inline void PSET_GetNextN( pset_t *ppset, portable_samplepair_t *pbuffer, int SampleCount, int op ) { portable_samplepair_t *pbf = pbuffer; prc_t *pprc; int count = ppset->cprcs;
switch ( ppset->type ) { default: case PSET_SIMPLE: { // x(n)--->P(0)--->y(n)
ppset->prcs[0].pfnGetNextN (ppset->prcs[0].pdata, pbf, SampleCount, op); return; } case PSET_LINEAR: {
// w0 w1 w2
// x(n)--->P(0)-->P(1)-->...P(count-1)--->y(n)
// w0 w1 w2 w3 w4 w5
// x(n)--->P(0)-->P(1)-->P(2)-->P(3)-->P(4)-->y(n)
// call batch processors in sequence - no internal state for batch processing
// point to first processor
pprc = &ppset->prcs[0];
for (int i = 0; i < count; i++) { pprc->pfnGetNextN (pprc->pdata, pbf, SampleCount, op); pprc++; }
return; } } }
// Get next sample from this preset. called once for every sample in buffer
// ppset is pointer to preset
// x is input sample
inline int PSET_GetNext ( pset_t *ppset, int x ) {
// pset_simple and pset_linear have no internal state:
// this is REQUIRED for all presets that have a batch getnextN equivalent!
if ( ppset->type == PSET_SIMPLE ) { // x(n)--->P(0)--->y(n)
return ppset->prcs[0].pfnGetNext (ppset->prcs[0].pdata, x); } prc_t *pprc; int count = ppset->cprcs; if ( ppset->type == PSET_LINEAR ) { int y = x;
// w0 w1 w2
// x(n)--->P(0)-->P(1)-->...P(count-1)--->y(n)
// w0 w1 w2 w3 w4 w5
// x(n)--->P(0)-->P(1)-->P(2)-->P(3)-->P(4)-->y(n)
// call processors in reverse order, from count to 1
//for (int i = count; i > 0; i--, pprc--)
// w[i] = pprc->pfnGetNext (pprc->pdata, w[i-1]);
// return w[count];
// point to first processor, update sequentially, no state preserved
pprc = &ppset->prcs[0];
switch (count) { default: case 5: y = pprc->pfnGetNext (pprc->pdata, y); pprc++; case 4: y = pprc->pfnGetNext (pprc->pdata, y); pprc++; case 3: y = pprc->pfnGetNext (pprc->pdata, y); pprc++; case 2: y = pprc->pfnGetNext (pprc->pdata, y); pprc++; case 1: case 0: y = pprc->pfnGetNext (pprc->pdata, y); }
return y; }
// all other preset types have internal state:
// initialize 0'th element of state array
int *w = ppset->w; w[0] = x;
switch ( ppset->type ) { default: case PSET_PARALLEL2: { // w0 w1 w3
// x(n)--->P(0)-->(+)-->y(n)
// ^
// w0 w2 |
// x(n)--->P(1)-----
pprc = &ppset->prcs[0];
w[3] = w[1] + w[2];
w[1] = pprc->pfnGetNext( pprc->pdata, w[0] ); pprc++; w[2] = pprc->pfnGetNext( pprc->pdata, w[0] );
return w[3]; }
case PSET_PARALLEL4: { // w0 w1 w2 w5
// x(n)--->P(0)-->P(1)-->(+)-->y(n)
// ^
// w0 w3 w4 |
// x(n)--->P(2)-->P(3)-----
pprc = &ppset->prcs[0];
w[5] = w[2] + w[4];
w[2] = pprc[1].pfnGetNext( pprc[1].pdata, w[1] ); w[4] = pprc[3].pfnGetNext( pprc[3].pdata, w[3] );
w[1] = pprc[0].pfnGetNext( pprc[0].pdata, w[0] ); w[3] = pprc[2].pfnGetNext( pprc[2].pdata, w[0] );
return w[5]; }
case PSET_PARALLEL5: { // w0 w1 w2 w5 w6
// x(n)--->P(0)-->P(1)-->(+)--P(4)-->y(n)
// ^
// w0 w3 w4 |
// x(n)--->P(2)-->P(3)-----
pprc = &ppset->prcs[0];
w[5] = w[2] + w[4];
w[2] = pprc[1].pfnGetNext( pprc[1].pdata, w[1] ); w[4] = pprc[3].pfnGetNext( pprc[3].pdata, w[3] );
w[1] = pprc[0].pfnGetNext( pprc[0].pdata, w[0] ); w[3] = pprc[2].pfnGetNext( pprc[2].pdata, w[0] );
return pprc[4].pfnGetNext( pprc[4].pdata, w[5] ); }
case PSET_FEEDBACK: { // w0 w1 w2 w3 w4 w7
// x(n)-P(0)--(+)-->P(1)-->P(2)-->---->y(n)
// ^ |
// | w6 w5 v
// -----P(4)<--P(3)--
pprc = &ppset->prcs[0]; // start with adders
w[2] = w[1] + w[6];
// evaluate in reverse order
w[6] = pprc[4].pfnGetNext( pprc[4].pdata, w[5] ); w[5] = pprc[3].pfnGetNext( pprc[3].pdata, w[4] );
w[4] = pprc[2].pfnGetNext( pprc[2].pdata, w[3] ); w[3] = pprc[1].pfnGetNext( pprc[1].pdata, w[2] ); w[1] = pprc[0].pfnGetNext( pprc[0].pdata, w[0] );
return w[4]; } case PSET_FEEDBACK3: { // w0 w1 w2
// x(n)---(+)-->P(0)--------->y(n)
// ^ |
// | w4 w3 v
// -----P(2)<--P(1)--
pprc = &ppset->prcs[0]; // start with adders
w[1] = w[0] + w[4];
// evaluate in reverse order
w[4] = pprc[2].pfnGetNext( pprc[2].pdata, w[3] ); w[3] = pprc[1].pfnGetNext( pprc[1].pdata, w[2] ); w[2] = pprc[0].pfnGetNext( pprc[0].pdata, w[1] ); return w[2]; } case PSET_FEEDBACK4: { // w0 w1 w2 w5
// x(n)---(+)-->P(0)-------->P(3)--->y(n)
// ^ |
// | w4 w3 v
// ---P(2)<--P(1)--
pprc = &ppset->prcs[0]; // start with adders
w[1] = w[0] + w[4];
// evaluate in reverse order
w[5] = pprc[3].pfnGetNext( pprc[3].pdata, w[2] ); w[4] = pprc[2].pfnGetNext( pprc[2].pdata, w[3] ); w[3] = pprc[1].pfnGetNext( pprc[1].pdata, w[2] ); w[2] = pprc[0].pfnGetNext( pprc[0].pdata, w[1] ); return w[2]; } case PSET_MOD: { // w0 w1 w3 w4
// x(n)------>P(1)--P(2)--P(3)--->y(n)
// w0 w2 ^
// x(n)------>P(0)....:
pprc = &ppset->prcs[0];
w[4] = pprc[3].pfnGetNext( pprc[3].pdata, w[3] );
w[3] = pprc[2].pfnGetNext( pprc[2].pdata, w[1] );
// modulate processor 2
pprc[2].pfnMod( pprc[2].pdata, ((float)w[2] / (float)PMAX));
// get modulator output
w[2] = pprc[0].pfnGetNext( pprc[0].pdata, w[0] );
w[1] = pprc[1].pfnGetNext( pprc[1].pdata, w[0] );
return w[4]; } case PSET_MOD2: { // w0 w2
// x(n)---------P(1)-->y(n)
// w0 w1 ^
// x(n)-->P(0)....:
pprc = &ppset->prcs[0];
// modulate processor 1
pprc[1].pfnMod( pprc[1].pdata, ((float)w[1] / (float)PMAX));
// get modulator output
w[1] = pprc[0].pfnGetNext( pprc[0].pdata, w[0] );
w[2] = pprc[1].pfnGetNext( pprc[1].pdata, w[0] );
return w[2];
} case PSET_MOD3: { // w0 w2 w3
// x(n)----------P(1)-->P(2)-->y(n)
// w0 w1 ^
// x(n)-->P(0).....:
pprc = &ppset->prcs[0];
w[3] = pprc[2].pfnGetNext( pprc[2].pdata, w[2] );
// modulate processor 1
pprc[1].pfnMod( pprc[1].pdata, ((float)w[1] / (float)PMAX));
// get modulator output
w[1] = pprc[0].pfnGetNext( pprc[0].pdata, w[0] );
w[2] = pprc[1].pfnGetNext( pprc[1].pdata, w[0] );
return w[2]; } } }
/////////////
// DSP system
/////////////
// Main interface
// Whenever the preset # changes on any of these processors, the old processor is faded out, new is faded in.
// dsp_chan is optionally set when a sound is played - a preset is sent with the start_static/dynamic sound.
//
// sound1---->dsp_chan--> -------------(+)---->dsp_water--->dsp_player--->out
// sound2---->dsp_chan--> | |
// sound3---------------> ----dsp_room---
// | |
// --dsp_indirect-
// dsp_room - set this cvar to a preset # to change the room dsp. room fx are more prevalent farther from player.
// use: when player moves into a new room, all sounds played in room take on its reverberant character
// dsp_water - set this cvar (once) to a preset # for serial underwater sound.
// use: when player goes under water, all sounds pass through this dsp (such as low pass filter)
// dsp_player - set this cvar to a preset # to cause all sounds to run through the effect (serial, in-line).
// use: player is deafened, player fires special weapon, player is hit by special weapon.
// dsp_facingaway- set this cvar to a preset # appropriate for sounds which are played facing away from player (weapon,voice)
//
// dsp_spatial - set by system to create modulated spatial delays for left/right/front/back ears - delay value
// modulates by distance to nearest l/r surface in world
// Dsp presets
ConVar dsp_room ("dsp_room", "0", FCVAR_DEMO ); // room dsp preset - sounds more distant from player (1ch)
ConVar dsp_water ("dsp_water", "14", FCVAR_DEMO ); // "14" underwater dsp preset - sound when underwater (1-2ch)
ConVar dsp_player ("dsp_player", "0", FCVAR_DEMO | FCVAR_SERVER_CAN_EXECUTE ); // dsp on player - sound when player hit by special device (1-2ch)
ConVar dsp_facingaway ("dsp_facingaway", "0", FCVAR_DEMO ); // "30" sounds that face away from player (weapons, voice) (1-4ch)
ConVar dsp_speaker ("dsp_speaker", "50", FCVAR_DEMO ); // "50" small distorted speaker sound (1ch)
ConVar dsp_spatial ("dsp_spatial", "40", FCVAR_DEMO ); // spatial delays for l/r front/rear ears
ConVar dsp_automatic ("dsp_automatic", "0", FCVAR_DEMO ); // automatic room type detection. if non zero, replaces dsp_room
int ipset_room_prev; int ipset_water_prev; int ipset_player_prev; int ipset_facingaway_prev; int ipset_speaker_prev; int ipset_spatial_prev; int ipset_automatic_prev;
// legacy room_type support
ConVar dsp_room_type ( "room_type", "0", FCVAR_DEMO ); int ipset_room_typeprev;
// DSP processors
int idsp_room; int idsp_water; int idsp_player; int idsp_facingaway; int idsp_speaker; int idsp_spatial; int idsp_automatic;
ConVar dsp_off ("dsp_off", "0", FCVAR_CHEAT | FCVAR_ALLOWED_IN_COMPETITIVE ); // set to 1 to disable all dsp processing
ConVar dsp_slow_cpu ("dsp_slow_cpu", "0", FCVAR_ARCHIVE|FCVAR_DEMO ); // set to 1 if cpu bound - ie: does not process dsp_room fx
ConVar snd_profile ("snd_profile", "0", FCVAR_DEMO ); // 1 - profile dsp, 2 - mix, 3 - load sound, 4 - all sound
ConVar dsp_volume ("dsp_volume", "1.0", FCVAR_ARCHIVE|FCVAR_DEMO ); // 0.0 - 2.0; master dsp volume control
ConVar dsp_vol_5ch ("dsp_vol_5ch", "0.5", FCVAR_DEMO ); // 0.0 - 1.0; attenuate master dsp volume for 5ch surround
ConVar dsp_vol_4ch ("dsp_vol_4ch", "0.5", FCVAR_DEMO ); // 0.0 - 1.0; attenuate master dsp volume for 4ch surround
ConVar dsp_vol_2ch ("dsp_vol_2ch", "1.0", FCVAR_DEMO ); // 0.0 - 1.0; attenuate master dsp volume for 2ch surround
ConVar dsp_enhance_stereo("dsp_enhance_stereo", "0", FCVAR_ARCHIVE ); // 1) use dsp_spatial delays on all reverb channels
// DSP preset executor
#define CDSPS 32 // max number dsp executors active
#define DSPCHANMAX 5 // max number of channels dsp can process (allocs a separte processor for each chan)
struct dsp_t { bool fused; int cchan; // 1-5 channels, ie: mono, FrontLeft, FrontRight, RearLeft, RearRight, FrontCenter
pset_t *ppset[DSPCHANMAX]; // current preset (1-5 channels)
int ipset; // current ipreset
pset_t *ppsetprev[DSPCHANMAX]; // previous preset (1-5 channels)
int ipsetprev; // previous ipreset
float xfade; // crossfade time between previous preset and new
float xfade_default; // default xfade value, set in DSP_Alloc
bool bexpfade; // true if exponential crossfade
int ipsetsav_oneshot; // previous preset before one-shot preset was set
rmp_t xramp; // crossfade ramp
};
dsp_t dsps[CDSPS];
void DSP_Init( int idsp ) { dsp_t *pdsp;
Assert( idsp < CDSPS );
if (idsp < 0 || idsp >= CDSPS) return; pdsp = &dsps[idsp];
Q_memset( pdsp, 0, sizeof (dsp_t) ); }
void DSP_Free( int idsp ) { dsp_t *pdsp;
Assert( idsp < CDSPS );
if (idsp < 0 || idsp >= CDSPS) return; pdsp = &dsps[idsp];
for (int i = 0; i < pdsp->cchan; i++) { if ( pdsp->ppset[i] ) PSET_Free( pdsp->ppset[i] ); if ( pdsp->ppsetprev[i] ) PSET_Free( pdsp->ppsetprev[i] ); }
Q_memset( pdsp, 0, sizeof (dsp_t) ); }
// Init all dsp processors - called once, during engine startup
void DSP_InitAll ( bool bLoadPresetFile ) { // only load template file on engine startup
if ( bLoadPresetFile ) DSP_LoadPresetFile();
// order is important, don't rearange.
FLT_InitAll(); DLY_InitAll(); RVA_InitAll(); LFOWAV_InitAll(); LFO_InitAll(); CRS_InitAll(); PTC_InitAll(); ENV_InitAll(); EFO_InitAll(); MDY_InitAll(); AMP_InitAll();
PSET_InitAll();
for (int idsp = 0; idsp < CDSPS; idsp++) DSP_Init( idsp ); }
// free all resources associated with dsp - called once, during engine shutdown
void DSP_FreeAll (void) { // order is important, don't rearange.
for (int idsp = 0; idsp < CDSPS; idsp++) DSP_Free( idsp );
AMP_FreeAll(); MDY_FreeAll(); EFO_FreeAll(); ENV_FreeAll(); PTC_FreeAll(); CRS_FreeAll(); LFO_FreeAll(); LFOWAV_FreeAll(); RVA_FreeAll(); DLY_FreeAll(); FLT_FreeAll(); }
// allocate a new dsp processor chain, kill the old processor. Called during dsp init only.
// ipset is new preset
// xfade is crossfade time when switching between presets (milliseconds)
// cchan is how many simultaneous preset channels to allocate (1-4)
// return index to new dsp
int DSP_Alloc( int ipset, float xfade, int cchan ) { dsp_t *pdsp; int i; int idsp; int cchans = clamp( cchan, 1, DSPCHANMAX);
// find free slot
for ( idsp = 0; idsp < CDSPS; idsp++ ) { if ( !dsps[idsp].fused ) break; }
if ( idsp >= CDSPS ) return -1;
pdsp = &dsps[idsp];
DSP_Init ( idsp ); pdsp->fused = true;
pdsp->cchan = cchans;
// allocate a preset processor for each channel
pdsp->ipset = ipset; pdsp->ipsetprev = 0; pdsp->ipsetsav_oneshot = 0;
for (i = 0; i < pdsp->cchan; i++) { pdsp->ppset[i] = PSET_Alloc ( ipset ); pdsp->ppsetprev[i] = NULL; }
// set up crossfade time in seconds
pdsp->xfade = xfade / 1000.0; pdsp->xfade_default = pdsp->xfade;
RMP_SetEnd(&pdsp->xramp);
return idsp; }
// call modulation function of specified processor within dsp preset
// idsp - dsp preset
// channel - channel 1-5 (l,r,rl,rr,fc)
// iproc - which processor to change (normally 0)
// value - new parameter value for processor
// NOTE: routine returns with no result or error if any parameter is invalid.
void DSP_ChangePresetValue( int idsp, int channel, int iproc, float value ) {
dsp_t *pdsp; pset_t *ppset; // preset
prc_Mod_t pfnMod; // modulation function
if (idsp < 0 || idsp >= CDSPS) return; if (channel >= DSPCHANMAX) return; if (iproc >= CPSET_PRCS) return;
// get ptr to processor preset
pdsp = &dsps[idsp];
// assert that this dsp processor has enough separate channels
Assert(channel <= pdsp->cchan);
ppset = pdsp->ppset[channel]; if (!ppset) return;
// get ptr to modulation function
pfnMod = ppset->prcs[iproc].pfnMod;
if (!pfnMod) return;
// call modulation function with new value
pfnMod (ppset->prcs[iproc].pdata, value); }
#define DSP_AUTOMATIC 1 // corresponds to Generic preset
// if dsp_room == DSP_AUTOMATIC, then use dsp_automatic value for dsp
// any subsequent reset of dsp_room will disable automatic room detection.
// return true if automatic room detection is enabled
bool DSP_CheckDspAutoEnabled( void ) { return (dsp_room.GetInt() == DSP_AUTOMATIC); }
// set dsp_automatic preset, used in place of dsp_room when automatic room detection enabled
void DSP_SetDspAuto( int dsp_preset ) { // set dsp_preset into dsp_automatic
dsp_automatic.SetValue( dsp_preset ); }
// wrapper on dsp_room GetInt so that dsp_automatic can override
int dsp_room_GetInt ( void ) { // if dsp_automatic is not enabled, get room
if (! DSP_CheckDspAutoEnabled()) return dsp_room.GetInt();
// automatic room detection is on, get dsp_automatic instead of dsp_room
return dsp_automatic.GetInt(); }
// wrapper on idsp_room preset so that idsp_automatic can override
int Get_idsp_room ( void ) {
// if dsp_automatic is not enabled, get room
if ( !DSP_CheckDspAutoEnabled()) return idsp_room;
// automatic room detection is on, return dsp_automatic preset instead of dsp_room preset
return idsp_automatic; }
// free previous preset if not 0
inline void DSP_FreePrevPreset( dsp_t *pdsp ) { // free previous presets if non-null - ie: rapid change of preset just kills old without xfade
if ( pdsp->ipsetprev ) { for (int i = 0; i < pdsp->cchan; i++) { if ( pdsp->ppsetprev[i] ) { PSET_Free( pdsp->ppsetprev[i] ); pdsp->ppsetprev[i] = NULL; } }
pdsp->ipsetprev = 0; }
}
extern ConVar dsp_mix_min; extern ConVar dsp_mix_max; extern ConVar dsp_db_min; extern ConVar dsp_db_mixdrop;
// alloc new preset if different from current
// xfade from prev to new preset
// free previous preset, copy current into previous, set up xfade from previous to new
void DSP_SetPreset( int idsp, int ipsetnew) { dsp_t *pdsp; pset_t *ppsetnew[DSPCHANMAX];
Assert (idsp >= 0 && idsp < CDSPS);
pdsp = &dsps[idsp];
// validate new preset range
if ( ipsetnew >= g_cpsettemplates || ipsetnew < 0 ) return;
// ignore if new preset is same as current preset
if ( ipsetnew == pdsp->ipset ) return;
// alloc new presets (each channel is a duplicate preset)
Assert (pdsp->cchan <= DSPCHANMAX);
for (int i = 0; i < pdsp->cchan; i++) { ppsetnew[i] = PSET_Alloc ( ipsetnew ); if ( !ppsetnew[i] ) { DevMsg("WARNING: DSP preset failed to allocate.\n"); return; } }
Assert (pdsp);
// free PREVIOUS previous preset if not 0
DSP_FreePrevPreset( pdsp );
for (int i = 0; i < pdsp->cchan; i++) { // current becomes previous
pdsp->ppsetprev[i] = pdsp->ppset[i]; // new becomes current
pdsp->ppset[i] = ppsetnew[i]; } pdsp->ipsetprev = pdsp->ipset; pdsp->ipset = ipsetnew;
if ( idsp == idsp_room || idsp == idsp_automatic ) { // set up new dsp mix min & max, db_min & db_drop params so that new channels get new mix values
// NOTE: only new sounds will get the new mix min/max values set in their dspmix param
// NOTE: so - no crossfade is needed between dspmix and dspmix prev, but this also means
// NOTE: that currently playing ambients will not see changes to dspmix at all.
float mix_min = pdsp->ppset[0]->mix_min; float mix_max = pdsp->ppset[0]->mix_max; float db_min = pdsp->ppset[0]->db_min; float db_mixdrop = pdsp->ppset[0]->db_mixdrop;
dsp_mix_min.SetValue( mix_min ); dsp_mix_max.SetValue( mix_max ); dsp_db_min.SetValue( db_min ); dsp_db_mixdrop.SetValue( db_mixdrop ); }
RMP_SetEnd( &pdsp->xramp ); // make sure previous dsp preset has data
Assert (pdsp->ppsetprev[0]);
// shouldn't be crossfading if current dsp preset == previous dsp preset
Assert (pdsp->ipset != pdsp->ipsetprev);
// if new preset is one-shot, keep previous preset to restore when one-shot times out
// but: don't restore previous one-shots!
pdsp->ipsetsav_oneshot = 0;
if ( PSET_IsOneShot( pdsp->ppset[0] ) && !PSET_IsOneShot( pdsp->ppsetprev[0] ) ) pdsp->ipsetsav_oneshot = pdsp->ipsetprev; // get new xfade time from previous preset (ie: fade out time). if 0 use default. if < 0, use exponential xfade
if ( fabs(pdsp->ppsetprev[0]->fade) > 0.0 ) { pdsp->xfade = fabs(pdsp->ppsetprev[0]->fade); pdsp->bexpfade = pdsp->ppsetprev[0]->fade < 0 ? 1 : 0; } else { // no previous preset - use defauts, set in DSP_Alloc
pdsp->xfade = pdsp->xfade_default; pdsp->bexpfade = false; }
RMP_Init( &(pdsp->xramp), pdsp->xfade, 0, PMAX, false ); }
#define DSP_AUTO_BASE 60 // presets 60-100 in g_psettemplates are reserved as autocreated presets
#define DSP_CAUTO_PRESETS 40 // must be same as DAS_CNODES!!!
// construct a dsp preset based on provided parameters,
// preset is constructed within g_psettemplates[] array.
// return preset #
// parameter batch
struct auto_params_t { // passed in params
bool bskyabove; // true if sky is mostly above player
int width; // max width of room in inches
int length; // max length of room in inches (length always > width)
int height; // max height of room in inches
float fdiffusion; // diffusion of room 0..1.0
float freflectivity; // average reflectivity of all surfaces in room 0..1.0
float surface_refl[6]; // reflectivity for left,right,front,back,ceiling,floor surfaces 0.0 for open surface (sky or no hit)
// derived params
int shape; // ADSP_ROOM, etc 0...4
int size; // ADSP_SIZE_SMALL, etc 0...3
int len; // ADSP_LENGTH_SHORT, etc 0...3
int wid; // ADSP_WIDTH_NARROW, etc 0...3
int ht; // ADSP_HEIGHT_LOW, etc 0...3
int reflectivity; // ADSP_DULL, etc 0..3
int diffusion; // ADSP_EMPTY, etc 0...3
};
// select type 1..5 based on params
// 1:simple reverb
// 2:diffusor + reverb
// 3:diffusor + delay + reverb
// 4:simple delay
// 5:diffusor + delay
#define AROOM_SMALL (10.0 * 12.0) // small room
#define AROOM_MEDIUM (20.0 * 12.0) // medium room
#define AROOM_LARGE (40.0 * 12.0) // large room
#define AROOM_HUGE (100.0 * 12.0) // huge room
#define AROOM_GIGANTIC (200.0 * 12.0) // gigantic room
#define AROOM_DUCT_WIDTH (4.0 * 12.0) // max width for duct
#define AROOM_DUCT_HEIGHT (6.0 * 12.0)
#define AROOM_HALL_WIDTH (8.0 * 12.0) // max width for hall
#define AROOM_HALL_HEIGHT (16.0 * 12.0) // max height for hall
#define AROOM_TUNNEL_WIDTH (20.0 * 12.0) // max width for tunnel
#define AROOM_TUNNEL_HEIGHT (30.0 * 12.0) // max height for tunnel
#define AROOM_STREET_WIDTH (12.0 * 12.0) // min width for street
#define AROOM_SHORT_LENGTH (12.0 * 12.0) // max length for short hall
#define AROOM_MEDIUM_LENGTH (24.0 * 12.0) // min length for medium hall
#define AROOM_LONG_LENGTH (48.0 * 12.0) // min length for long hall
#define AROOM_VLONG_LENGTH (96.0 * 12.0) // min length for very long hall
#define AROOM_XLONG_LENGTH (192.0 * 12.0) // min length for huge hall
#define AROOM_LOW_HEIGHT (4.0 * 12.0) // short ceiling
#define AROOM_MEDIUM_HEIGHT (128) // medium ceiling
#define AROOM_TALL_HEIGHT (18.0 * 12.0) // tall ceiling
#define AROOM_VTALL_HEIGHT (32.0 * 12.0) // very tall ceiling
#define AROOM_XTALL_HEIGHT (64.0 * 12.0) // huge tall ceiling
#define AROOM_NARROW_WIDTH (6.0 * 12.0) // narrow width
#define AROOM_MEDIUM_WIDTH (12.0 * 12.0) // medium width
#define AROOM_WIDE_WIDTH (24.0 * 12.0) // wide width
#define AROOM_VWIDE_WIDTH (48.0 * 12.0) // very wide
#define AROOM_XWIDE_WIDTH (96.0 * 12.0) // huge width
#define BETWEEN(a,b,c) ( ((a) > (b)) && ((a) <= (c)) )
#define ADSP_IsShaft(pa) (pa->height > (3.0 * pa->length))
#define ADSP_IsRoom(pa) (pa->length <= (2.5 * pa->width))
#define ADSP_IsHall(pa) ((pa->length > (2.5 * pa->width)) && (BETWEEN(pa->width, AROOM_DUCT_WIDTH, AROOM_HALL_WIDTH)))
#define ADSP_IsTunnel(pa) ((pa->length > (4.0 * pa->width)) && (pa->width > AROOM_HALL_WIDTH))
#define ADSP_IsDuct(pa) ((pa->length > (4.0 * pa->width)) && (pa->width <= AROOM_DUCT_WIDTH))
#define ADSP_IsCourtyard(pa) (pa->length <= (2.5 * pa->width))
#define ADSP_IsAlley(pa) ((pa->length > (2.5 * pa->width)) && (pa->width <= AROOM_STREET_WIDTH))
#define ADSP_IsStreet(pa) ((pa->length > (2.5 * pa->width)) && (pa->width > AROOM_STREET_WIDTH))
#define ADSP_IsSmallRoom(pa) (pa->length <= AROOM_SMALL)
#define ADSP_IsMediumRoom(pa) ((BETWEEN(pa->length, AROOM_SMALL, AROOM_MEDIUM)) ) // && (BETWEEN(pa->width, AROOM_SMALL, AROOM_MEDIUM)))
#define ADSP_IsLargeRoom(pa) (BETWEEN(pa->length, AROOM_MEDIUM, AROOM_LARGE) ) // && BETWEEN(pa->width, AROOM_MEDIUM, AROOM_LARGE))
#define ADSP_IsHugeRoom(pa) (BETWEEN(pa->length, AROOM_LARGE, AROOM_HUGE) ) // && BETWEEN(pa->width, AROOM_LARGE, AROOM_HUGE))
#define ADSP_IsGiganticRoom(pa) ((pa->length > AROOM_HUGE) ) // && (pa->width > AROOM_HUGE))
#define ADSP_IsShortLength(pa) (pa->length <= AROOM_SHORT_LENGTH)
#define ADSP_IsMediumLength(pa) (BETWEEN(pa->length, AROOM_SHORT_LENGTH, AROOM_MEDIUM_LENGTH))
#define ADSP_IsLongLength(pa) (BETWEEN(pa->length, AROOM_MEDIUM_LENGTH, AROOM_LONG_LENGTH))
#define ADSP_IsVLongLength(pa) (BETWEEN(pa->length, AROOM_LONG_LENGTH, AROOM_VLONG_LENGTH))
#define ADSP_IsXLongLength(pa) (pa->length > AROOM_VLONG_LENGTH)
#define ADSP_IsLowHeight(pa) (pa->height <= AROOM_LOW_HEIGHT)
#define ADSP_IsMediumHeight(pa) (BETWEEN(pa->height, AROOM_LOW_HEIGHT, AROOM_MEDIUM_HEIGHT))
#define ADSP_IsTallHeight(pa) (BETWEEN(pa->height, AROOM_MEDIUM_HEIGHT, AROOM_TALL_HEIGHT))
#define ADSP_IsVTallHeight(pa) (BETWEEN(pa->height, AROOM_TALL_HEIGHT, AROOM_VTALL_HEIGHT))
#define ADSP_IsXTallHeight(pa) (pa->height > AROOM_VTALL_HEIGHT)
#define ADSP_IsNarrowWidth(pa) (pa->width <= AROOM_NARROW_WIDTH)
#define ADSP_IsMediumWidth(pa) (BETWEEN(pa->width, AROOM_NARROW_WIDTH, AROOM_MEDIUM_WIDTH))
#define ADSP_IsWideWidth(pa) (BETWEEN(pa->width, AROOM_MEDIUM_WIDTH, AROOM_WIDE_WIDTH))
#define ADSP_IsVWideWidth(pa) (BETWEEN(pa->width, AROOM_WIDE_WIDTH, AROOM_VWIDE_WIDTH))
#define ADSP_IsXWideWidth(pa) (pa->width > AROOM_VWIDE_WIDTH)
#define ADSP_IsInside(pa) (!(pa->bskyabove))
// room diffusion
#define ADSP_EMPTY 0
#define ADSP_SPARSE 1
#define ADSP_CLUTTERED 2
#define ADSP_FULL 3
#define ADSP_DIFFUSION_MAX 4
#define AROOM_DIF_EMPTY 0.01 // 1% of space by volume is other objects
#define AROOM_DIF_SPARSE 0.1 // 10% "
#define AROOM_DIF_CLUTTERED 0.3 // 30% "
#define AROOM_DIF_FULL 0.5 // 50% "
#define ADSP_IsEmpty(pa) (pa->fdiffusion <= AROOM_DIF_EMPTY)
#define ADSP_IsSparse(pa) (BETWEEN(pa->fdiffusion, AROOM_DIF_EMPTY, AROOM_DIF_SPARSE))
#define ADSP_IsCluttered(pa) (BETWEEN(pa->fdiffusion, AROOM_DIF_SPARSE, AROOM_DIF_CLUTTERED))
#define ADSP_IsFull(pa) (pa->fdiffusion > AROOM_DIF_CLUTTERED)
#define ADSP_IsDiffuse(pa) (pa->diffusion > ADSP_SPARSE)
// room acoustic reflectivity
// tile 0.3 * 3.3 = 0.99
// metal 0.25 * 3.3 = 0.83
// concrete,rock,brick,glass,gravel 0.2 * 3.3 = 0.66
// metal panel/vent, wood, water 0.1 * 3.3 = 0.33
// carpet,sand,snow,dirt 0.01 * 3.3 = 0.03
#define ADSP_DULL 0
#define ADSP_FLAT 1
#define ADSP_REFLECTIVE 2
#define ADSP_BRIGHT 3
#define ADSP_REFLECTIVITY_MAX 4
#define AROOM_REF_DULL 0.04
#define AROOM_REF_FLAT 0.50
#define AROOM_REF_REFLECTIVE 0.80
#define AROOM_REF_BRIGHT 0.99
#define ADSP_IsDull(pa) (pa->freflectivity <= AROOM_REF_DULL)
#define ADSP_IsFlat(pa) (BETWEEN(pa->freflectivity, AROOM_REF_DULL, AROOM_REF_FLAT))
#define ADSP_IsReflective(pa) (BETWEEN(pa->freflectivity, AROOM_REF_FLAT, AROOM_REF_REFLECTIVE))
#define ADSP_IsBright(pa) (pa->freflectivity > AROOM_REF_REFLECTIVE)
#define ADSP_IsRefl(pa) (pa->reflectivity > ADSP_FLAT)
// room shapes
#define ADSP_ROOM 0
#define ADSP_DUCT 1
#define ADSP_HALL 2
#define ADSP_TUNNEL 3
#define ADSP_STREET 4
#define ADSP_ALLEY 5
#define ADSP_COURTYARD 6
#define ADSP_OPEN_SPACE 7 // NOTE: 7..10 must remain in order !!!
#define ADSP_OPEN_WALL 8
#define ADSP_OPEN_STREET 9
#define ADSP_OPEN_COURTYARD 10
// room sizes
#define ADSP_SIZE_SMALL 0 // NOTE: must remain 0..4!!!
#define ADSP_SIZE_MEDIUM 1
#define ADSP_SIZE_LARGE 2
#define ADSP_SIZE_HUGE 3
#define ADSP_SIZE_GIGANTIC 4
#define ADSP_SIZE_MAX 5
#define ADSP_LENGTH_SHORT 0
#define ADSP_LENGTH_MEDIUM 1
#define ADSP_LENGTH_LONG 2
#define ADSP_LENGTH_VLONG 3
#define ADSP_LENGTH_XLONG 4
#define ADSP_LENGTH_MAX 5
#define ADSP_WIDTH_NARROW 0
#define ADSP_WIDTH_MEDIUM 1
#define ADSP_WIDTH_WIDE 2
#define ADSP_WIDTH_VWIDE 3
#define ADSP_WIDTH_XWIDE 4
#define ADSP_WIDTH_MAX 5
#define ADSP_HEIGHT_LOW 0
#define ADSP_HEIGTH_MEDIUM 1
#define ADSP_HEIGHT_TALL 2
#define ADSP_HEIGHT_VTALL 3
#define ADSP_HEIGHT_XTALL 4
#define ADSP_HEIGHT_MAX 5
// convert numeric size params to #defined size params
void ADSP_GetSize( auto_params_t *pa ) { pa->size = ((ADSP_IsSmallRoom(pa) ? 1 : 0) * ADSP_SIZE_SMALL) + ((ADSP_IsMediumRoom(pa) ? 1 : 0) * ADSP_SIZE_MEDIUM) + ((ADSP_IsLargeRoom(pa) ? 1 : 0) * ADSP_SIZE_LARGE) + ((ADSP_IsHugeRoom(pa) ? 1 : 0) * ADSP_SIZE_HUGE) + ((ADSP_IsGiganticRoom(pa) ? 1 : 0) * ADSP_SIZE_GIGANTIC);
pa->len = ((ADSP_IsShortLength(pa) ? 1 : 0) * ADSP_LENGTH_SHORT) + ((ADSP_IsMediumLength(pa) ? 1 : 0) * ADSP_LENGTH_MEDIUM) + ((ADSP_IsLongLength(pa) ? 1 : 0) * ADSP_LENGTH_LONG) + ((ADSP_IsVLongLength(pa) ? 1 : 0) * ADSP_LENGTH_VLONG) + ((ADSP_IsXLongLength(pa) ? 1 : 0) * ADSP_LENGTH_XLONG); pa->wid = ((ADSP_IsNarrowWidth(pa) ? 1 : 0) * ADSP_WIDTH_NARROW) + ((ADSP_IsMediumWidth(pa) ? 1 : 0) * ADSP_WIDTH_MEDIUM) + ((ADSP_IsWideWidth(pa) ? 1 : 0) * ADSP_WIDTH_WIDE) + ((ADSP_IsVWideWidth(pa) ? 1 : 0) * ADSP_WIDTH_VWIDE) + ((ADSP_IsXWideWidth(pa) ? 1 : 0) * ADSP_WIDTH_XWIDE);
pa->ht = ((ADSP_IsLowHeight(pa) ? 1 : 0) * ADSP_HEIGHT_LOW) + ((ADSP_IsMediumHeight(pa) ? 1 : 0) * ADSP_HEIGTH_MEDIUM) + ((ADSP_IsTallHeight(pa) ? 1 : 0) * ADSP_HEIGHT_TALL) + ((ADSP_IsVTallHeight(pa) ? 1 : 0) * ADSP_HEIGHT_VTALL) + ((ADSP_IsXTallHeight(pa) ? 1 : 0) * ADSP_HEIGHT_XTALL);
pa->reflectivity = ((ADSP_IsDull(pa) ? 1 : 0) * ADSP_DULL) + ((ADSP_IsFlat(pa) ? 1 : 0) * ADSP_FLAT) + ((ADSP_IsReflective(pa) ? 1 : 0) * ADSP_REFLECTIVE) + ((ADSP_IsBright(pa) ? 1 : 0) * ADSP_BRIGHT);
pa->diffusion = ((ADSP_IsEmpty(pa) ? 1 : 0) * ADSP_EMPTY) + ((ADSP_IsSparse(pa) ? 1 : 0) * ADSP_SPARSE) + ((ADSP_IsCluttered(pa) ? 1 : 0) * ADSP_CLUTTERED) + ((ADSP_IsFull(pa) ? 1 : 0) * ADSP_FULL);
Assert(pa->size < ADSP_SIZE_MAX); Assert(pa->len < ADSP_LENGTH_MAX); Assert(pa->wid < ADSP_WIDTH_MAX); Assert(pa->ht < ADSP_HEIGHT_MAX); Assert(pa->reflectivity < ADSP_REFLECTIVITY_MAX); Assert(pa->diffusion < ADSP_DIFFUSION_MAX);
if ( pa->shape != ADSP_COURTYARD && pa->shape != ADSP_OPEN_COURTYARD ) { // fix up size for streets, alleys, halls, ducts, tunnelsy
if (pa->shape == ADSP_STREET || pa->shape == ADSP_ALLEY ) pa->size = pa->wid; else pa->size = (pa->len + pa->wid) / 2;
}
}
void ADSP_GetOutsideSize( auto_params_t *pa ) { ADSP_GetSize( pa ); }
// return # of sides that had max length or sky hits (out of 6 sides).
int ADSP_COpenSides( auto_params_t *pa ) { int count = 0;
// only look at left,right,front,back walls - ignore floor, ceiling
for (int i = 0; i < 4; i++) { if (pa->surface_refl[i] == 0.0) count++; }
return count; }
// given auto params, return shape and size of room
void ADSP_GetAutoShape( auto_params_t *pa ) {
// INSIDE:
// shapes: duct, hall, tunnel, shaft (vertical duct, hall or tunnel)
// sizes: short->long, narrow->wide, low->tall
// shapes: room
// sizes: small->large, low->tall
// OUTSIDE:
// shapes: street, alley
// sizes: short->long, narrow->wide
// shapes: courtyard
// sizes: small->large
// shapes: open_space, wall, open_street, open_corner, open_courtyard
// sizes: open, narrow->wide
bool bshaft = false; int t;
if (ADSP_IsInside(pa)) { if (ADSP_IsShaft(pa)) { // temp swap height and length
bshaft = true; t = pa->height; pa->height = pa->length; pa->length = t; if (das_debug.GetInt() > 1) DevMsg("VERTICAL SHAFT Detected \n"); }
// get shape
if (ADSP_IsDuct(pa)) { pa->shape = ADSP_DUCT; ADSP_GetSize( pa ); if (das_debug.GetInt() > 1) DevMsg("DUCT Detected \n"); goto autoshape_exit; }
if (ADSP_IsHall(pa)) { // get size
pa->shape = ADSP_HALL; ADSP_GetSize( pa );
if (das_debug.GetInt() > 1) DevMsg("HALL Detected \n");
goto autoshape_exit; }
if (ADSP_IsTunnel(pa)) { // get size
pa->shape = ADSP_TUNNEL; ADSP_GetSize( pa );
if (das_debug.GetInt() > 1) DevMsg("TUNNEL Detected \n");
goto autoshape_exit; } // default
// (ADSP_IsRoom(pa))
{ // get size
pa->shape = ADSP_ROOM; ADSP_GetSize( pa );
if (das_debug.GetInt() > 1) DevMsg("ROOM Detected \n");
goto autoshape_exit; } }
// outside:
if (ADSP_COpenSides(pa) > 0) // side hit sky, or side has max length
{ // get shape - courtyard, street, wall or open space
// 10..7
pa->shape = ADSP_OPEN_COURTYARD - (ADSP_COpenSides(pa) - 1); ADSP_GetOutsideSize( pa ); if (das_debug.GetInt() > 1) DevMsg("OPEN SIDED OUTDOOR AREA Detected \n");
goto autoshape_exit; }
// all sides closed:
// get shape - closed street or alley or courtyard
if (ADSP_IsCourtyard(pa)) { pa->shape = ADSP_COURTYARD; ADSP_GetOutsideSize( pa ); if (das_debug.GetInt() > 1) DevMsg("OUTSIDE COURTYARD Detected \n");
goto autoshape_exit; }
if (ADSP_IsAlley(pa)) { pa->shape = ADSP_ALLEY; ADSP_GetOutsideSize( pa );
if (das_debug.GetInt() > 1) DevMsg("OUTSIDE ALLEY Detected \n"); goto autoshape_exit; } // default to 'street' if sides are closed
// if (ADSP_IsStreet(pa))
{ pa->shape = ADSP_STREET; ADSP_GetOutsideSize( pa ); if (das_debug.GetInt() > 1) DevMsg("OUTSIDE STREET Detected \n"); goto autoshape_exit; } autoshape_exit:
// swap height & length if needed
if (bshaft) { t = pa->height; pa->height = pa->length; pa->length = t; } }
int MapReflectivityToDLYCutoff[] = { 1000, // DULL
2000, // FLAT
4000, // REFLECTIVE
6000 // BRIGHT
};
float MapSizeToDLYFeedback[] = { 0.9, // 0.6, // SMALL
0.8, // 0.5, // MEDIUM
0.7, // 0.4, // LARGE
0.6, // 0.3, // HUGE
0.5, // 0.2, // GIGANTIC
};
void ADSP_SetupAutoDelay( prc_t *pprc_dly, auto_params_t *pa ) { // shapes:
// inside: duct, long hall, long tunnel, large room
// outside: open courtyard, street wall, space
// outside: closed courtyard, alley, street
// size 0..4
// len 0..3
// wid 0..3
// reflectivity: 0..3
// diffusion 0..3
// dtype: delay type DLY_PLAIN, DLY_LOWPASS, DLY_ALLPASS
// delay: delay in milliseconds (room max size in feet)
// feedback: feedback 0-1.0
// gain: final gain of output stage, 0-1.0
int size = pa->length * 2.0; if (pa->shape == ADSP_ALLEY || pa->shape == ADSP_STREET || pa->shape == ADSP_OPEN_STREET) size = pa->width * 2.0;
pprc_dly->type = PRC_DLY;
pprc_dly->prm[dly_idtype] = DLY_LOWPASS; // delay with feedback
pprc_dly->prm[dly_idelay] = clamp((size / 12.0), 5.0, 500.0); pprc_dly->prm[dly_ifeedback] = MapSizeToDLYFeedback[pa->len];
// reduce gain based on distance reflection travels
// float g = 1.0 - ( clamp(pprc_dly->prm[dly_idelay], 10.0, 1000.0) / (1000.0 - 10.0) );
// pprc_dly->prm[dly_igain] = g;
pprc_dly->prm[dly_iftype] = FLT_LP; if (ADSP_IsInside(pa)) pprc_dly->prm[dly_icutoff] = MapReflectivityToDLYCutoff[pa->reflectivity]; else pprc_dly->prm[dly_icutoff] = (int)((float)(MapReflectivityToDLYCutoff[pa->reflectivity]) * 0.75);
pprc_dly->prm[dly_iqwidth] = 0;
pprc_dly->prm[dly_iquality] = QUA_LO;
float l = clamp((pa->length * 2.0 / 12.0), 14.0, 500.0); float w = clamp((pa->width * 2.0 / 12.0), 14.0, 500.0);
// convert to multitap delay
pprc_dly->prm[dly_idtype] = DLY_LOWPASS_4TAP;
pprc_dly->prm[dly_idelay] = l; pprc_dly->prm[dly_itap1] = w; pprc_dly->prm[dly_itap2] = l; // max(7, l * 0.7 );
pprc_dly->prm[dly_itap3] = l; // max(7, w * 0.7 );
pprc_dly->prm[dly_igain] = 1.0; }
int MapReflectivityToRVACutoff[] = { 1000, // DULL
2000, // FLAT
4000, // REFLECTIVE
6000 // BRIGHT
};
float MapSizeToRVANumDelays[] = { 3, // SMALL 3 reverbs
6, // MEDIUM 6 reverbs
6, // LARGE 6 reverbs
9, // HUGE 9 reverbs
12, // GIGANTIC 12 reverbs
};
float MapSizeToRVAFeedback[] = { 0.75, // SMALL
0.8, // MEDIUM
0.9, // LARGE
0.95, // HUGE
0.98, // GIGANTIC
};
void ADSP_SetupAutoReverb( prc_t *pprc_rva, auto_params_t *pa ) { // shape: hall, tunnel or room
// size 0..4
// reflectivity: 0..3
// diffusion 0..3
// size: 0-2.0 scales nominal delay parameters (18 to 47 ms * scale = delay)
// numdelays: 0-12 controls # of parallel or series delays
// decay: 0-2.0 scales feedback parameters (.7 to .9 * scale/2.0 = feedback)
// fparallel: if true, filters are built into delays, otherwise filter output only
// fmoddly: if true, all delays are modulating delays
float gain = 1.0;
pprc_rva->type = PRC_RVA;
pprc_rva->prm[rva_size_max] = 50.0; pprc_rva->prm[rva_size_min] = 30.0; if (ADSP_IsRoom(pa)) pprc_rva->prm[rva_inumdelays] = MapSizeToRVANumDelays[pa->size]; else pprc_rva->prm[rva_inumdelays] = MapSizeToRVANumDelays[pa->len];
pprc_rva->prm[rva_ifeedback] = 0.9; pprc_rva->prm[rva_icutoff] = MapReflectivityToRVACutoff[pa->reflectivity]; pprc_rva->prm[rva_ifparallel] = 1; pprc_rva->prm[rva_imoddly] = ADSP_IsEmpty(pa) ? 0 : 4; pprc_rva->prm[rva_imodrate] = 3.48;
pprc_rva->prm[rva_iftaps] = 0; // 0.1 // use extra delay taps to increase density
pprc_rva->prm[rva_width] = clamp( ((float)(pa->width) / 12.0), 6.0, 500.0); // in feet
pprc_rva->prm[rva_depth] = clamp( ((float)(pa->length) / 12.0), 6.0, 500.0); pprc_rva->prm[rva_height] = clamp( ((float)(pa->height) / 12.0), 6.0, 500.0);
// room
pprc_rva->prm[rva_fbwidth] = 0.9; // MapSizeToRVAFeedback[pa->size]; // larger size = more feedback
pprc_rva->prm[rva_fbdepth] = 0.9; // MapSizeToRVAFeedback[pa->size];
pprc_rva->prm[rva_fbheight] = 0.5; // MapSizeToRVAFeedback[pa->size];
// feedback is based on size of room:
if (ADSP_IsInside(pa)) { if (pa->shape == ADSP_HALL) { pprc_rva->prm[rva_fbwidth] = 0.7; //MapSizeToRVAFeedback[pa->wid];
pprc_rva->prm[rva_fbdepth] = -0.5; //MapSizeToRVAFeedback[pa->len];
pprc_rva->prm[rva_fbheight] = 0.3; //MapSizeToRVAFeedback[pa->ht];
}
if (pa->shape == ADSP_TUNNEL) { pprc_rva->prm[rva_fbwidth] = 0.9; pprc_rva->prm[rva_fbdepth] = -0.8; // fixed pre-delay, no feedback
pprc_rva->prm[rva_fbheight] = 0.3; } } else { if (pa->shape == ADSP_ALLEY) { pprc_rva->prm[rva_fbwidth] = 0.9; pprc_rva->prm[rva_fbdepth] = -0.8; // fixed pre-delay, no feedback
pprc_rva->prm[rva_fbheight] = 0.0; } }
if (!ADSP_IsInside(pa)) pprc_rva->prm[rva_fbheight] = 0.0; pprc_rva->prm[rva_igain] = gain; }
// diffusor templates for auto create
// size: 0-1.0 scales all delays (13ms to 41ms * scale = delay)
// numdelays: 0-4.0 controls # of series delays
// decay: 0-1.0 scales all feedback parameters
// prctype size #dly feedback
#if 0
#define PRC_DFRA_S {PRC_DFR, {0.5, 2, 0.10}, NULL,NULL,NULL,NULL,NULL} // S room
#define PRC_DFRA_M {PRC_DFR, {0.75, 2, 0.12}, NULL,NULL,NULL,NULL,NULL} // M room
#define PRC_DFRA_L {PRC_DFR, {1.0, 3, 0.13}, NULL,NULL,NULL,NULL,NULL} // L room
#define PRC_DFRA_VL {PRC_DFR, {1.0, 3, 0.15}, NULL,NULL,NULL,NULL,NULL} // VL room
prc_t g_prc_dfr_auto[] = {PRC_DFRA_S, PRC_DFRA_M, PRC_DFRA_L, PRC_DFRA_VL, PRC_DFRA_VL};
//$BUGBUGBUG: I think this should be sizeof(prc_t), not sizeof(pset_t)...
#define CDFRTEMPLATES (sizeof(g_prc_dfr_auto)/sizeof(pset_t)) // number of diffusor templates
// copy diffusor template from preset list, based on room size
void ADSP_SetupAutoDiffusor( prc_t *pprc_dfr, auto_params_t *pa ) { int i = clamp(pa->size, 0, (int)CDFRTEMPLATES - 1);
// copy diffusor preset based on size
*pprc_dfr = g_prc_dfr_auto[i]; } #endif
// return index to processor given processor type and preset
// skips N processors of similar type
// returns -1 if type not found
int ADSP_FindProc( pset_t *ppset, int proc_type, int skip ) { int skipcount = skip;
for (int i = 0; i < ppset->cprcs; i++) { // look for match on processor type
if ( ppset->prcs[i].type == proc_type ) { // skip first N procs of similar type,
// return index to processor
if (!skipcount) return i;
skipcount--; }
}
return -1; }
// interpolate parameter:
// pnew - target preset
// pmin - preset with parameter with min value
// pmax - preset with parameter with max value
// proc_type - type of processor to look for ie: PRC_RVA or PRC_DLY
// skipprocs - skip n processors of type
// iparam - which parameter within processor to interpolate
// index -
// index_max: use index/index_max as interpolater between pmin param and pmax param
// if bexp is true, interpolate exponentially as (index/index_max)^2
// NOTE: returns with no result if processor type is not found in all presets.
void ADSP_InterpParam( pset_t *pnew, pset_t *pmin, pset_t *pmax, int proc_type, int skipprocs, int iparam, int index, int index_max, bool bexp ) { // find processor index in pnew
int iproc_new = ADSP_FindProc( pnew, proc_type, skipprocs); int iproc_min = ADSP_FindProc( pmin, proc_type, skipprocs); int iproc_max = ADSP_FindProc( pmax, proc_type, skipprocs);
// make sure processor type found in all presets
if ( iproc_new < 0 || iproc_min < 0 || iproc_max < 0 ) return; float findex = (float)index/(float)index_max; float vmin = pmin->prcs[iproc_min].prm[iparam]; float vmax = pmax->prcs[iproc_max].prm[iparam]; float vinterp;
// interpolate
if (!bexp) vinterp = vmin + (vmax - vmin) * findex; else vinterp = vmin + (vmax - vmin) * findex * findex;
pnew->prcs[iproc_new].prm[iparam] = vinterp;
return; }
// directly set parameter
void ADSP_SetParam( pset_t *pnew, int proc_type, int skipprocs, int iparam, float value ) { int iproc_new = ADSP_FindProc( pnew, proc_type, skipprocs);
if (iproc_new >= 0) pnew->prcs[iproc_new].prm[iparam] = value; }
// directly set parameter if min or max is negative
void ADSP_SetParamIfNegative( pset_t *pnew, pset_t *pmin, pset_t *pmax, int proc_type, int skipprocs, int iparam, int index, int index_max, bool bexp, float value ) { // find processor index in pnew
int iproc_new = ADSP_FindProc( pnew, proc_type, skipprocs); int iproc_min = ADSP_FindProc( pmin, proc_type, skipprocs); int iproc_max = ADSP_FindProc( pmax, proc_type, skipprocs);
// make sure processor type found in all presets
if ( iproc_new < 0 || iproc_min < 0 || iproc_max < 0 ) return; float vmin = pmin->prcs[iproc_min].prm[iparam]; float vmax = pmax->prcs[iproc_max].prm[iparam];
if ( vmin < 0.0 || vmax < 0.0 ) ADSP_SetParam( pnew, proc_type, skipprocs, iparam, value ); else ADSP_InterpParam( pnew, pmin, pmax, proc_type, skipprocs, iparam, index, index_max, bexp);
return; }
// given min and max preset and auto parameters, create new preset
// NOTE: the # and type of processors making up pmin and pmax presets must be identical!
void ADSP_InterpolatePreset( pset_t *pnew, pset_t *pmin, pset_t *pmax, auto_params_t *pa, int iskip ) { int i;
// if size > mid size, then copy basic processors from MAX preset,
// otherwise, copy from MIN preset
if ( !iskip ) { // only copy on 1st call
if ( pa->size > ADSP_SIZE_MEDIUM ) { *pnew = *pmax; } else { *pnew = *pmin; } }
// DFR
// interpolate all DFR params on size
for (i = 0; i < dfr_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_DFR, iskip, i, pa->size, ADSP_SIZE_MAX , 0);
// RVA
// interpolate size_max, size_min, feedback, #delays, moddly, imodrate, based on ap size
ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_ifeedback, pa->size, ADSP_SIZE_MAX, 0); ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_size_min, pa->size, ADSP_SIZE_MAX, 1); ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_size_max, pa->size, ADSP_SIZE_MAX, 1); ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_igain, pa->size, ADSP_SIZE_MAX, 0); ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_inumdelays, pa->size, ADSP_SIZE_MAX , 0); ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_imoddly, pa->size, ADSP_SIZE_MAX , 0); ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_imodrate, pa->size, ADSP_SIZE_MAX , 0);
// interpolate width,depth,height based on ap width length & height - exponential interpolation
// if pmin or pmax parameters are < 0, directly set value from w/l/h
float w = clamp( ((float)(pa->width) / 12.0), 6.0, 500.0); // in feet
float l = clamp( ((float)(pa->length) / 12.0), 6.0, 500.0); float h = clamp( ((float)(pa->height) / 12.0), 6.0, 500.0);
ADSP_SetParamIfNegative( pnew, pmin, pmax, PRC_RVA, iskip, rva_width, pa->wid, ADSP_WIDTH_MAX, 1, w); ADSP_SetParamIfNegative( pnew, pmin, pmax, PRC_RVA, iskip, rva_depth, pa->len, ADSP_LENGTH_MAX, 1, l); ADSP_SetParamIfNegative( pnew, pmin, pmax, PRC_RVA, iskip, rva_height, pa->ht, ADSP_HEIGHT_MAX, 1, h);
// interpolate w/d/h feedback based on ap w/d/f
ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_fbwidth, pa->wid, ADSP_WIDTH_MAX , 0); ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_fbdepth, pa->len, ADSP_LENGTH_MAX , 0); ADSP_InterpParam( pnew, pmin, pmax, PRC_RVA, iskip, rva_fbheight, pa->ht, ADSP_HEIGHT_MAX , 0);
// interpolate cutoff based on ap reflectivity
// NOTE: cutoff goes from max to min! ie: small bright - large dull
ADSP_InterpParam( pnew, pmax, pmin, PRC_RVA, iskip, rva_icutoff, pa->reflectivity, ADSP_REFLECTIVITY_MAX , 0);
// don't interpolate: fparallel, ftaps
// DLY
// directly set delay value from pa->length if pmin or pmax value is < 0
l = clamp((pa->length * 2.0 / 12.0), 14.0, 500.0); w = clamp((pa->width * 2.0 / 12.0), 14.0, 500.0);
ADSP_SetParamIfNegative( pnew, pmin, pmax, PRC_DLY, iskip, dly_idelay, pa->len, ADSP_LENGTH_MAX, 1, l);
// interpolate feedback, gain, based on max size (length)
ADSP_InterpParam( pnew, pmin, pmax, PRC_DLY, iskip, dly_ifeedback, pa->len, ADSP_LENGTH_MAX , 0); ADSP_InterpParam( pnew, pmin, pmax, PRC_DLY, iskip, dly_igain, pa->len, ADSP_LENGTH_MAX , 0);
// directly set tap value from pa->width if pmin or pmax value is < 0
ADSP_SetParamIfNegative( pnew, pmin, pmax, PRC_DLY, iskip, dly_itap1, pa->len, ADSP_LENGTH_MAX, 1, w); ADSP_SetParamIfNegative( pnew, pmin, pmax, PRC_DLY, iskip, dly_itap2, pa->len, ADSP_LENGTH_MAX, 1, l); ADSP_SetParamIfNegative( pnew, pmin, pmax, PRC_DLY, iskip, dly_itap3, pa->len, ADSP_LENGTH_MAX, 1, l); // interpolate cutoff and qwidth based on reflectivity NOTE: this can affect gain!
// NOTE: cutoff goes from max to min! ie: small bright - large dull
ADSP_InterpParam( pnew, pmax, pmin, PRC_DLY, iskip, dly_icutoff, pa->len, ADSP_LENGTH_MAX , 0); ADSP_InterpParam( pnew, pmax, pmin, PRC_DLY, iskip, dly_iqwidth, pa->len, ADSP_LENGTH_MAX , 0); // interpolate all other parameters for all other processor types based on size
// PRC_MDY, PRC_AMP, PRC_FLT, PTC, CRS, ENV, EFO, LFO
for (i = 0; i < mdy_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_MDY, iskip, i, pa->len, ADSP_LENGTH_MAX , 0); for (i = 0; i < amp_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_AMP, iskip, i, pa->size, ADSP_SIZE_MAX , 0);
for (i = 0; i < flt_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_FLT, iskip, i, pa->size, ADSP_SIZE_MAX , 0);
for (i = 0; i < ptc_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_PTC, iskip, i, pa->size, ADSP_SIZE_MAX , 0);
for (i = 0; i < crs_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_CRS, iskip, i, pa->size, ADSP_SIZE_MAX , 0);
for (i = 0; i < env_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_ENV, iskip, i, pa->size, ADSP_SIZE_MAX , 0);
for (i = 0; i < efo_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_EFO, iskip, i, pa->size, ADSP_SIZE_MAX , 0);
for (i = 0; i < lfo_cparam; i++) ADSP_InterpParam( pnew, pmin, pmax, PRC_LFO, iskip, i, pa->size, ADSP_SIZE_MAX , 0);
}
// these convars store the index to the first preset for each shape type in dsp_presets.txt
ConVar adsp_room_min ("adsp_room_min", "102"); ConVar adsp_duct_min ("adsp_duct_min", "106"); ConVar adsp_hall_min ("adsp_hall_min", "110"); ConVar adsp_tunnel_min ("adsp_tunnel_min", "114"); ConVar adsp_street_min ("adsp_street_min", "118"); ConVar adsp_alley_min ("adsp_alley_min", "122"); ConVar adsp_courtyard_min ("adsp_courtyard_min", "126"); ConVar adsp_openspace_min ("adsp_openspace_min", "130"); ConVar adsp_openwall_min ("adsp_openwall_min", "130"); ConVar adsp_openstreet_min ("adsp_openstreet_min", "118"); ConVar adsp_opencourtyard_min ("adsp_opencourtyard_min", "126");
// given room parameters, construct and return a dsp preset representing the room.
// bskyabove, width, length, height, fdiffusion, freflectivity are all passed-in room parameters
// psurf_refl is a passed-in array of reflectivity values for 6 surfaces
// inode is the location within g_psettemplates[] that the dsp preset will be constructed (inode = dsp preset#)
// cnode should always = DSP_CAUTO_PRESETS
// returns idsp preset.
int DSP_ConstructPreset( bool bskyabove, int width, int length, int height, float fdiffusion, float freflectivity, float *psurf_refl, int inode, int cnodes ) { auto_params_t ap; auto_params_t *pa; pset_t new_pset; // preset
pset_t pset_min; pset_t pset_max;
int ipreset; int ipset_min; int ipset_max;
if (inode >= DSP_CAUTO_PRESETS) { Assert(false); // check DAS_CNODES == DSP_CAUTO_PRESETS!!!
return 0; } // fill parameter struct
ap.bskyabove = bskyabove; ap.width = width; ap.length = length; ap.height = height; ap.fdiffusion = fdiffusion; ap.freflectivity = freflectivity;
for (int i = 0; i < 6; i++) ap.surface_refl[i] = psurf_refl[i];
if (ap.bskyabove) ap.surface_refl[4] = 0.0; // select shape, size based on params
ADSP_GetAutoShape( &ap ); // set up min/max presets based on shape
switch ( ap.shape ) { default: case ADSP_ROOM: ipset_min = adsp_room_min.GetInt(); break; case ADSP_DUCT: ipset_min = adsp_duct_min.GetInt(); break; case ADSP_HALL: ipset_min = adsp_hall_min.GetInt(); break; case ADSP_TUNNEL: ipset_min = adsp_tunnel_min.GetInt(); break; case ADSP_STREET: ipset_min = adsp_street_min.GetInt(); break; case ADSP_ALLEY: ipset_min = adsp_alley_min.GetInt(); break; case ADSP_COURTYARD: ipset_min = adsp_courtyard_min.GetInt(); break; case ADSP_OPEN_SPACE: ipset_min = adsp_openspace_min.GetInt(); break; case ADSP_OPEN_WALL: ipset_min = adsp_openwall_min.GetInt(); break; case ADSP_OPEN_STREET: ipset_min = adsp_openstreet_min.GetInt(); break; case ADSP_OPEN_COURTYARD: ipset_min = adsp_opencourtyard_min.GetInt(); break; }
// presets in dsp_presets.txt are ordered as:
// <shape><empty><min>
// <shape><empty><max>
// <shape><diffuse><min>
// <shape><diffuse><max>
pa = ≈ if ( ADSP_IsDiffuse(pa) ) ipset_min += 2;
ipset_max = ipset_min + 1;
pset_min = g_psettemplates[ipset_min]; pset_max = g_psettemplates[ipset_max];
// given min and max preset and auto parameters, create new preset
// interpolate between 1st instances of each processor type (ie: PRC_DLY) appearing in preset
ADSP_InterpolatePreset( &new_pset, &pset_min, &pset_max, &ap, 0 );
// interpolate between 2nd instances of each processor type (ie: PRC_DLY) appearing in preset
ADSP_InterpolatePreset( &new_pset, &pset_min, &pset_max, &ap, 1 );
// copy constructed preset back into node's template location
ipreset = DSP_AUTO_BASE + inode;
g_psettemplates[ipreset] = new_pset;
return ipreset; }
///////////////////////////////////////
// Helpers: called only from DSP_Process
///////////////////////////////////////
// return true if batch processing version of preset exists
inline bool FBatchPreset( pset_t *ppset ) { switch (ppset->type) { case PSET_LINEAR: return true; case PSET_SIMPLE: return true; default: return false; } }
// Helper: called only from DSP_Process
// mix front stereo buffer to mono buffer, apply dsp fx
inline void DSP_ProcessStereoToMono(dsp_t *pdsp, portable_samplepair_t *pbfront, portable_samplepair_t *pbrear, int sampleCount, bool bcrossfading ) { portable_samplepair_t *pbf = pbfront; // pointer to buffer of front stereo samples to process
int count = sampleCount; int av; int x;
if ( !bcrossfading ) { if ( !pdsp->ipset ) return;
if ( FBatchPreset(pdsp->ppset[0])) { // convert Stereo to Mono in place, then batch process fx: perf KDB
// front->left + front->right / 2 into front->left, front->right duplicated.
while ( count-- ) { pbf->left = (pbf->left + pbf->right) >> 1; pbf++; } // process left (mono), duplicate output into right
PSET_GetNextN( pdsp->ppset[0], pbfront, sampleCount, OP_LEFT_DUPLICATE); } else { // avg left and right -> mono fx -> duplcate out left and right
while ( count-- ) { av = ( ( pbf->left + pbf->right ) >> 1 ); x = PSET_GetNext( pdsp->ppset[0], av ); x = CLIP_DSP( x ); pbf->left = pbf->right = x; pbf++; } } return; }
// crossfading to current preset from previous preset
if ( bcrossfading ) { int r; int fl; int fr; int flp; int frp; int xf_fl; int xf_fr; bool bexp = pdsp->bexpfade; bool bfadetostereo = (pdsp->ipset == 0); bool bfadefromstereo = (pdsp->ipsetprev == 0);
Assert ( !(bfadetostereo && bfadefromstereo) ); // don't call if ipset & ipsetprev both 0!
if ( bfadetostereo || bfadefromstereo ) { // special case if fading to or from preset 0, stereo passthrough
while ( count-- ) { av = ( ( pbf->left + pbf->right ) >> 1 );
// get current preset values
if ( pdsp->ipset ) { fl = fr = PSET_GetNext( pdsp->ppset[0], av ); } else { fl = pbf->left; fr = pbf->right; } // get previous preset values
if ( pdsp->ipsetprev ) { frp = flp = PSET_GetNext( pdsp->ppsetprev[0], av ); } else { flp = pbf->left; frp = pbf->right; } fl = CLIP_DSP(fl); fr = CLIP_DSP(fr); flp = CLIP_DSP(flp); frp = CLIP_DSP(frp);
// get current ramp value
r = RMP_GetNext( &pdsp->xramp );
// crossfade from previous to current preset
if (!bexp) { xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE(fr, frp, r); // crossfade front left previous to front left
} else { xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE_EXP(fr, frp, r); // crossfade front left previous to front left
}
pbf->left = xf_fl; // crossfaded front left, duplicate in right channel
pbf->right = xf_fr; pbf++; }
return; }
// crossfade mono to mono preset
while ( count-- ) { av = ( ( pbf->left + pbf->right ) >> 1 );
// get current preset values
fl = PSET_GetNext( pdsp->ppset[0], av ); // get previous preset values
flp = PSET_GetNext( pdsp->ppsetprev[0], av ); fl = CLIP_DSP(fl); flp = CLIP_DSP(flp); // get current ramp value
r = RMP_GetNext( &pdsp->xramp );
// crossfade from previous to current preset
if (!bexp) xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
else xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
pbf->left = xf_fl; // crossfaded front left, duplicate in right channel
pbf->right = xf_fl; pbf++; } } }
// Helper: called only from DSP_Process
// DSP_Process stereo in to stereo out (if more than 2 procs, ignore them)
inline void DSP_ProcessStereoToStereo(dsp_t *pdsp, portable_samplepair_t *pbfront, portable_samplepair_t *pbrear, int sampleCount, bool bcrossfading ) { portable_samplepair_t *pbf = pbfront; // pointer to buffer of front stereo samples to process
int count = sampleCount; int fl, fr;
if ( !bcrossfading ) {
if ( !pdsp->ipset ) return;
if ( FBatchPreset(pdsp->ppset[0]) && FBatchPreset(pdsp->ppset[1]) ) {
// process left & right
PSET_GetNextN( pdsp->ppset[0], pbfront, sampleCount, OP_LEFT ); PSET_GetNextN( pdsp->ppset[1], pbfront, sampleCount, OP_RIGHT ); } else { // left -> left fx, right -> right fx
while ( count-- ) { fl = PSET_GetNext( pdsp->ppset[0], pbf->left ); fr = PSET_GetNext( pdsp->ppset[1], pbf->right );
fl = CLIP_DSP( fl ); fr = CLIP_DSP( fr );
pbf->left = fl; pbf->right = fr; pbf++; } } return; }
// crossfading to current preset from previous preset
if ( bcrossfading ) { int r; int flp, frp; int xf_fl, xf_fr; bool bexp = pdsp->bexpfade;
while ( count-- ) { // get current preset values
fl = PSET_GetNext( pdsp->ppset[0], pbf->left ); fr = PSET_GetNext( pdsp->ppset[1], pbf->right ); // get previous preset values
flp = PSET_GetNext( pdsp->ppsetprev[0], pbf->left ); frp = PSET_GetNext( pdsp->ppsetprev[1], pbf->right ); // get current ramp value
r = RMP_GetNext( &pdsp->xramp );
fl = CLIP_DSP( fl ); fr = CLIP_DSP( fr ); flp = CLIP_DSP( flp ); frp = CLIP_DSP( frp );
// crossfade from previous to current preset
if (!bexp) { xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE(fr, frp, r); } else { xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE_EXP(fr, frp, r); } pbf->left = xf_fl; // crossfaded front left
pbf->right = xf_fr; pbf++; } } }
// Helper: called only from DSP_Process
// DSP_Process quad in to mono out (front left = front right)
inline void DSP_ProcessQuadToMono(dsp_t *pdsp, portable_samplepair_t *pbfront, portable_samplepair_t *pbrear, int sampleCount, bool bcrossfading ) { portable_samplepair_t *pbf = pbfront; // pointer to buffer of front stereo samples to process
portable_samplepair_t *pbr = pbrear; // pointer to buffer of rear stereo samples to process
int count = sampleCount; int x; int av;
if ( !bcrossfading ) { if ( !pdsp->ipset ) return;
if ( FBatchPreset(pdsp->ppset[0]) ) {
// convert Quad to Mono in place, then batch process fx: perf KDB
// left front + rear -> left, right front + rear -> right
while ( count-- ) { pbf->left = ((pbf->left + pbf->right + pbr->left + pbr->right) >> 2); pbf++; pbr++; } // process left (mono), duplicate into right
PSET_GetNextN( pdsp->ppset[0], pbfront, sampleCount, OP_LEFT_DUPLICATE);
// copy processed front to rear
count = sampleCount; pbf = pbfront; pbr = pbrear;
while ( count-- ) { pbr->left = pbf->left; pbr->right = pbf->right; pbf++; pbr++; }
} else { // avg fl,fr,rl,rr into mono fx, duplicate on all channels
while ( count-- ) { av = ((pbf->left + pbf->right + pbr->left + pbr->right) >> 2); x = PSET_GetNext( pdsp->ppset[0], av ); x = CLIP_DSP( x ); pbr->left = pbr->right = pbf->left = pbf->right = x; pbf++; pbr++; } } return; }
if ( bcrossfading ) { int r; int fl, fr, rl, rr; int flp, frp, rlp, rrp; int xf_fl, xf_fr, xf_rl, xf_rr; bool bexp = pdsp->bexpfade; bool bfadetoquad = (pdsp->ipset == 0); bool bfadefromquad = (pdsp->ipsetprev == 0);
if ( bfadetoquad || bfadefromquad ) { // special case if previous or current preset is 0 (quad passthrough)
while ( count-- ) { av = ((pbf->left + pbf->right + pbr->left + pbr->right) >> 2);
// get current preset values
// current preset is 0, which implies fading to passthrough quad output
// need to fade from mono to quad
if ( pdsp->ipset ) { rl = rr = fl = fr = PSET_GetNext( pdsp->ppset[0], av ); } else { fl = pbf->left; fr = pbf->right; rl = pbr->left; rr = pbr->right; } // get previous preset values
if ( pdsp->ipsetprev ) { rrp = rlp = frp = flp = PSET_GetNext( pdsp->ppsetprev[0], av ); } else { flp = pbf->left; frp = pbf->right; rlp = pbr->left; rrp = pbr->right; } fl = CLIP_DSP(fl); fr = CLIP_DSP(fr); flp = CLIP_DSP(flp); frp = CLIP_DSP(frp); rl = CLIP_DSP(rl); rr = CLIP_DSP(rr); rlp = CLIP_DSP(rlp); rrp = CLIP_DSP(rrp);
// get current ramp value
r = RMP_GetNext( &pdsp->xramp );
// crossfade from previous to current preset
if (!bexp) { xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE(fr, frp, r); // crossfade front left previous to front left
xf_rl = XFADE(rl, rlp, r); // crossfade front left previous to front left
xf_rr = XFADE(rr, rrp, r); // crossfade front left previous to front left
} else { xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE_EXP(fr, frp, r); // crossfade front left previous to front left
xf_rl = XFADE_EXP(rl, rlp, r); // crossfade front left previous to front left
xf_rr = XFADE_EXP(rr, rrp, r); // crossfade front left previous to front left
}
pbf->left = xf_fl; pbf->right = xf_fr; pbr->left = xf_rl; pbr->right = xf_rr;
pbf++; pbr++; }
return; }
while ( count-- ) { av = ((pbf->left + pbf->right + pbr->left + pbr->right) >> 2);
// get current preset values
fl = PSET_GetNext( pdsp->ppset[0], av ); // get previous preset values
flp = PSET_GetNext( pdsp->ppsetprev[0], av ); // get current ramp value
r = RMP_GetNext( &pdsp->xramp );
fl = CLIP_DSP( fl ); flp = CLIP_DSP( flp );
// crossfade from previous to current preset
if (!bexp) xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
else xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
pbf->left = xf_fl; // crossfaded front left, duplicated to all channels
pbf->right = xf_fl; pbr->left = xf_fl; pbr->right = xf_fl;
pbf++; pbr++; } } }
// Helper: called only from DSP_Process
// DSP_Process quad in to stereo out (preserve stereo spatialization, throw away front/rear)
inline void DSP_ProcessQuadToStereo(dsp_t *pdsp, portable_samplepair_t *pbfront, portable_samplepair_t *pbrear, int sampleCount, bool bcrossfading ) { portable_samplepair_t *pbf = pbfront; // pointer to buffer of front stereo samples to process
portable_samplepair_t *pbr = pbrear; // pointer to buffer of rear stereo samples to process
int count = sampleCount; int fl, fr; if ( !bcrossfading ) { if ( !pdsp->ipset ) return;
if ( FBatchPreset(pdsp->ppset[0]) && FBatchPreset(pdsp->ppset[1]) ) {
// convert Quad to Stereo in place, then batch process fx: perf KDB
// left front + rear -> left, right front + rear -> right
while ( count-- ) { pbf->left = (pbf->left + pbr->left) >> 1; pbf->right = (pbf->right + pbr->right) >> 1; pbf++; pbr++; } // process left & right
PSET_GetNextN( pdsp->ppset[0], pbfront, sampleCount, OP_LEFT); PSET_GetNextN( pdsp->ppset[1], pbfront, sampleCount, OP_RIGHT );
// copy processed front to rear
count = sampleCount; pbf = pbfront; pbr = pbrear;
while ( count-- ) { pbr->left = pbf->left; pbr->right = pbf->right; pbf++; pbr++; }
} else { // left front + rear -> left fx, right front + rear -> right fx
while ( count-- ) { fl = PSET_GetNext( pdsp->ppset[0], (pbf->left + pbr->left) >> 1); fr = PSET_GetNext( pdsp->ppset[1], (pbf->right + pbr->right) >> 1); fl = CLIP_DSP( fl ); fr = CLIP_DSP( fr );
pbr->left = pbf->left = fl; pbr->right = pbf->right = fr; pbf++; pbr++; } } return; }
// crossfading to current preset from previous preset
if ( bcrossfading ) { int r; int rl, rr; int flp, frp, rlp, rrp; int xf_fl, xf_fr, xf_rl, xf_rr; int avl, avr; bool bexp = pdsp->bexpfade; bool bfadetoquad = (pdsp->ipset == 0); bool bfadefromquad = (pdsp->ipsetprev == 0);
if ( bfadetoquad || bfadefromquad ) { // special case if previous or current preset is 0 (quad passthrough)
while ( count-- ) { avl = (pbf->left + pbr->left) >> 1; avr = (pbf->right + pbr->right) >> 1;
// get current preset values
// current preset is 0, which implies fading to passthrough quad output
// need to fade from stereo to quad
if ( pdsp->ipset ) { rl = fl = PSET_GetNext( pdsp->ppset[0], avl ); rr = fr = PSET_GetNext( pdsp->ppset[0], avr ); } else { fl = pbf->left; fr = pbf->right; rl = pbr->left; rr = pbr->right; } // get previous preset values
if ( pdsp->ipsetprev ) { rlp = flp = PSET_GetNext( pdsp->ppsetprev[0], avl ); rrp = frp = PSET_GetNext( pdsp->ppsetprev[0], avr ); } else { flp = pbf->left; frp = pbf->right; rlp = pbr->left; rrp = pbr->right; } fl = CLIP_DSP(fl); fr = CLIP_DSP(fr); flp = CLIP_DSP(flp); frp = CLIP_DSP(frp); rl = CLIP_DSP(rl); rr = CLIP_DSP(rr); rlp = CLIP_DSP(rlp); rrp = CLIP_DSP(rrp);
// get current ramp value
r = RMP_GetNext( &pdsp->xramp );
// crossfade from previous to current preset
if (!bexp) { xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE(fr, frp, r); // crossfade front left previous to front left
xf_rl = XFADE(rl, rlp, r); // crossfade front left previous to front left
xf_rr = XFADE(rr, rrp, r); // crossfade front left previous to front left
} else { xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE_EXP(fr, frp, r); // crossfade front left previous to front left
xf_rl = XFADE_EXP(rl, rlp, r); // crossfade front left previous to front left
xf_rr = XFADE_EXP(rr, rrp, r); // crossfade front left previous to front left
}
pbf->left = xf_fl; pbf->right = xf_fr; pbr->left = xf_rl; pbr->right = xf_rr;
pbf++; pbr++; }
return; }
while ( count-- ) { avl = (pbf->left + pbr->left) >> 1; avr = (pbf->right + pbr->right) >> 1;
// get current preset values
fl = PSET_GetNext( pdsp->ppset[0], avl ); fr = PSET_GetNext( pdsp->ppset[1], avr ); // get previous preset values
flp = PSET_GetNext( pdsp->ppsetprev[0], avl ); frp = PSET_GetNext( pdsp->ppsetprev[1], avr );
fl = CLIP_DSP( fl ); fr = CLIP_DSP( fr ); // get previous preset values
flp = CLIP_DSP( flp ); frp = CLIP_DSP( frp );
// get current ramp value
r = RMP_GetNext( &pdsp->xramp );
// crossfade from previous to current preset
if (!bexp) { xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE(fr, frp, r); } else { xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE_EXP(fr, frp, r); } pbf->left = xf_fl; // crossfaded front left
pbf->right = xf_fr;
pbr->left = xf_fl; // duplicate front channel to rear channel
pbr->right = xf_fr;
pbf++; pbr++; } } }
// Helper: called only from DSP_Process
// DSP_Process quad in to quad out
inline void DSP_ProcessQuadToQuad(dsp_t *pdsp, portable_samplepair_t *pbfront, portable_samplepair_t *pbrear, int sampleCount, bool bcrossfading ) { portable_samplepair_t *pbf = pbfront; // pointer to buffer of front stereo samples to process
portable_samplepair_t *pbr = pbrear; // pointer to buffer of rear stereo samples to process
int count = sampleCount; int fl, fr, rl, rr;
if ( !bcrossfading ) { if ( !pdsp->ipset ) return;
// each channel gets its own processor
if ( FBatchPreset(pdsp->ppset[0]) && FBatchPreset(pdsp->ppset[1]) && FBatchPreset(pdsp->ppset[2]) && FBatchPreset(pdsp->ppset[3])) { // batch process fx front & rear, left & right: perf KDB
PSET_GetNextN( pdsp->ppset[0], pbfront, sampleCount, OP_LEFT); PSET_GetNextN( pdsp->ppset[1], pbfront, sampleCount, OP_RIGHT ); PSET_GetNextN( pdsp->ppset[2], pbrear, sampleCount, OP_LEFT ); PSET_GetNextN( pdsp->ppset[3], pbrear, sampleCount, OP_RIGHT ); } else { while ( count-- ) { fl = PSET_GetNext( pdsp->ppset[0], pbf->left ); fr = PSET_GetNext( pdsp->ppset[1], pbf->right ); rl = PSET_GetNext( pdsp->ppset[2], pbr->left ); rr = PSET_GetNext( pdsp->ppset[3], pbr->right ); pbf->left = CLIP_DSP( fl ); pbf->right = CLIP_DSP( fr ); pbr->left = CLIP_DSP( rl ); pbr->right = CLIP_DSP( rr );
pbf++; pbr++; } } return; }
// crossfading to current preset from previous preset
if ( bcrossfading ) { int r; int flp, frp, rlp, rrp; int xf_fl, xf_fr, xf_rl, xf_rr; bool bexp = pdsp->bexpfade;
while ( count-- ) { // get current preset values
fl = PSET_GetNext( pdsp->ppset[0], pbf->left ); fr = PSET_GetNext( pdsp->ppset[1], pbf->right ); rl = PSET_GetNext( pdsp->ppset[2], pbr->left ); rr = PSET_GetNext( pdsp->ppset[3], pbr->right );
// get previous preset values
flp = PSET_GetNext( pdsp->ppsetprev[0], pbf->left ); frp = PSET_GetNext( pdsp->ppsetprev[1], pbf->right ); rlp = PSET_GetNext( pdsp->ppsetprev[2], pbr->left ); rrp = PSET_GetNext( pdsp->ppsetprev[3], pbr->right );
// get current ramp value
r = RMP_GetNext( &pdsp->xramp );
// crossfade from previous to current preset
if (!bexp) { xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE(fr, frp, r); xf_rl = XFADE(rl, rlp, r); xf_rr = XFADE(rr, rrp, r); } else { xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE_EXP(fr, frp, r); xf_rl = XFADE_EXP(rl, rlp, r); xf_rr = XFADE_EXP(rr, rrp, r); } pbf->left = CLIP_DSP(xf_fl); // crossfaded front left
pbf->right = CLIP_DSP(xf_fr); pbr->left = CLIP_DSP(xf_rl); pbr->right = CLIP_DSP(xf_rr);
pbf++; pbr++; } } }
// Helper: called only from DSP_Process
// DSP_Process quad + center in to mono out (front left = front right)
inline void DSP_Process5To1(dsp_t *pdsp, portable_samplepair_t *pbfront, portable_samplepair_t *pbrear, portable_samplepair_t *pbcenter, int sampleCount, bool bcrossfading ) { portable_samplepair_t *pbf = pbfront; // pointer to buffer of front stereo samples to process
portable_samplepair_t *pbr = pbrear; // pointer to buffer of rear stereo samples to process
portable_samplepair_t *pbc = pbcenter; // pointer to buffer of center mono samples to process
int count = sampleCount; int x; int av;
if ( !bcrossfading ) { if ( !pdsp->ipset ) return;
if ( FBatchPreset(pdsp->ppset[0]) ) {
// convert Quad + Center to Mono in place, then batch process fx: perf KDB
// left front + rear -> left, right front + rear -> right
while ( count-- ) { // pbf->left = ((pbf->left + pbf->right + pbr->left + pbr->right + pbc->left) / 5);
av = (pbf->left + pbf->right + pbr->left + pbr->right + pbc->left) * 51; // 51/255 = 1/5
av >>= 8; pbf->left = av; pbf++; pbr++; pbc++; } // process left (mono), duplicate into right
PSET_GetNextN( pdsp->ppset[0], pbfront, sampleCount, OP_LEFT_DUPLICATE);
// copy processed front to rear & center
count = sampleCount; pbf = pbfront; pbr = pbrear; pbc = pbcenter;
while ( count-- ) { pbr->left = pbf->left; pbr->right = pbf->right; pbc->left = pbf->left; pbf++; pbr++; pbc++; }
} else { // avg fl,fr,rl,rr,fc into mono fx, duplicate on all channels
while ( count-- ) { // av = ((pbf->left + pbf->right + pbr->left + pbr->right + pbc->left) / 5);
av = (pbf->left + pbf->right + pbr->left + pbr->right + pbc->left) * 51; // 51/255 = 1/5
av >>= 8; x = PSET_GetNext( pdsp->ppset[0], av ); x = CLIP_DSP( x ); pbr->left = pbr->right = pbf->left = pbf->right = pbc->left = x; pbf++; pbr++; pbc++; } } return; }
if ( bcrossfading ) { int r; int fl, fr, rl, rr, fc; int flp, frp, rlp, rrp, fcp; int xf_fl, xf_fr, xf_rl, xf_rr, xf_fc; bool bexp = pdsp->bexpfade; bool bfadetoquad = (pdsp->ipset == 0); bool bfadefromquad = (pdsp->ipsetprev == 0);
if ( bfadetoquad || bfadefromquad ) { // special case if previous or current preset is 0 (quad passthrough)
while ( count-- ) { // av = ((pbf->left + pbf->right + pbr->left + pbr->right) >> 2);
av = (pbf->left + pbf->right + pbr->left + pbr->right + pbc->left) * 51; // 51/255 = 1/5
av >>= 8;
// get current preset values
// current preset is 0, which implies fading to passthrough quad output
// need to fade from mono to quad
if ( pdsp->ipset ) { fc = rl = rr = fl = fr = PSET_GetNext( pdsp->ppset[0], av ); } else { fl = pbf->left; fr = pbf->right; rl = pbr->left; rr = pbr->right; fc = pbc->left; } // get previous preset values
if ( pdsp->ipsetprev ) { fcp = rrp = rlp = frp = flp = PSET_GetNext( pdsp->ppsetprev[0], av ); } else { flp = pbf->left; frp = pbf->right; rlp = pbr->left; rrp = pbr->right; fcp = pbc->left; } fl = CLIP_DSP(fl); fr = CLIP_DSP(fr); flp = CLIP_DSP(flp); frp = CLIP_DSP(frp); rl = CLIP_DSP(rl); rr = CLIP_DSP(rr); rlp = CLIP_DSP(rlp); rrp = CLIP_DSP(rrp); fc = CLIP_DSP(fc); fcp = CLIP_DSP(fcp);
// get current ramp value
r = RMP_GetNext( &pdsp->xramp );
// crossfade from previous to current preset
if (!bexp) { xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE(fr, frp, r); // crossfade front left previous to front left
xf_rl = XFADE(rl, rlp, r); // crossfade front left previous to front left
xf_rr = XFADE(rr, rrp, r); // crossfade front left previous to front left
xf_fc = XFADE(fc, fcp, r); // crossfade front left previous to front left
} else { xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE_EXP(fr, frp, r); // crossfade front left previous to front left
xf_rl = XFADE_EXP(rl, rlp, r); // crossfade front left previous to front left
xf_rr = XFADE_EXP(rr, rrp, r); // crossfade front left previous to front left
xf_fc = XFADE_EXP(fc, fcp, r); // crossfade front left previous to front left
}
pbf->left = xf_fl; pbf->right = xf_fr; pbr->left = xf_rl; pbr->right = xf_rr; pbc->left = xf_fc;
pbf++; pbr++; pbc++; }
return; }
while ( count-- ) { // av = ((pbf->left + pbf->right + pbr->left + pbr->right) >> 2);
av = (pbf->left + pbf->right + pbr->left + pbr->right + pbc->left) * 51; // 51/255 = 1/5
av >>= 8;
// get current preset values
fl = PSET_GetNext( pdsp->ppset[0], av ); // get previous preset values
flp = PSET_GetNext( pdsp->ppsetprev[0], av ); // get current ramp value
r = RMP_GetNext( &pdsp->xramp );
fl = CLIP_DSP( fl ); flp = CLIP_DSP( flp );
// crossfade from previous to current preset
if (!bexp) xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
else xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
pbf->left = xf_fl; // crossfaded front left, duplicated to all channels
pbf->right = xf_fl; pbr->left = xf_fl; pbr->right = xf_fl; pbc->left = xf_fl;
pbf++; pbr++; pbc++; } } }
// Helper: called only from DSP_Process
// DSP_Process quad + center in to quad + center out
inline void DSP_Process5To5(dsp_t *pdsp, portable_samplepair_t *pbfront, portable_samplepair_t *pbrear, portable_samplepair_t *pbcenter, int sampleCount, bool bcrossfading ) { portable_samplepair_t *pbf = pbfront; // pointer to buffer of front stereo samples to process
portable_samplepair_t *pbr = pbrear; // pointer to buffer of rear stereo samples to process
portable_samplepair_t *pbc = pbcenter; // pointer to buffer of center mono samples to process
int count = sampleCount; int fl, fr, rl, rr, fc;
if ( !bcrossfading ) { if ( !pdsp->ipset ) return;
// each channel gets its own processor
if ( FBatchPreset(pdsp->ppset[0]) && FBatchPreset(pdsp->ppset[1]) && FBatchPreset(pdsp->ppset[2]) && FBatchPreset(pdsp->ppset[3])) { // batch process fx front & rear, left & right: perf KDB
PSET_GetNextN( pdsp->ppset[0], pbfront, sampleCount, OP_LEFT); PSET_GetNextN( pdsp->ppset[1], pbfront, sampleCount, OP_RIGHT ); PSET_GetNextN( pdsp->ppset[2], pbrear, sampleCount, OP_LEFT ); PSET_GetNextN( pdsp->ppset[3], pbrear, sampleCount, OP_RIGHT ); PSET_GetNextN( pdsp->ppset[4], pbcenter, sampleCount, OP_LEFT ); } else { while ( count-- ) { fl = PSET_GetNext( pdsp->ppset[0], pbf->left ); fr = PSET_GetNext( pdsp->ppset[1], pbf->right ); rl = PSET_GetNext( pdsp->ppset[2], pbr->left ); rr = PSET_GetNext( pdsp->ppset[3], pbr->right ); fc = PSET_GetNext( pdsp->ppset[4], pbc->left ); pbf->left = CLIP_DSP( fl ); pbf->right = CLIP_DSP( fr ); pbr->left = CLIP_DSP( rl ); pbr->right = CLIP_DSP( rr ); pbc->left = CLIP_DSP( fc );
pbf++; pbr++; pbc++; } } return; }
// crossfading to current preset from previous preset
if ( bcrossfading ) { int r; int flp, frp, rlp, rrp, fcp; int xf_fl, xf_fr, xf_rl, xf_rr, xf_fc; bool bexp = pdsp->bexpfade;
while ( count-- ) { // get current preset values
fl = PSET_GetNext( pdsp->ppset[0], pbf->left ); fr = PSET_GetNext( pdsp->ppset[1], pbf->right ); rl = PSET_GetNext( pdsp->ppset[2], pbr->left ); rr = PSET_GetNext( pdsp->ppset[3], pbr->right ); fc = PSET_GetNext( pdsp->ppset[4], pbc->left );
// get previous preset values
flp = PSET_GetNext( pdsp->ppsetprev[0], pbf->left ); frp = PSET_GetNext( pdsp->ppsetprev[1], pbf->right ); rlp = PSET_GetNext( pdsp->ppsetprev[2], pbr->left ); rrp = PSET_GetNext( pdsp->ppsetprev[3], pbr->right ); fcp = PSET_GetNext( pdsp->ppsetprev[4], pbc->left );
// get current ramp value
r = RMP_GetNext( &pdsp->xramp );
// crossfade from previous to current preset
if (!bexp) { xf_fl = XFADE(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE(fr, frp, r); xf_rl = XFADE(rl, rlp, r); xf_rr = XFADE(rr, rrp, r); xf_fc = XFADE(fc, fcp, r); } else { xf_fl = XFADE_EXP(fl, flp, r); // crossfade front left previous to front left
xf_fr = XFADE_EXP(fr, frp, r); xf_rl = XFADE_EXP(rl, rlp, r); xf_rr = XFADE_EXP(rr, rrp, r); xf_fc = XFADE_EXP(fc, fcp, r); } pbf->left = CLIP_DSP(xf_fl); // crossfaded front left
pbf->right = CLIP_DSP(xf_fr); pbr->left = CLIP_DSP(xf_rl); pbr->right = CLIP_DSP(xf_rr); pbc->left = CLIP_DSP(xf_fc);
pbf++; pbr++; pbc++; } } }
// This is an evil hack, but we need to restore the old presets after letting the sound system update for a few frames, so we just
// "defer" the restore until the top of the next call to CheckNewDspPresets. I put in a bit of warning in case we ever have code
// outside of this time period modifying any of the dsp convars. It doesn't seem to be an issue just save/loading between levels
static bool g_bNeedPresetRestore = false;
//-----------------------------------------------------------------------------
// Purpose:
//-----------------------------------------------------------------------------
struct PreserveDSP_t { ConVar *cvar; float oldvalue; };
static PreserveDSP_t g_PreserveDSP[] = { { &dsp_room }, { &dsp_water }, { &dsp_player }, { &dsp_facingaway }, { &dsp_speaker }, { &dsp_spatial }, { &dsp_automatic } };
//-----------------------------------------------------------------------------
// Purpose: Called at the top of CheckNewDspPresets to restore ConVars to real values
//-----------------------------------------------------------------------------
void DSP_CheckRestorePresets() { if ( !g_bNeedPresetRestore ) return;
g_bNeedPresetRestore = false;
int i; int c = ARRAYSIZE( g_PreserveDSP );
// Restore
for ( i = 0 ; i < c; ++i ) { PreserveDSP_t& slot = g_PreserveDSP[ i ];
ConVar *cv = slot.cvar; Assert( cv ); if ( cv->GetFloat() != 0.0f ) { // NOTE: dsp_speaker is being (correctly) save/restored by maps, which would trigger this warning
//Warning( "DSP_CheckRestorePresets: Value of %s was changed between DSP_ClearState and CheckNewDspPresets, not restoring to old value\n", cv->GetName() );
continue; } cv->SetValue( slot.oldvalue ); }
// reinit all dsp processors (only load preset file on engine init, however)
AllocDsps( false );
// flush dsp automatic nodes
g_bdas_init_nodes = 0; g_bdas_room_init = 0; }
//-----------------------------------------------------------------------------
// Purpose:
//-----------------------------------------------------------------------------
void DSP_ClearState() { // if we already cleared dsp state, and a restore is pending,
// don't clear again
if ( g_bNeedPresetRestore ) return;
// always save a cleared dsp automatic value to force reset of all adsp code
dsp_automatic.SetValue(0);
// Tracker 7155: YWB: This is a pretty ugly hack to zero out all of the dsp convars and bootstrap the dsp system into using them for a few frames
int i; int c = ARRAYSIZE( g_PreserveDSP );
for ( i = 0 ; i < c; ++i ) { PreserveDSP_t& slot = g_PreserveDSP[ i ];
ConVar *cv = slot.cvar; Assert( cv ); slot.oldvalue = cv->GetFloat(); cv->SetValue( 0 ); } // force all dsp presets to end crossfades, end one-shot presets, & release and reset all resources
// immediately.
FreeDsps( false ); // free all dsp states, but don't discard preset templates
// This forces the ConVars which we set to zero above to be reloaded to their old values at the time we issue the CheckNewDspPresets
// command. This seems to happen early enough in level changes were we don't appear to be trying to stomp real settings...
g_bNeedPresetRestore = true; }
// return true if dsp's preset is one-shot and it has expired
bool DSP_HasExpired( int idsp ) { dsp_t *pdsp;
Assert( idsp < CDSPS );
if (idsp < 0 || idsp >= CDSPS) return false;
pdsp = &dsps[idsp];
// if first preset has expired, dsp has expired
if ( PSET_IsOneShot( pdsp->ppset[0] ) ) return PSET_HasExpired( pdsp->ppset[0] ); else return false; }
// returns true if dsp is crossfading from previous dsp preset
bool DSP_IsCrossfading( int idsp ) { dsp_t *pdsp;
Assert( idsp < CDSPS );
if (idsp < 0 || idsp >= CDSPS) return false;
pdsp = &dsps[idsp];
return !RMP_HitEnd( &pdsp->xramp );
}
// returns previous preset # before oneshot preset was set
int DSP_OneShotPrevious( int idsp ) { dsp_t *pdsp; int idsp_prev;
Assert( idsp < CDSPS );
if (idsp < 0 || idsp >= CDSPS) return 0;
pdsp = &dsps[idsp]; idsp_prev = pdsp->ipsetsav_oneshot;
return idsp_prev; }
// given idsp (processor index), return true if
// both current and previous presets are 0 for this processor
bool DSP_PresetIsOff( int idsp ) { dsp_t *pdsp;
if (idsp < 0 || idsp >= CDSPS) return true;
Assert ( idsp < CDSPS ); // make sure idsp is valid
pdsp = &dsps[idsp];
// if current and previous preset 0, return - preset 0 is 'off'
return ( !pdsp->ipset && !pdsp->ipsetprev ); }
// returns true if dsp is off for room effects
bool DSP_RoomDSPIsOff() { return DSP_PresetIsOff( Get_idsp_room() ); }
// Main DSP processing routine:
// process samples in buffers using pdsp processor
// continue crossfade between 2 dsp processors if crossfading on switch
// pfront - front stereo buffer to process
// prear - rear stereo buffer to process (may be NULL)
// pcenter - front center mono buffer (may be NULL)
// sampleCount - number of samples in pbuf to process
// This routine also maps the # processing channels in the pdsp to the number of channels
// supplied. ie: if the pdsp has 4 channels and pbfront and pbrear are both non-null, the channels
// map 1:1 through the processors.
void DSP_Process( int idsp, portable_samplepair_t *pbfront, portable_samplepair_t *pbrear, portable_samplepair_t *pbcenter, int sampleCount ) { bool bcrossfading; int cchan_in; // input channels (2,4 or 5)
int cprocs; // output cannels (1, 2 or 4)
dsp_t *pdsp;
if (idsp < 0 || idsp >= CDSPS) return;
// Don't pull dsp data in if player is not connected (during load/level change)
if ( !g_pSoundServices->IsConnected() ) return;
Assert ( idsp < CDSPS ); // make sure idsp is valid
pdsp = &dsps[idsp];
Assert (pbfront);
// return right away if fx processing is turned off
if ( dsp_off.GetInt() ) return;
// if current and previous preset 0, return - preset 0 is 'off'
if ( !pdsp->ipset && !pdsp->ipsetprev ) return;
if ( sampleCount < 0 ) return;
bcrossfading = !RMP_HitEnd( &pdsp->xramp );
// if not crossfading, and previous channel is not null, free previous
if ( !bcrossfading ) DSP_FreePrevPreset( pdsp );
// if current and previous preset 0 (ie: just freed previous), return - preset 0 is 'off'
if ( !pdsp->ipset && !pdsp->ipsetprev ) return;
cchan_in = (pbrear ? 4 : 2) + (pbcenter ? 1 : 0); cprocs = pdsp->cchan;
Assert(cchan_in == 2 || cchan_in == 4 || cchan_in == 5 );
// if oneshot preset, update the duration counter (only update front left counter)
PSET_UpdateDuration( pdsp->ppset[0], sampleCount );
// NOTE: when mixing between different channel sizes,
// always AVERAGE down to fewer channels and DUPLICATE up more channels.
// The following routines always process cchan_in channels.
// ie: QuadToMono still updates 4 values in buffer
// DSP_Process stereo in to mono out (ie: left and right are averaged)
if ( cchan_in == 2 && cprocs == 1) { DSP_ProcessStereoToMono( pdsp, pbfront, pbrear, sampleCount, bcrossfading ); return; }
// DSP_Process stereo in to stereo out (if more than 2 procs, ignore them)
if ( cchan_in == 2 && cprocs >= 2) { DSP_ProcessStereoToStereo( pdsp, pbfront, pbrear, sampleCount, bcrossfading ); return; }
// DSP_Process quad in to mono out
if ( cchan_in == 4 && cprocs == 1) { DSP_ProcessQuadToMono( pdsp, pbfront, pbrear, sampleCount, bcrossfading ); return; }
// DSP_Process quad in to stereo out (preserve stereo spatialization, loose front/rear)
if ( cchan_in == 4 && cprocs == 2) { DSP_ProcessQuadToStereo( pdsp, pbfront, pbrear, sampleCount, bcrossfading ); return; }
// DSP_Process quad in to quad out
if ( cchan_in == 4 && cprocs == 4) { DSP_ProcessQuadToQuad( pdsp, pbfront, pbrear, sampleCount, bcrossfading ); return; }
// DSP_Process quad + center in to mono out
if ( cchan_in == 5 && cprocs == 1) { DSP_Process5To1( pdsp, pbfront, pbrear, pbcenter, sampleCount, bcrossfading ); return; }
if ( cchan_in == 5 && cprocs == 2) { // undone: not used in AllocDsps
Assert(false); //DSP_Process5to2( pdsp, pbfront, pbrear, pbcenter, sampleCount, bcrossfading );
return; }
if ( cchan_in == 5 && cprocs == 4) { // undone: not used in AllocDsps
Assert(false); //DSP_Process5to4( pdsp, pbfront, pbrear, pbcenter, sampleCount, bcrossfading );
return; }
// DSP_Process quad + center in to quad + center out
if ( cchan_in == 5 && cprocs == 5) { DSP_Process5To5( pdsp, pbfront, pbrear, pbcenter, sampleCount, bcrossfading ); return; }
}
// DSP helpers
// free all dsp processors
void FreeDsps( bool bReleaseTemplateMemory ) {
DSP_Free(idsp_room); DSP_Free(idsp_water); DSP_Free(idsp_player); DSP_Free(idsp_facingaway); DSP_Free(idsp_speaker); DSP_Free(idsp_spatial); DSP_Free(idsp_automatic);
idsp_room = 0; idsp_water = 0; idsp_player = 0; idsp_facingaway = 0; idsp_speaker = 0; idsp_spatial = 0; idsp_automatic = 0;
for ( int i = SOUND_BUFFER_SPECIAL_START; i < g_paintBuffers.Count(); ++i ) { paintbuffer_t *pSpecialBuffer = MIX_GetPPaintFromIPaint( i ); if ( pSpecialBuffer->nSpecialDSP != 0 ) { DSP_Free( pSpecialBuffer->idsp_specialdsp ); pSpecialBuffer->idsp_specialdsp = 0; pSpecialBuffer->nPrevSpecialDSP = 0; pSpecialBuffer->nSpecialDSP = 0; } }
DSP_FreeAll(); // only unlock and free psettemplate memory on engine shutdown
if ( bReleaseTemplateMemory ) DSP_ReleaseMemory(); }
// alloc dsp processors, load dsp preset array from file on engine init only
bool AllocDsps( bool bLoadPresetFile ) { int csurround = (g_AudioDevice->IsSurround() ? 2: 0); // surround channels to allocate
int ccenter = (g_AudioDevice->IsSurroundCenter() ? 1 : 0); // center channels to allocate
DSP_InitAll( bLoadPresetFile );
idsp_room = -1; idsp_water = -1; idsp_player = -1; idsp_facingaway = -1; idsp_speaker = -1; idsp_spatial = -1; idsp_automatic = -1;
// alloc dsp room channel (mono, stereo if dsp_stereo is 1)
// dsp room is mono, 300ms default fade time
idsp_room = DSP_Alloc( dsp_room.GetInt(), 200, 1 );
// dsp automatic overrides dsp_room, if dsp_room set to DSP_AUTOMATIC (1)
idsp_automatic = DSP_Alloc( dsp_automatic.GetInt(), 200, 1 ) ;
// alloc stereo or quad series processors for player or water
// water and player presets are mono
idsp_water = DSP_Alloc( dsp_water.GetInt(), 100, 1 ); idsp_player = DSP_Alloc( dsp_player.GetInt(), 100, 1 );
// alloc facing away filters (stereo, quad or 5ch)
idsp_facingaway = DSP_Alloc( dsp_facingaway.GetInt(), 100, 2 + csurround + ccenter );
// alloc speaker preset (mono)
idsp_speaker = DSP_Alloc( dsp_speaker.GetInt(), 300, 1 );
// alloc spatial preset (2-5 chan)
idsp_spatial = DSP_Alloc( dsp_spatial.GetInt(), 300, 2 + csurround + ccenter );
// init prev values
ipset_room_prev = dsp_room.GetInt(); ipset_water_prev = dsp_water.GetInt(); ipset_player_prev = dsp_player.GetInt(); ipset_facingaway_prev = dsp_facingaway.GetInt(); ipset_room_typeprev = dsp_room_type.GetInt(); ipset_speaker_prev = dsp_speaker.GetInt(); ipset_spatial_prev = dsp_spatial.GetInt(); ipset_automatic_prev = dsp_automatic.GetInt();
if (idsp_room < 0 || idsp_water < 0 || idsp_player < 0 || idsp_facingaway < 0 || idsp_speaker < 0 || idsp_spatial < 0 || idsp_automatic < 0) { DevMsg ("WARNING: DSP processor failed to initialize! \n" );
FreeDsps( true ); return false; } return true; }
// count number of dsp presets specified in preset file
// counts outer {} pairs, ignoring inner {} pairs.
int DSP_CountFilePresets( const char *pstart ) { int cpresets = 0; bool binpreset = false; bool blookleft = false;
while ( 1 ) { pstart = COM_Parse( pstart ); if ( strlen(com_token) <= 0) break;
if ( com_token[0] == '{' ) // left paren
{ if (!binpreset) { cpresets++; // found preset:
blookleft = true; // look for another left
binpreset = true; } else { blookleft = false; // inside preset: next, look for matching right paren
}
continue; }
if ( com_token[0] == '}' ) // right paren
{ if (binpreset) { if (!blookleft) // looking for right paren
{ blookleft = true; // found it, now look for another left
} else { // expected inner left paren, found outer right - end of preset definition
binpreset = false; blookleft = true; } } else { // error - unexpected } paren
DevMsg("PARSE ERROR!!! dsp_presets.txt: unexpected '}' \n"); continue; } }
}
return cpresets; }
struct dsp_stringmap_t { char sz[33]; int i; };
// token map for dsp_preset.txt
dsp_stringmap_t gdsp_stringmap[] = { // PROCESSOR TYPE:
{"NULL", PRC_NULL}, {"DLY", PRC_DLY}, {"RVA", PRC_RVA}, {"FLT", PRC_FLT}, {"CRS", PRC_CRS}, {"PTC", PRC_PTC}, {"ENV", PRC_ENV}, {"LFO", PRC_LFO}, {"EFO", PRC_EFO}, {"MDY", PRC_MDY}, {"DFR", PRC_DFR}, {"AMP", PRC_AMP},
// FILTER TYPE:
{"LP", FLT_LP}, {"HP", FLT_HP}, {"BP", FLT_BP},
// FILTER QUALITY:
{"LO", QUA_LO}, {"MED", QUA_MED}, {"HI", QUA_HI}, {"VHI", QUA_VHI},
// DELAY TYPE:
{"PLAIN", DLY_PLAIN}, {"ALLPASS", DLY_ALLPASS}, {"LOWPASS", DLY_LOWPASS}, {"DLINEAR", DLY_LINEAR}, {"FLINEAR", DLY_FLINEAR}, {"LOWPASS_4TAP",DLY_LOWPASS_4TAP}, {"PLAIN_4TAP", DLY_PLAIN_4TAP},
// LFO TYPE:
{"SIN", LFO_SIN}, {"TRI", LFO_TRI}, {"SQR", LFO_SQR}, {"SAW", LFO_SAW}, {"RND", LFO_RND}, {"LOG_IN", LFO_LOG_IN}, {"LOG_OUT", LFO_LOG_OUT}, {"LIN_IN", LFO_LIN_IN}, {"LIN_OUT", LFO_LIN_OUT},
// ENVELOPE TYPE:
{"LIN", ENV_LIN}, {"EXP", ENV_EXP},
// PRESET CONFIGURATION TYPE:
{"SIMPLE", PSET_SIMPLE}, {"LINEAR", PSET_LINEAR}, {"PARALLEL2", PSET_PARALLEL2}, {"PARALLEL4", PSET_PARALLEL4}, {"PARALLEL5", PSET_PARALLEL5}, {"FEEDBACK", PSET_FEEDBACK}, {"FEEDBACK3", PSET_FEEDBACK3}, {"FEEDBACK4", PSET_FEEDBACK4}, {"MOD1", PSET_MOD}, {"MOD2", PSET_MOD2}, {"MOD3", PSET_MOD3} };
int gcdsp_stringmap = sizeof(gdsp_stringmap) / sizeof (dsp_stringmap_t);
#define isnumber(c) (c == '+' || c == '-' || c == '0' || c == '1' || c == '2' || c == '3' || c == '4' || c == '5' || c == '6' || c == '7'|| c == '8' || c == '9')\
// given ptr to null term. string, return integer or float value from g_dsp_stringmap
float DSP_LookupStringToken( char *psz, int ipset ) { int i; float fipset = (float)ipset;
if (isnumber(psz[0])) return atof(psz);
for (i = 0; i < gcdsp_stringmap; i++) { if (!strcmpi(gdsp_stringmap[i].sz, psz)) return gdsp_stringmap[i].i; }
// not found
DevMsg("DSP PARSE ERROR! token not found in dsp_presets.txt. Preset: %3.0f \n", fipset ); return 0; }
// load dsp preset file, parse presets into g_psettemplate array
// format for each preset:
// { <preset #> <preset type> <#processors> <gain> { <processor type> <param0>...<param15> } {...} {...} }
#define CHAR_LEFT_PAREN '{'
#define CHAR_RIGHT_PAREN '}'
// free preset template memory
void DSP_ReleaseMemory( void ) { if (g_psettemplates) { delete[] g_psettemplates; g_psettemplates = NULL; } }
bool DSP_LoadPresetFile( void ) { char szFile[ MAX_OSPATH ]; char *pbuffer; const char *pstart; bool bResult = false; int cpresets; int ipreset; int itype; int cproc; float mix_min; float mix_max; float db_min; float db_mixdrop; int j; bool fdone; float duration; float fadeout;
Q_snprintf( szFile, sizeof( szFile ), "scripts/dsp_presets.txt" );
MEM_ALLOC_CREDIT();
CUtlBuffer buf;
if ( !g_pFullFileSystem->ReadFile( szFile, "GAME", buf ) ) { Error( "DSP_LoadPresetFile: unable to open '%s'\n", szFile ); return bResult; } pbuffer = (char *)buf.PeekGet(); // Use malloc - free at end of this routine
pstart = pbuffer;
// figure out how many presets we're loading - count outer parens.
cpresets = DSP_CountFilePresets( pstart ); g_cpsettemplates = cpresets;
g_psettemplates = new pset_t[cpresets]; if (!g_psettemplates) { Warning( "DSP Preset Loader: Out of memory.\n"); goto load_exit; } memset (g_psettemplates, 0, cpresets * sizeof(pset_t));
// parse presets into g_psettemplates array
pstart = pbuffer;
// for each preset...
for ( j = 0; j < cpresets; j++) { // check for end of file or next CHAR_LEFT_PAREN
while (1) { pstart = COM_Parse( pstart ); if ( strlen(com_token) <= 0) break;
if ( com_token[0] != CHAR_LEFT_PAREN ) continue;
break; }
// found start of a new preset definition
// get preset #, type, cprocessors, gain
pstart = COM_Parse( pstart ); ipreset = atoi( com_token );
pstart = COM_Parse( pstart ); itype = (int)DSP_LookupStringToken( com_token , ipreset);
pstart = COM_Parse( pstart ); mix_min = atof( com_token );
pstart = COM_Parse( pstart ); mix_max = atof( com_token );
pstart = COM_Parse( pstart ); duration = atof( com_token );
pstart = COM_Parse( pstart ); fadeout = atof( com_token ); pstart = COM_Parse( pstart ); db_min = atof( com_token );
pstart = COM_Parse( pstart ); db_mixdrop = atof( com_token );
g_psettemplates[ipreset].fused = true; g_psettemplates[ipreset].mix_min = mix_min; g_psettemplates[ipreset].mix_max = mix_max; g_psettemplates[ipreset].duration = duration; g_psettemplates[ipreset].fade = fadeout; g_psettemplates[ipreset].db_min = db_min; g_psettemplates[ipreset].db_mixdrop = db_mixdrop; // parse each processor for this preset
fdone = false; cproc = 0;
while (1) { // find CHAR_LEFT_PAREN - start of new processor
while (1) { pstart = COM_Parse( pstart );
if ( strlen(com_token) <= 0) break;
if (com_token[0] == CHAR_LEFT_PAREN) break;
if (com_token[0] == CHAR_RIGHT_PAREN) { // if found right paren, no more processors: done with this preset
fdone = true; break; } } if ( fdone ) break;
// get processor type
pstart = COM_Parse( pstart ); g_psettemplates[ipreset].prcs[cproc].type = (int)DSP_LookupStringToken( com_token, ipreset );
// get param 0..n or stop when hit closing CHAR_RIGHT_PAREN
int ip = 0;
while (1) { pstart = COM_Parse( pstart );
if ( strlen(com_token) <= 0) break;
if ( com_token[0] == CHAR_RIGHT_PAREN ) break;
g_psettemplates[ipreset].prcs[cproc].prm[ip++] = DSP_LookupStringToken( com_token, ipreset ); // cap at max params
ip = min(ip, CPRCPARAMS); } cproc++; if (cproc > CPSET_PRCS) DevMsg("DSP PARSE ERROR!!! dsp_presets.txt: missing } or too many processors in preset #: %d \n", ipreset); cproc = min(cproc, CPSET_PRCS); // don't overflow # procs
} // if cproc == 1, type is always SIMPLE
if ( cproc == 1) itype = PSET_SIMPLE;
g_psettemplates[ipreset].type = itype; g_psettemplates[ipreset].cprcs = cproc; }
bResult = true;
load_exit: return bResult; }
//-----------------------------------------------------------------------------
// Purpose: Called by client on level shutdown to clear ear ringing dsp effects
// could be extended to other stuff
//-----------------------------------------------------------------------------
void DSP_FastReset( int dspType ) { int c = ARRAYSIZE( g_PreserveDSP );
// Restore
for ( int i = 0 ; i < c; ++i ) { PreserveDSP_t& slot = g_PreserveDSP[ i ];
if ( slot.cvar == &dsp_player ) { slot.oldvalue = dspType; return; } } }
// Helper to check for change in preset of any of 4 processors
// if switching to a new preset, alloc new preset, simulate both presets in DSP_Process & xfade,
// called a few times per frame.
void CheckNewDspPresets( void ) { bool b_slow_cpu = dsp_slow_cpu.GetInt() == 0 ? false : true;
DSP_CheckRestorePresets();
// room fx are on only if cpu is not slow
int iroom = b_slow_cpu ? 0 : dsp_room.GetInt() ; int ifacingaway = b_slow_cpu ? 0 : dsp_facingaway.GetInt(); int iroomtype = b_slow_cpu ? 0 : dsp_room_type.GetInt(); int ispatial = b_slow_cpu ? 0 : dsp_spatial.GetInt(); int iautomatic = b_slow_cpu ? 0 : dsp_automatic.GetInt();
// always use dsp to process these
int iwater = dsp_water.GetInt(); int iplayer = dsp_player.GetInt(); int ispeaker = dsp_speaker.GetInt();
// check for expired one-shot presets on player and room.
// Only check if a) no new preset has been set and b) not crossfading from previous preset (ie; previous is null)
if ( iplayer == ipset_player_prev && !DSP_IsCrossfading( idsp_player ) ) { if ( DSP_HasExpired ( idsp_player ) ) { iplayer = DSP_OneShotPrevious( idsp_player); // preset has expired - revert to previous preset before one-shot
dsp_player.SetValue(iplayer); } }
if ( iroom == ipset_room_prev && !DSP_IsCrossfading( idsp_room ) ) { if ( DSP_HasExpired ( idsp_room ) ) { iroom = DSP_OneShotPrevious( idsp_room ); // preset has expired - revert to previous preset before one-shot
dsp_room.SetValue(iroom); } }
// legacy code support for "room_type" Cvar
if ( iroomtype != ipset_room_typeprev ) { // force dsp_room = room_type
ipset_room_typeprev = iroomtype; dsp_room.SetValue(iroomtype); }
// NOTE: don't change presets if currently crossfading from a previous preset
if ( iroom != ipset_room_prev && !DSP_IsCrossfading( idsp_room) ) { DSP_SetPreset( idsp_room, iroom ); ipset_room_prev = iroom;
// force room_type = dsp_room
dsp_room_type.SetValue(iroom); ipset_room_typeprev = iroom; }
if ( iwater != ipset_water_prev && !DSP_IsCrossfading( idsp_water) ) { DSP_SetPreset( idsp_water, iwater ); ipset_water_prev = iwater; }
if ( iplayer != ipset_player_prev && !DSP_IsCrossfading( idsp_player)) { DSP_SetPreset( idsp_player, iplayer ); ipset_player_prev = iplayer; }
if ( ifacingaway != ipset_facingaway_prev && !DSP_IsCrossfading( idsp_facingaway) ) { DSP_SetPreset( idsp_facingaway, ifacingaway ); ipset_facingaway_prev = ifacingaway; }
if ( ispeaker != ipset_speaker_prev && !DSP_IsCrossfading( idsp_speaker) ) { DSP_SetPreset( idsp_speaker, ispeaker ); ipset_speaker_prev = ispeaker; }
if ( ispatial != ipset_spatial_prev && !DSP_IsCrossfading( idsp_spatial) ) { DSP_SetPreset( idsp_spatial, ispatial ); ipset_spatial_prev = ispatial; }
if ( iautomatic != ipset_automatic_prev && !DSP_IsCrossfading( idsp_automatic) ) { DSP_SetPreset( idsp_automatic, iautomatic ); ipset_automatic_prev = iautomatic; }
for ( int i = SOUND_BUFFER_SPECIAL_START; i < g_paintBuffers.Count(); ++i ) { paintbuffer_t *pSpecialBuffer = MIX_GetPPaintFromIPaint( i ); if ( pSpecialBuffer->nSpecialDSP != pSpecialBuffer->nPrevSpecialDSP && !DSP_IsCrossfading( pSpecialBuffer->idsp_specialdsp ) ) { DSP_SetPreset( pSpecialBuffer->idsp_specialdsp, pSpecialBuffer->nSpecialDSP ); pSpecialBuffer->nPrevSpecialDSP = pSpecialBuffer->nSpecialDSP; } } }
// create idsp_room preset from set of values, reload the preset.
// modifies psettemplates in place.
// ipreset is the preset # ie: 40
// iproc is the processor to modify within the preset (typically 0)
// pvalues is an array of floating point parameters
// cparams is the # of elements in pvalues
// USED FOR DEBUG ONLY.
void DSP_DEBUGSetParams(int ipreset, int iproc, float *pvalues, int cparams) { pset_t new_pset; // preset
int cparam = clamp (cparams, 0, CPRCPARAMS); prc_t *pprct;
// copy template preset from template array
new_pset = g_psettemplates[ipreset];
// get iproc processor
pprct = &(new_pset.prcs[iproc]);
// copy parameters in to processor
for (int i = 0; i < cparam; i++) { pprct->prm[i] = pvalues[i]; }
// copy constructed preset back into template location
g_psettemplates[ipreset] = new_pset;
// setup new preset
dsp_room.SetValue( 0 );
CheckNewDspPresets();
dsp_room.SetValue( ipreset );
CheckNewDspPresets(); }
// reload entire preset file, reset all current dsp presets
// NOTE: this is debug code only. It doesn't do all mem free work correctly!
void DSP_DEBUGReloadPresetFile( void ) { int iroom = dsp_room.GetInt(); int iwater = dsp_water.GetInt(); int iplayer = dsp_player.GetInt(); // int ifacingaway = dsp_facingaway.GetInt();
// int iroomtype = dsp_room_type.GetInt();
int ispeaker = dsp_speaker.GetInt(); int ispatial = dsp_spatial.GetInt(); // int iautomatic = dsp_automatic.GetInt();
// reload template array
DSP_ReleaseMemory();
DSP_LoadPresetFile();
// force presets to reload
dsp_room.SetValue( 0 ); dsp_water.SetValue( 0 ); dsp_player.SetValue( 0 ); //dsp_facingaway.SetValue( 0 );
//dsp_room_type.SetValue( 0 );
dsp_speaker.SetValue( 0 ); dsp_spatial.SetValue( 0 ); //dsp_automatic.SetValue( 0 );
CUtlVector< int > specialDSPs; for ( int i = SOUND_BUFFER_SPECIAL_START; i < g_paintBuffers.Count(); ++i ) { paintbuffer_t *pSpecialBuffer = MIX_GetPPaintFromIPaint( i );
specialDSPs.AddToTail( pSpecialBuffer->nSpecialDSP ); pSpecialBuffer->nSpecialDSP = 0; } CheckNewDspPresets();
dsp_room.SetValue( iroom ); dsp_water.SetValue( iwater ); dsp_player.SetValue( iplayer ); //dsp_facingaway.SetValue( ifacingaway );
//dsp_room_type.SetValue( iroomtype );
dsp_speaker.SetValue( ispeaker ); dsp_spatial.SetValue( ispatial ); //dsp_automatic.SetValue( iautomatic );
int nSpecialDSPNum = 0; for ( int i = SOUND_BUFFER_SPECIAL_START; i < g_paintBuffers.Count(); ++i ) { paintbuffer_t *pSpecialBuffer = MIX_GetPPaintFromIPaint( i );
pSpecialBuffer->nSpecialDSP = specialDSPs[ nSpecialDSPNum ]; nSpecialDSPNum++; }
CheckNewDspPresets();
// flush dsp automatic nodes
g_bdas_init_nodes = 0; g_bdas_room_init = 0; }
// UNDONE: stock reverb presets:
// carpet hallway
// tile hallway
// wood hallway
// metal hallway
// train tunnel
// sewer main tunnel
// concrete access tunnel
// cave tunnel
// sand floor cave tunnel
// metal duct shaft
// elevator shaft
// large elevator shaft
// parking garage
// aircraft hangar
// cathedral
// train station
// small cavern
// large cavern
// huge cavern
// watery cavern
// long, low cavern
// wood warehouse
// metal warehouse
// concrete warehouse
// small closet room
// medium drywall room
// medium wood room
// medium metal room
// elevator
// small metal room
// medium metal room
// large metal room
// huge metal room
// small metal room dense
// medium metal room dense
// large metal room dense
// huge metal room dense
// small concrete room
// medium concrete room
// large concrete room
// huge concrete room
// small concrete room dense
// medium concrete room dense
// large concrete room dense
// huge concrete room dense
// soundproof room
// carpet lobby
// swimming pool
// open park
// open courtyard
// wide parkinglot
// narrow street
// wide street, short buildings
// wide street, tall buildings
// narrow canyon
// wide canyon
// huge canyon
// small valley
// wide valley
// wreckage & rubble
// small building cluster
// wide open plain
// high vista
// alien interior small
// alien interior medium
// alien interior large
// alien interior huge
// special fx presets:
// alien citadel
// teleport aftershock (these presets all ADSR timeout and reset the dsp_* to 0)
// on target teleport
// off target teleport
// death fade
// beam stasis
// scatterbrain
// pulse only
// slomo
// hypersensitive
// supershocker
// physwhacked
// forcefieldfry
// juiced
// zoomed in
// crabbed
// barnacle gut
// bad transmission
////////////////////////
// Dynamics processing
////////////////////////
// compressor defines
#define COMP_MAX_AMP 32767 // abs max amplitude
#define COMP_THRESH 20000 // start compressing at this threshold
// compress input value - smoothly limit output y to -32767 <= y <= 32767
// UNDONE: not tested or used
inline int S_Compress( int xin ) {
return CLIP( xin >> 2 ); // DEBUG - disabled
float Yn, Xn, Cn, Fn; float C0 = 20000; // threshold
float p = .3; // compression ratio
float g = 1; // gain after compression
Xn = (float)xin;
// Compressor formula:
// Cn = l*Cn-1 + (1-l)*|Xn| // peak detector with memory
// f(Cn) = (Cn/C0)^(p-1) for Cn > C0 // gain function above threshold
// f(Cn) = 1 for C <= C0 // unity gain below threshold
// Yn = f(Cn) * Xn // compressor output
// UNDONE: curves discontinuous at threshold, causes distortion, try catmul-rom
//float l = .5; // compressor memory
//Cn = l * (*pCnPrev) + (1 - l) * fabs((float)xin);
//*pCnPrev = Cn;
Cn = fabs((float)xin);
if (Cn < C0) Fn = 1; else Fn = powf((Cn / C0),(p - 1)); Yn = Fn * Xn * g; //if (Cn > 0)
// Msg("%d -> %d\n", xin, (int)Yn); // DEBUG
//if (fabs(Yn) > 32767)
// Yn = Yn; // DEBUG
return (CLIP((int)Yn)); }
|