[Commits] [svn:einsteintoolkit] GRHydro/trunk/ (Rev. 426)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/10/25 11:37 AM

Modified:
 /trunk/
  schedule.ccl
 /trunk/src/
  GRHydro_Boundaries.F90

Log:
 GRHydro: Don't sync Lorentz factor and pressure since they are always computed in Con2Prim and are not needed for initial guesses.
 
 From: Christian Reisswig <reisswig at scriwalker.(none)>

File Changes:

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

File [modified]: GRHydro_Boundaries.F90
Delta lines: +12 -8
===================================================================
--- trunk/src/GRHydro_Boundaries.F90	2012-10-25 16:36:57 UTC (rev 425)
+++ trunk/src/GRHydro_Boundaries.F90	2012-10-25 16:37:00 UTC (rev 426)
@@ -513,13 +513,15 @@
 
   if (CCTK_EQUALS(bound,"flat")) then
     
+    ! Lorentz factor is not needed since it is computed in Con2Prim which runs after any sync!
+    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !      "HydroBase::w_lorentz", "Flat")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-          "HydroBase::w_lorentz", "Flat")
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::rho", "Flat")
+    ! Pressure is not needed since it is computed in Con2Prim which runs after any sync!
+    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !      "HydroBase::press", "Flat")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-          "HydroBase::press", "Flat")
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::eps", "Flat")
     if (general_coordinates .ne. 0) then
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
@@ -559,13 +561,15 @@
   endif
 
   if (CCTK_EQUALS(bound,"none")) then      
+    ! Lorentz factor is not needed since it is computed in Con2Prim which runs after any sync!
+    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !      "HydroBase::w_lorentz", "None")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-          "HydroBase::w_lorentz", "None")
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::rho", "None")
+    ! Pressure is not needed since it is computed in Con2Prim which runs after any sync!
+    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !      "HydroBase::press", "None")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-          "HydroBase::press", "None")
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::eps", "None")
     if (general_coordinates .ne. 0) then
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &

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

File [modified]: schedule.ccl
Delta lines: +2 -2
===================================================================
--- trunk/schedule.ccl	2012-10-25 16:36:57 UTC (rev 425)
+++ trunk/schedule.ccl	2012-10-25 16:37:00 UTC (rev 426)
@@ -1150,8 +1150,8 @@
       {
 	 LANG: Fortran
 	 OPTIONS: LEVEL
-         SYNC: HydroBase::w_lorentz
-         SYNC: HydroBase::press
+#         SYNC: HydroBase::w_lorentz   # don't need that since computed in Con2Prim which runs in postregrid!
+#         SYNC: HydroBase::press       # don't need that since computed in Con2Prim which runs in postregrid!
 #         SYNC: HydroBase::vel
 #         SYNC: HydroBase::Bvec
          SYNC: hydrobase::entropy