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

[knarf at cct.lsu.edu]

User: knarf
Date: 2011/12/13 10:20 PM

Modified:
 /trunk/m/
  WeylScal4.m

Log:
 Enable-remaining-Psis-and-add-curvIr

File Changes:

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

File [modified]: WeylScal4.m
Delta lines: +7 -4
===================================================================
--- trunk/m/WeylScal4.m	2011-12-12 20:19:58 UTC (rev 90)
+++ trunk/m/WeylScal4.m	2011-12-14 04:20:21 UTC (rev 91)
@@ -35,7 +35,7 @@
   R, gamma,g, gInv, k, ltet, n, rm, im, rmbar, imbar, tn, va, vb, vc,
   wa, wb, wc, ea, eb, ec,
   Ro, Rojo, R4p, Psi0r, Psi0i, Psi1r, Psi1i, Psi2r, Psi2i, Psi3r,
-  Psi3i, Psi4r,Psi4i
+  Psi3i, Psi4r,Psi4i, curvIr
 }];
 
 (* Psi0,2,4 behave as (pseudo)scalars *)
@@ -69,7 +69,7 @@
 
 (* Cactus group definitions *)
 
-scalars = { (*Psi0r, Psi0i, Psi1r, Psi1i, Psi2r, Psi2i, Psi3r, Psi3i, *) Psi4r,Psi4i};
+scalars = {Psi0r, Psi0i, Psi1r, Psi1i, Psi2r, Psi2i, Psi3r, Psi3i, Psi4r,Psi4i, curvIr};
 
 scalarGroups = Map[CreateGroupFromTensor, scalars];
 
@@ -211,7 +211,7 @@
 
 	Psi4i -> R4p[la,lb,lc,ld] n[ua] n[uc] ( - rm[ub] im[ud] - im[ub] rm[ud] )
 		 + 2 Ro[la,lb,lc] n[ub] nn ( - rm[ua] im[uc] - im[ua] rm[uc] )
-		 + Rojo[la,lb] nn nn ( - rm[ua] im[ub] - im[ua] rm[ub] )  (*,
+		 + Rojo[la,lb] nn nn ( - rm[ua] im[ub] - im[ua] rm[ub] )  ,
 
 
 	Psi3r -> R4p[la,lb,lc,ld] ltet[ua] n[ub] rm[uc] n[ud]
@@ -252,7 +252,10 @@
 
 	Psi0i -> R4p[la,lb,lc,ld] ltet[ua] ltet[uc] (rm[ub] im[ud] + im[ub] rm[ud])
 		 + 2 Ro[la,lb,lc] nn ltet[ub] (rm[ua] im[uc] + im[ua] rm[uc])
-		 + Rojo[la,lb] nn nn (rm[ua] im[ub] + im[ua] rm[ub]) *)
+		 + Rojo[la,lb] nn nn (rm[ua] im[ub] + im[ua] rm[ub]),
+
+        curvIr -> ComplexExpand[Re[3 (Psi2r+I Psi2i)^2 - 4 (Psi1r+I Psi1i) (Psi3r + I Psi3i) + (Psi4r + I Psi4i) (Psi0r + I Psi0i)]]
+
   }
 };