diff options
author | Peter Diener <diener@diener-3.lsu.edu> | 2008-07-01 14:21:21 -0500 |
---|---|---|
committer | Peter Diener <diener@diener-3.lsu.edu> | 2008-07-01 14:21:21 -0500 |
commit | 650e98e017aeecbff39cf9bdb0691aeb5f3e13ab (patch) | |
tree | 3deb750f996876321818833a487c31d3de1b8c81 /m | |
parent | e4cf3ce9b6b4bf007ca1134db90b1beb94bd672a (diff) | |
parent | 9c1c9c787b183dba5ef41fbfb9028dbbf8b986e5 (diff) |
Merge commit 'origin/master'
Diffstat (limited to 'm')
-rw-r--r-- | m/McLachlan_BSSN.m | 56 |
1 files changed, 10 insertions, 46 deletions
diff --git a/m/McLachlan_BSSN.m b/m/McLachlan_BSSN.m index bc6abea..7f32d26 100644 --- a/m/McLachlan_BSSN.m +++ b/m/McLachlan_BSSN.m @@ -271,8 +271,7 @@ convertFromADMBaseGammaCalc = convertToADMBaseCalc = { Name -> BSSN <> "_convertToADMBase", - Schedule -> {"IN MoL_PostStep AFTER (" <> BSSN <> "_ApplyBCs " <> BSSN <> "_enforce)"}, - ConditionalOnKeyword -> {"evolution_method", BSSN}, + Schedule -> {"IN ML_BSSN_convertToADMBaseGroup"}, Where -> Interior, Shorthands -> {e4phi, g[la,lb], K[la,lb]}, Equations -> @@ -302,48 +301,7 @@ convertToADMBaseCalc = + LapseAdvectionCoeff beta[ua] PD[alpha,la], dtbetax -> + ShiftGammaCoeff B1 + ShiftAdvectionCoeff beta[ub] PD[beta[ua],lb], - dtbetay -> + ShiftGammaCoeff B2 - + ShiftAdvectionCoeff beta[ub] PD[beta[ua],lb], - dtbetaz -> + ShiftGammaCoeff B3 - + ShiftAdvectionCoeff beta[ub] PD[beta[ua],lb] - } -}; - -convertToADMBaseLapseShiftCalc = -{ - Name -> BSSN <> "_convertToADMBase", - Schedule -> {"IN MoL_PostStep AFTER (" <> BSSN <> "_ApplyBCs " <> BSSN <> "_enforce)"}, - ConditionalOnKeyword -> {"evolution_method", BSSN}, - Where -> Interior, - 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 *) - dtalp -> - harmonicF alpha^harmonicN - ((1 - LapseAdvectionCoeff) A + LapseAdvectionCoeff trK) - + LapseAdvectionCoeff beta[ua] PD[alpha,la], - dtbetax -> + ShiftGammaCoeff B1 - + ShiftAdvectionCoeff beta[ub] PD[beta[ua],lb], - dtbetay -> + ShiftGammaCoeff B2 + dtbetay -> + ShiftGammaCoeff B2 + ShiftAdvectionCoeff beta[ub] PD[beta[ua],lb], dtbetaz -> + ShiftGammaCoeff B3 + ShiftAdvectionCoeff beta[ub] PD[beta[ua],lb] @@ -353,7 +311,7 @@ convertToADMBaseLapseShiftCalc = boundaryCalcADMBase = { Name -> BSSN <> "_ADMBaseBoundary", - Schedule -> {"IN MoL_PostStep AFTER " <> BSSN <> "_convertToADMBase"}, + Schedule -> {"IN ML_BSSN_convertToADMBaseGroup AFTER " <> BSSN <> "_convertToADMBase"}, ConditionalOnKeyword -> {"my_boundary_condition", "Minkowski"}, Where -> BoundaryWithGhosts, Equations -> @@ -567,7 +525,6 @@ enforceCalc = { Name -> BSSN <> "_enforce", Schedule -> {"IN MoL_PostStep BEFORE " <> BSSN <> "_BoundConds"}, - ConditionalOnKeyword -> {"evolution_method", BSSN}, Shorthands -> {detgt, gtu[ua,ub], trAt}, Equations -> { @@ -802,6 +759,13 @@ keywordParameters = (* Description -> "ddd", *) AllowedValues -> {"none", "Minkowski"}, Default -> "none" + }, + { + Name -> "calculate_ADMBase_variables_at", + Visibility -> "restricted", + (* Description -> "ddd", *) + AllowedValues -> {"MoL_PostStep", "CCTK_EVOL", "CCTK_ANALYSIS"}, + Default -> "MoL_PostStep" } }; |