[Commits] [svn:einsteintoolkit] IllinoisGRMHD/trunk/ (Rev. 22)

[zachetie at gmail.com]

User: zetienne
Date: 2014/03/25 11:52 AM

Modified:
 /trunk/
  interface.ccl, schedule.ccl

Log:
 Remove some dead lines from schedule.ccl. Also, add warnings to interface.ccl. Hold off on using IllinoisGRMHD until semi-staggered prolongation/restriction operators are added to the latest Carpet.

File Changes:

Directory: /trunk/
==================

File [modified]: interface.ccl
Delta lines: +4 -0
===================================================================
--- trunk/interface.ccl	2014-03-25 16:29:20 UTC (rev 21)
+++ trunk/interface.ccl	2014-03-25 16:52:35 UTC (rev 22)
@@ -16,17 +16,20 @@
 
 # These variables are semi-staggered:
 # Ax is defined on the semi-staggered grid (i,j+1/2,k+1/2)
+# WARNING: WILL NOT WORK PROPERLY WITHOUT SEMI-STAGGERED PROLONGATION/RESTRICTION:
 #cctk_real em_Ax type = GF Timelevels=3 tags='Prolongation="STAGGER011"'
 cctk_real em_Ax type = GF Timelevels=3
 {
   Ax
 } "x-component of the vector potential, evolved when constrained_transport_scheme==3"
 # Ay is defined on the semi-staggered grid (i+1/2,j,k+1/2)
+# WARNING: WILL NOT WORK PROPERLY WITHOUT SEMI-STAGGERED PROLONGATION/RESTRICTION:
 #cctk_real em_Ay type = GF Timelevels=3 tags='Prolongation="STAGGER101"'
 cctk_real em_Ay type = GF Timelevels=3
 {
   Ay
 } "y-component of the vector potential, evolved when constrained_transport_scheme==3"
+# WARNING: WILL NOT WORK PROPERLY WITHOUT SEMI-STAGGERED PROLONGATION/RESTRICTION:
 # Az is defined on the semi-staggered grid (i+1/2,j+1/2,k)
 #cctk_real em_Az type = GF Timelevels=3 tags='Prolongation="STAGGER110"'
 cctk_real em_Az type = GF Timelevels=3
@@ -35,6 +38,7 @@
 } "z-component of the vector potential, evolved when constrained_transport_scheme==3"
 
 # psi6phi is defined on the staggered grid (i+1/2,j+1/2,k+1/2)
+# WARNING: WILL NOT WORK PROPERLY WITHOUT SEMI-STAGGERED PROLONGATION/RESTRICTION:
 #cctk_real em_psi6phi type = GF Timelevels=3 tags='Prolongation="STAGGER111"'
 cctk_real em_psi6phi type = GF Timelevels=3
 {

File [modified]: schedule.ccl
Delta lines: +1 -21
===================================================================
--- trunk/schedule.ccl	2014-03-25 16:29:20 UTC (rev 21)
+++ trunk/schedule.ccl	2014-03-25 16:52:35 UTC (rev 22)
@@ -100,7 +100,7 @@
 # Note however that B^i depends on derivatives of A_{\mu}, so a SYNC is necessary on B^i.
 schedule compute_B_and_Bstagger_from_A in HydroBase_Boundaries after outer_boundaries_on_A_mu
 {
-  # This is strictly a processor sync, as prolongation is disabled for all primitives & B^i's.
+  # This is strictly a processor sync, as prolongation is disabled for all B^i's.
   SYNC: grmhd_primitives_Bi
   LANG: C
 } "Compute B and B_stagger from A"
@@ -119,23 +119,3 @@
   LANG: C
 } "Apply outflow-only, flat BCs on {P,rho_b,vx,vy,vz}. Outflow only == velocities pointed inward from outer boundary are set to zero."
 ##########################################################
-
-##########################################################
-# CON2PRIM STUFF
-# ( 1) Compute B^i from A_i everywhere, synchronize B^i
-# ( 2) Call con2prim to get primitives on all points (really only need interior, but who cares)
-
-## Easiest primitives to solve for: B^i
-#schedule compute_B_and_Bstagger_from_A in HydroBase_Con2Prim
-#{
-## This is strictly a processor sync, as prolongation is disabled for all primitives & B^i's.
-#  SYNC: grmhd_primitives_Bi
-#  LANG: C
-#} "Compute B and B_stagger from A"
-#
-## Nontrivial primitives solve, for P,rho_b,vx,vy,vz:
-#schedule conserv_to_prims in HydroBase_Con2Prim after compute_B_and_Bstagger_from_A
-#{
-#  LANG: C
-#} "Compute primitive variables from conservatives. This is non-trivial, requiring a Newton-Raphson root-finder."
-##########################################################