[Commits] [svn:einsteintoolkit] WeylScal4/trunk/m/ (Rev. 121)

[ian.hinder at aei.mpg.de]

User: hinder
Date: 2012/08/11 12:24 PM

Modified:
 /trunk/m/
  WeylScal4.m

Log:
 Add simple Riemann symmetries
 
 This cuts Krancing time by approximately a factor of 2.  Needs Kranc
 commit 4707e0c3f03cc52d82f88735e602c52cf74f0cbe.

File Changes:

Directory: /trunk/m/
====================

File [modified]: WeylScal4.m
Delta lines: +14 -0
===================================================================
--- trunk/m/WeylScal4.m	2012-08-11 17:24:42 UTC (rev 120)
+++ trunk/m/WeylScal4.m	2012-08-11 17:24:50 UTC (rev 121)
@@ -59,6 +59,20 @@
 
 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]]];