mirror/kissfft
1
0
Fork 0
mirror of https://github.com/mborgerding/kissfft.git synced 2025-05-27 21:20:27 -04:00
Code Issues Packages Projects Releases Wiki Activity

fixed error measuremenet function

Browse Source
This commit is contained in:
Mark Borgerding 2004-04-03 00:26:30 +00:00
parent dd689e1661
commit 2de422432e
1 changed files with 4 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
test/mk_test.py
View File

@ -70,15 +70,14 @@ double snr_compare( kiss_fft_cpx * test_vec_out,kiss_fft_cpx * testbuf, int n)
{ {
int k; int k;
double sigpow,noisepow,err,snr,scale=0; double sigpow,noisepow,err,snr,scale=0;
kiss_fft_cpx err;
sigpow = noisepow = .000000000000000000000000000001; sigpow = noisepow = .000000000000000000000000000001;
for (k=0;k<n;++k) { for (k=0;k<n;++k) {
sigpow += test_vec_out[k].r * test_vec_out[k].i + sigpow += test_vec_out[k].r * test_vec_out[k].r +
test_vec_out[k].i * test_vec_out[k].i; test_vec_out[k].i * test_vec_out[k].i;
err = test_vec_out[k].r - testbuf[k].r; C_SUB(err,test_vec_out[k],testbuf[k].r);
noisepow += err * err; noisepow += err.r * err.r + err.i + err.i;
err = test_vec_out[k].i - testbuf[k].i;
noisepow += err * err;
if (test_vec_out[k].r) if (test_vec_out[k].r)
scale += testbuf[k].r / test_vec_out[k].r; scale += testbuf[k].r / test_vec_out[k].r;