2 ! PR tree-optimization/30092
3 ! This caused once an ICE due to internal tree changes
6 integer, parameter :: N
= 30
7 real, dimension(N
) :: rho
, pre
, cs
12 call EOS(N
, rho
, pre
, cs
, gamma
)
13 if (abs(CS(1) - sqrt(gamma
*pre(1)/rho(1))) > epsilon(cs
)) &
16 SUBROUTINE EOS(NODES
, DENS
, PRES
, CS
, CGAMMA
)
20 REAL, DIMENSION(NODES
) :: DENS
, PRES
, CS
21 REAL, PARAMETER :: RGAS
= 8.314
22 CS(:NODES
) = SQRT(CGAMMA
*PRES(:NODES
)/DENS(:NODES
))