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@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