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

[cott at tapir.caltech.edu]

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