57 lines
1.6 KiB
Fortran
57 lines
1.6 KiB
Fortran
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
|
|
*
|
|
* %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|