Possible g77 bug

[James Sternberg <jsternbe at atomic1 dot phys dot utk dot edu>]

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