# HG changeset patch
# User vitor
# Date 1225400657 0
# Node ID 1c07635d73341bd0160bacc2f12558dbd109cc35
# Parent  103e41d627815fa7e92e0969f69fa731f7d6a6a2
Add ff_dot_productf() to celp_math.{c,h}
Part of the QCELP patch by Kenan Gillet, kenan.gillet gmail com

diff -r 103e41d62781 -r 1c07635d7334 celp_math.c
--- a/celp_math.c	Thu Oct 30 21:00:01 2008 +0000
+++ b/celp_math.c	Thu Oct 30 21:04:17 2008 +0000
@@ -195,3 +195,14 @@
 
     return (power_int << 15) + value;
 }
+
+float ff_dot_productf(const float* a, const float* b, int length)
+{
+    float sum = 0;
+    int i;
+
+    for(i=0; i<length; i++)
+        sum += a[i] * b[i];
+
+    return sum;
+}
diff -r 103e41d62781 -r 1c07635d7334 celp_math.h
--- a/celp_math.h	Thu Oct 30 21:00:01 2008 +0000
+++ b/celp_math.h	Thu Oct 30 21:04:17 2008 +0000
@@ -83,4 +83,14 @@
     else           return value <<  offset;
 }
 
+/**
+ * returns the dot product.
+ * @param a input data array
+ * @param b input data array
+ * @param length number of elements
+ *
+ * @return dot product = sum of elementwise products
+ */
+extern float ff_dot_productf(const float* a, const float* b, int length);
+
 #endif /* AVCODEC_CELP_MATH_H */