diff options
| author | Erik Schnetter <schnetter@gmail.com> | 2012-05-02 09:06:16 -0400 |
|---|---|---|
| committer | Erik Schnetter <schnetter@gmail.com> | 2012-05-02 09:06:16 -0400 |
| commit | 620e3b49672b16ad59089b4c8388a3099228614f (patch) | |
| tree | e8a22fc11704ef58ccc84026a3e5bce41cd53a77 /m | |
| parent | 7c501f7b1a37a66326cdf8140b7c91246187156b (diff) | |
Replace pi by Pi
Diffstat (limited to 'm')
| -rw-r--r-- | m/McLachlan_BSSN_Peter.m | 10 |
1 files changed, 4 insertions, 6 deletions
diff --git a/m/McLachlan_BSSN_Peter.m b/m/McLachlan_BSSN_Peter.m index c5b5b7d..52a887f 100644 --- a/m/McLachlan_BSSN_Peter.m +++ b/m/McLachlan_BSSN_Peter.m @@ -196,8 +196,6 @@ ddetgtExpr[la_] = Sum [D[Det[MatrixOfComponents[gt[la, lb]]], X] PD[X, la],
{X, Union[Flatten[MatrixOfComponents[gt[la, lb]]]]}];
-pi = N[Pi,40];
-
(******************************************************************************)
(* Groups *)
(******************************************************************************)
@@ -531,14 +529,14 @@ evolCalc = - Xtn[uj] PD[beta[ui],lj]
+ (2/3) Xtn[ui] PD[beta[uj],lj]
(* Equation (4.28) in Baumgarte & Shapiro (Phys. Rept. 376 (2003) 41-131) *)
- + addMatter (- 16 pi alpha gtu[ui,uj] S[lj]),
+ + addMatter (- 16 Pi alpha gtu[ui,uj] S[lj]),
(* PRD 62, 044034 (2000), eqn. (11) *)
dot[trK] -> - gu[ua,ub] CD[alpha,la,lb]
+ alpha (Atm[ua,lb] Atm[ub,la] + (1/3) trK^2)
+ beta[ua] PDu[trK,la]
(* Equation (4.21) in Baumgarte & Shapiro (Phys. Rept. 376 (2003) 41-131) *)
- + addMatter (4 pi alpha (rho + trS)),
+ + addMatter (4 Pi alpha (rho + trS)),
(* PRD 62, 044034 (2000), eqn. (12) *)
(* TODO: use Hamiltonian constraint to make tracefree *)
@@ -757,7 +755,7 @@ constraintsCalc = (* H -> trR - Km[ua,lb] Km[ub,la] + trK^2, *)
(* PRD 67, 084023 (2003), eqn. (19) *)
- H -> trR - Atm[ua,lb] Atm[ub,la] + (2/3) trK^2 - addMatter 16 pi rho,
+ H -> trR - Atm[ua,lb] Atm[ub,la] + (2/3) trK^2 - addMatter 16 Pi rho,
(* gK[la,lb,lc] -> CD[K[la,lb],lc], *)
(* gK[la,lb,lc] -> + 4 e4phi PD[phi,lc] At[la,lb] + e4phi CD[At[la,lb],lc]
@@ -767,7 +765,7 @@ constraintsCalc = M[li] -> + gtu[uj,uk] (CDt[At[li,lj],lk] + 6 At[li,lj] PD[phi,lk])
- (2/3) PD[trK,li]
- - addMatter 8 pi S[li],
+ - addMatter 8 Pi S[li],
(* TODO: use PRD 67, 084023 (2003), eqn. (20) *)
(* det gamma-tilde *)
|