From 2de422432ee0ebc9455bd374b59f921e0b0f824c Mon Sep 17 00:00:00 2001
From: Mark Borgerding <mark@borgerding.net>
Date: Sat, 3 Apr 2004 00:26:30 +0000
Subject: [PATCH] fixed error measuremenet function

---
 test/mk_test.py | 9 ++++-----
 1 file changed, 4 insertions(+), 5 deletions(-)

diff --git a/test/mk_test.py b/test/mk_test.py
index 7481c06..4271e8c 100755
--- a/test/mk_test.py
+++ b/test/mk_test.py
@@ -70,15 +70,14 @@ double snr_compare( kiss_fft_cpx * test_vec_out,kiss_fft_cpx * testbuf, int n)
 {
     int k;
     double sigpow,noisepow,err,snr,scale=0;
+    kiss_fft_cpx err;
     sigpow = noisepow = .000000000000000000000000000001; 
 
     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;
-        err = test_vec_out[k].r - testbuf[k].r;
-        noisepow += err * err;
-        err = test_vec_out[k].i - testbuf[k].i;
-        noisepow += err * err;
+        C_SUB(err,test_vec_out[k],testbuf[k].r);
+        noisepow += err.r * err.r + err.i + err.i;
 
         if (test_vec_out[k].r)
             scale += testbuf[k].r / test_vec_out[k].r;