diff --git a/applications/test/vector/vectorTest.C b/applications/test/vector/vectorTest.C
index f48f65d711..65e995ae4a 100644
--- a/applications/test/vector/vectorTest.C
+++ b/applications/test/vector/vectorTest.C
@@ -27,7 +27,11 @@ int main()
     Info<< magSqr(d) << endl;
 
     vector d2(0.5, 0.51, -0.5);
-    Info<< cmptMax(d2) << " " << cmptSum(d2) << " " << cmptMag(d2) << endl;
+    Info<< cmptMax(d2) << " "
+        << cmptSum(d2) << " "
+        << cmptProduct(d2) << " "
+        << cmptMag(d2)
+        << endl;
     Info<< min(d, d2) << endl;
     return 0;
 }
diff --git a/src/OpenFOAM/primitives/VectorSpace/VectorSpaceI.H b/src/OpenFOAM/primitives/VectorSpace/VectorSpaceI.H
index 1eff60d99b..d655990e21 100644
--- a/src/OpenFOAM/primitives/VectorSpace/VectorSpaceI.H
+++ b/src/OpenFOAM/primitives/VectorSpace/VectorSpaceI.H
@@ -406,6 +406,17 @@ inline Cmpt cmptAv
     return cmptSum(vs)/nCmpt;
 }
 
+template<class Form, class Cmpt, int nCmpt>
+inline Cmpt cmptProduct
+(
+    const VectorSpace<Form, Cmpt, nCmpt>& vs
+)
+{
+    Cmpt product = vs.v_[0];
+    VectorSpaceOps<nCmpt,1>::SeqOp(product, vs, multiplyEqOp<Cmpt>());
+    return product;
+}
+
 
 template<class Form, class Cmpt, int nCmpt>
 inline Form cmptMag