[Bug fortran/31120] ICE with integer_exponentiation_1.f90 and -ffast-math

["dominiq at lps dot ens dot fr" <gcc-bugzilla at gcc dot gnu dot org>]

------- Comment #4 from dominiq at lps dot ens dot fr  2007-03-14 13:58 -------
Subject: Re:  ICE with integer_exponentiation_1.f90 and -ffast-math

> And Dominique, I would appreciate if you could test the patch on ppc-darwin7.

I'll do it tonight, but before could you test the following code:

integer(4)  :: i
integer(8)  :: j
real(4)     :: a
real(8)     :: b, c, d
a = 1.0 + epsilon(a)
b = 1.0 + epsilon(b)
print *, a**huge(i), a**(-huge(i)), a**(-huge(i)-1)
print '(3(1PG30.17))', b**huge(i), b**(-huge(i)), b**(-huge(i)-1)
c = b**huge(i)
d= b**(-huge(i))
print '(3(1PG30.17))', c*d, d/b, exp(-huge(i)*log(b))
d=1/b
do i = 1, 31
  d = d*d
end do
print '(2(1PG30.17))', d, b*d
!print *, b**huge(j), b**(-huge(j)), b**(-huge(j)-1)
end

If I uncomment the last print, I get:

[karma] f90/bug% gfc test_pow.f90
Out of stack space.
Try running 'limit stacksize unlimited' in the shell to raise its limit.

Is this also fixed by your patch?


-- 


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=31120