Add simple Riemann symmetries

This cuts Krancing time by approximately a factor of 2.  Needs Kranc
commit 4707e0c3f03cc52d82f88735e602c52cf74f0cbe.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/WeylScal4/trunk@121 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843

1 files changed, 14 insertions, 0 deletions

diff --git a/m/WeylScal4.m b/m/WeylScal4.m
index 053b123..338c708 100644
--- a/m/WeylScal4.m
+++ b/m/WeylScal4.m
@@ -59,6 +59,20 @@ Map[AssertSymmetricDecreasing,
 
 AssertSymmetricDecreasing[gamma[ua,lb,lc], lb, lc];
 
+Tensor[R, i_, j_, k_, l_] /; i > j := -R[j, i, k, l];
+Tensor[R, i_, j_, k_, l_] /; i == j := 0;
+Tensor[R, i_, j_, k_, l_] /; k > l := -R[i, j, l, k];
+Tensor[R, i_, j_, k_, l_] /; k == l := 0;
+Tensor[R, i_, j_, k_, l_] /; i > k := R[k, l, i, j];
+Tensor[R, i_, j_, k_, l_] /; i == k && j > l := R[k, l, i, j];
+
+Tensor[R4p, i_, j_, k_, l_] /; i > j := -R4p[j, i, k, l];
+Tensor[R4p, i_, j_, k_, l_] /; i == j := 0;
+Tensor[R4p, i_, j_, k_, l_] /; k > l := -R4p[i, j, l, k];
+Tensor[R4p, i_, j_, k_, l_] /; k == l := 0;
+Tensor[R4p, i_, j_, k_, l_] /; i > k := R4p[k, l, i, j];
+Tensor[R4p, i_, j_, k_, l_] /; i == k && j > l := R4p[k, l, i, j];
+
 (* Determinants of the metrics in terms of their components
    (Mathematica symbolic expressions) *)
 gDet = Det[MatrixOfComponents[g[la,lb]]];