From 087e9552ecda35df6885be1a98810238ad81f342 Mon Sep 17 00:00:00 2001
From: drowe67 <drowe67@01035d8c-6547-0410-b346-abe4f91aad63>
Date: Tue, 16 Nov 2010 09:35:38 +0000
Subject: [PATCH] support for HPF

git-svn-id: https://svn.code.sf.net/p/freetel/code@236 01035d8c-6547-0410-b346-abe4f91aad63
---
 codec2/src/c2sim.c | 4 +++-
 codec2/src/sine.h  | 4 +++-
 2 files changed, 6 insertions(+), 2 deletions(-)

diff --git a/codec2/src/c2sim.c b/codec2/src/c2sim.c
index 4f809eb8..3fb03b56 100644
--- a/codec2/src/c2sim.c
+++ b/codec2/src/c2sim.c
@@ -125,6 +125,7 @@ int main(int argc, char *argv[])
   float ak_interp[LPC_MAX];
 
   void *nlp_states;
+  float hpf_states[2];
 
   for(i=0; i<M; i++)
       Sn[i] = 1.0;
@@ -146,6 +147,7 @@ int main(int argc, char *argv[])
       prev_lsps[i] = i*PI/(LPC_ORD+1);
   }
   e = prev_e = 1;
+  hpf_states[0] = hpf_states[1] = 0.0;
 
   nlp_states = nlp_create();
 
@@ -240,7 +242,7 @@ int main(int argc, char *argv[])
     for(i=0; i<M-N; i++)
       Sn[i] = Sn[i+N];
     for(i=0; i<N; i++)
-      Sn[i+M-N] = buf[i];
+	Sn[i+M-N] = hpf((float)buf[i], hpf_states);
  
     /* Estimate pitch */
 
diff --git a/codec2/src/sine.h b/codec2/src/sine.h
index 88eee37f..ae578bf7 100644
--- a/codec2/src/sine.h
+++ b/codec2/src/sine.h
@@ -32,10 +32,12 @@
 #include "comp.h"
 
 void make_analysis_window(float w[], COMP W[]);
+float hpf(float x, float states[]);
 void dft_speech(COMP Sw[], float Sn[], float w[]);
 void two_stage_pitch_refinement(MODEL *model, COMP Sw[]);
 void estimate_amplitudes(MODEL *model, COMP Sw[], COMP W[]);
-float est_voicing_mbe(MODEL *model, COMP Sw[], COMP W[], COMP Sw_[],COMP Ew[], float prev_Wo);
+float est_voicing_mbe(MODEL *model, COMP Sw[], COMP W[], COMP Sw_[],COMP Ew[], 
+		      float prev_Wo);
 void make_synthesis_window(float Pn[]);
 void synthesise(float Sn_[], MODEL *model, float Pn[], int shift);
 
-- 
2.25.1