shithub: opus

Download patch

ref: 4508ebd58562e23980157a3e954085bf31640870
parent: 44830b044069aa5870906087aae3ccf78d5749f5
author: Jean-Marc Valin <[email protected]>
date: Sun Feb 24 18:48:15 EST 2008 

real fft snr comparison testcase includes DC


--- a/tests/real-fft-test.c
+++ b/tests/real-fft-test.c
@@ -26,12 +26,15 @@
     int k;
     double sigpow=1e-10, noisepow=1e-10, err,snr;
 
-    for (k=1;k<n;++k) {
+    vec1[0].i = vec1[n].r;
+    for (k=0;k<n;++k) {
         sigpow += (double)vec1[k].r * (double)vec1[k].r + 
                   (double)vec1[k].i * (double)vec1[k].i;
         err = (double)vec1[k].r - (double)vec2[2*k];
+        /*printf ("%f %f\n", (double)vec1[k].r, (double)vec2[2*k]);*/
         noisepow += err * err;
         err = (double)vec1[k].i - (double)vec2[2*k+1];
+        /*printf ("%f %f\n", (double)vec1[k].i, (double)vec2[2*k+1]);*/
         noisepow += err * err;
 
     }
@@ -49,7 +52,7 @@
     int k;
     double sigpow=1e-10, noisepow=1e-10, err,snr;
 
-    for (k=1;k<n;++k) {
+    for (k=0;k<n;++k) {
         sigpow += (double)vec1[k] * (double)vec1[k];
         err = (double)vec1[k] - (double)vec2[k];
         noisepow += err * err;
@@ -82,8 +85,6 @@
     kiss_fft_scalar rout[NFFT+2];
     kiss_fft_scalar zero;
     memset(&zero,0,sizeof(zero) ); // ugly way of setting short,int,float,double, or __m128 to zero
-
-    srand(time(0));
 
     for (i=0;i<NFFT;++i) {
         rin[i] = rand_scalar();