diff options
Diffstat (limited to 'src')
-rw-r--r-- | src/GRHydro_Con2Prim.F90 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/GRHydro_Con2Prim.F90 b/src/GRHydro_Con2Prim.F90 index 958d033..9f5e511 100644 --- a/src/GRHydro_Con2Prim.F90 +++ b/src/GRHydro_Con2Prim.F90 @@ -556,7 +556,7 @@ subroutine Con2Prim_pt(handle, dens, sx, sy, sz, tau, rho, velx, vely, & tmp = (utau + pnew + udens)**2 - s2 plow = max(pmin, sqrt(s2) - utau - udens) - if (pnew .lt. plow .or. tmp .le. 0.0d0 .or. mustbisect .eq. .true.) then + if (pnew .lt. plow .or. tmp .le. 0.0d0 .or. mustbisect) then ! Ok, Newton-Raphson ended up finding something unphysical. ! Let's try to find our root via bisection (which converges slower but is more robust) @@ -1039,7 +1039,7 @@ subroutine Con2Prim_pt_hot(cctk_iteration, ii,jj,kk,handle, dens, & tmp = (utau + pnew + udens)**2 - s2 plow = max(pminl, sqrt(s2) - utau - udens) - if (pnew .lt. plow .or. tmp .le. 0.0d0 .or. mustbisect .eq. .true.) then + if (pnew .lt. plow .or. tmp .le. 0.0d0 .or. mustbisect) then ! Ok, Newton-Raphson ended up finding something unphysical. ! Let's try to find our root via bisection (which converges slower but is more robust) |