Remove old autovect-branch by moving to "dead" directory.
[official-gcc.git] / old-autovect-branch / gcc / testsuite / gfortran.dg / g77 / 19990305-0.f
blob056d2b7a3a12d5d7619a6ce93d85caec26aef480
1 c { dg-do compile }
2 * Date: Fri, 5 Mar 1999 00:35:44 -0500 (EST)
3 * From: Denes Molnar <molnard@phys.columbia.edu>
4 * To: fortran@gnu.org
5 * Subject: f771 gets fatal signal 6
6 * Content-Type: TEXT/PLAIN; charset=US-ASCII
7 * X-UIDL: 8d81e9cbdcc96209c6e9b298d966ba7f
9 * Hi,
12 * Comiling object from the source code below WORKS FINE with
13 * 'g77 -o hwuci2 -c hwuci2.F'
14 * but FAILS with fatal signal 6
15 * 'g77 -o hwuci2 -O -c hwuci2.F'
17 * Any explanations?
19 * I am running GNU Fortran 0.5.23 with GCC 2.8.1 (glibc1).
22 * Denes Molnar
24 * %%%%%%%%%%%%%%%%%%%%%%%%%
25 * %the source:
26 * %%%%%%%%%%%%%%%%%%%%%%%%%
28 CDECK ID>, HWUCI2.
29 *CMZ :- -23/08/94 13.22.29 by Mike Seymour
30 *-- Author : Ulrich Baur & Nigel Glover, adapted by Ian Knowles
31 C-----------------------------------------------------------------------
32 FUNCTION HWUCI2(A,B,Y0)
33 C-----------------------------------------------------------------------
34 C Integral LOG(A-EPSI-BY(1-Y))/(Y-Y0)
35 C-----------------------------------------------------------------------
36 IMPLICIT NONE
37 complex(kind=8) HWUCI2,HWULI2,EPSI,Y1,Y2,Z1,Z2,Z3,Z4
38 DOUBLE PRECISION A,B,Y0,ZERO,ONE,FOUR,HALF
39 EXTERNAL HWULI2
40 COMMON/SMALL/EPSI
41 PARAMETER (ZERO=0.D0, ONE =1.D0, FOUR= 4.D0, HALF=0.5D0)
42 IF(B.EQ.ZERO)THEN
43 HWUCI2=CMPLX(ZERO,ZERO)
44 ELSE
45 Y1=HALF*(ONE+SQRT(ONE-FOUR*(A+EPSI)/B))
46 Y2=ONE-Y1
47 Z1=Y0/(Y0-Y1)
48 Z2=(Y0-ONE)/(Y0-Y1)
49 Z3=Y0/(Y0-Y2)
50 Z4=(Y0-ONE)/(Y0-Y2)
51 HWUCI2=HWULI2(Z1)-HWULI2(Z2)+HWULI2(Z3)-HWULI2(Z4)
52 ENDIF
53 RETURN
54 END
56 * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%