Subversion Repositories openrisc

c { dg-do compile }

* Date: Fri, 5 Mar 1999 00:35:44 -0500 (EST)

* From: Denes Molnar <molnard@phys.columbia.edu>

* To: fortran@gnu.org

* Subject: f771 gets fatal signal 6

* Content-Type: TEXT/PLAIN; charset=US-ASCII

* X-UIDL: 8d81e9cbdcc96209c6e9b298d966ba7f

*

* Hi,

*

*

* Comiling object from the source code below WORKS FINE with

*   'g77 -o hwuci2 -c hwuci2.F'

* but FAILS with fatal signal 6

*   'g77 -o hwuci2 -O -c hwuci2.F'

*

* Any explanations?

*

* I am running GNU Fortran 0.5.23 with GCC 2.8.1 (glibc1).

*

*

* Denes Molnar

*

* %%%%%%%%%%%%%%%%%%%%%%%%%

* %the source:

* %%%%%%%%%%%%%%%%%%%%%%%%%

*

CDECK  ID>, HWUCI2.

*CMZ :-        -23/08/94  13.22.29  by  Mike Seymour

*-- Author :    Ulrich Baur & Nigel Glover, adapted by Ian Knowles

C-----------------------------------------------------------------------

      FUNCTION HWUCI2(A,B,Y0)

C-----------------------------------------------------------------------

C     Integral  LOG(A-EPSI-BY(1-Y))/(Y-Y0)

C-----------------------------------------------------------------------

      IMPLICIT NONE

      complex(kind=8) HWUCI2,HWULI2,EPSI,Y1,Y2,Z1,Z2,Z3,Z4

      DOUBLE PRECISION A,B,Y0,ZERO,ONE,FOUR,HALF

      EXTERNAL HWULI2

      COMMON/SMALL/EPSI

      PARAMETER (ZERO=0.D0, ONE =1.D0, FOUR= 4.D0, HALF=0.5D0)

      IF(B.EQ.ZERO)THEN

         HWUCI2=CMPLX(ZERO,ZERO)

      ELSE

         Y1=HALF*(ONE+SQRT(ONE-FOUR*(A+EPSI)/B))

         Y2=ONE-Y1

         Z1=Y0/(Y0-Y1)

         Z2=(Y0-ONE)/(Y0-Y1)

         Z3=Y0/(Y0-Y2)

         Z4=(Y0-ONE)/(Y0-Y2)

         HWUCI2=HWULI2(Z1)-HWULI2(Z2)+HWULI2(Z3)-HWULI2(Z4)

      ENDIF

      RETURN

      END

*

* %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%