symmTensor: Add support for the innerSqr function

which takes the inner product of a symmTensor with itself and returns a symmTensor

applications/test/tensor/Test-tensor.C

@@ -24,9 +24,11 @@ int main()
 
     tensor t6(1,0,-4,0,5,4,-4,4,3);
     //tensor t6(1,2,0,2,5,0,0,0,0);
+
     Info<< "tensor " << t6 << endl;
     vector e = eigenValues(t6);
     Info<< "eigenvalues " << e << endl;
+
     tensor ev = eigenVectors(t6);
     Info<< "eigenvectors " << ev << endl;
 
@@ -58,6 +60,17 @@ int main()
     Info<< "Check for dot product of symmetric tensors "
         << (st1 & st2) << endl;
 
+    Info<< "Check for inner sqr of a symmetric tensor "
+        << innerSqr(st1) << " " << innerSqr(st1) - (st1 & st1) << endl;
+
+    Info<< "Check for symmetric part of dot product of symmetric tensors "
+        << twoSymm(st1&st2) - ((st1&st2) + (st2&st1)) << endl;
+
+    tensor sk1 = skew(t6);
+    tensor sk2 = skew(t7);
+    Info<< "Check for symmetric part of dot product of skew tensors "
+        << twoSymm(sk1&sk2) - ((sk1&sk2) - (sk2&sk1)) << endl;
+
     vector v1(1, 2, 3);
 
     Info<< sqr(v1) << endl;
@@ -184,6 +197,5 @@ int main()
     Debug(U_triple);
     */
 
-
     return 0;
 }