diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-10-27 04:16:35 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-10-27 04:16:35 +0000 |
commit | e937d71b0b60c2e1d59332d2dfa650b840f06ebc (patch) | |
tree | 9341dfd1e39efe4b899dc948393540a96ddf9745 /src/GRHydro_Con2PrimHot.F90 | |
parent | 777d3d5417c9833eed108f298698ac934b7dd65e (diff) |
GRHydro: replace Fortan pointers with Cray pointers
Intel 12 seems to like these better. Otherwise GRHydro_Con2PrimM.F90
takes very long to compile.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@580 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Con2PrimHot.F90')
-rw-r--r-- | src/GRHydro_Con2PrimHot.F90 | 40 |
1 files changed, 18 insertions, 22 deletions
diff --git a/src/GRHydro_Con2PrimHot.F90 b/src/GRHydro_Con2PrimHot.F90 index 07180dc..4cd07d1 100644 --- a/src/GRHydro_Con2PrimHot.F90 +++ b/src/GRHydro_Con2PrimHot.F90 @@ -9,12 +9,6 @@ subroutine Conservative2PrimitiveHot(CCTK_ARGUMENTS) implicit none - - ! save memory when MP is not used - ! TARGET as to be before DECLARE_CCTK_ARGUMENTS for gcc 4.1 - TARGET gaa, gab, gac, gbb, gbc, gcc - TARGET gxx, gxy, gxz, gyy, gyz, gzz - TARGET lvel, vel DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS @@ -32,8 +26,10 @@ subroutine Conservative2PrimitiveHot(CCTK_ARGUMENTS) ! save memory when MP is not used CCTK_INT :: GRHydro_UseGeneralCoordinates - CCTK_REAL, DIMENSION(:,:,:), POINTER :: g11, g12, g13, g22, g23, g33 - CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup + CCTK_REAL, DIMENSION(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)) :: g11, g12, g13, g22, g23, g33 + pointer (pg11,g11), (pg12,g12), (pg13,g13), (pg22,g22), (pg23,g23), (pg33,g33) + CCTK_REAL, DIMENSION(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3),3) :: vup + pointer (pvup,vup) ! begin EOS Omni vars integer :: n,keytemp,anyerr,keyerr(1), reset_to_atmo @@ -47,21 +43,21 @@ subroutine Conservative2PrimitiveHot(CCTK_ARGUMENTS) ! save memory when MP is not used if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then - g11 => gaa - g12 => gab - g13 => gac - g22 => gbb - g23 => gbc - g33 => gcc - vup => lvel + pg11 = loc(gaa) + pg12 = loc(gab) + pg13 = loc(gac) + pg22 = loc(gbb) + pg23 = loc(gbc) + pg33 = loc(gcc) + pvup = loc(lvel) else - g11 => gxx - g12 => gxy - g13 => gxz - g22 => gyy - g23 => gyz - g33 => gzz - vup => vel + pg11 = loc(gxx) + pg12 = loc(gxy) + pg13 = loc(gxz) + pg22 = loc(gyy) + pg23 = loc(gyz) + pg33 = loc(gzz) + pvup = loc(vel) end if nx = cctk_lsh(1) |