aboutsummaryrefslogtreecommitdiff
path: root/m/McLachlan_ADMQuantities.m
diff options
context:
space:
mode:
authorErik Schnetter <schnetter@gmail.com>2012-01-02 11:58:25 -0500
committerErik Schnetter <schnetter@gmail.com>2012-01-02 11:58:25 -0500
commit71894c54a1d778df2b27277c004157aa6be93da7 (patch)
tree9f30a20a4b8853960f660ee5b2ba4ce53238204b /m/McLachlan_ADMQuantities.m
parent0e41ae0288c7830f38d3306d14fe8a66861d15aa (diff)
Don't replace Pi by its numeric equivalent -- Kranc does this already
Diffstat (limited to 'm/McLachlan_ADMQuantities.m')
-rw-r--r--m/McLachlan_ADMQuantities.m10
1 files changed, 4 insertions, 6 deletions
diff --git a/m/McLachlan_ADMQuantities.m b/m/McLachlan_ADMQuantities.m
index 003a9f8..2f78d56 100644
--- a/m/McLachlan_ADMQuantities.m
+++ b/m/McLachlan_ADMQuantities.m
@@ -149,8 +149,6 @@ T11=eTxx; T12=eTxy; T22=eTyy; T13=eTxz; T23=eTyz; T33=eTzz;
detgtExpr = Det [MatrixOfComponents [gt[la,lb]]];
-pi = N[Pi,40];
-
(******************************************************************************)
(* Groups *)
(******************************************************************************)
@@ -257,16 +255,16 @@ ADMQuantitiesCalc =
(* ADM quantities *)
(* See PRD 66, 084026 (2002) *)
- Madm -> 1/(16 pi)
- (+ ephi^5 (+ 16 pi addMatter rho
+ Madm -> 1/(16 Pi)
+ (+ ephi^5 (+ 16 Pi addMatter rho
+ Atm[ua,lb] Atm[ub,la]
- 2/3 trK^2)
- gtu[ua,ub] Gt[uc,la,ld] Gtlu[lb,lc,ud]
+ (1 - ephi) trRt),
- Jadm[li] -> 1/(16 pi) Eps[li,lj,uk] ephi^6
+ Jadm[li] -> 1/(16 Pi) Eps[li,lj,uk] ephi^6
(+ 2 Atm[uj,lk]
- + 16 pi x[uj] S[lk]
+ + 16 Pi x[uj] S[lk]
+ 4/3 x[uj] PD[trK,lk]
- x[uj] dgtu[ul,um,lk] At[ll,lm])
}