diff options
author | Barry Wardell <barry.wardell@gmail.com> | 2012-04-27 10:59:10 +0100 |
---|---|---|
committer | Barry Wardell <barry.wardell@gmail.com> | 2012-04-27 10:59:10 +0100 |
commit | f71ebd54d52ffe0b4755ae64fcec82c48a1d76c7 (patch) | |
tree | bd516b6e2110d62fa5cf136b4b0fec593a3b9fd9 /m | |
parent | b5b0de9ec6f63188057c3b8b4bdb020e903738ae (diff) | |
parent | af04ce2e3d96d19b94dd4b416a6617b649c6da34 (diff) |
Merge branch 'master' into CCZ4
Conflicts:
m/McLachlan_BSSN.m
Diffstat (limited to 'm')
-rw-r--r-- | m/McLachlanW.m | 12 | ||||
-rw-r--r-- | m/McLachlan_ADM.m | 2 | ||||
-rw-r--r-- | m/McLachlan_ADMConstraints.m | 8 | ||||
-rw-r--r-- | m/McLachlan_ADMQuantities.m | 10 | ||||
-rw-r--r-- | m/McLachlan_BSSN.m | 37 | ||||
-rw-r--r-- | m/McLachlan_WeylScalars.m | 6 |
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 *) (******************************************************************************) |