aboutsummaryrefslogtreecommitdiff
path: root/m
diff options
context:
space:
mode:
authorBarry Wardell <barry.wardell@gmail.com>2012-04-27 10:59:10 +0100
committerBarry Wardell <barry.wardell@gmail.com>2012-04-27 10:59:10 +0100
commitf71ebd54d52ffe0b4755ae64fcec82c48a1d76c7 (patch)
treebd516b6e2110d62fa5cf136b4b0fec593a3b9fd9 /m
parentb5b0de9ec6f63188057c3b8b4bdb020e903738ae (diff)
parentaf04ce2e3d96d19b94dd4b416a6617b649c6da34 (diff)
Merge branch 'master' into CCZ4
Conflicts: m/McLachlan_BSSN.m
Diffstat (limited to 'm')
-rw-r--r--m/McLachlanW.m12
-rw-r--r--m/McLachlan_ADM.m2
-rw-r--r--m/McLachlan_ADMConstraints.m8
-rw-r--r--m/McLachlan_ADMQuantities.m10
-rw-r--r--m/McLachlan_BSSN.m37
-rw-r--r--m/McLachlan_WeylScalars.m6
6 files changed, 30 insertions, 45 deletions
diff --git a/m/McLachlanW.m b/m/McLachlanW.m
index 4883dde..e8a329b 100644
--- a/m/McLachlanW.m
+++ b/m/McLachlanW.m
@@ -107,8 +107,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 *)
(******************************************************************************)
@@ -481,14 +479,14 @@ evolCalcBSSNW =
- 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)
+ ( betam[ua] PDm[trK,la] + betap[ua] PDp[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 *)
@@ -501,7 +499,7 @@ evolCalcBSSNW =
+ At[la,lc] PD[beta[uc],lb] + At[lb,lc] PD[beta[uc],la]
- (2/3) At[la,lb] PD[beta[uc],lc]
(* Equation (4.23) in Baumgarte & Shapiro (Phys. Rept. 376 (2003) 41-131) *)
- + addMatter (- W2 alpha 8 pi
+ + addMatter (- W2 alpha 8 Pi
(T[la,lb] - (1/3) g[la,lb] trS)),
(* dot[alpha] -> - harmonicF alpha^harmonicN trK, *)
@@ -682,7 +680,7 @@ constraintsCalcBSSNW =
(* 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]
@@ -692,7 +690,7 @@ constraintsCalcBSSNW =
M[li] -> + gtu[uj,uk] (CDt[At[li,lj],lk] + 6 At[li,lj] pdphi[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 *)
diff --git a/m/McLachlan_ADM.m b/m/McLachlan_ADM.m
index 92f0590..94ab2ae 100644
--- a/m/McLachlan_ADM.m
+++ b/m/McLachlan_ADM.m
@@ -82,8 +82,6 @@ ddetgExpr[la_] =
Sum [D[Det[MatrixOfComponents[g[la, lb]]], X] PD[X, la],
{X, Union[Flatten[MatrixOfComponents[g[la, lb]]]]}];
-pi = N[Pi,40];
-
(******************************************************************************)
(* Groups *)
(******************************************************************************)
diff --git a/m/McLachlan_ADMConstraints.m b/m/McLachlan_ADMConstraints.m
index 19346e7..eaddae6 100644
--- a/m/McLachlan_ADMConstraints.m
+++ b/m/McLachlan_ADMConstraints.m
@@ -82,14 +82,10 @@ ddetgExpr[la_] =
Sum [D[Det[MatrixOfComponents[g[la, lb]]], X] PD[X, la],
{X, Union[Flatten[MatrixOfComponents[g[la, lb]]]]}];
-pi = N[Pi,40];
-
(******************************************************************************)
(* Groups *)
(******************************************************************************)
-SetGroupTimelevels[g_,tl_] = Join[g, {Timelevels -> tl}];
-
evolvedGroups = {};
evaluatedGroups =
{SetGroupName [CreateGroupFromTensor [H ], prefix <> "Ham"],
@@ -163,9 +159,9 @@ ADMConstraintsCalc =
(* ADM constraints *)
H -> + trR - Km[ua,lb] Km[ub,la] + trK^2
- - addMatter 16 pi rho,
+ - addMatter 16 Pi rho,
M[la] -> + gu[ub,uc] (CD[K[lc,la], lb] - CD[K[lc,lb], la])
- - addMatter 8 pi S[la]
+ - addMatter 8 Pi S[la]
}
};
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])
}
diff --git a/m/McLachlan_BSSN.m b/m/McLachlan_BSSN.m
index 30cda03..6943a3e 100644
--- a/m/McLachlan_BSSN.m
+++ b/m/McLachlan_BSSN.m
@@ -159,8 +159,6 @@ T11=eTxx; T12=eTxy; T22=eTyy; T13=eTxz; T23=eTyz; T33=eTzz;
(* Expressions *)
(******************************************************************************)
-pi = N[Pi,40];
-
detgExpr = Det [MatrixOfComponents [g [la,lb]]];
ddetgExpr[la_] =
Sum [D[Det[MatrixOfComponents[g[la, lb]]], X] PD[X, la],
@@ -171,7 +169,7 @@ ddetgtExpr[la_] =
Sum [D[Det[MatrixOfComponents[gt[la, lb]]], X] PD[X, la],
{X, Union[Flatten[MatrixOfComponents[gt[la, lb]]]]}];
-etaExpr = Min [SpatialBetaDriverRadius / r, 1];
+etaExpr = SpatialBetaDriverRadius / Max [r, SpatialBetaDriverRadius];
thetaExpr = Min [Exp [1 - r / SpatialShiftGammaCoeffRadius], 1];
@@ -180,8 +178,6 @@ thetaExpr = Min [Exp [1 - r / SpatialShiftGammaCoeffRadius], 1];
(* Groups *)
(******************************************************************************)
-SetGroupTimelevels[g_,tl_] = Join[g, {Timelevels -> tl}];
-
evolvedGroups =
{SetGroupName [CreateGroupFromTensor [phi ], prefix <> "log_confac"],
SetGroupName [CreateGroupFromTensor [gt[la,lb]], prefix <> "metric" ],
@@ -642,7 +638,7 @@ evolCalc =
+ 2 gtu[ui,uj] (alpha PD[Tet,lj] - Tet PD[alpha,lj])
- 2 alpha e4phi dampk1 Zet[ui]
(* 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]),
dot[Xt[ui]] -> dotXt[ui],
(* Equation for Theta *)
@@ -665,7 +661,7 @@ evolCalc =
+ 2 dotTet + 2 PD[alpha,la] Zet[ua]
+ dampk1 alpha Tet - dampk1k2 alpha Tet
(* Equation (4.21) in Baumgarte & Shapiro (Phys. Rept. 376 (2003) 41-131) *)
- + addMatter (4 pi alpha (rho + trS)),
+ + addMatter (4 Pi alpha (rho + trS)),
dot[trK] -> dottrK,
(* PRD 62, 044034 (2000), eqn. (12) *)
@@ -683,7 +679,7 @@ evolCalc =
(* damping term in trA, alternative to trA cleaning *)
(* - (dampA/3) alpha At[lc,ld] gtu[uc,ud] gt[la,lb] *)
(* Equation (4.23) in Baumgarte & Shapiro (Phys. Rept. 376 (2003) 41-131) *)
- + addMatter (- em4phi alpha 8 pi
+ + addMatter (- em4phi alpha 8 Pi
(T[la,lb] - (1/3) g[la,lb] trS)),
(* dot[alpha] -> - harmonicF alpha^harmonicN trK, *)
@@ -700,10 +696,9 @@ evolCalc =
*)
dot[alpha] -> - harmonicF alpha^harmonicN
(+ LapseACoeff A
- + (1 - LapseACoeff) (trK - 2 Tet)),
+ + (1 - LapseACoeff) (trK + AlphaDriver (alpha - 1) - 2 Tet)),
dot[A] -> + LapseACoeff (dottrK - 2 dotTet - AlphaDriver A),
-
eta -> etaExpr,
theta -> thetaExpr,
@@ -756,16 +751,22 @@ advectCalc =
dot[At[la,lb]] -> dot[At[la,lb]] + Upwind[beta[uc], At[la,lb], lc],
- dot[alpha] -> dot[alpha] + LapseAdvectionCoeff Upwind[beta[ua], alpha, la],
+ dot[alpha] -> dot[alpha]
+ + LapseAdvectionCoeff Upwind[beta[ua], alpha, la],
- dot[A] -> dot[A] + LapseAdvectionCoeff Upwind[beta[ua], A, la],
+ dot[A] -> dot[A]
+ + LapseACoeff (
+ + LapseAdvectionCoeff Upwind[beta[ua], A, la]
+ + (1 - LapseAdvectionCoeff) Upwind[beta[ua], trK, la]),
- dot[beta[ua]] -> dot[beta[ua]] + ShiftAdvectionCoeff Upwind[beta[ub], beta[ua], lb],
+ dot[beta[ua]] -> dot[beta[ua]]
+ + ShiftAdvectionCoeff Upwind[beta[ub], beta[ua], lb],
dot[B[ua]] -> dot[B[ua]]
- + ShiftBCoeff Upwind[beta[uj], Xt[ua], lj] (* take care *)
- + ShiftAdvectionCoeff Upwind[beta[ub], B[ua], lb]
- - ShiftAdvectionCoeff Upwind[beta[ub], Xt[ua], lb]
+ + ShiftBCoeff (
+ + ShiftAdvectionCoeff Upwind[beta[ub], B[ua], lb]
+ + ((1 - ShiftAdvectionCoeff)
+ Upwind[beta[ub], Xt[ua], lb]))
}
};
@@ -1082,7 +1083,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]
@@ -1092,7 +1093,7 @@ constraintsCalc =
M[li] -> + gtu[uj,uk] (CDt[At[li,lj],lk] + 6 At[li,lj] cdphi[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 *)
diff --git a/m/McLachlan_WeylScalars.m b/m/McLachlan_WeylScalars.m
index 5e8d805..df10897 100644
--- a/m/McLachlan_WeylScalars.m
+++ b/m/McLachlan_WeylScalars.m
@@ -105,12 +105,6 @@ Map [DefineTensor,
dtbetax, dtbetay, dtbetaz}];
(******************************************************************************)
-(* Expressions *)
-(******************************************************************************)
-
-pi = N[Pi,40];
-
-(******************************************************************************)
(* Groups *)
(******************************************************************************)