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