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_FluxSplit.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_FluxSplit.F90')
-rw-r--r-- | src/GRHydro_FluxSplit.F90 | 40 |
1 files changed, 9 insertions, 31 deletions
diff --git a/src/GRHydro_FluxSplit.F90 b/src/GRHydro_FluxSplit.F90 index d54e446..1c701e8 100644 --- a/src/GRHydro_FluxSplit.F90 +++ b/src/GRHydro_FluxSplit.F90 @@ -73,11 +73,7 @@ subroutine GRHydro_FSAlpha(CCTK_ARGUMENTS) gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),& gyz(i,j,k),gzz(i,j,k)) - if (shift_state .eq. 0) then - beta = 0.d0 - else - beta = betax(i,j,k) - end if + beta = betax(i,j,k) call eigenvalues(GRHydro_eos_handle, & rho (i,j,k), & @@ -119,11 +115,7 @@ subroutine GRHydro_FSAlpha(CCTK_ARGUMENTS) gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),& gyz(i,j,k),gzz(i,j,k)) - if (shift_state .eq. 0) then - beta = 0.d0 - else - beta = betay(i,j,k) - end if + beta = betay(i,j,k) call eigenvalues(GRHydro_eos_handle, & rho (i,j,k), & @@ -165,11 +157,7 @@ subroutine GRHydro_FSAlpha(CCTK_ARGUMENTS) gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),& gyz(i,j,k),gzz(i,j,k)) - if (shift_state .eq. 0) then - beta = 0.d0 - else - beta = betaz(i,j,k) - end if + beta = betaz(i,j,k) call eigenvalues(GRHydro_eos_handle, & rho (i,j,k), & @@ -310,11 +298,8 @@ subroutine GRHydro_SplitFlux(CCTK_ARGUMENTS) do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil + 2 do j = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil + 2 - if (shift_state .eq. 0) then - dummy = 0.d0 - else - dummy = betax(:,j,k) - end if + dummy = betax(:,j,k) + do i = 1, cctk_lsh(1) det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\ gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) @@ -347,11 +332,8 @@ subroutine GRHydro_SplitFlux(CCTK_ARGUMENTS) do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil + 2 do i = GRHydro_stencil, cctk_lsh(1) - GRHydro_stencil + 2 - if (shift_state .eq. 0) then - dummy = 0.d0 - else - dummy = betay(i,:,k) - end if + dummy = betay(i,:,k) + do j = 1, cctk_lsh(2) det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\ gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) @@ -384,12 +366,8 @@ subroutine GRHydro_SplitFlux(CCTK_ARGUMENTS) do j = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil + 2 do i = GRHydro_stencil, cctk_lsh(1) - GRHydro_stencil + 2 - if (shift_state .eq. 0) then - dummy = 0.d0 - else - dummy = betaz(i,j,:) - end if - + dummy = betaz(i,j,:) + do k = 1, cctk_lsh(3) det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\ gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) |