gcc/testsuite/gfortran.dg/g77/19990305-0.f

   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
   8 *
   9 * Hi,
  10 *
  11 *
  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'
  16 *
  17 * Any explanations?
  18 *
  19 * I am running GNU Fortran 0.5.23 with GCC 2.8.1 (glibc1).
  20 *
  21 *
  22 * Denes Molnar
  23 *
  24 * %%%%%%%%%%%%%%%%%%%%%%%%%
  25 * %the source:
  26 * %%%%%%%%%%%%%%%%%%%%%%%%%
  27 *
  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
  55 *
  56 * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%