diff options
| author | Ian Hinder <ian.hinder@aei.mpg.de> | 2011-06-09 21:23:21 +0200 |
|---|---|---|
| committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2011-06-09 21:23:21 +0200 |
| commit | ec32fb2e34d6a2e8b7534195c16bbc2008835ec3 (patch) | |
| tree | 4f00f842b8fd614ce2600d86378a6051cc81f162 /m | |
| parent | 6172261aa96876c3fd77121c6dacd2804d0aff7a (diff) | |
McLachlan_BSSN.m: Apply dissipation consistently
Previously, the dissipation term for Xt^i was being added both in the Xt^i RHS and in the B^i RHS, due to the way the temporary variable for dot[X^i} was used in the Gamma driver shift.
Diffstat (limited to 'm')
| -rw-r--r-- | m/McLachlan_BSSN.m | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/m/McLachlan_BSSN.m b/m/McLachlan_BSSN.m index e2adac1..724e1c7 100644 --- a/m/McLachlan_BSSN.m +++ b/m/McLachlan_BSSN.m @@ -700,13 +700,13 @@ evolCalc = + gtu[uj,ul] PD[beta[ui],lj,ll] + (1/3) gtu[ui,uj] PD[beta[ul],lj,ll] + Upwind[beta[uj], Xt[ui], lj] - + epsdiss[uj] PDdiss[Xt[ui],lj] - 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]), - dot[Xt[ui]] -> dotXt[ui], - + dot[Xt[ui]] -> dotXt[ui] + + epsdiss[uj] PDdiss[Xt[ui],lj], + (* PRD 62, 044034 (2000), eqn. (11) *) dottrK -> - em4phi ( gtu[ua,ub] ( PD[alpha,la,lb] + 2 cdphi[la] PD[alpha,lb] ) |