array_2.f90

来自「用于进行gcc测试」· F90 代码 · 共 25 行

! { 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