# HG changeset patch
# User vitor
# Date 1225400737 0
# Node ID 4c95f44c4c23cedbbe396e4be3a9cff77faf1073
# Parent  1c07635d73341bd0160bacc2f12558dbd109cc35
Add a LPC filter
Part of the QCELP patch by Kenan Gillet, kenan.gillet gmail com

diff -r 1c07635d7334 -r 4c95f44c4c23 celp_filters.c
--- a/celp_filters.c	Thu Oct 30 21:04:17 2008 +0000
+++ b/celp_filters.c	Thu Oct 30 21:05:37 2008 +0000
@@ -84,3 +84,24 @@
 
     return 0;
 }
+
+void ff_celp_lp_synthesis_filterf(
+        float *out,
+        const float* filter_coeffs,
+        const float* in,
+        int buffer_length,
+        int filter_length)
+{
+    int i,n;
+
+    // These two lines are to avoid a -1 subtraction in the main loop
+    filter_length++;
+    filter_coeffs--;
+
+    for(n=0; n<buffer_length; n++)
+    {
+        out[n] = in[n];
+        for(i=1; i<filter_length; i++)
+            out[n] += filter_coeffs[i] * out[n-i];
+    }
+}
diff -r 1c07635d7334 -r 4c95f44c4c23 celp_filters.h
--- a/celp_filters.h	Thu Oct 30 21:04:17 2008 +0000
+++ b/celp_filters.h	Thu Oct 30 21:05:37 2008 +0000
@@ -69,4 +69,28 @@
         int stop_on_overflow,
         int rounder);
 
+/**
+ * LP synthesis filter.
+ * @param out [out] pointer to output buffer
+ *        - the array out[-filter_length, -1] must
+ *        contain the previous result of this filter
+ * @param filter_coeffs filter coefficients.
+ * @param in input signal
+ * @param buffer_length amount of data to process
+ * @param filter_length filter length (10 for 10th order LP filter)
+ *
+ * @return 1 if overflow occurred, 0 - otherwise
+ *
+ * @note Output buffer must contain 10 samples of past
+ *       speech data before pointer.
+ *
+ * Routine applies 1/A(z) filter to given speech data.
+ */
+void ff_celp_lp_synthesis_filterf(
+        float *out,
+        const float* filter_coeffs,
+        const float* in,
+        int buffer_length,
+        int filter_length);
+
 #endif /* AVCODEC_CELP_FILTERS_H */