aboutsummaryrefslogtreecommitdiff
path: root/ML_BSSN/schedule.ccl
diff options
context:
space:
mode:
Diffstat (limited to 'ML_BSSN/schedule.ccl')
-rw-r--r--ML_BSSN/schedule.ccl240
1 files changed, 137 insertions, 103 deletions
diff --git a/ML_BSSN/schedule.ccl b/ML_BSSN/schedule.ccl
index 2aa1baf..4fe7352 100644
--- a/ML_BSSN/schedule.ccl
+++ b/ML_BSSN/schedule.ccl
@@ -17,7 +17,7 @@ STORAGE: ML_Kdot[timelevels]
STORAGE: ML_Xtdot[timelevels]
-STORAGE: ML_phidot[timelevels]
+STORAGE: ML_beta_term[timelevels]
STORAGE: ML_shift[timelevels]
@@ -243,13 +243,6 @@ schedule ML_BSSN_RHS IN ML_BSSN_evolCalcGroup
WRITES: ML_BSSN::gt22rhs(Interior)
WRITES: ML_BSSN::gt23rhs(Interior)
WRITES: ML_BSSN::gt33rhs(Interior)
- WRITES: ML_BSSN::Kdot11(Interior)
- WRITES: ML_BSSN::Kdot12(Interior)
- WRITES: ML_BSSN::Kdot13(Interior)
- WRITES: ML_BSSN::Kdot22(Interior)
- WRITES: ML_BSSN::Kdot23(Interior)
- WRITES: ML_BSSN::Kdot33(Interior)
- WRITES: ML_BSSN::phidot(Interior)
WRITES: ML_BSSN::phirhs(Interior)
WRITES: ML_BSSN::term1(Interior)
WRITES: ML_BSSN::term2(Interior)
@@ -258,19 +251,60 @@ schedule ML_BSSN_RHS IN ML_BSSN_evolCalcGroup
WRITES: ML_BSSN::Xt1rhs(Interior)
WRITES: ML_BSSN::Xt2rhs(Interior)
WRITES: ML_BSSN::Xt3rhs(Interior)
- WRITES: ML_BSSN::Xtdot1(Interior)
- WRITES: ML_BSSN::Xtdot2(Interior)
- WRITES: ML_BSSN::Xtdot3(Interior)
} "ML_BSSN_RHS"
-schedule ML_BSSN_lapse_evol IN ML_BSSN_evolCalcGroup
+ schedule eval_timediff IN MoL_PreStep BEFORE quasimaximal_slicing_axi_solve
+ #schedule eval_timediff IN ML_BSSN_evolCalcGroup BEFORE quasimaximal_slicing_axi_solve
+ {
+ LANG: C
+ READS: grid::r(Everywhere)
+ READS: ML_BSSN::At11(Everywhere)
+ READS: ML_BSSN::At12(Everywhere)
+ READS: ML_BSSN::At13(Everywhere)
+ READS: ML_BSSN::At22(Everywhere)
+ READS: ML_BSSN::At23(Everywhere)
+ READS: ML_BSSN::At33(Everywhere)
+ READS: ML_BSSN::B1(Everywhere)
+ READS: ML_BSSN::B2(Everywhere)
+ READS: ML_BSSN::B3(Everywhere)
+ READS: ML_BSSN::Xt1(Everywhere)
+ READS: ML_BSSN::Xt2(Everywhere)
+ READS: ML_BSSN::Xt3(Everywhere)
+ READS: ML_BSSN::alpha(Everywhere)
+ READS: ML_BSSN::phi(Everywhere)
+ READS: ML_BSSN::gt11(Everywhere)
+ READS: ML_BSSN::gt12(Everywhere)
+ READS: ML_BSSN::gt13(Everywhere)
+ READS: ML_BSSN::gt22(Everywhere)
+ READS: ML_BSSN::gt23(Everywhere)
+ READS: ML_BSSN::gt33(Everywhere)
+ READS: ML_BSSN::beta1(Everywhere)
+ READS: ML_BSSN::beta2(Everywhere)
+ READS: ML_BSSN::beta3(Everywhere)
+ READS: ML_BSSN::trK(Everywhere)
+ WRITES: ML_BSSN::Kdot11(Interior)
+ WRITES: ML_BSSN::Kdot12(Interior)
+ WRITES: ML_BSSN::Kdot13(Interior)
+ WRITES: ML_BSSN::Kdot22(Interior)
+ WRITES: ML_BSSN::Kdot23(Interior)
+ WRITES: ML_BSSN::Kdot33(Interior)
+ WRITES: ML_BSSN::beta_term(Interior)
+ WRITES: ML_BSSN::Xtdot1(Interior)
+ WRITES: ML_BSSN::Xtdot2(Interior)
+ WRITES: ML_BSSN::Xtdot3(Interior)
+ } "eval_timediff"
+
+if (CCTK_EQUALS(lapse_evolution_method, "ML_BSSN"))
{
- LANG: C
- READS: QuasiMaximalSlicing::W(Everywhere)
- READS: ML_BSSN::alpha(Everywhere)
- READS: ML_BSSN::trK(Everywhere)
- WRITES: ML_BSSN::alpharhs(Interior)
-} "ML_BSSN_lapse_evol"
+ schedule ML_BSSN_lapse_evol IN ML_BSSN_evolCalcGroup
+ {
+ LANG: C
+ READS: QuasiMaximalSlicing::W(Everywhere)
+ READS: ML_BSSN::alpha(Everywhere)
+ READS: ML_BSSN::trK(Everywhere)
+ WRITES: ML_BSSN::alpharhs(Interior)
+ } "ML_BSSN_lapse_evol"
+}
if (CCTK_EQUALS(apply_dissipation, "always"))
@@ -579,91 +613,91 @@ schedule ML_BSSN_convertToADMBase IN ML_BSSN_convertToADMBaseGroup
} "ML_BSSN_convertToADMBase"
-if (CCTK_EQUALS(dt_lapse_shift_method, "correct"))
-{
- schedule ML_BSSN_convertToADMBaseDtLapseShift IN ML_BSSN_convertToADMBaseGroup
- {
- LANG: C
- SYNC: ADMBase::dtshift
- READS: grid::r(Everywhere)
- READS: ML_BSSN::B1(Everywhere)
- READS: ML_BSSN::B2(Everywhere)
- READS: ML_BSSN::B3(Everywhere)
- READS: ML_BSSN::Xt1(Everywhere)
- READS: ML_BSSN::Xt2(Everywhere)
- READS: ML_BSSN::Xt3(Everywhere)
- READS: ML_BSSN::alpha(Everywhere)
- READS: ML_BSSN::phi(Everywhere)
- READS: ML_BSSN::gt11(Everywhere)
- READS: ML_BSSN::gt12(Everywhere)
- READS: ML_BSSN::gt13(Everywhere)
- READS: ML_BSSN::gt22(Everywhere)
- READS: ML_BSSN::gt23(Everywhere)
- READS: ML_BSSN::gt33(Everywhere)
- READS: ML_BSSN::beta1(Everywhere)
- READS: ML_BSSN::beta2(Everywhere)
- READS: ML_BSSN::beta3(Everywhere)
- WRITES: ADMBase::dtbetax(Interior)
- WRITES: ADMBase::dtbetay(Interior)
- WRITES: ADMBase::dtbetaz(Interior)
- } "ML_BSSN_convertToADMBaseDtLapseShift"
-}
-
-
-if (CCTK_EQUALS(dt_lapse_shift_method, "correct"))
-{
- schedule ML_BSSN_convertToADMBaseDtLapseShiftBoundary IN ML_BSSN_convertToADMBaseGroup
- {
- LANG: C
- READS: grid::r(Everywhere)
- READS: ML_BSSN::B1(Everywhere)
- READS: ML_BSSN::B2(Everywhere)
- READS: ML_BSSN::B3(Everywhere)
- READS: ML_BSSN::Xt1(Everywhere)
- READS: ML_BSSN::Xt2(Everywhere)
- READS: ML_BSSN::Xt3(Everywhere)
- READS: ML_BSSN::gt11(Everywhere)
- READS: ML_BSSN::gt12(Everywhere)
- READS: ML_BSSN::gt13(Everywhere)
- READS: ML_BSSN::gt22(Everywhere)
- READS: ML_BSSN::gt23(Everywhere)
- READS: ML_BSSN::gt33(Everywhere)
- READS: ML_BSSN::beta1(Everywhere)
- READS: ML_BSSN::beta2(Everywhere)
- READS: ML_BSSN::beta3(Everywhere)
- WRITES: ADMBase::dtbetax(ERROR(BoundaryWithGhosts))
- WRITES: ADMBase::dtbetay(ERROR(BoundaryWithGhosts))
- WRITES: ADMBase::dtbetaz(ERROR(BoundaryWithGhosts))
- } "ML_BSSN_convertToADMBaseDtLapseShiftBoundary"
-}
-
-
-if (CCTK_EQUALS(dt_lapse_shift_method, "noLapseShiftAdvection"))
-{
- schedule ML_BSSN_convertToADMBaseFakeDtLapseShift IN ML_BSSN_convertToADMBaseGroup
- {
- LANG: C
- READS: grid::r(Everywhere)
- READS: ML_BSSN::B1(Everywhere)
- READS: ML_BSSN::B2(Everywhere)
- READS: ML_BSSN::B3(Everywhere)
- READS: ML_BSSN::Xt1(Everywhere)
- READS: ML_BSSN::Xt2(Everywhere)
- READS: ML_BSSN::Xt3(Everywhere)
- READS: ML_BSSN::gt11(Everywhere)
- READS: ML_BSSN::gt12(Everywhere)
- READS: ML_BSSN::gt13(Everywhere)
- READS: ML_BSSN::gt22(Everywhere)
- READS: ML_BSSN::gt23(Everywhere)
- READS: ML_BSSN::gt33(Everywhere)
- READS: ML_BSSN::beta1(Everywhere)
- READS: ML_BSSN::beta2(Everywhere)
- READS: ML_BSSN::beta3(Everywhere)
- WRITES: ADMBase::dtbetax(Everywhere)
- WRITES: ADMBase::dtbetay(Everywhere)
- WRITES: ADMBase::dtbetaz(Everywhere)
- } "ML_BSSN_convertToADMBaseFakeDtLapseShift"
-}
+#if (CCTK_EQUALS(dt_lapse_shift_method, "correct"))
+#{
+# schedule ML_BSSN_convertToADMBaseDtLapseShift IN ML_BSSN_convertToADMBaseGroup
+# {
+# LANG: C
+# SYNC: ADMBase::dtshift
+# READS: grid::r(Everywhere)
+# READS: ML_BSSN::B1(Everywhere)
+# READS: ML_BSSN::B2(Everywhere)
+# READS: ML_BSSN::B3(Everywhere)
+# READS: ML_BSSN::Xt1(Everywhere)
+# READS: ML_BSSN::Xt2(Everywhere)
+# READS: ML_BSSN::Xt3(Everywhere)
+# READS: ML_BSSN::alpha(Everywhere)
+# READS: ML_BSSN::phi(Everywhere)
+# READS: ML_BSSN::gt11(Everywhere)
+# READS: ML_BSSN::gt12(Everywhere)
+# READS: ML_BSSN::gt13(Everywhere)
+# READS: ML_BSSN::gt22(Everywhere)
+# READS: ML_BSSN::gt23(Everywhere)
+# READS: ML_BSSN::gt33(Everywhere)
+# READS: ML_BSSN::beta1(Everywhere)
+# READS: ML_BSSN::beta2(Everywhere)
+# READS: ML_BSSN::beta3(Everywhere)
+# WRITES: ADMBase::dtbetax(Interior)
+# WRITES: ADMBase::dtbetay(Interior)
+# WRITES: ADMBase::dtbetaz(Interior)
+# } "ML_BSSN_convertToADMBaseDtLapseShift"
+#}
+
+
+#if (CCTK_EQUALS(dt_lapse_shift_method, "correct"))
+#{
+# schedule ML_BSSN_convertToADMBaseDtLapseShiftBoundary IN ML_BSSN_convertToADMBaseGroup
+# {
+# LANG: C
+# READS: grid::r(Everywhere)
+# READS: ML_BSSN::B1(Everywhere)
+# READS: ML_BSSN::B2(Everywhere)
+# READS: ML_BSSN::B3(Everywhere)
+# READS: ML_BSSN::Xt1(Everywhere)
+# READS: ML_BSSN::Xt2(Everywhere)
+# READS: ML_BSSN::Xt3(Everywhere)
+# READS: ML_BSSN::gt11(Everywhere)
+# READS: ML_BSSN::gt12(Everywhere)
+# READS: ML_BSSN::gt13(Everywhere)
+# READS: ML_BSSN::gt22(Everywhere)
+# READS: ML_BSSN::gt23(Everywhere)
+# READS: ML_BSSN::gt33(Everywhere)
+# READS: ML_BSSN::beta1(Everywhere)
+# READS: ML_BSSN::beta2(Everywhere)
+# READS: ML_BSSN::beta3(Everywhere)
+# WRITES: ADMBase::dtbetax(ERROR(BoundaryWithGhosts))
+# WRITES: ADMBase::dtbetay(ERROR(BoundaryWithGhosts))
+# WRITES: ADMBase::dtbetaz(ERROR(BoundaryWithGhosts))
+# } "ML_BSSN_convertToADMBaseDtLapseShiftBoundary"
+#}
+#
+#
+#if (CCTK_EQUALS(dt_lapse_shift_method, "noLapseShiftAdvection"))
+#{
+# schedule ML_BSSN_convertToADMBaseFakeDtLapseShift IN ML_BSSN_convertToADMBaseGroup
+# {
+# LANG: C
+# READS: grid::r(Everywhere)
+# READS: ML_BSSN::B1(Everywhere)
+# READS: ML_BSSN::B2(Everywhere)
+# READS: ML_BSSN::B3(Everywhere)
+# READS: ML_BSSN::Xt1(Everywhere)
+# READS: ML_BSSN::Xt2(Everywhere)
+# READS: ML_BSSN::Xt3(Everywhere)
+# READS: ML_BSSN::gt11(Everywhere)
+# READS: ML_BSSN::gt12(Everywhere)
+# READS: ML_BSSN::gt13(Everywhere)
+# READS: ML_BSSN::gt22(Everywhere)
+# READS: ML_BSSN::gt23(Everywhere)
+# READS: ML_BSSN::gt33(Everywhere)
+# READS: ML_BSSN::beta1(Everywhere)
+# READS: ML_BSSN::beta2(Everywhere)
+# READS: ML_BSSN::beta3(Everywhere)
+# WRITES: ADMBase::dtbetax(Everywhere)
+# WRITES: ADMBase::dtbetay(Everywhere)
+# WRITES: ADMBase::dtbetaz(Everywhere)
+# } "ML_BSSN_convertToADMBaseFakeDtLapseShift"
+#}
schedule group ML_BSSN_constraints_group in MoL_PseudoEvolution after MoL_PostStep
{
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-06 08:48:15 +0000