diff --git a/Numeric/Numeric.h b/Numeric/Numeric.h
index 58590b1a52447e29083f1d308641c732fe0f30fa..520481d175d1a6badb5d7b04940edce3ccade19a 100644
--- a/Numeric/Numeric.h
+++ b/Numeric/Numeric.h
@@ -25,6 +25,11 @@ inline double crossProd(double a[3], double b[3], int i)
   int i2 = (i+2) % 3;
   return a[i1]*b[i2] - a[i2]*b[i1];
 }
+inline double scalProd(double a[3], double b[3])
+{
+  return a[0]*b[0] + a[1]*b[1] + a[2]*b[2];
+}
+
 inline void prodve(double a[3], double b[3], double c[3])
 {
   c[2] = a[0] * b[1] - a[1] * b[0];