#include <stdbool.h>
#include <stdio.h>
+#ifdef COMPILER_MSVC
+#include <float.h>
+#define isfinite_local(val) (bool)_finite((double)val)
+#else
+#define isfinite_local isfinite
+#endif
+
#include "lv2/lv2plug.in/ns/lv2core/lv2.h"
#ifdef LV2_EXTENDED
self->s[0] = self->s[1] = 0.0;
}
+static void linear_svf_protect(struct linear_svf *self)
+{
+ if (!isfinite_local (self->s[0]) || !isfinite_local (self->s[1])) {
+ linear_svf_reset (self);
+ }
+}
+
typedef struct {
float* f0[BANDS];
float* g[BANDS];
static void linear_svf_set_peq(struct linear_svf *self, float gdb, float sample_rate, float cutoff, float bandwidth)
{
double f0 = (double)cutoff;
- double q = (double)pow(2.0, 1.0 / bandwidth) / (pow(2.0, bandwidth) - 1.0);
+ double q = (double)pow(2.0, 0.5 * bandwidth) / (pow(2.0, bandwidth) - 1.0);
double sr = (double)sample_rate;
double A = pow(10.0, gdb/40.0);
if (!is_eq(aeq->v_f0[i], *aeq->f0[i], 0.1)) {
aeq->v_f0[i] += tau * (*aeq->f0[i] - aeq->v_f0[i]);
changed = true;
- } else {
- aeq->v_f0[i] = *aeq->f0[i];
}
if (*aeq->filtog[i] <= 0 || *aeq->enable <= 0) {
if (!is_eq(aeq->v_g[i], 0.f, 0.05)) {
aeq->v_g[i] += tau * (0.0 - aeq->v_g[i]);
changed = true;
- } else {
- aeq->v_g[i] = 0.0;
}
} else {
if (!is_eq(aeq->v_g[i], *aeq->g[i], 0.05)) {
aeq->v_g[i] += tau * (*aeq->g[i] - aeq->v_g[i]);
changed = true;
- } else {
- aeq->v_g[i] = *aeq->g[i];
}
}
if (!is_eq(aeq->v_bw[i], *aeq->bw[i], 0.001)) {
aeq->v_bw[i] += tau * (*aeq->bw[i] - aeq->v_bw[i]);
changed = true;
- } else {
- aeq->v_bw[i] = *aeq->bw[i];
}
}
offset += block;
}
+ for (uint32_t j = 0; j < BANDS; j++) {
+ linear_svf_protect(&aeq->v_filter[j]);
+ }
+
#ifdef LV2_EXTENDED
if (aeq->need_expose && aeq->queue_draw) {
aeq->need_expose = false;