From: drowe67 <drowe67@01035d8c-6547-0410-b346-abe4f91aad63>
Date: Sat, 17 Nov 2012 01:27:55 +0000 (+0000)
Subject: prevent oaccsional NaN from snr_est due to log(0) errors
X-Git-Url: http://git.whiteaudio.com/gitweb/?a=commitdiff_plain;h=7994eab96779479b7d4190bccb617f6f48311aaf;p=freetel-svn-tracking.git

prevent oaccsional NaN from snr_est due to log(0) errors

git-svn-id: https://svn.code.sf.net/p/freetel/code@994 01035d8c-6547-0410-b346-abe4f91aad63
---

diff --git a/codec2-dev/src/fdmdv.c b/codec2-dev/src/fdmdv.c
index 312358bc..b16d8018 100644
--- a/codec2-dev/src/fdmdv.c
+++ b/codec2-dev/src/fdmdv.c
@@ -1274,7 +1274,7 @@ float calc_snr(float sig_est[], float noise_est[])
     S = 0.0;
     for(c=0; c<NC+1; c++)
 	S += pow(sig_est[c], 2.0);
-    SdB = 10.0*log10(S);
+    SdB = 10.0*log10(S+1E-12);
     
     /* Average noise mag across all carriers and square to get an
        average noise power.  This is an estimate of the noise power in
@@ -1286,7 +1286,7 @@ float calc_snr(float sig_est[], float noise_est[])
 	mean += noise_est[c];
     mean /= (NC+1);
     N50 = pow(mean, 2.0);
-    N50dB = 10.0*log10(N50);
+    N50dB = 10.0*log10(N50+1E-12);
 
     /* Now multiply by (3000 Hz)/(50 Hz) to find the total noise power
        in 3000 Hz */