GRHydro: add OpenMP parallelization to IntialAtmosphereReset

From: Roland Haas <roland.haas@physics.gatech.edu>

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@418 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 7 insertions, 0 deletions

diff --git a/src/GRHydro_UpdateMask.F90 b/src/GRHydro_UpdateMask.F90
index b8c3208..b5a84a5 100644
--- a/src/GRHydro_UpdateMask.F90
+++ b/src/GRHydro_UpdateMask.F90
@@ -201,6 +201,7 @@ subroutine GRHydro_InitAtmosMask(CCTK_ARGUMENTS)
 
   CCTK_INT :: i,j,k
 
+  !$OMP PARALLEL DO PRIVATE(i,j,k)
   do k = 1, cctk_lsh(3)
     do j = 1, cctk_lsh(2)
       do i = 1, cctk_lsh(1)
@@ -211,6 +212,7 @@ subroutine GRHydro_InitAtmosMask(CCTK_ARGUMENTS)
       end do
     end do
   end do
+  !$OMP END PARALLEL DO
 
   call CCTK_INFO("Setting up the atmosphere mask: points with rho<rho_min are set to be atmosphere")
   
@@ -445,6 +447,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS)
     rho_min = rho_min * initial_atmosphere_factor
   endif
 
+!$OMP PARALLEL DO PRIVATE(det,keytemp,i,j,k,anyerr,keyerr)
   do k = 1, cctk_lsh(3)
     do j = 1, cctk_lsh(2)
       do i = 1, cctk_lsh(1)
@@ -477,6 +480,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS)
                   temperature(i,j,k),y_e(i,j,k))
              y_e_con(i,j,k) = dens(i,j,k) * y_e(i,j,k)
           else
+             keytemp = 0
              call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,&
                   rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr)
              call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,&
@@ -518,6 +522,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS)
                   temperature_p(i,j,k),y_e_p(i,j,k))
              y_e_con_p(i,j,k) = dens_p(i,j,k) * y_e_p(i,j,k)
             else
+               keytemp = 0
                call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,&
                     rho_p(i,j,k),eps_p(i,j,k),xtemp,xye,press_p(i,j,k),keyerr,anyerr)
                call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,&
@@ -559,6 +564,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS)
                   temperature_p_p(i,j,k),y_e_p_p(i,j,k))
              y_e_con_p_p(i,j,k) = dens_p_p(i,j,k) * y_e_p_p(i,j,k)
             else
+               keytemp = 0
                call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,&
                     rho_p_p(i,j,k),eps_p_p(i,j,k),xtemp,xye,press_p_p(i,j,k),keyerr,anyerr)
                call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,&
@@ -576,6 +582,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS)
       end do
     end do
   end do
+!$OMP END PARALLEL DO
 
 end subroutine GRHydro_InitialAtmosphereReset