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