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