diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-03-28 01:46:34 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-03-28 01:46:34 +0000 |
| commit | b0cdccb34f58a51ee5357cbcd2174b847952244c (patch) | |
| tree | eaf258046d724987f4503112a798701859675fab | |
| parent | c0f1312a52b05b79e3a4d57277b25def6702b062 (diff) | |
GRHydro: reduce size of critical section in eigenbalues routine
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@494 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_Eigenproblem.F90 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/GRHydro_Eigenproblem.F90 b/src/GRHydro_Eigenproblem.F90 index 1a69403..7fe8bdb 100644 --- a/src/GRHydro_Eigenproblem.F90 +++ b/src/GRHydro_Eigenproblem.F90 @@ -86,16 +86,16 @@ subroutine eigenvalues(handle,rho,velx,vely,velz,eps, & (1.0d0 + eps + press/rho) if(cs2.lt.0.0d0) then - !$OMP CRITICAL if (abs(cs2) .gt. 1.0d-4) then + !$OMP CRITICAL write(warnline,'(a60,6g16.7)') 'abs(cs2), rho, dpdrho, press*dpdeps/rho**2, eps, press/rho: ', abs(cs2), rho, dpdrho, press * dpdeps / (rho**2), eps, press/rho call CCTK_WARN(1,warnline) call CCTK_WARN(1,"cs2 < 0! Check speed of sound calculation!") + !$OMP END CRITICAL cs2 = 0.0d0 else cs2 = 0.0d0 endif - !$OMP END CRITICAL endif |