diff options
-rw-r--r-- | m/WeylScal4.m | 11 |
1 files 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)]] + } }; |