2 * Date: Fri
, 5 Mar
1999 00:35:44 -0500 (EST
)
3 * From
: Denes Molnar
<molnard@phys
.columbia
.edu
>
5 * Subject
: f771 gets fatal signal
6
6 * Content
-Type
: TEXT
/PLAIN
; charset
=US
-ASCII
7 * X
-UIDL
: 8d81e9cbdcc96209c6e9b298d966ba7f
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'
19 * I am running GNU Fortran
0.5.23 with GCC
2.8.1 (glibc1
).
24 * %%%%%%%%%%%%%%%%%%%%%%%%%
26 * %%%%%%%%%%%%%%%%%%%%%%%%%
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-----------------------------------------------------------------------
37 complex(kind
=8) HWUCI2
,HWULI2
,EPSI
,Y1
,Y2
,Z1
,Z2
,Z3
,Z4
38 DOUBLE PRECISION A
,B
,Y0
,ZERO
,ONE
,FOUR
,HALF
41 PARAMETER (ZERO
=0.D0
, ONE
=1.D0
, FOUR
= 4.D0
, HALF
=0.5D0
)
43 HWUCI2
=CMPLX
(ZERO
,ZERO
)
45 Y1
=HALF*
(ONE
+SQRT
(ONE
-FOUR*
(A
+EPSI
)/B
))
51 HWUCI2
=HWULI2
(Z1
)-HWULI2
(Z2
)+HWULI2
(Z3
)-HWULI2
(Z4
)
56 * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%