[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 581)
schnetter at cct.lsu.edu
schnetter at cct.lsu.edu
Mon Oct 28 11:35:35 CDT 2013
User: eschnett
Date: 2013/10/28 11:35 AM
Modified:
/trunk/src/
GRHydro_HLLE.F90
Log:
Correct array indices in etaXYZ
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_HLLE.F90
Delta lines: +3 -3
===================================================================
--- trunk/src/GRHydro_HLLE.F90 2013-10-27 04:16:35 UTC (rev 580)
+++ trunk/src/GRHydro_HLLE.F90 2013-10-28 16:35:34 UTC (rev 581)
@@ -33,15 +33,15 @@
! eta across face in x-direction
#define etaX(i,j,k, v, c) \
- (0.5d0*(abs(v(i+1,j,k,0)-v(i,j,k,0)) + abs(c(i+1,j,k)-c(i,j,k))))
+ (0.5d0*(abs(v((i)+1,j,k,1)-v(i,j,k,1)) + abs(c((i)+1,j,k)-c(i,j,k))))
! eta across face in x-direction
#define etaY(i,j,k, v, c) \
- (0.5d0*(abs(v(i,j+1,k,1)-v(i,j,k,1)) + abs(c(i,j+1,k)-c(i,j,k))))
+ (0.5d0*(abs(v(i,(j)+1,k,2)-v(i,j,k,2)) + abs(c(i,(j)+1,k)-c(i,j,k))))
! eta across face in x-direction
#define etaZ(i,j,k, v, c) \
- (0.5d0*(abs(v(i,j,k+1,2)-v(i,j,k,2)) + abs(c(i,j,k+1)-c(i,j,k))))
+ (0.5d0*(abs(v(i,j,(k)+1,3)-v(i,j,k,3)) + abs(c(i,j,(k)+1)-c(i,j,k))))
More information about the Commits
mailing list