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

20 Dec 2010

User: cott
Date: 2010/12/21 12:53 AM
Modified:
 /trunk/src/
  GRHydro_Prim2Con.F90
Log:
 * improve thread safety
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_Prim2Con.F90
Delta lines: +8 -1
===================================================================

--- trunk/src/GRHydro_Prim2Con.F90	2010-12-21 06:52:27 UTC (rev 192)
+++ trunk/src/GRHydro_Prim2Con.F90	2010-12-21 06:53:01 UTC (rev 193)
@@ -49,6 +49,9 @@
   CCTK_REAL :: xtemp
if(evolve_temper.ne.1) then
+     !$OMP PARALLEL DO PRIVATE(i, j, avg_detl, avg_detr,&
+     !$OMP                      gxxl,gxyl,gxzl,gyyl,gyzl,gzzl, &
+     !$OMP                      gxxr,gxyr,gxzr,gyyr,gyzr,gzzr)
      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
@@ -88,7 +91,11 @@
            end do
         end do
      end do
+     !$OMP END PARALLEL DO
   else
+     !$OMP PARALLEL DO PRIVATE(i, j, avg_detl, avg_detr, xtemp,&
+     !$OMP                      gxxl,gxyl,gxzl,gyyl,gyzl,gzzl, &
+     !$OMP                      gxxr,gxyr,gxzr,gyyr,gyzr,gzzr)
      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
@@ -147,7 +154,7 @@
            end do
         end do
      end do
-
+     !$OMP END PARALLEL DO
   endif
end subroutine primitive2conservative

    