ENH: cmptProduct function, i.e. cmptProduct(vector(a, b, c)) = a*b*c.

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;
 }

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