3 ! { dg-options "-Ofast" }
5 integer, parameter :: r8
= selected_real_kind(13), i4
= kind(1)
6 integer (i4
), parameter :: a
= 400, b
= 2
7 real (r8
), parameter, dimension(b
) :: c
= (/ .001_r8
, 10.00_r8
/)
8 real (r8
) :: d
, e
, f
, g
, h
9 real (r8
), parameter :: j
&
10 = 10**(log10(c(1))-(log10(c(b
))-log10(c(1)))/real(a
))
14 f
= (log10(e
)-log10(d
))/real(a
)