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))