blob: 982daa75a915f7fce42ad81bf31baaed1de6939f [file] [log] [blame]
! { dg-do run }
! PR tree-optimization/30092
! This caused once an ICE due to internal tree changes
program test
implicit none
integer, parameter :: N = 30
real, dimension(N) :: rho, pre, cs
real :: gamma
gamma = 2.1314
rho = 5.0
pre = 3.0
call EOS(N, rho, pre, cs, gamma)
if (abs(CS(1) - sqrt(gamma*pre(1)/rho(1))) > epsilon(cs)) &
STOP 1
contains
SUBROUTINE EOS(NODES, DENS, PRES, CS, CGAMMA)
IMPLICIT NONE
INTEGER NODES
REAL CGAMMA
REAL, DIMENSION(NODES) :: DENS, PRES, CS
REAL, PARAMETER :: RGAS = 8.314
CS(:NODES) = SQRT(CGAMMA*PRES(:NODES)/DENS(:NODES))
END SUBROUTINE EOS
end program test