From 00055ac08315061e3fd7a85487896e6297dd7726 Mon Sep 17 00:00:00 2001 From: kwagner Date: Mon, 22 Feb 2010 12:37:13 +0000 Subject: [PATCH] Define some constants as double precision objects, in order to remove compiler warnings. git-svn-id: https://lxsd.femto-st.fr/svn/fvn@70 b657c933-2333-4658-acf2-d3c7c2708721 --- fvn_fnlib/beskes.f | 2 +- fvn_fnlib/dbskes.f | 2 +- fvn_fnlib/z8lgmc.f | 4 ++-- fvn_fnlib/z9lgmc.f | 24 ++++++++++++------------ fvn_fnlib/zasin.f | 4 ++-- fvn_fnlib/zatanh.f | 2 +- fvn_fnlib/zbeta.f | 2 +- fvn_fnlib/zcosh.f | 2 +- fvn_fnlib/zcot.f | 2 +- fvn_fnlib/zexprl.f | 2 +- fvn_fnlib/zlngam.f | 2 +- fvn_fnlib/zlnrel.f | 2 +- fvn_fnlib/zpsi.f | 32 ++++++++++++++++---------------- fvn_fnlib/zsinh.f | 2 +- fvn_fnlib/ztan.f | 2 +- fvn_fnlib/ztanh.f | 2 +- 16 files changed, 44 insertions(+), 44 deletions(-) diff --git a/fvn_fnlib/beskes.f b/fvn_fnlib/beskes.f index 973e515..e0c5ce6 100644 --- a/fvn_fnlib/beskes.f +++ b/fvn_fnlib/beskes.f @@ -1,6 +1,6 @@ subroutine beskes (xnu, x, nin, bke) c july 1980 edition. w. fullerton, c3, los alamos scientific lab. - dimension bke(1) + dimension bke(*) external r1mach data alnbig / 0. / c diff --git a/fvn_fnlib/dbskes.f b/fvn_fnlib/dbskes.f index f614d4a..26ae701 100644 --- a/fvn_fnlib/dbskes.f +++ b/fvn_fnlib/dbskes.f @@ -1,6 +1,6 @@ subroutine dbskes (xnu, x, nin, bke) c july 1980 edition. w. fullerton, c3, los alamos scientific lab. - double precision xnu, x, bke(1), bknu1, v, vincr, vend, alnbig, + double precision xnu, x, bke(*), bknu1, v, vincr, vend, alnbig, 1 d1mach, dlog external d1mach data alnbig / 0.d0 / diff --git a/fvn_fnlib/z8lgmc.f b/fvn_fnlib/z8lgmc.f index d429145..e2d426e 100644 --- a/fvn_fnlib/z8lgmc.f +++ b/fvn_fnlib/z8lgmc.f @@ -9,14 +9,14 @@ c is small. c when real(z) is negative, c8lgmc merely returns a correction which c may be wrong by a multiple of 2*pi*i. c - complex(kind(1.d0)) zin, z, corr, z9lgmc, z0lgmc, zexp, zlnrel + complex(kind(1.d0)) zin, z, corr, z9lgmc, z0lgmc, zlnrel real(kind(1.d0)) d1mach,pi,bound,sqeps,eps,x,y,absz,test integer nterm,n,i,irold,ir c complex(kind(1.d0)) tmp_arg external z0lgmc, z9lgmc, zlnrel, 1 d1mach data pi / 3.1415926535 8979324d0 / - data bound, sqeps, eps / 3*0.0 / + data bound, sqeps, eps / 3*0.0d0 / c if (bound.ne.0.0) go to 10 nterm = -0.30*log(d1mach(3)) diff --git a/fvn_fnlib/z9lgmc.f b/fvn_fnlib/z9lgmc.f index 608a76f..611a0aa 100644 --- a/fvn_fnlib/z9lgmc.f +++ b/fvn_fnlib/z9lgmc.f @@ -13,19 +13,19 @@ c external d1mach c dimension bern(11) - data bern( 1) / .08333333333 3333333e0 / - data bern( 2) / -.002777777777 7777778e0 / - data bern( 3) / .0007936507936 5079365e0 / - data bern( 4) / -.0005952380952 3809524e0 / - data bern( 5) / .0008417508417 5084175e0 / - data bern( 6) / -.001917526917 5269175e0 / - data bern( 7) / .006410256410 2564103e0 / - data bern( 8) / -.02955065359 4771242e0 / - data bern( 9) / .1796443723 6883057e0 / - data bern(10) / -1.392432216 9059011e0 / - data bern(11) / 13.40286404 4168392e0 / + data bern( 1) / .08333333333 3333333d0 / + data bern( 2) / -.002777777777 7777778d0 / + data bern( 3) / .0007936507936 5079365d0 / + data bern( 4) / -.0005952380952 3809524d0 / + data bern( 5) / .0008417508417 5084175d0 / + data bern( 6) / -.001917526917 5269175d0 / + data bern( 7) / .006410256410 2564103d0 / + data bern( 8) / -.02955065359 4771242d0 / + data bern( 9) / .1796443723 6883057d0 / + data bern(10) / -1.392432216 9059011d0 / + data bern(11) / 13.40286404 4168392d0 / c - data nterm, bound, xbig, xmax / 0, 3*0.0 / + data nterm, bound, xbig, xmax / 0d0, 3*0.0d0 / c if (nterm.ne.0) go to 10 c diff --git a/fvn_fnlib/zasin.f b/fvn_fnlib/zasin.f index f900a57..dcb6319 100644 --- a/fvn_fnlib/zasin.f +++ b/fvn_fnlib/zasin.f @@ -12,8 +12,8 @@ c data pi2 /1.5707963267 9489661923d0/ data pi /3.1415926535 8979324d0/ - data ci /(0.,1.)/ - data nterms, rmin / 0, 0.0 / + data ci /(0.d0,1.d0)/ + data nterms, rmin / 0, 0.0d0 / c if (nterms.ne.0) go to 10 c nterms = alog(eps)/alog(rmax) where rmax = 0.1 diff --git a/fvn_fnlib/zatanh.f b/fvn_fnlib/zatanh.f index 0e61676..f619bb3 100644 --- a/fvn_fnlib/zatanh.f +++ b/fvn_fnlib/zatanh.f @@ -3,7 +3,7 @@ c april 1977 version. w. fullerton, c3, los alamos scientific lab. complex(kind(1.d0)) z, ci, zatan external zatan - data ci /(0.,1.)/ + data ci /(0.d0,1.d0)/ c zatanh = -ci*zatan(ci*z) c diff --git a/fvn_fnlib/zbeta.f b/fvn_fnlib/zbeta.f index fc5d8e6..c94b878 100644 --- a/fvn_fnlib/zbeta.f +++ b/fvn_fnlib/zbeta.f @@ -4,7 +4,7 @@ c july 1977 edition. w. fullerton, c3, los alamos scientific lab. complex(kind(1.d0)) a, b, zgamma, zlbeta real(kind(1.d0)) xmax,xmin external zgamma, zlbeta - data xmax / 0.0 / + data xmax / 0.0d0 / c if (xmax.eq.0.0) call d9gaml (xmin, xmax) c diff --git a/fvn_fnlib/zcosh.f b/fvn_fnlib/zcosh.f index 5c7e7d7..5b13723 100644 --- a/fvn_fnlib/zcosh.f +++ b/fvn_fnlib/zcosh.f @@ -2,7 +2,7 @@ implicit none c april 1977 version. w. fullerton, c3, los alamos scientific lab. complex(kind(1.d0)) z, ci - data ci /(0.,1.)/ + data ci /(0.d0,1.d0)/ c zcosh = cos (ci*z) c diff --git a/fvn_fnlib/zcot.f b/fvn_fnlib/zcot.f index fe3ee50..b975df7 100644 --- a/fvn_fnlib/zcot.f +++ b/fvn_fnlib/zcot.f @@ -7,7 +7,7 @@ c march 1979 edition. w. fullerton, c3, los alamos scientific lab. real(kind(1.d0)) x,y,x2,y2,sn2x,den integer irold,irold2 external d1mach - data eps, xmax, ylarge, ybig, rmin, ymin / 5*0.0, 1.5 / + data eps, xmax, ylarge, ybig, rmin, ymin / 5*0.0d0, 1.5d0 / c if (eps.ne.0.0) go to 10 eps = d1mach(4) diff --git a/fvn_fnlib/zexprl.f b/fvn_fnlib/zexprl.f index dd69c4c..79c1121 100644 --- a/fvn_fnlib/zexprl.f +++ b/fvn_fnlib/zexprl.f @@ -12,7 +12,7 @@ c real(kind(1.d0)) sqeps,r,xn,xln,rbnd,d1mach,alneps external d1mach integer nterms,irold,irold2,i - data nterms, rbnd, sqeps / 0, 2*0.0 / + data nterms, rbnd, sqeps / 0d0, 2*0.0d0 / c if (nterms.ne.0) go to 10 alneps = log(d1mach(3)) diff --git a/fvn_fnlib/zlngam.f b/fvn_fnlib/zlngam.f index deff484..949dca8 100644 --- a/fvn_fnlib/zlngam.f +++ b/fvn_fnlib/zlngam.f @@ -11,7 +11,7 @@ c z in the vicinity of 1 and 2. data pi / 3.1415926535 8979324d0 / data sq2pil / 0.9189385332 0467274d0 / c - data bound, dxrel, rmax / 3*0.0 / + data bound, dxrel, rmax / 3*0.0d0 / c if (bound.ne.0.) go to 10 n = -0.30*log(d1mach(3)) diff --git a/fvn_fnlib/zlnrel.f b/fvn_fnlib/zlnrel.f index a838537..e0bf412 100644 --- a/fvn_fnlib/zlnrel.f +++ b/fvn_fnlib/zlnrel.f @@ -13,7 +13,7 @@ c complex(kind(1.d0)) z real(kind(1.d0)) dlnrel,d1mach,zarg,sqeps,x,rho external dlnrel, zarg, d1mach - data sqeps /0.0/ + data sqeps /0.0d0/ c if (sqeps.eq.0.) sqeps = sqrt (d1mach(4)) c diff --git a/fvn_fnlib/zpsi.f b/fvn_fnlib/zpsi.f index fcb3cd6..7f5837a 100644 --- a/fvn_fnlib/zpsi.f +++ b/fvn_fnlib/zpsi.f @@ -8,22 +8,22 @@ c may 1978 edition. w. fullerton, c3, los alamos scientific lab. integer nterm,ndx,n,i external zcot, d1mach c - data bern( 1) / .8333333333 3333333 e-1 / - data bern( 2) / -.8333333333 3333333 e-2 / - data bern( 3) / .3968253968 2539683 e-2 / - data bern( 4) / -.4166666666 6666667 e-2 / - data bern( 5) / .7575757575 7575758 e-2 / - data bern( 6) / -.2109279609 2796093 e-1 / - data bern( 7) / .8333333333 3333333 e-1 / - data bern( 8) / -.4432598039 2156863 e0 / - data bern( 9) / .3053954330 2701197 e1 / - data bern(10) / -.2645621212 1212121 e2 / - data bern(11) / .2814601449 2753623 e3 / - data bern(12) / -.3454885393 7728938 e4 / - data bern(13) / .5482758333 3333333 e5 / -c - data pi / 3.141592653 589793 e0 / - data nterm, bound, dxrel, rmin, rbig / 0, 4*0.0 / + data bern( 1) / .8333333333 3333333 d-1 / + data bern( 2) / -.8333333333 3333333 d-2 / + data bern( 3) / .3968253968 2539683 d-2 / + data bern( 4) / -.4166666666 6666667 d-2 / + data bern( 5) / .7575757575 7575758 d-2 / + data bern( 6) / -.2109279609 2796093 d-1 / + data bern( 7) / .8333333333 3333333 d-1 / + data bern( 8) / -.4432598039 2156863 d0 / + data bern( 9) / .3053954330 2701197 d1 / + data bern(10) / -.2645621212 1212121 d2 / + data bern(11) / .2814601449 2753623 d3 / + data bern(12) / -.3454885393 7728938 d4 / + data bern(13) / .5482758333 3333333 d5 / +c + data pi / 3.141592653 589793 d0 / + data nterm, bound, dxrel, rmin, rbig / 0d0, 4*0.0d0 / c if (nterm.ne.0) go to 10 nterm = -0.30*log(d1mach(3)) diff --git a/fvn_fnlib/zsinh.f b/fvn_fnlib/zsinh.f index 20833a7..3978e18 100644 --- a/fvn_fnlib/zsinh.f +++ b/fvn_fnlib/zsinh.f @@ -2,7 +2,7 @@ implicit none c april 1977 version. w. fullerton, c3, los alamos scientific lab. complex(kind(1.d0)) z, ci - data ci /(0.,1.)/ + data ci /(0.d0,1.d0)/ c zsinh = -ci*sin(ci*z) c diff --git a/fvn_fnlib/ztan.f b/fvn_fnlib/ztan.f index 7fb5d02..617a751 100644 --- a/fvn_fnlib/ztan.f +++ b/fvn_fnlib/ztan.f @@ -6,7 +6,7 @@ c march 1979 edition. w. fullerton, c3, los alamos scientific lab. real(kind(1.d0)) x,y,x2,y2 integer irold,irold2 external d1mach - data eps, xmax, ylarge, ybig, ymin / 4*0.0, 1.50 / + data eps, xmax, ylarge, ybig, ymin / 4*0.0d0, 1.50d0 / c if (eps.ne.0.0) go to 10 eps = d1mach(4) diff --git a/fvn_fnlib/ztanh.f b/fvn_fnlib/ztanh.f index 6cba088..75c89e1 100644 --- a/fvn_fnlib/ztanh.f +++ b/fvn_fnlib/ztanh.f @@ -2,7 +2,7 @@ c april 1977 version. w. fullerton, c3, los alamos scientific lab. complex(kind(1.d0)) z, ci, ztan external ztan - data ci /(0.,1.)/ + data ci /(0.d0,1.d0)/ c ztanh = -ci*ztan(ci*z) c -- 2.16.4