OSDN Git Service

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