diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2014-04-15 19:48:49 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2014-04-15 19:48:49 +0000 |
| commit | c3651bf964ce2be13f8cac49af7f646240787f8e (patch) | |
| tree | 65f1defc104b8234a8094445d95d99a9712c5a9f | |
| parent | 265d03d0aab057f5305949dc2ebcc058a982860e (diff) | |
GRHydro: correct vel slicing in optimized C++ PPM reconstruction
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@611 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_PPMReconstruct_drv_opt.cc | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/GRHydro_PPMReconstruct_drv_opt.cc b/src/GRHydro_PPMReconstruct_drv_opt.cc index 7629019..3f5610f 100644 --- a/src/GRHydro_PPMReconstruct_drv_opt.cc +++ b/src/GRHydro_PPMReconstruct_drv_opt.cc @@ -281,8 +281,8 @@ void CCTK_FNAME(GRHydro_PPMReconstruct_drv_opt)(cGH const * const restrict & res oned_slice(nx,nstart,1,rho,rho1d); oned_slice(nx,nstart, 1,vup,velx1d); - oned_slice(nx,nstart+nxyz, 1,vup,velx1d); - oned_slice(nx,nstart+2*nxyz,1,vup,velx1d); + oned_slice(nx,nstart+nxyz, 1,vup,vely1d); + oned_slice(nx,nstart+2*nxyz,1,vup,velz1d); oned_slice(nx,nstart,1,eps,eps1d); oned_slice(nx,nstart,1,press,press1d); @@ -371,8 +371,8 @@ void CCTK_FNAME(GRHydro_PPMReconstruct_drv_opt)(cGH const * const restrict & res oned_slice(ny,nstart,nx,rho,rho1d); oned_slice(ny,nstart, nx,vup,velx1d); - oned_slice(ny,nstart+nxyz, nx,vup,velx1d); - oned_slice(ny,nstart+2*nxyz,nx,vup,velx1d); + oned_slice(ny,nstart+nxyz, nx,vup,vely1d); + oned_slice(ny,nstart+2*nxyz,nx,vup,velz1d); oned_slice(ny,nstart,nx,eps,eps1d); oned_slice(ny,nstart,nx,press,press1d); @@ -463,8 +463,8 @@ void CCTK_FNAME(GRHydro_PPMReconstruct_drv_opt)(cGH const * const restrict & res oned_slice(nz,nstart,nxy,rho,rho1d); oned_slice(nz,nstart, nxy,vup,velx1d); - oned_slice(nz,nstart+nxyz, nxy,vup,velx1d); - oned_slice(nz,nstart+2*nxyz,nxy,vup,velx1d); + oned_slice(nz,nstart+nxyz, nxy,vup,vely1d); + oned_slice(nz,nstart+2*nxyz,nxy,vup,velz1d); oned_slice(nz,nstart,nxy,eps,eps1d); oned_slice(nz,nstart,nxy,press,press1d); |