diff options
author | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-08-02 21:32:00 +0000 |
---|---|---|
committer | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-08-02 21:32:00 +0000 |
commit | 415ab7729def9fad475a4f40e7906ebeb49a29ee (patch) | |
tree | cdf3ba087941dd4909fc382ec8610f4f748c68f2 /src/GRHydro_RoeSolver.F90 | |
parent | 4b2e68805d07a963e081b8d4284ee1d6e621efb9 (diff) |
* Optimize: remove support for shift_state = 0 (except for shock tubes and
Cowling calculations of spherically symmetric objects, there is no reason
not to have storage for the shift).
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@259 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_RoeSolver.F90')
-rw-r--r-- | src/GRHydro_RoeSolver.F90 | 22 |
1 files changed, 9 insertions, 13 deletions
diff --git a/src/GRHydro_RoeSolver.F90 b/src/GRHydro_RoeSolver.F90 index 26f9999..18d2f4d 100644 --- a/src/GRHydro_RoeSolver.F90 +++ b/src/GRHydro_RoeSolver.F90 @@ -107,21 +107,17 @@ subroutine GRHydro_RoeSolve(CCTK_ARGUMENTS) !!$ Set metric terms at interface - if (shift_state .ne. 0) then - if (flux_direction == 1) then - avg_beta = 0.5d0 * (betax(i+xoffset,j+yoffset,k+zoffset) + & - betax(i,j,k)) - else if (flux_direction == 2) then - avg_beta = 0.5d0 * (betay(i+xoffset,j+yoffset,k+zoffset) + & + if (flux_direction == 1) then + avg_beta = 0.5d0 * (betax(i+xoffset,j+yoffset,k+zoffset) + & + betax(i,j,k)) + else if (flux_direction == 2) then + avg_beta = 0.5d0 * (betay(i+xoffset,j+yoffset,k+zoffset) + & betay(i,j,k)) - else if (flux_direction == 3) then - avg_beta = 0.5d0 * (betaz(i+xoffset,j+yoffset,k+zoffset) + & - betaz(i,j,k)) - else - call CCTK_WARN(0, "Flux direction not x,y,z") - end if + else if (flux_direction == 3) then + avg_beta = 0.5d0 * (betaz(i+xoffset,j+yoffset,k+zoffset) + & + betaz(i,j,k)) else - avg_beta = 0.d0 + call CCTK_WARN(0, "Flux direction not x,y,z") end if avg_alp = 0.5 * (alp(i,j,k) + alp(i+xoffset,j+yoffset,k+zoffset)) |