From b0cdccb34f58a51ee5357cbcd2174b847952244c Mon Sep 17 00:00:00 2001 From: rhaas Date: Thu, 28 Mar 2013 01:46:34 +0000 Subject: GRHydro: reduce size of critical section in eigenbalues routine From: Roland Haas git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@494 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Eigenproblem.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'src/GRHydro_Eigenproblem.F90') 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 -- cgit v1.2.3