diff options
| author | Peter Diener <diener@linux-hn3d.site> | 2009-04-27 10:36:52 -0500 |
|---|---|---|
| committer | Peter Diener <diener@linux-hn3d.site> | 2009-04-27 10:36:52 -0500 |
| commit | 8e0b4a3127e6b19307297fbc4c27f776363fd55e (patch) | |
| tree | 99c430a126460f619d147ab461ecc62451fe9313 /m | |
| parent | d4ba15953909c5f0a72abf6bd41e8cb169fc737a (diff) | |
| parent | b880dce3c928b74552836fe26211d1e8f89b74c1 (diff) | |
Merge branch 'HEAD' of ssh://relativity.phys.lsu.edu/home/perturbed/gitroot/McLachlan
Conflicts:
ML_BSSN_MP/schedule.ccl
m/McLachlan_BSSN.m
Diffstat (limited to 'm')
| -rw-r--r-- | m/McLachlan_BSSN.m | 42 |
1 files changed, 41 insertions, 1 deletions
diff --git a/m/McLachlan_BSSN.m b/m/McLachlan_BSSN.m index e2765b3..c21df13 100644 --- a/m/McLachlan_BSSN.m +++ b/m/McLachlan_BSSN.m @@ -391,6 +391,9 @@ convertToADMBaseCalc = } }; +(* TODO: Need to apply a boundary condition to the ADM variables in + all cases! *) +(* boundaryCalcADMBase = { Name -> BSSN <> "_ADMBaseBoundary", @@ -421,6 +424,43 @@ boundaryCalcADMBase = dtbetaz -> 0 } }; +*) + +boundaryCalcADMBase = +{ + Name -> BSSN <> "_ADMBaseBoundary", + Schedule -> {"IN " <> BSSN <> "_convertToADMBaseGroup AFTER " <> BSSN <> "_convertToADMBase"}, + Where -> BoundaryWithGhosts, + Shorthands -> {e4phi, g[la,lb], K[la,lb]}, + Equations -> + { + e4phi -> Exp [4 phi], + g[la,lb] -> e4phi gt[la,lb], + gxx -> g11, + gxy -> g12, + gxz -> g13, + gyy -> g22, + gyz -> g23, + gzz -> g33, + K[la,lb] -> e4phi At[la,lb] + (1/3) g[la,lb] trK, + kxx -> K11, + kxy -> K12, + kxz -> K13, + kyy -> K22, + kyz -> K23, + kzz -> K33, + alp -> alpha, + betax -> beta1, + betay -> beta2, + betaz -> beta3, + (* see RHS, but omit derivatives near the boundary *) + dtalp -> - harmonicF alpha^harmonicN + ((1 - LapseAdvectionCoeff) A + LapseAdvectionCoeff trK), + dtbetax -> + ShiftGammaCoeff B1, + dtbetay -> + ShiftGammaCoeff B2, + dtbetaz -> + ShiftGammaCoeff B3 + } +}; (******************************************************************************) (* Evolution equations *) @@ -545,7 +585,7 @@ evolCalc = + addMatter (4 pi alpha (rho + trS)), (* PRD 62, 044034 (2000), eqn. (12) *) - (* TODO: use Hamiltonian constraint to make tracefree *) + (* TODO: Should we use the Hamiltonian constraint to make Rij tracefree? *) Ats[la,lb] -> - CD[alpha,la,lb] + alpha R[la,lb], trAts -> gu[ua,ub] Ats[la,lb], dot[At[la,lb]] -> + em4phi (+ Ats[la,lb] - (1/3) g[la,lb] trAts ) |