Skip to content

Commit

Permalink
Use tabs
Browse files Browse the repository at this point in the history
  • Loading branch information
@sakertooth
sakertooth committed Mar 12, 2024
1 parent f7a723f commit 07afd8a
Showing 1 changed file with 156 additions and 156 deletions.
312 changes: 156 additions & 156 deletions plugins/ReverbSC/revsc.c
View file
Original file line number Diff line number Diff line change
Expand Up @@ -34,14 +34,14 @@
/* reverbParams[n][3] = random seed (0 - 32767) */

static const SPFLOAT reverbParams[8][4] = {
{ (2473.0 / DEFAULT_SRATE), 0.0010, 3.100, 1966.0 },
{ (2767.0 / DEFAULT_SRATE), 0.0011, 3.500, 29491.0 },
{ (3217.0 / DEFAULT_SRATE), 0.0017, 1.110, 22937.0 },
{ (3557.0 / DEFAULT_SRATE), 0.0006, 3.973, 9830.0 },
{ (3907.0 / DEFAULT_SRATE), 0.0010, 2.341, 20643.0 },
{ (4127.0 / DEFAULT_SRATE), 0.0011, 1.897, 22937.0 },
{ (2143.0 / DEFAULT_SRATE), 0.0017, 0.891, 29491.0 },
{ (1933.0 / DEFAULT_SRATE), 0.0006, 3.221, 14417.0 }
{ (2473.0 / DEFAULT_SRATE), 0.0010, 3.100, 1966.0 },
{ (2767.0 / DEFAULT_SRATE), 0.0011, 3.500, 29491.0 },
{ (3217.0 / DEFAULT_SRATE), 0.0017, 1.110, 22937.0 },
{ (3557.0 / DEFAULT_SRATE), 0.0006, 3.973, 9830.0 },
{ (3907.0 / DEFAULT_SRATE), 0.0010, 2.341, 20643.0 },
{ (4127.0 / DEFAULT_SRATE), 0.0011, 1.897, 22937.0 },
{ (2143.0 / DEFAULT_SRATE), 0.0017, 0.891, 29491.0 },
{ (1933.0 / DEFAULT_SRATE), 0.0006, 3.221, 14417.0 }
};

static int delay_line_max_samples(SPFLOAT sr, SPFLOAT iPitchMod, int n);
Expand All @@ -50,113 +50,113 @@ static int delay_line_bytes_alloc(SPFLOAT sr, SPFLOAT iPitchMod, int n);
static const SPFLOAT outputGain = 0.35;
static const SPFLOAT jpScale = 0.25;
int sp_revsc_create(sp_revsc **p){
*p = malloc(sizeof(sp_revsc));
return SP_OK;
*p = malloc(sizeof(sp_revsc));
return SP_OK;
}

int sp_revsc_init(sp_data *sp, sp_revsc *p)
{
p->iSampleRate = sp->sr;
p->sampleRate = sp->sr;
p->feedback = 0.97;
p->lpfreq = 10000;
p->iPitchMod = 1;
p->iSkipInit = 0;
p->dampFact = 1.0;
p->prv_LPFreq = 0.0;
p->initDone = 1;
p->iSampleRate = sp->sr;
p->sampleRate = sp->sr;
p->feedback = 0.97;
p->lpfreq = 10000;
p->iPitchMod = 1;
p->iSkipInit = 0;
p->dampFact = 1.0;
p->prv_LPFreq = 0.0;
p->initDone = 1;
int i = 0, nBytes = 0;
for(i = 0; i < 8; i++){
nBytes += delay_line_bytes_alloc(sp->sr, 1, i);
}
sp_auxdata_alloc(&p->aux, nBytes);
nBytes = 0;
for (i = 0; i < 8; i++) {
p->delayLines[i].buf = (SPFLOAT*) (((char*) p->aux.ptr) + nBytes);
init_delay_line(p, &p->delayLines[i], i);
nBytes += delay_line_bytes_alloc(sp->sr, 1, i);
}

return SP_OK;
nBytes += delay_line_bytes_alloc(sp->sr, 1, i);
}
sp_auxdata_alloc(&p->aux, nBytes);
nBytes = 0;
for (i = 0; i < 8; i++) {
p->delayLines[i].buf = (SPFLOAT*) (((char*) p->aux.ptr) + nBytes);
init_delay_line(p, &p->delayLines[i], i);
nBytes += delay_line_bytes_alloc(sp->sr, 1, i);
}

return SP_OK;
}


int sp_revsc_destroy(sp_revsc **p)
{
sp_revsc *pp = *p;
sp_auxdata_free(&pp->aux);
free(*p);
return SP_OK;
sp_revsc *pp = *p;
sp_auxdata_free(&pp->aux);
free(*p);
return SP_OK;
}

static int delay_line_max_samples(SPFLOAT sr, SPFLOAT iPitchMod, int n)
{
SPFLOAT maxDel = 0.0f;

maxDel = reverbParams[n][0];
maxDel += (reverbParams[n][1] * (SPFLOAT) iPitchMod * 1.125);
return (int) (maxDel * sr + 16.5);
maxDel += (reverbParams[n][1] * (SPFLOAT) iPitchMod * 1.125);
return (int) (maxDel * sr + 16.5);
}

static int delay_line_bytes_alloc(SPFLOAT sr, SPFLOAT iPitchMod, int n)
{
int nBytes = 0;
int nBytes = 0;

nBytes += (delay_line_max_samples(sr, iPitchMod, n) * (int) sizeof(SPFLOAT));
return nBytes;
nBytes += (delay_line_max_samples(sr, iPitchMod, n) * (int) sizeof(SPFLOAT));
return nBytes;
}

static void next_random_lineseg(sp_revsc *p, sp_revsc_dl *lp, int n)
{
SPFLOAT prvDel = 0.0f, nxtDel = 0.0f, phs_incVal = 0.0f;

/* update random seed */
if (lp->seedVal < 0)
lp->seedVal += 0x10000;
lp->seedVal = (lp->seedVal * 15625 + 1) & 0xFFFF;
if (lp->seedVal >= 0x8000)
lp->seedVal -= 0x10000;
/* length of next segment in samples */
lp->randLine_cnt = (int) ((p->sampleRate / reverbParams[n][2]) + 0.5);
prvDel = (SPFLOAT) lp->writePos;
prvDel -= ((SPFLOAT) lp->readPos
+ ((SPFLOAT) lp->readPosFrac / (SPFLOAT) DELAYPOS_SCALE));
while (prvDel < 0.0)
prvDel += lp->bufferSize;
prvDel = prvDel / p->sampleRate; /* previous delay time in seconds */
nxtDel = (SPFLOAT) lp->seedVal * reverbParams[n][1] / 32768.0;
/* next delay time in seconds */
nxtDel = reverbParams[n][0] + (nxtDel * (SPFLOAT) p->iPitchMod);
/* calculate phase increment per sample */
phs_incVal = (prvDel - nxtDel) / (SPFLOAT) lp->randLine_cnt;
phs_incVal = phs_incVal * p->sampleRate + 1.0;
lp->readPosFrac_inc = (int) (phs_incVal * DELAYPOS_SCALE + 0.5);
if (lp->seedVal < 0)
lp->seedVal += 0x10000;
lp->seedVal = (lp->seedVal * 15625 + 1) & 0xFFFF;
if (lp->seedVal >= 0x8000)
lp->seedVal -= 0x10000;
/* length of next segment in samples */
lp->randLine_cnt = (int) ((p->sampleRate / reverbParams[n][2]) + 0.5);
prvDel = (SPFLOAT) lp->writePos;
prvDel -= ((SPFLOAT) lp->readPos
+ ((SPFLOAT) lp->readPosFrac / (SPFLOAT) DELAYPOS_SCALE));
while (prvDel < 0.0)
prvDel += lp->bufferSize;
prvDel = prvDel / p->sampleRate; /* previous delay time in seconds */
nxtDel = (SPFLOAT) lp->seedVal * reverbParams[n][1] / 32768.0;
/* next delay time in seconds */
nxtDel = reverbParams[n][0] + (nxtDel * (SPFLOAT) p->iPitchMod);
/* calculate phase increment per sample */
phs_incVal = (prvDel - nxtDel) / (SPFLOAT) lp->randLine_cnt;
phs_incVal = phs_incVal * p->sampleRate + 1.0;
lp->readPosFrac_inc = (int) (phs_incVal * DELAYPOS_SCALE + 0.5);
}

static int init_delay_line(sp_revsc *p, sp_revsc_dl *lp, int n)
{
SPFLOAT readPos = 0.0f;
/* int i; */

/* calculate length of delay line */
lp->bufferSize = delay_line_max_samples(p->sampleRate, 1, n);
lp->dummy = 0;
lp->writePos = 0;
/* set random seed */
lp->seedVal = (int) (reverbParams[n][3] + 0.5);
/* set initial delay time */
readPos = (SPFLOAT) lp->seedVal * reverbParams[n][1] / 32768;
readPos = reverbParams[n][0] + (readPos * (SPFLOAT) p->iPitchMod);
readPos = (SPFLOAT) lp->bufferSize - (readPos * p->sampleRate);
lp->readPos = (int) readPos;
readPos = (readPos - (SPFLOAT) lp->readPos) * (SPFLOAT) DELAYPOS_SCALE;
lp->readPosFrac = (int) (readPos + 0.5);
/* initialise first random line segment */
next_random_lineseg(p, lp, n);
/* clear delay line to zero */
lp->filterState = 0.0;
memset(lp->buf, 0, sizeof(SPFLOAT) * lp->bufferSize);
return SP_OK;
/* calculate length of delay line */
lp->bufferSize = delay_line_max_samples(p->sampleRate, 1, n);
lp->dummy = 0;
lp->writePos = 0;
/* set random seed */
lp->seedVal = (int) (reverbParams[n][3] + 0.5);
/* set initial delay time */
readPos = (SPFLOAT) lp->seedVal * reverbParams[n][1] / 32768;
readPos = reverbParams[n][0] + (readPos * (SPFLOAT) p->iPitchMod);
readPos = (SPFLOAT) lp->bufferSize - (readPos * p->sampleRate);
lp->readPos = (int) readPos;
readPos = (readPos - (SPFLOAT) lp->readPos) * (SPFLOAT) DELAYPOS_SCALE;
lp->readPosFrac = (int) (readPos + 0.5);
/* initialise first random line segment */
next_random_lineseg(p, lp, n);
/* clear delay line to zero */
lp->filterState = 0.0;
memset(lp->buf, 0, sizeof(SPFLOAT) * lp->bufferSize);
return SP_OK;
}


Expand All @@ -170,107 +170,107 @@ int sp_revsc_compute(sp_data *sp, sp_revsc *p, SPFLOAT *in1, SPFLOAT *in2, SPFLO
int bufferSize = 0; /* Local copy */
SPFLOAT dampFact = p->dampFact;

if (p->initDone <= 0) return SP_NOT_OK;
if (p->initDone <= 0) return SP_NOT_OK;

/* calculate tone filter coefficient if frequency changed */
/* calculate tone filter coefficient if frequency changed */

if (p->lpfreq != p->prv_LPFreq) {
p->prv_LPFreq = p->lpfreq;
dampFact = 2.0 - cos(p->prv_LPFreq * (2 * M_PI) / p->sampleRate);
dampFact = p->dampFact = dampFact - sqrt(dampFact * dampFact - 1.0);
}
if (p->lpfreq != p->prv_LPFreq) {
p->prv_LPFreq = p->lpfreq;
dampFact = 2.0 - cos(p->prv_LPFreq * (2 * M_PI) / p->sampleRate);
dampFact = p->dampFact = dampFact - sqrt(dampFact * dampFact - 1.0);
}

/* calculate "resultant junction pressure" and mix to input signals */
/* calculate "resultant junction pressure" and mix to input signals */

ainL = aoutL = aoutR = 0.0;
for (n = 0; n < 8; n++) {
ainL += p->delayLines[n].filterState;
}
ainL *= jpScale;
ainR = ainL + *in2;
ainL = ainL + *in1;
ainL = aoutL = aoutR = 0.0;
for (n = 0; n < 8; n++) {
ainL += p->delayLines[n].filterState;
}
ainL *= jpScale;
ainR = ainL + *in2;
ainL = ainL + *in1;

/* loop through all delay lines */
/* loop through all delay lines */

for (n = 0; n < 8; n++) {
lp = &p->delayLines[n];
bufferSize = lp->bufferSize;
for (n = 0; n < 8; n++) {
lp = &p->delayLines[n];
bufferSize = lp->bufferSize;

/* send input signal and feedback to delay line */
/* send input signal and feedback to delay line */

lp->buf[lp->writePos] = (SPFLOAT) ((n & 1 ? ainR : ainL)
- lp->filterState);
if (++lp->writePos >= bufferSize) {
lp->writePos -= bufferSize;
}
lp->buf[lp->writePos] = (SPFLOAT) ((n & 1 ? ainR : ainL)
- lp->filterState);
if (++lp->writePos >= bufferSize) {
lp->writePos -= bufferSize;
}

/* read from delay line with cubic interpolation */
/* read from delay line with cubic interpolation */

if (lp->readPosFrac >= DELAYPOS_SCALE) {
lp->readPos += (lp->readPosFrac >> DELAYPOS_SHIFT);
lp->readPosFrac &= DELAYPOS_MASK;
}
if (lp->readPos >= bufferSize)
lp->readPos -= bufferSize;
readPos = lp->readPos;
frac = (SPFLOAT) lp->readPosFrac * (1.0 / (SPFLOAT) DELAYPOS_SCALE);
if (lp->readPosFrac >= DELAYPOS_SCALE) {
lp->readPos += (lp->readPosFrac >> DELAYPOS_SHIFT);
lp->readPosFrac &= DELAYPOS_MASK;
}
if (lp->readPos >= bufferSize)
lp->readPos -= bufferSize;
readPos = lp->readPos;
frac = (SPFLOAT) lp->readPosFrac * (1.0 / (SPFLOAT) DELAYPOS_SCALE);

/* calculate interpolation coefficients */
/* calculate interpolation coefficients */

a2 = frac * frac; a2 -= 1.0; a2 *= (1.0 / 6.0);
a1 = frac; a1 += 1.0; a1 *= 0.5; am1 = a1 - 1.0;
a0 = 3.0 * a2; a1 -= a0; am1 -= a2; a0 -= frac;
a2 = frac * frac; a2 -= 1.0; a2 *= (1.0 / 6.0);
a1 = frac; a1 += 1.0; a1 *= 0.5; am1 = a1 - 1.0;
a0 = 3.0 * a2; a1 -= a0; am1 -= a2; a0 -= frac;

/* read four samples for interpolation */
/* read four samples for interpolation */

if (readPos > 0 && readPos < (bufferSize - 2)) {
vm1 = (SPFLOAT) (lp->buf[readPos - 1]);
v0 = (SPFLOAT) (lp->buf[readPos]);
v1 = (SPFLOAT) (lp->buf[readPos + 1]);
v2 = (SPFLOAT) (lp->buf[readPos + 2]);
}
else {
if (readPos > 0 && readPos < (bufferSize - 2)) {
vm1 = (SPFLOAT) (lp->buf[readPos - 1]);
v0 = (SPFLOAT) (lp->buf[readPos]);
v1 = (SPFLOAT) (lp->buf[readPos + 1]);
v2 = (SPFLOAT) (lp->buf[readPos + 2]);
}
else {

/* at buffer wrap-around, need to check index */
/* at buffer wrap-around, need to check index */

if (--readPos < 0) readPos += bufferSize;
vm1 = (SPFLOAT) lp->buf[readPos];
if (++readPos >= bufferSize) readPos -= bufferSize;
v0 = (SPFLOAT) lp->buf[readPos];
if (++readPos >= bufferSize) readPos -= bufferSize;
v1 = (SPFLOAT) lp->buf[readPos];
if (++readPos >= bufferSize) readPos -= bufferSize;
v2 = (SPFLOAT) lp->buf[readPos];
}
v0 = (am1 * vm1 + a0 * v0 + a1 * v1 + a2 * v2) * frac + v0;
if (--readPos < 0) readPos += bufferSize;
vm1 = (SPFLOAT) lp->buf[readPos];
if (++readPos >= bufferSize) readPos -= bufferSize;
v0 = (SPFLOAT) lp->buf[readPos];
if (++readPos >= bufferSize) readPos -= bufferSize;
v1 = (SPFLOAT) lp->buf[readPos];
if (++readPos >= bufferSize) readPos -= bufferSize;
v2 = (SPFLOAT) lp->buf[readPos];
}
v0 = (am1 * vm1 + a0 * v0 + a1 * v1 + a2 * v2) * frac + v0;

/* update buffer read position */
/* update buffer read position */

lp->readPosFrac += lp->readPosFrac_inc;
lp->readPosFrac += lp->readPosFrac_inc;

/* apply feedback gain and lowpass filter */
/* apply feedback gain and lowpass filter */

v0 *= (SPFLOAT) p->feedback;
v0 = (lp->filterState - v0) * dampFact + v0;
lp->filterState = v0;
v0 *= (SPFLOAT) p->feedback;
v0 = (lp->filterState - v0) * dampFact + v0;
lp->filterState = v0;

/* mix to output */
/* mix to output */

if (n & 1) {
aoutR += v0;
}else{
aoutL += v0;
}
if (n & 1) {
aoutR += v0;
}else{
aoutL += v0;
}

/* start next random line segment if current one has reached endpoint */
/* start next random line segment if current one has reached endpoint */

if (--(lp->randLine_cnt) <= 0) {
next_random_lineseg(p, lp, n);
}
}
/* someday, use aoutR for multimono out */
if (--(lp->randLine_cnt) <= 0) {
next_random_lineseg(p, lp, n);
}
}
/* someday, use aoutR for multimono out */

*out1 = aoutL * outputGain;
*out2 = aoutR * outputGain;
return SP_OK;
*out1 = aoutL * outputGain;
*out2 = aoutR * outputGain;
return SP_OK;
}

0 comments on commit 07afd8a

Please sign in to comment.