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

[roland.haas at physics.gatech.edu]

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