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

Modified: /trunk/src/ GRHydro_CalcUpdate.F90

Log: GRHydro: fix OpenMP data race condition in UpdateCalculation for constraint transport

File Changes:

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

File [modified]: GRHydro_CalcUpdate.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_CalcUpdate.F90 2012-05-13 19:12:48 UTC (rev 325) +++ trunk/src/GRHydro_CalcUpdate.F90 2012-05-13 19:13:31 UTC (rev 326) @@ -47,7 +47,7 @@

if (use_weighted_fluxes == 0) then

- !$OMP PARALLEL DO PRIVATE(i,j,k,itracer,alp_l,alp_r,Bvec_l,Bvec_r) + !$OMP PARALLEL DO PRIVATE(i,j,k,itracer,alp_l,alp_r,alp_tmp,Bvec_l,Bvec_r) do k = GRHydro_stencil + 1 - transport_constraints, cctk_lsh(3) - GRHydro_stencil ! we need to compute Evec on all faces/edges where the fluxes are defined do j = GRHydro_stencil + 1 - transport_constraints, cctk_lsh(2) - GRHydro_stencil do i = GRHydro_stencil + 1 - transport_constraints, cctk_lsh(1) - GRHydro_stencil