Mirrors/VILLASnode
1
0
Fork 0
mirror of https://git.rwth-aachen.de/acs/public/villas/node/ synced 2025-03-09 00:00:00 +01:00
Code Issues Releases Wiki Activity

dp: remove stability tweak

Browse source
This commit is contained in:
Steffen Vogel 2019-03-15 17:31:11 +01:00
parent 3ca1ae89b8
commit 5eb006bed9
1 changed files with 1 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
lib/hooks/dp.c
View file

@ -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)));