From 5eb006bed93b75e545381c713b8f9bff35d7fe20 Mon Sep 17 00:00:00 2001 From: Steffen Vogel Date: Fri, 15 Mar 2019 17:31:11 +0100 Subject: [PATCH] dp: remove stability tweak --- lib/hooks/dp.c | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/lib/hooks/dp.c b/lib/hooks/dp.c index cf1a954ed..f64789454 100644 --- a/lib/hooks/dp.c +++ b/lib/hooks/dp.c @@ -56,7 +56,6 @@ struct dp { static void dp_step(struct dp *d, double *in, float complex *out) { int n = d->window.steps; - double r = 0.9999999999; double complex om, corr; double newest = *in; double oldest = window_update(&d->window, newest); @@ -65,8 +64,7 @@ static void dp_step(struct dp *d, double *in, float complex *out) om = 2.0 * M_PI * J * d->fharmonics[i] / n; /* Recursive update */ - //d->coeffs[i] = cexp(om) * (d->coeffs[i] + (newest - oldest)); - d->coeffs[i] = d->coeffs[i] * r * cexp(om) - powf(r, n) * oldest + newest; + d->coeffs[i] = cexp(om) * (d->coeffs[i] + (newest - oldest)); /* Correction for stationary phasor */ corr = cexp(-om * (d->t - (d->window.steps + 1)));