25 lines
687 B
Fortran
25 lines
687 B
Fortran
! { 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)) &
|
|
call abort()
|
|
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
|