From e6e4a81ce90ed50c834e7ab4566e092113bff2fc Mon Sep 17 00:00:00 2001 From: knarf Date: Wed, 14 Dec 2011 04:20:21 +0000 Subject: Enable-remaining-Psis-and-add-curvIr git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/WeylScal4/trunk@91 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843 --- m/WeylScal4.m | 11 +++++++---- 1 file changed, 7 insertions(+), 4 deletions(-) diff --git a/m/WeylScal4.m b/m/WeylScal4.m index ba54b6d..7a7f86b 100644 --- a/m/WeylScal4.m +++ b/m/WeylScal4.m @@ -35,7 +35,7 @@ Map[DefineTensor, 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 @@ gDet = Det[MatrixOfComponents[g[la,lb]]]; (* 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 @@ PsisCalc[fdOrder_, PD_] := 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 @@ PsisCalc[fdOrder_, PD_] := 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)]] + } }; -- cgit v1.2.3