diff --git a/libavcodec/celp_math.c b/libavcodec/celp_math.c
index 9fa7dfc70d..df05917d20 100644
--- a/libavcodec/celp_math.c
+++ b/libavcodec/celp_math.c
@@ -195,3 +195,14 @@ int ff_log2(uint32_t value)
 
     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 --git a/libavcodec/celp_math.h b/libavcodec/celp_math.h
index 029b5310af..ce0726fba1 100644
--- a/libavcodec/celp_math.h
+++ b/libavcodec/celp_math.h
@@ -83,4 +83,14 @@ static inline int bidir_sal(int value, int offset)
     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 */