User: rhaas Date: 2012/08/27 02:19 PM

Modified: /trunk/src/ GRHydro_HLLEM.F90

Log: GRHydro: Adding back OpenMP support to HLLEM.F90

From: Philipp Moesta pmoesta@tapir.caltech.edu

File Changes:

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

File [modified]: GRHydro_HLLEM.F90 Delta lines: +12 -0 =================================================================== --- trunk/src/GRHydro_HLLEM.F90 2012-08-27 19:19:29 UTC (rev 413) +++ trunk/src/GRHydro_HLLEM.F90 2012-08-27 19:19:32 UTC (rev 414) @@ -123,6 +123,17 @@

! constraint transport needs to be able to average fluxes in the directions ! other that flux_direction + + !$OMP PARALLEL DO PRIVATE(k,j,i,f1,lamminus,lamplus,cons_p,cons_m,fplus,fminus,qdiff,psidcf,psidcp,psidcm,prim_p,prim_m,& + !$OMP avg_betax,avg_betay,avg_betaz,avg_beta,avg_alp,& + !$OMP gxxh,gxyh,gxzh,gyyh,gyzh,gzzh,avg_det,sdet,uxxh,uxyh,uxzh,uyyh,uyzh,uzzh,& + !$OMP vxtp,vxtm,vytp,vytm,vztp,vztm,& + !$OMP velxlowp,velxlowm,Bvecxlowp,Bvecxlowm,& + !$OMP velylowp,velylowm,Bvecylowp,Bvecylowm,& + !$OMP velzlowp,velzlowm,Bveczlowp,Bveczlowm,& + !$OMP bdotvp,bdotvm,b2p,b2m,v2p,v2m,wp,wm,& + !$OMP bxlowp,bxlowm,bylowp,bylowm,bzlowp,bzlowm,& + !$OMP rhoenth_p,rhoenth_m,ab0p,ab0m,vA2p,vA2m,pressstarp,pressstarm,usendh,psidcdiff,psidcfp,psidcfm,charmin,charmax,chartop,charpm,m) do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil + transport_constraints*(1-zoffset) do j = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil + transport_constraints*(1-yoffset) do i = GRHydro_stencil, cctk_lsh(1) - GRHydro_stencil + transport_constraints*(1-xoffset) @@ -608,6 +619,7 @@ end do end do end do + !$OMP END PARALLEL DO #undef faulty_gxx #undef faulty_gxy #undef faulty_gxz