This is the mail archive of the
gcc-bugs@gcc.gnu.org
mailing list for the GCC project.
Possible g77 bug
- To: egcs-bugs at cygnus dot com
- Subject: Possible g77 bug
- From: James Sternberg <jsternbe at atomic1 dot phys dot utk dot edu>
- Date: Tue, 19 May 1998 10:34:11 -0400 (EDT)
One of my co-workers gave me the following piece of code which he was
trying to compile on our alpha-linux box using egcs-1.0.3. When we try to
compile it, the compilation process never reaches the end. No errors or
warnings are given. I have tried compiling with egcs-1.0.2 on the same
machine with the same result. It *does* compile properly with egcs-1.0.2
on an MIPS machine.
If someone else is using egcs on alpha-linux, then perhaps they can see if
they get the same result. If it is just something screwed up on my end, I
am sorry for the inconvenience.
-James Sternberg
----------------------CUT HERE-------------------------------
Complex*16 function Darwin (a,x)
Complex*16 a,x
Complex*16 x2, Xb, Xb3, Gamma
Complex*16 cbsqrt
external Gamma
Real*8 tpif
parameter (tpif = 1.583233487086160d0 ) ! (2\pi)^(1/4)
Real*8 d91, d92, d93, d94, d121
parameter (d91=-7.d0/5760.d0, d92 = -7.d0/320.d0)
parameter (d93 = -49.d0/320.d0, d94 = 31.d0/12.d0)
parameter (d121 = 153.d0/8.d0)
x2 = x*x
Xb = cdsqrt(x2+4.d0*a)
Xb3=1.d0/Xb**3
Darwin = tpif/cbsqrt(Xb*Gamma(a+0.5d0)) *
* cdexp (-0.25d0*x*Xb - a*cdlog((x+Xb)/(2.d0*cdsqrt(a))) +
+ (-(x2/48.d0+0.5d0*a)*x/a +
+ (0.75*x2-2.d0*a +
+ (-(((((d91*x2+d92*a)*x2+d93*a*a)*x2+d94*a**3)*x2-
-19.d0*a**4)*x)/
/ a**3 +
+ ((d121*x2-186.d0*a)*x2+80.d0*a*a)*Xb3)
& *Xb3)*Xb3)*Xb3)
return
end