[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 198)

[schnetter at cct.lsu.edu]

User: eschnett
Date: 2010/12/22 08:47 PM

Modified:
 /trunk/src/
  GRHydro_Prim2Con.F90

Log:
 Loop in most prim2con routines over all grid points, not only the
 evolved ones. In this way, the conserved variables are defined on all
 grid points.

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Prim2Con.F90
Delta lines: +18 -18
===================================================================
--- trunk/src/GRHydro_Prim2Con.F90	2010-12-22 06:52:38 UTC (rev 197)
+++ trunk/src/GRHydro_Prim2Con.F90	2010-12-23 02:47:41 UTC (rev 198)
@@ -318,9 +318,9 @@
 
   if(evolve_temper.ne.1) then
      !$OMP PARALLEL DO PRIVATE(i, j, det)
-     do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1
-        do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1
-           do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1
+     do k = 1,cctk_lsh(3)
+        do j = 1,cctk_lsh(2)
+           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))
@@ -337,9 +337,9 @@
      !$OMP END PARALLEL DO
   else
      !$OMP PARALLEL DO PRIVATE(i, j, det)
-     do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1
-        do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1
-           do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1
+     do k = 1,cctk_lsh(3)
+        do j = 1,cctk_lsh(2)
+           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))
@@ -360,9 +360,9 @@
 
   if(evolve_Y_e.ne.0) then
      !$OMP PARALLEL DO PRIVATE(i, j)
-     do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1
-        do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1
-           do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1
+     do k = 1,cctk_lsh(3)
+        do j = 1,cctk_lsh(2)
+           do i = 1,cctk_lsh(1)
               Y_e_con(i,j,k) = Y_e(i,j,k) * dens(i,j,k)
            enddo
         enddo
@@ -551,9 +551,9 @@
   CCTK_REAL :: det
   
   !$OMP PARALLEL DO PRIVATE(i, j, det)
-  do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1
-    do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1
-      do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1
+     do k = 1,cctk_lsh(3)
+        do j = 1,cctk_lsh(2)
+           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))
@@ -571,9 +571,9 @@
 
   if(evolve_Y_e.ne.0) then
      !$OMP PARALLEL DO PRIVATE(i, j)
-     do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1
-        do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1
-           do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1
+     do k = 1,cctk_lsh(3)
+        do j = 1,cctk_lsh(2)
+           do i = 1,cctk_lsh(1)
               Y_e_con(i,j,k) = Y_e(i,j,k) * dens(i,j,k)
            enddo
         enddo
@@ -764,9 +764,9 @@
   
   ierr = EOS_SetGFs(cctkGH, EOS_Prim2ConCellsCall)
 
-  do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil + 1
-    do j = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil + 1
-      do i = GRHydro_stencil, cctk_lsh(1) - GRHydro_stencil + 1
+  do k = 1,cctk_lsh(3)
+    do j = 1,cctk_lsh(2)
+      do i = 1,cctk_lsh(1)
         
         gxxpt = gxx(i,j,k)
         gxypt = gxy(i,j,k)