1! PR middle-end/82004
2! { dg-do run }
3! { dg-options "-Ofast" }
4
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))
11
12  d = c(1)
13  e = c(b)
14  f = (log10(e)-log10(d))/real(a)
15  g = log10(d) - f
16  h = 10**(g)
17  if (h.ne.j) stop 1
18end
19