ref: ad283e1403a93be876d71ffcc7ff76152c92d453
dir: /tangled/bigverb.c/
#line 93 "bigverb.org"
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#define SK_BIGVERB_PRIV
#include "bigverb.h"
#line 56 "bigverb.org"
struct bigverb_paramset {
int delay; /* in samples, 44.1 kHz */
int drift; /* 1/10 milliseconds */
int randfreq; /* Hertz * 1000 */
int seed;
};
static const struct bigverb_paramset params[8] = {
{0x09a9, 0x0a, 0xc1c, 0x07ae},
{0x0acf, 0x0b, 0xdac, 0x7333},
{0x0c91, 0x11, 0x456, 0x5999},
{0x0de5, 0x06, 0xf85, 0x2666},
{0x0f43, 0x0a, 0x925, 0x50a3},
{0x101f, 0x0b, 0x769, 0x5999},
{0x085f, 0x11, 0x37b, 0x7333},
{0x078d, 0x06, 0xc95, 0x3851}
};
#line 93 "bigverb.org"
#line 851 "bigverb.org"
#define FRACSCALE 0x10000000
#line 861 "bigverb.org"
#define FRACMASK 0xFFFFFFF
#line 869 "bigverb.org"
#define FRACNBITS 28
#line 93 "bigverb.org"
#line 434 "bigverb.org"
static int get_delay_size(const struct bigverb_paramset *p, int sr);
#line 511 "bigverb.org"
static void delay_init(sk_bigverb_delay *d,
const struct bigverb_paramset *p,
SKFLT *buf,
size_t sz,
int sr);
#line 619 "bigverb.org"
static SKFLT delay_compute(sk_bigverb_delay *d,
SKFLT in,
SKFLT fdbk,
SKFLT filt,
int sr);
#line 896 "bigverb.org"
static void generate_next_line(sk_bigverb_delay *d, int sr);
#line 93 "bigverb.org"
#line 148 "bigverb.org"
sk_bigverb * sk_bigverb_new(int sr)
{
sk_bigverb *bv;
bv = calloc(1, sizeof(sk_bigverb));
bv->sr = sr;
#line 207 "bigverb.org"
sk_bigverb_size(bv, 0.93);
#line 233 "bigverb.org"
sk_bigverb_cutoff(bv, 10000.0);
#line 252 "bigverb.org"
bv->pcutoff = -1;
#line 272 "bigverb.org"
bv->filt = 1.0;
#line 402 "bigverb.org"
bv->buf = NULL;
#line 148 "bigverb.org"
#line 407 "bigverb.org"
{
unsigned long total_size;
int i;
SKFLT *buf;
total_size = 0;
buf = NULL;
#line 449 "bigverb.org"
for (i = 0; i < 8; i++) {
total_size += get_delay_size(¶ms[i], sr);
}
#line 407 "bigverb.org"
#line 460 "bigverb.org"
buf = calloc(1, sizeof(SKFLT) * total_size);
bv->buf = buf;
#line 407 "bigverb.org"
#line 476 "bigverb.org"
{
unsigned long bufpos;
bufpos = 0;
for (i = 0; i < 8; i++) {
unsigned int sz;
sz = get_delay_size(¶ms[i], sr);
delay_init(&bv->delay[i], ¶ms[i],
&buf[bufpos], sz, sr);
bufpos += sz;
}
}
#line 417 "bigverb.org"
}
#line 157 "bigverb.org"
return bv;
}
#line 172 "bigverb.org"
void sk_bigverb_del(sk_bigverb *bv)
{
#line 466 "bigverb.org"
free(bv->buf);
#line 175 "bigverb.org"
free(bv);
bv = NULL;
}
#line 212 "bigverb.org"
void sk_bigverb_size(sk_bigverb *bv, SKFLT size)
{
bv->size = size;
}
#line 257 "bigverb.org"
void sk_bigverb_cutoff(sk_bigverb *bv, SKFLT cutoff)
{
bv->cutoff = cutoff;
}
#line 296 "bigverb.org"
void sk_bigverb_tick(sk_bigverb *bv,
SKFLT inL, SKFLT inR,
SKFLT *outL, SKFLT *outR)
{
/* TODO: implement */
SKFLT lsum, rsum;
lsum = 0;
rsum = 0;
#line 329 "bigverb.org"
if (bv->pcutoff != bv->cutoff) {
bv->pcutoff = bv->cutoff;
bv->filt = 2.0 - cos(bv->pcutoff * 2 * M_PI / bv->sr);
bv->filt = bv->filt - sqrt(bv->filt * bv->filt - 1.0);
}
#line 296 "bigverb.org"
#line 344 "bigverb.org"
{
int i;
SKFLT jp;
jp = 0;
for (i = 0; i < 8; i++) {
jp += bv->delay[i].y;
}
jp *= 0.25;
inL = jp + inL;
inR = jp + inR;
}
#line 296 "bigverb.org"
#line 370 "bigverb.org"
{
int i;
for (i = 0; i < 8; i++) {
if (i & 1) {
rsum += delay_compute(&bv->delay[i],
inR,
bv->size,
bv->filt,
bv->sr);
} else {
lsum += delay_compute(&bv->delay[i],
inL,
bv->size,
bv->filt,
bv->sr);
}
}
}
rsum *= 0.35f;
lsum *= 0.35f;
#line 309 "bigverb.org"
*outL = lsum;
*outR = rsum;
}
#line 439 "bigverb.org"
static int get_delay_size(const struct bigverb_paramset *p, int sr)
{
SKFLT sz;
sz = (SKFLT)p->delay/44100 + (p->drift * 0.0001) * 1.125;
return floor(16 + sz*sr);
}
#line 520 "bigverb.org"
static void delay_init(sk_bigverb_delay *d,
const struct bigverb_paramset *p,
SKFLT *buf,
size_t sz,
int sr)
{
SKFLT readpos;
#line 541 "bigverb.org"
d->buf = buf;
d->sz = sz;
#line 554 "bigverb.org"
d->wpos = 0;
#line 578 "bigverb.org"
d->rng = p->seed;
#line 592 "bigverb.org"
readpos = ((SKFLT)p->delay / 44100);
readpos += d->rng * (p->drift * 0.0001) / 32768.0;
readpos = sz - (readpos * sr);
d->irpos = floor(readpos);
d->frpos = floor((readpos - d->irpos) * FRACSCALE);
#line 578 "bigverb.org"
#line 592 "bigverb.org"
#line 886 "bigverb.org"
d->inc = 0;
d->counter = 0;
#line 935 "bigverb.org"
d->maxcount = round((sr / ((SKFLT)p->randfreq * 0.001)));
#line 973 "bigverb.org"
d->dels = p->delay / 44100.0;
#line 983 "bigverb.org"
d->drift = p->drift;
#line 604 "bigverb.org"
generate_next_line(d, sr);
#line 1034 "bigverb.org"
d->y = 0.0;
#line 528 "bigverb.org"
}
#line 628 "bigverb.org"
static SKFLT delay_compute(sk_bigverb_delay *del,
SKFLT in,
SKFLT fdbk,
SKFLT filt,
int sr)
{
SKFLT out;
SKFLT frac_norm;
SKFLT a, b, c, d;
SKFLT s[4];
out = 0;
#line 666 "bigverb.org"
del->buf[del->wpos] = in - del->y;
#line 628 "bigverb.org"
#line 674 "bigverb.org"
del->wpos++;
if (del->wpos >= del->sz) del->wpos -= del->sz;
#line 628 "bigverb.org"
#line 686 "bigverb.org"
if (del->frpos >= FRACSCALE) {
del->irpos += del->frpos >> FRACNBITS;
del->frpos &= FRACMASK;
}
#line 628 "bigverb.org"
#line 696 "bigverb.org"
if (del->irpos >= del->sz) del->irpos -= del->sz;
#line 628 "bigverb.org"
#line 705 "bigverb.org"
frac_norm = del->frpos / (SKFLT)FRACSCALE;
#line 628 "bigverb.org"
#line 718 "bigverb.org"
{
SKFLT tmp[2];
d = ((frac_norm * frac_norm) - 1) / 6.0;
tmp[0] = ((frac_norm + 1.0) * 0.5);
tmp[1] = 3.0 * d;
a = tmp[0] - 1.0 - d;
c = tmp[0] - tmp[1];
b = tmp[1] - frac_norm;
}
#line 628 "bigverb.org"
#line 737 "bigverb.org"
{
int n;
SKFLT *x;
n = del->irpos;
x = del->buf;
if (n > 0 && n < (del->sz - 2)) {
s[0] = x[n - 1];
s[1] = x[n];
s[2] = x[n + 1];
s[3] = x[n + 2];
} else {
int k;
n--;
if (n < 0) n += del->sz;
s[0] = x[n];
for (k = 0; k < 3; k++) {
n++;
if (n >= del->sz) n -= del->sz;
s[k + 1] = x[n];
}
}
}
#line 628 "bigverb.org"
#line 772 "bigverb.org"
out = (a*s[0] + b*s[1] + c*s[2] + d*s[3]) * frac_norm + s[1];
#line 628 "bigverb.org"
#line 780 "bigverb.org"
del->frpos += del->inc;
#line 628 "bigverb.org"
#line 790 "bigverb.org"
out *= fdbk;
out += (del->y - out) * filt;
del->y = out;
#line 628 "bigverb.org"
#line 817 "bigverb.org"
del->counter--;
if (del->counter <= 0) {
generate_next_line(del, sr);
}
#line 650 "bigverb.org"
return out;
}
#line 1008 "bigverb.org"
static void generate_next_line(sk_bigverb_delay *d, int sr)
{
SKFLT curdel;
SKFLT nxtdel;
SKFLT inc;
#line 915 "bigverb.org"
if (d->rng < 0) d->rng += 0x10000;
/* 5^6 = 15625 */
d->rng = (1 + d->rng * 0x3d09);
d->rng &= 0xFFFF;
if (d->rng >= 0x8000) d->rng -= 0x10000;
#line 1008 "bigverb.org"
#line 940 "bigverb.org"
d->counter = d->maxcount;
#line 1008 "bigverb.org"
#line 950 "bigverb.org"
curdel = d->wpos -
(d->irpos + (d->frpos/(SKFLT)FRACSCALE));
while (curdel < 0) curdel += d->sz;
curdel /= sr;
#line 961 "bigverb.org"
nxtdel = (d->rng * (d->drift * 0.0001) / 32768.0) + d->dels;
#line 1008 "bigverb.org"
#line 994 "bigverb.org"
inc = ((curdel - nxtdel) / (SKFLT)d->counter)*sr;
inc += 1;
#line 1008 "bigverb.org"
#line 1003 "bigverb.org"
d->inc = floor(inc * FRACSCALE);
#line 1018 "bigverb.org"
}
#line 93 "bigverb.org"