diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-03-28 01:46:46 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-03-28 01:46:46 +0000 |
commit | 4715ab5ae5a87ed2ab3ee23673536a98b040e2f2 (patch) | |
tree | 5415b193b9fc1c9c72cf47f85221f40b03263fdf /src/GRHydro_HLLEM.F90 | |
parent | 91171b59c4f68833e46319bf919ce6fb2cf5be77 (diff) |
GRHydro: change CCTK_WARN(0,...) to CCTK_ERROR
this should tell the compiler that the function never returns and help
optimizing the code by indicating eg slow branches of if statements
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@497 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_HLLEM.F90')
-rw-r--r-- | src/GRHydro_HLLEM.F90 | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/src/GRHydro_HLLEM.F90 b/src/GRHydro_HLLEM.F90 index 780367c..2b025b4 100644 --- a/src/GRHydro_HLLEM.F90 +++ b/src/GRHydro_HLLEM.F90 @@ -119,7 +119,7 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) call SpaceMask_GetStateBits(trivial, "Hydro_RiemannProblemZ", & &"trivial") else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if ! constraint transport needs to be able to average fluxes in the directions @@ -229,7 +229,7 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) else if (flux_direction == 3) then avg_beta=avg_betaz else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if avg_alp = 0.5 * (alp(i,j,k) + alp(i+xoffset,j+yoffset,k+zoffset)) @@ -337,7 +337,7 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) endif else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if else !!! The end of this branch is right at the bottom of the routine @@ -349,7 +349,7 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) else if (flux_direction == 3) then usendh = uzzh else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if !!$ Calculate the jumps in the conserved variables @@ -541,7 +541,7 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) endif else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if @@ -779,7 +779,7 @@ subroutine GRHydro_HLLE_TracerM(CCTK_ARGUMENTS) call SpaceMask_GetStateBits(trivial, "Hydro_RiemannProblemZ", & &"trivial") else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil @@ -840,7 +840,7 @@ subroutine GRHydro_HLLE_TracerM(CCTK_ARGUMENTS) avg_beta = 0.5d0 * (beta3(i+xoffset,j+yoffset,k+zoffset) + & beta3(i,j,k)) else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if avg_alp = 0.5 * (alp(i,j,k) + alp(i+xoffset,j+yoffset,k+zoffset)) @@ -865,7 +865,7 @@ subroutine GRHydro_HLLE_TracerM(CCTK_ARGUMENTS) else if (flux_direction == 3) then usendh = uzzh else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if !!$ b^2 = (1-v^2)B^2+(B dot v)^2 @@ -929,7 +929,7 @@ subroutine GRHydro_HLLE_TracerM(CCTK_ARGUMENTS) fminus(:) = (velzminus(i+xoffset,j+yoffset,k+zoffset) - avg_beta / avg_alp) * & cons_tracerminus(i+xoffset,j+yoffset,k+zoffset,:) else - call CCTK_WARN(0, "Flux direction not x,y,z") + call CCTK_ERROR("Flux direction not x,y,z") end if !!$ Find minimum and maximum wavespeeds |