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

[roland.haas at physics.gatech.edu]

User: rhaas
Date: 2012/05/13 02:13 PM

Modified:
 /trunk/src/
  GRHydro_CalcUpdate.F90

Log:
 GRHydro: OpenMP parallelize constraint transport loop

File Changes:

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

File [modified]: GRHydro_CalcUpdate.F90
Delta lines: +2 -0
===================================================================
--- trunk/src/GRHydro_CalcUpdate.F90	2012-05-13 19:13:31 UTC (rev 326)
+++ trunk/src/GRHydro_CalcUpdate.F90	2012-05-13 19:13:54 UTC (rev 327)
@@ -293,6 +293,7 @@
   if (transport_constraints.ne.0 .and. flux_direction.eq.1) then ! HACK: x direction is last
     ! FIXME: I think one could wrap all of this into a single do loop and remove the
     !        Evec storage
+    !$OMP PARALLEL DO PRIVATE(i,j,k)
     do k = GRHydro_stencil + 1, cctk_lsh(3) - GRHydro_stencil
       do j = GRHydro_stencil + 1, cctk_lsh(2) - GRHydro_stencil
         do i = GRHydro_stencil + 1, cctk_lsh(1) - GRHydro_stencil
@@ -311,6 +312,7 @@
         enddo
       enddo
     enddo
+    !$OMP END PARALLEL DO
   end if
   
   return