Patch use-current-dcdflib for stats Bug #58394
Patch version 2016-01-20 23:53 UTC
Return to Bug #58394 |
Download this patch
Patch Revisions:
Developer: y.uchiyama.1015@gmail.com
diff --git a/dcdflib.c b/dcdflib.c
index b7f5900..7f6ea99 100644
--- a/dcdflib.c
+++ b/dcdflib.c
@@ -22,13 +22,13 @@
*/
double algdiv(double *a,double *b)
{
-double c0 = .833333333333333e-01;
-double c1 = -.277777777760991e-02;
-double c2 = .793650666825390e-03;
-double c3 = -.595202931351870e-03;
-double c4 = .837308034031215e-03;
-double c5 = -.165322962780713e-02;
-double algdiv,c,d,h,s11,s3,s5,s7,s9,t,u,v,w,x,x2,T1;
+static double c0 = .833333333333333e-01;
+static double c1 = -.277777777760991e-02;
+static double c2 = .793650666825390e-03;
+static double c3 = -.595202931351870e-03;
+static double c4 = .837308034031215e-03;
+static double c5 = -.165322962780713e-02;
+static double algdiv,c,d,h,s11,s3,s5,s7,s9,t,u,v,w,x,x2,T1;
/*
..
.. Executable Statements ..
@@ -109,25 +109,25 @@ double alngam(double *x)
*/
{
#define hln2pi 0.91893853320467274178e0
-double coef[5] = {
+static double coef[5] = {
0.83333333333333023564e-1,-0.27777777768818808e-2,0.79365006754279e-3,
-0.594997310889e-3,0.8065880899e-3
};
-double scoefd[4] = {
+static double scoefd[4] = {
0.62003838007126989331e2,0.9822521104713994894e1,-0.8906016659497461257e1,
0.1000000000000000000e1
};
-double scoefn[9] = {
+static double scoefn[9] = {
0.62003838007127258804e2,0.36036772530024836321e2,0.20782472531792126786e2,
0.6338067999387272343e1,0.215994312846059073e1,0.3980671310203570498e0,
0.1093115956710439502e0,0.92381945590275995e-2,0.29737866448101651e-2
};
-int K1 = 9;
-int K3 = 4;
-int K5 = 5;
-double alngam,offset,prod,xx;
-int i,n;
-double T2,T4,T6;
+static int K1 = 9;
+static int K3 = 4;
+static int K5 = 5;
+static double alngam,offset,prod,xx;
+static int i,n;
+static double T2,T4,T6;
/*
..
.. Executable Statements ..
@@ -192,13 +192,13 @@ double alnrel(double *a)
-----------------------------------------------------------------------
*/
{
-double p1 = -.129418923021993e+01;
-double p2 = .405303492862024e+00;
-double p3 = -.178874546012214e-01;
-double q1 = -.162752256355323e+01;
-double q2 = .747811014037616e+00;
-double q3 = -.845104217945565e-01;
-double alnrel,t,t2,w,x;
+static double p1 = -.129418923021993e+01;
+static double p2 = .405303492862024e+00;
+static double p3 = -.178874546012214e-01;
+static double q1 = -.162752256355323e+01;
+static double q2 = .747811014037616e+00;
+static double q3 = -.845104217945565e-01;
+static double alnrel,t,t2,w,x;
/*
..
.. Executable Statements ..
@@ -223,8 +223,8 @@ double apser(double *a,double *b,double *x,double *eps)
-----------------------------------------------------------------------
*/
{
-double g = .577215664901533e0;
-double apser,aj,bx,c,j,s,t,tol;
+static double g = .577215664901533e0;
+static double apser,aj,bx,c,j,s,t,tol;
/*
..
.. Executable Statements ..
@@ -259,9 +259,9 @@ double basym(double *a,double *b,double *lambda,double *eps)
-----------------------------------------------------------------------
*/
{
-double e0 = 1.12837916709551e0;
-double e1 = .353553390593274e0;
-int num = 20;
+static double e0 = 1.12837916709551e0;
+static double e1 = .353553390593274e0;
+static int num = 20;
/*
------------------------
****** NUM IS THE MAXIMUM VALUE THAT N CAN TAKE IN THE DO LOOP
@@ -272,11 +272,11 @@ int num = 20;
E1 = 2**(-3/2)
------------------------
*/
-int K3 = 1;
-double basym,bsum,dsum,f,h,h2,hn,j0,j1,r,r0,r1,s,sum,t,t0,t1,u,w,w0,z,z0,
+static int K3 = 1;
+static double basym,bsum,dsum,f,h,h2,hn,j0,j1,r,r0,r1,s,sum,t,t0,t1,u,w,w0,z,z0,
z2,zn,znm1;
-int i,im1,imj,j,m,mm1,mmj,n,np1;
-double a0[21],b0[21],c[21],d[21],T1,T2;
+static int i,im1,imj,j,m,mm1,mmj,n,np1;
+static double a0[21],b0[21],c[21],d[21],T1,T2;
/*
..
.. Executable Statements ..
@@ -368,13 +368,13 @@ double bcorr(double *a0,double *b0)
-----------------------------------------------------------------------
*/
{
-double c0 = .833333333333333e-01;
-double c1 = -.277777777760991e-02;
-double c2 = .793650666825390e-03;
-double c3 = -.595202931351870e-03;
-double c4 = .837308034031215e-03;
-double c5 = -.165322962780713e-02;
-double bcorr,a,b,c,h,s11,s3,s5,s7,s9,t,w,x,x2;
+static double c0 = .833333333333333e-01;
+static double c1 = -.277777777760991e-02;
+static double c2 = .793650666825390e-03;
+static double c3 = -.595202931351870e-03;
+static double c4 = .837308034031215e-03;
+static double c5 = -.165322962780713e-02;
+static double bcorr,a,b,c,h,s11,s3,s5,s7,s9,t,w,x,x2;
/*
..
.. Executable Statements ..
@@ -415,10 +415,10 @@ double betaln(double *a0,double *b0)
--------------------------
*/
{
-double e = .918938533204673e0;
-double betaln,a,b,c,h,u,v,w,z;
-int i,n;
-double T1;
+static double e = .918938533204673e0;
+static double betaln,a,b,c,h,u,v,w,z;
+static int i,n;
+static double T1;
/*
..
.. Executable Statements ..
@@ -521,7 +521,7 @@ double bfrac(double *a,double *b,double *x,double *y,double *lambda,
-----------------------------------------------------------------------
*/
{
-double bfrac,alpha,an,anp1,beta,bn,bnp1,c,c0,c1,e,n,p,r,r0,s,t,w,yp1;
+static double bfrac,alpha,an,anp1,beta,bn,bnp1,c,c0,c1,e,n,p,r,r0,s,t,w,yp1;
/*
..
.. Executable Statements ..
@@ -590,9 +590,9 @@ void bgrat(double *a,double *b,double *x,double *y,double *w,
-----------------------------------------------------------------------
*/
{
-double bm1,bp2n,cn,coef,dj,j,l,lnx,n2,nu,p,q,r,s,sum,t,t2,u,v,z;
-int i,n,nm1;
-double c[30],d[30],T1;
+static double bm1,bp2n,cn,coef,dj,j,l,lnx,n2,nu,p,q,r,s,sum,t,t2,u,v,z;
+static int i,n,nm1;
+static double c[30],d[30],T1;
/*
..
.. Executable Statements ..
@@ -669,8 +669,8 @@ double bpser(double *a,double *b,double *x,double *eps)
-----------------------------------------------------------------------
*/
{
-double bpser,a0,apb,b0,c,n,sum,t,tol,u,w,z;
-int i,m;
+static double bpser,a0,apb,b0,c,n,sum,t,tol,u,w,z;
+static int i,m;
/*
..
.. Executable Statements ..
@@ -796,10 +796,10 @@ void bratio(double *a,double *b,double *x,double *y,double *w,
-----------------------------------------------------------------------
*/
{
-int K1 = 1;
-double a0,b0,eps,lambda,t,x0,y0,z;
-int ierr1,ind,n;
-double T2,T3,T4,T5;
+static int K1 = 1;
+static double a0,b0,eps,lambda,t,x0,y0,z;
+static int ierr1,ind,n;
+static double T2,T3,T4,T5;
/*
..
.. Executable Statements ..
@@ -1003,15 +1003,15 @@ double brcmp1(int *mu,double *a,double *b,double *x,double *y)
-----------------------------------------------------------------------
*/
{
-double Const = .398942280401433e0;
-double brcmp1,a0,apb,b0,c,e,h,lambda,lnx,lny,t,u,v,x0,y0,z;
-int i,n;
+static double Const = .398942280401433e0;
+static double brcmp1,a0,apb,b0,c,e,h,lambda,lnx,lny,t,u,v,x0,y0,z;
+static int i,n;
/*
-----------------
CONST = 1/SQRT(2*PI)
-----------------
*/
-double T1,T2,T3,T4;
+static double T1,T2,T3,T4;
/*
..
.. Executable Statements ..
@@ -1141,15 +1141,15 @@ double brcomp(double *a,double *b,double *x,double *y)
-----------------------------------------------------------------------
*/
{
-double Const = .398942280401433e0;
-double brcomp,a0,apb,b0,c,e,h,lambda,lnx,lny,t,u,v,x0,y0,z;
-int i,n;
+static double Const = .398942280401433e0;
+static double brcomp,a0,apb,b0,c,e,h,lambda,lnx,lny,t,u,v,x0,y0,z;
+static int i,n;
/*
-----------------
CONST = 1/SQRT(2*PI)
-----------------
*/
-double T1,T2;
+static double T1,T2;
/*
..
.. Executable Statements ..
@@ -1280,10 +1280,10 @@ double bup(double *a,double *b,double *x,double *y,int *n,double *eps)
-----------------------------------------------------------------------
*/
{
-int K1 = 1;
-int K2 = 0;
-double bup,ap1,apb,d,l,r,t,w;
-int i,k,kp1,mu,nm1;
+static int K1 = 1;
+static int K2 = 0;
+static double bup,ap1,apb,d,l,r,t,w;
+static int i,k,kp1,mu,nm1;
/*
..
.. Executable Statements ..
@@ -1451,14 +1451,14 @@ void cdfbet(int *which,double *p,double *q,double *x,double *y,
#define zero 1.0e-100
#define inf 1.0e100
#define one 1.0e0
-int K1 = 1;
-double K2 = 0.0e0;
-double K3 = 1.0e0;
-double K8 = 0.5e0;
-double K9 = 5.0e0;
-double fx,xhi,xlo,cum,ccum,xy,pq;
-unsigned long qhi,qleft,qporq;
-double T4,T5,T6,T7,T10,T11,T12,T13,T14,T15;
+static int K1 = 1;
+static double K2 = 0.0e0;
+static double K3 = 1.0e0;
+static double K8 = 0.5e0;
+static double K9 = 5.0e0;
+static double fx,xhi,xlo,cum,ccum,xy,pq;
+static unsigned long qhi,qleft,qporq;
+static double T4,T5,T6,T7,T10,T11,T12,T13,T14,T15;
/*
..
.. Executable Statements ..
@@ -1817,14 +1817,14 @@ void cdfbin(int *which,double *p,double *q,double *s,double *xn,
#define zero 1.0e-100
#define inf 1.0e100
#define one 1.0e0
-int K1 = 1;
-double K2 = 0.0e0;
-double K3 = 0.5e0;
-double K4 = 5.0e0;
-double K11 = 1.0e0;
-double fx,xhi,xlo,cum,ccum,pq,prompr;
-unsigned long qhi,qleft,qporq;
-double T5,T6,T7,T8,T9,T10,T12,T13;
+static int K1 = 1;
+static double K2 = 0.0e0;
+static double K3 = 0.5e0;
+static double K4 = 5.0e0;
+static double K11 = 1.0e0;
+static double fx,xhi,xlo,cum,ccum,pq,prompr;
+static unsigned long qhi,qleft,qporq;
+static double T5,T6,T7,T8,T9,T10,T12,T13;
/*
..
.. Executable Statements ..
@@ -2176,13 +2176,13 @@ void cdfchi(int *which,double *p,double *q,double *x,double *df,
#define atol 1.0e-50
#define zero 1.0e-100
#define inf 1.0e100
-int K1 = 1;
-double K2 = 0.0e0;
-double K4 = 0.5e0;
-double K5 = 5.0e0;
-double fx,cum,ccum,pq,porq;
-unsigned long qhi,qleft,qporq;
-double T3,T6,T7,T8,T9,T10,T11;
+static int K1 = 1;
+static double K2 = 0.0e0;
+static double K4 = 0.5e0;
+static double K5 = 5.0e0;
+static double fx,cum,ccum,pq,porq;
+static unsigned long qhi,qleft,qporq;
+static double T3,T6,T7,T8,T9,T10,T11;
/*
..
.. Executable Statements ..
@@ -2473,12 +2473,12 @@ void cdfchn(int *which,double *p,double *q,double *x,double *df,
#define zero 1.0e-100
#define one ( 1.0e0 - 1.0e-16 )
#define inf 1.0e100
-double K1 = 0.0e0;
-double K3 = 0.5e0;
-double K4 = 5.0e0;
-double fx,cum,ccum;
-unsigned long qhi,qleft;
-double T2,T5,T6,T7,T8,T9,T10,T11,T12,T13;
+static double K1 = 0.0e0;
+static double K3 = 0.5e0;
+static double K4 = 5.0e0;
+static double fx,cum,ccum;
+static unsigned long qhi,qleft;
+static double T2,T5,T6,T7,T8,T9,T10,T11,T12,T13;
/*
..
.. Executable Statements ..
@@ -2741,13 +2741,13 @@ void cdff(int *which,double *p,double *q,double *f,double *dfn,
#define atol 1.0e-50
#define zero 1.0e-100
#define inf 1.0e100
-int K1 = 1;
-double K2 = 0.0e0;
-double K4 = 0.5e0;
-double K5 = 5.0e0;
-double pq,fx,cum,ccum;
-unsigned long qhi,qleft,qporq;
-double T3,T6,T7,T8,T9,T10,T11,T12,T13,T14,T15;
+static int K1 = 1;
+static double K2 = 0.0e0;
+static double K4 = 0.5e0;
+static double K5 = 5.0e0;
+static double pq,fx,cum,ccum;
+static unsigned long qhi,qleft,qporq;
+static double T3,T6,T7,T8,T9,T10,T11,T12,T13,T14,T15;
/*
..
.. Executable Statements ..
@@ -3072,12 +3072,12 @@ void cdffnc(int *which,double *p,double *q,double *f,double *dfn,
#define zero 1.0e-100
#define one ( 1.0e0 - 1.0e-16 )
#define inf 1.0e100
-double K1 = 0.0e0;
-double K3 = 0.5e0;
-double K4 = 5.0e0;
-double fx,cum,ccum;
-unsigned long qhi,qleft;
-double T2,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,T15,T16,T17;
+static double K1 = 0.0e0;
+static double K3 = 0.5e0;
+static double K4 = 5.0e0;
+static double fx,cum,ccum;
+static unsigned long qhi,qleft;
+static double T2,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,T15,T16,T17;
/*
..
.. Executable Statements ..
@@ -3387,13 +3387,13 @@ void cdfgam(int *which,double *p,double *q,double *x,double *shape,
#define atol 1.0e-50
#define zero 1.0e-100
#define inf 1.0e100
-int K1 = 1;
-double K5 = 0.5e0;
-double K6 = 5.0e0;
-double xx,fx,xscale,cum,ccum,pq,porq;
-int ierr;
-unsigned long qhi,qleft,qporq;
-double T2,T3,T4,T7,T8,T9;
+static int K1 = 1;
+static double K5 = 0.5e0;
+static double K6 = 5.0e0;
+static double xx,fx,xscale,cum,ccum,pq,porq;
+static int ierr;
+static unsigned long qhi,qleft,qporq;
+static double T2,T3,T4,T7,T8,T9;
/*
..
.. Executable Statements ..
@@ -3689,14 +3689,14 @@ void cdfnbn(int *which,double *p,double *q,double *s,double *xn,
#define atol 1.0e-50
#define inf 1.0e100
#define one 1.0e0
-int K1 = 1;
-double K2 = 0.0e0;
-double K4 = 0.5e0;
-double K5 = 5.0e0;
-double K11 = 1.0e0;
-double fx,xhi,xlo,pq,prompr,cum,ccum;
-unsigned long qhi,qleft,qporq;
-double T3,T6,T7,T8,T9,T10,T12,T13;
+static int K1 = 1;
+static double K2 = 0.0e0;
+static double K4 = 0.5e0;
+static double K5 = 5.0e0;
+static double K11 = 1.0e0;
+static double fx,xhi,xlo,pq,prompr,cum,ccum;
+static unsigned long qhi,qleft,qporq;
+static double T3,T6,T7,T8,T9,T10,T12,T13;
/*
..
.. Executable Statements ..
@@ -4054,8 +4054,8 @@ void cdfnor(int *which,double *p,double *q,double *x,double *mean,
**********************************************************************/
{
-int K1 = 1;
-double z,pq;
+static int K1 = 1;
+static double z,pq;
/*
..
.. Executable Statements ..
@@ -4242,13 +4242,13 @@ void cdfpoi(int *which,double *p,double *q,double *s,double *xlam,
#define tol 1.0e-8
#define atol 1.0e-50
#define inf 1.0e100
-int K1 = 1;
-double K2 = 0.0e0;
-double K4 = 0.5e0;
-double K5 = 5.0e0;
-double fx,cum,ccum,pq;
-unsigned long qhi,qleft,qporq;
-double T3,T6,T7,T8,T9,T10;
+static int K1 = 1;
+static double K2 = 0.0e0;
+static double K4 = 0.5e0;
+static double K5 = 5.0e0;
+static double fx,cum,ccum,pq;
+static unsigned long qhi,qleft,qporq;
+static double T3,T6,T7,T8,T9,T10;
/*
..
.. Executable Statements ..
@@ -4500,12 +4500,12 @@ void cdft(int *which,double *p,double *q,double *t,double *df,
#define inf 1.0e100
#define rtinf 1.0e100
#define maxdf 1.0e10
-int K1 = 1;
-double K4 = 0.5e0;
-double K5 = 5.0e0;
-double fx,cum,ccum,pq;
-unsigned long qhi,qleft,qporq;
-double T2,T3,T6,T7,T8,T9,T10,T11;
+static int K1 = 1;
+static double K4 = 0.5e0;
+static double K5 = 5.0e0;
+static double fx,cum,ccum,pq;
+static unsigned long qhi,qleft,qporq;
+static double T2,T3,T6,T7,T8,T9,T10,T11;
/*
..
.. Executable Statements ..
@@ -4751,11 +4751,11 @@ void cdftnc(int *which,double *p,double *q,double *t,double *df,
#define zero 1.0e-100
#define one ( 1.0e0 - 1.0e-16 )
#define inf 1.0e100
-double K3 = 0.5e0;
-double K4 = 5.0e0;
-double ccum,cum,fx;
-unsigned long qhi,qleft;
-double T1,T2,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14;
+static double K3 = 0.5e0;
+static double K4 = 5.0e0;
+static double ccum,cum,fx;
+static unsigned long qhi,qleft;
+static double T1,T2,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14;
/*
..
.. Executable Statements ..
@@ -4942,7 +4942,7 @@ void cumbet(double *x,double *y,double *a,double *b,double *cum,
**********************************************************************
*/
{
-int ierr;
+static int ierr;
/*
..
.. Executable Statements ..
@@ -5013,7 +5013,7 @@ void cumbin(double *s,double *xn,double *pr,double *ompr,
**********************************************************************
*/
{
-double T1,T2;
+static double T1,T2;
/*
..
.. Executable Statements ..
@@ -5069,7 +5069,7 @@ void cumchi(double *x,double *df,double *cum,double *ccum)
**********************************************************************
*/
{
-double a,xx;
+static double a,xx;
/*
..
.. Executable Statements ..
@@ -5130,11 +5130,11 @@ void cumchn(double *x,double *df,double *pnonc,double *cum,
{
#define dg(i) (*df + 2.0e0 * (double)(i))
#define qsmall(xx) (int)(sum < 1.0e-20 || (xx) < eps * sum)
-double eps = 1.0e-5;
-double adj,centaj,centwt,chid2,dfd2,lcntaj,lcntwt,lfact,pcent,pterm,sum,
+static double eps = 1.0e-5;
+static double adj,centaj,centwt,chid2,dfd2,lcntaj,lcntwt,lfact,pcent,pterm,sum,
sumadj,term,wt,xnonc;
-int i,icent;
-double T1,T2,T3;
+static int i,icent;
+static double T1,T2,T3;
/*
..
.. Executable Statements ..
@@ -5310,9 +5310,9 @@ void cumf(double *f,double *dfn,double *dfd,double *cum,double *ccum)
{
#define half 0.5e0
#define done 1.0e0
-double dsum,prod,xx,yy;
-int ierr;
-double T1,T2;
+static double dsum,prod,xx,yy;
+static int ierr;
+static double T1,T2;
/*
..
.. Executable Statements ..
@@ -5410,11 +5410,11 @@ void cumfnc(double *f,double *dfn,double *dfd,double *pnonc,
#define qsmall(x) (int)(sum < 1.0e-20 || (x) < eps*sum)
#define half 0.5e0
#define done 1.0e0
-double eps = 1.0e-4;
-double dsum,dummy,prod,xx,yy,adn,aup,b,betdn,betup,centwt,dnterm,sum,
+static double eps = 1.0e-4;
+static double dsum,dummy,prod,xx,yy,adn,aup,b,betdn,betup,centwt,dnterm,sum,
upterm,xmult,xnonc;
-int i,icent,ierr;
-double T1,T2,T3,T4,T5,T6;
+static int i,icent,ierr;
+static double T1,T2,T3,T4,T5,T6;
/*
..
.. Executable Statements ..
@@ -5556,7 +5556,7 @@ void cumgam(double *x,double *a,double *cum,double *ccum)
**********************************************************************
*/
{
-int K1 = 0;
+static int K1 = 0;
/*
..
.. Executable Statements ..
@@ -5628,7 +5628,7 @@ void cumnbn(double *s,double *xn,double *pr,double *ompr,
**********************************************************************
*/
{
-double T1;
+static double T1;
/*
..
.. Executable Statements ..
@@ -5725,43 +5725,43 @@ void cumnor(double *arg,double *result,double *ccum)
------------------------------------------------------------------
*/
{
-double a[5] = {
+static double a[5] = {
2.2352520354606839287e00,1.6102823106855587881e02,1.0676894854603709582e03,
1.8154981253343561249e04,6.5682337918207449113e-2
};
-double b[4] = {
+static double b[4] = {
4.7202581904688241870e01,9.7609855173777669322e02,1.0260932208618978205e04,
4.5507789335026729956e04
};
-double c[9] = {
+static double c[9] = {
3.9894151208813466764e-1,8.8831497943883759412e00,9.3506656132177855979e01,
5.9727027639480026226e02,2.4945375852903726711e03,6.8481904505362823326e03,
1.1602651437647350124e04,9.8427148383839780218e03,1.0765576773720192317e-8
};
-double d[8] = {
+static double d[8] = {
2.2266688044328115691e01,2.3538790178262499861e02,1.5193775994075548050e03,
6.4855582982667607550e03,1.8615571640885098091e04,3.4900952721145977266e04,
3.8912003286093271411e04,1.9685429676859990727e04
};
-double half = 0.5e0;
-double p[6] = {
+static double half = 0.5e0;
+static double p[6] = {
2.1589853405795699e-1,1.274011611602473639e-1,2.2235277870649807e-2,
1.421619193227893466e-3,2.9112874951168792e-5,2.307344176494017303e-2
};
-double one = 1.0e0;
-double q[5] = {
+static double one = 1.0e0;
+static double q[5] = {
1.28426009614491121e00,4.68238212480865118e-1,6.59881378689285515e-2,
3.78239633202758244e-3,7.29751555083966205e-5
};
-double sixten = 1.60e0;
-double sqrpi = 3.9894228040143267794e-1;
-double thrsh = 0.66291e0;
-double root32 = 5.656854248e0;
-double zero = 0.0e0;
-int K1 = 1;
-int K2 = 2;
-int i;
-double del,eps,temp,x,xden,xnum,y,xsq,min;
+static double sixten = 1.60e0;
+static double sqrpi = 3.9894228040143267794e-1;
+static double thrsh = 0.66291e0;
+static double root32 = 5.656854248e0;
+static double zero = 0.0e0;
+static int K1 = 1;
+static int K2 = 2;
+static int i;
+static double del,eps,temp,x,xden,xnum,y,xsq,min;
/*
------------------------------------------------------------------
Machine dependent constants
@@ -5889,7 +5889,7 @@ void cumpoi(double *s,double *xlam,double *cum,double *ccum)
**********************************************************************
*/
{
-double chi,df;
+static double chi,df;
/*
..
.. Executable Statements ..
@@ -5939,8 +5939,8 @@ void cumt(double *t,double *df,double *cum,double *ccum)
**********************************************************************
*/
{
-double K2 = 0.5e0;
-double xx,a,oma,tt,yy,dfptt,T1;
+static double K2 = 0.5e0;
+static double xx,a,oma,tt,yy,dfptt,T1;
/*
..
.. Executable Statements ..
@@ -6009,12 +6009,12 @@ void cumtnc(double *t,double *df,double *pnonc,double *cum,
#define onep5 1.5e0
#define conv 1.0e-7
#define tiny 1.0e-10
-double alghdf,b,bb,bbcent,bcent,cent,d,dcent,dpnonc,dum1,dum2,e,ecent,
+static double alghdf,b,bb,bbcent,bcent,cent,d,dcent,dpnonc,dum1,dum2,e,ecent,
halfdf,lambda,lnomx,lnx,omx,pnonc2,s,scent,ss,sscent,t2,term,tt,twoi,x,xi,
xlnd,xlne;
-int ierr;
-unsigned long qrevs;
-double T1,T2,T3,T4,T5,T6,T7,T8,T9,T10;
+static int ierr;
+static unsigned long qrevs;
+static double T1,T2,T3,T4,T5,T6,T7,T8,T9,T10;
/*
..
.. Executable Statements ..
@@ -6222,8 +6222,8 @@ double devlpl(double a[],int *n,double *x)
**********************************************************************
*/
{
-double devlpl,term;
-int i;
+static double devlpl,term;
+static int i;
/*
..
.. Executable Statements ..
@@ -6279,9 +6279,9 @@ double dinvnr(double *p,double *q)
#define r2pi 0.3989422804014326e0
#define nhalf -0.5e0
#define dennor(x) (r2pi*exp(nhalf*(x)*(x)))
-double dinvnr,strtx,xcur,cum,ccum,pp,dx;
-int i;
-unsigned long qporq;
+static double dinvnr,strtx,xcur,cum,ccum,pp,dx;
+static int i;
+static unsigned long qporq;
/*
..
.. Executable Statements ..
@@ -6336,10 +6336,10 @@ static void E0000(int IENTRY,int *status,double *x,double *fx,
double *zrelto,double *zsmall,double *zstpmu)
{
#define qxmon(zx,zy,zz) (int)((zx) <= (zy) && (zy) <= (zz))
-double absstp,abstol,big,fbig,fsmall,relstp,reltol,small,step,stpmul,xhi,
+static double absstp,abstol,big,fbig,fsmall,relstp,reltol,small,step,stpmul,xhi,
xlb,xlo,xsave,xub,yy;
-int i99999;
-unsigned long qbdd,qcond,qdum1,qdum2,qincr,qlim,qok,qup;
+static int i99999;
+static unsigned long qbdd,qcond,qdum1,qdum2,qincr,qlim,qok,qup;
switch(IENTRY){case 0: goto DINVR; case 1: goto DSTINV;}
DINVR:
if(*status > 0) goto S310;
@@ -6704,20 +6704,20 @@ double dt1(double *p,double *q,double *df)
**********************************************************************
*/
{
-double coef[4][5] = {
+static double coef[4][5] = {
{1.0e0,1.0e0,0.0e0,0.0e0,0.0e0},
{3.0e0,16.0e0,5.0e0,0.0e0,0.0e0},
{-15.0e0,17.0e0,19.0e0,3.0e0,0.0e0},
{-945.0e0,-1920.0e0,1482.0e0,776.0e0,79.0e0}
};
-double denom[4] = {
+static double denom[4] = {
4.0e0,96.0e0,384.0e0,92160.0e0
};
-int ideg[4] = {
+static int ideg[4] = {
2,3,4,5
};
-double dt1,denpow,sum,term,x,xp,xx;
-int i;
+static double dt1,denpow,sum,term,x,xp,xx;
+static int i;
/*
..
.. Executable Statements ..
@@ -6747,9 +6747,9 @@ static void E0001(int IENTRY,int *status,double *x,double *fx,
double *zxhi,double *zxlo)
{
#define ftol(zx) (0.5e0*fifdmax1(abstol,reltol*fabs((zx))))
-double a,abstol,b,c,d,fa,fb,fc,fd,fda,fdb,m,mb,p,q,reltol,tol,w,xxhi,xxlo;
-int ext,i99999;
-unsigned long first,qrzero;
+static double a,abstol,b,c,d,fa,fb,fc,fd,fda,fdb,m,mb,p,q,reltol,tol,w,xxhi,xxlo;
+static int ext,i99999;
+static unsigned long first,qrzero;
switch(IENTRY){case 0: goto DZROR; case 1: goto DSTZR;}
DZROR:
if(*status > 0) goto S280;
@@ -7021,33 +7021,33 @@ double erf1(double *x)
-----------------------------------------------------------------------
*/
{
-double c = .564189583547756e0;
-double a[5] = {
+static double c = .564189583547756e0;
+static double a[5] = {
.771058495001320e-04,-.133733772997339e-02,.323076579225834e-01,
.479137145607681e-01,.128379167095513e+00
};
-double b[3] = {
+static double b[3] = {
.301048631703895e-02,.538971687740286e-01,.375795757275549e+00
};
-double p[8] = {
+static double p[8] = {
-1.36864857382717e-07,5.64195517478974e-01,7.21175825088309e+00,
4.31622272220567e+01,1.52989285046940e+02,3.39320816734344e+02,
4.51918953711873e+02,3.00459261020162e+02
};
-double q[8] = {
+static double q[8] = {
1.00000000000000e+00,1.27827273196294e+01,7.70001529352295e+01,
2.77585444743988e+02,6.38980264465631e+02,9.31354094850610e+02,
7.90950925327898e+02,3.00459260956983e+02
};
-double r[5] = {
+static double r[5] = {
2.10144126479064e+00,2.62370141675169e+01,2.13688200555087e+01,
4.65807828718470e+00,2.82094791773523e-01
};
-double s[4] = {
+static double s[4] = {
9.41537750555460e+01,1.87114811799590e+02,9.90191814623914e+01,
1.80124575948747e+01
};
-double erf1,ax,bot,t,top,x2;
+static double erf1,ax,bot,t,top,x2;
/*
..
.. Executable Statements ..
@@ -7092,34 +7092,34 @@ double erfc1(int *ind,double *x)
-----------------------------------------------------------------------
*/
{
-double c = .564189583547756e0;
-double a[5] = {
+static double c = .564189583547756e0;
+static double a[5] = {
.771058495001320e-04,-.133733772997339e-02,.323076579225834e-01,
.479137145607681e-01,.128379167095513e+00
};
-double b[3] = {
+static double b[3] = {
.301048631703895e-02,.538971687740286e-01,.375795757275549e+00
};
-double p[8] = {
+static double p[8] = {
-1.36864857382717e-07,5.64195517478974e-01,7.21175825088309e+00,
4.31622272220567e+01,1.52989285046940e+02,3.39320816734344e+02,
4.51918953711873e+02,3.00459261020162e+02
};
-double q[8] = {
+static double q[8] = {
1.00000000000000e+00,1.27827273196294e+01,7.70001529352295e+01,
2.77585444743988e+02,6.38980264465631e+02,9.31354094850610e+02,
7.90950925327898e+02,3.00459260956983e+02
};
-double r[5] = {
+static double r[5] = {
2.10144126479064e+00,2.62370141675169e+01,2.13688200555087e+01,
4.65807828718470e+00,2.82094791773523e-01
};
-double s[4] = {
+static double s[4] = {
9.41537750555460e+01,1.87114811799590e+02,9.90191814623914e+01,
1.80124575948747e+01
};
-int K1 = 1;
-double erfc1,ax,bot,e,t,top,w;
+static int K1 = 1;
+static double erfc1,ax,bot,e,t,top,w;
/*
..
.. Executable Statements ..
@@ -7195,7 +7195,7 @@ double esum(int *mu,double *x)
-----------------------------------------------------------------------
*/
{
-double esum,w;
+static double esum,w;
/*
..
.. Executable Statements ..
@@ -7230,11 +7230,11 @@ double exparg(int *l)
--------------------------------------------------------------------
*/
{
-int K1 = 4;
-int K2 = 9;
-int K3 = 10;
-double exparg,lnb;
-int b,m;
+static int K1 = 4;
+static int K2 = 9;
+static int K3 = 10;
+static double exparg,lnb;
+static int b,m;
/*
..
.. Executable Statements ..
@@ -7277,8 +7277,8 @@ double fpser(double *a,double *b,double *x,double *eps)
SET FPSER = X**A
*/
{
-int K1 = 1;
-double fpser,an,c,s,t,tol;
+static int K1 = 1;
+static double fpser,an,c,s,t,tol;
/*
..
.. Executable Statements ..
@@ -7314,23 +7314,23 @@ double gam1(double *a)
------------------------------------------------------------------
*/
{
-double s1 = .273076135303957e+00;
-double s2 = .559398236957378e-01;
-double p[7] = {
+static double s1 = .273076135303957e+00;
+static double s2 = .559398236957378e-01;
+static double p[7] = {
.577215664901533e+00,-.409078193005776e+00,-.230975380857675e+00,
.597275330452234e-01,.766968181649490e-02,-.514889771323592e-02,
.589597428611429e-03
};
-double q[5] = {
+static double q[5] = {
.100000000000000e+01,.427569613095214e+00,.158451672430138e+00,
.261132021441447e-01,.423244297896961e-02
};
-double r[9] = {
+static double r[9] = {
-.422784335098468e+00,-.771330383816272e+00,-.244757765222226e+00,
.118378989872749e+00,.930357293360349e-03,-.118290993445146e-01,
.223047661158249e-02,.266505979058923e-03,-.132674909766242e-03
};
-double gam1,bot,d,t,top,w,T1;
+static double gam1,bot,d,t,top,w,T1;
/*
..
.. Executable Statements ..
@@ -7423,40 +7423,40 @@ void gaminv(double *a,double *x,double *x0,double *p,double *q,
-------------------
*/
{
-double a0 = 3.31125922108741e0;
-double a1 = 11.6616720288968e0;
-double a2 = 4.28342155967104e0;
-double a3 = .213623493715853e0;
-double b1 = 6.61053765625462e0;
-double b2 = 6.40691597760039e0;
-double b3 = 1.27364489782223e0;
-double b4 = .036117081018842e0;
-double c = .577215664901533e0;
-double ln10 = 2.302585e0;
-double tol = 1.e-5;
-double amin[2] = {
+static double a0 = 3.31125922108741e0;
+static double a1 = 11.6616720288968e0;
+static double a2 = 4.28342155967104e0;
+static double a3 = .213623493715853e0;
+static double b1 = 6.61053765625462e0;
+static double b2 = 6.40691597760039e0;
+static double b3 = 1.27364489782223e0;
+static double b4 = .036117081018842e0;
+static double c = .577215664901533e0;
+static double ln10 = 2.302585e0;
+static double tol = 1.e-5;
+static double amin[2] = {
500.0e0,100.0e0
};
-double bmin[2] = {
+static double bmin[2] = {
1.e-28,1.e-13
};
-double dmin[2] = {
+static double dmin[2] = {
1.e-06,1.e-04
};
-double emin[2] = {
+static double emin[2] = {
2.e-03,6.e-03
};
-double eps0[2] = {
+static double eps0[2] = {
1.e-10,1.e-08
};
-int K1 = 1;
-int K2 = 2;
-int K3 = 3;
-int K8 = 0;
-double am1,amax,ap1,ap2,ap3,apn,b,c1,c2,c3,c4,c5,d,e,e2,eps,g,h,pn,qg,qn,
+static int K1 = 1;
+static int K2 = 2;
+static int K3 = 3;
+static int K8 = 0;
+static double am1,amax,ap1,ap2,ap3,apn,b,c1,c2,c3,c4,c5,d,e,e2,eps,g,h,pn,qg,qn,
r,rta,s,s2,sum,t,u,w,xmax,xmin,xn,y,z;
-int iop;
-double T4,T5,T6,T7,T9;
+static int iop;
+static double T4,T5,T6,T7,T9;
/*
..
.. Executable Statements ..
@@ -7742,16 +7742,16 @@ double gamln(double *a)
--------------------------
*/
{
-double c0 = .833333333333333e-01;
-double c1 = -.277777777760991e-02;
-double c2 = .793650666825390e-03;
-double c3 = -.595202931351870e-03;
-double c4 = .837308034031215e-03;
-double c5 = -.165322962780713e-02;
-double d = .418938533204673e0;
-double gamln,t,w;
-int i,n;
-double T1;
+static double c0 = .833333333333333e-01;
+static double c1 = -.277777777760991e-02;
+static double c2 = .793650666825390e-03;
+static double c3 = -.595202931351870e-03;
+static double c4 = .837308034031215e-03;
+static double c5 = -.165322962780713e-02;
+static double d = .418938533204673e0;
+static double gamln,t,w;
+static int i,n;
+static double T1;
/*
..
.. Executable Statements ..
@@ -7789,31 +7789,31 @@ double gamln1(double *a)
-----------------------------------------------------------------------
*/
{
-double p0 = .577215664901533e+00;
-double p1 = .844203922187225e+00;
-double p2 = -.168860593646662e+00;
-double p3 = -.780427615533591e+00;
-double p4 = -.402055799310489e+00;
-double p5 = -.673562214325671e-01;
-double p6 = -.271935708322958e-02;
-double q1 = .288743195473681e+01;
-double q2 = .312755088914843e+01;
-double q3 = .156875193295039e+01;
-double q4 = .361951990101499e+00;
-double q5 = .325038868253937e-01;
-double q6 = .667465618796164e-03;
-double r0 = .422784335098467e+00;
-double r1 = .848044614534529e+00;
-double r2 = .565221050691933e+00;
-double r3 = .156513060486551e+00;
-double r4 = .170502484022650e-01;
-double r5 = .497958207639485e-03;
-double s1 = .124313399877507e+01;
-double s2 = .548042109832463e+00;
-double s3 = .101552187439830e+00;
-double s4 = .713309612391000e-02;
-double s5 = .116165475989616e-03;
-double gamln1,w,x;
+static double p0 = .577215664901533e+00;
+static double p1 = .844203922187225e+00;
+static double p2 = -.168860593646662e+00;
+static double p3 = -.780427615533591e+00;
+static double p4 = -.402055799310489e+00;
+static double p5 = -.673562214325671e-01;
+static double p6 = -.271935708322958e-02;
+static double q1 = .288743195473681e+01;
+static double q2 = .312755088914843e+01;
+static double q3 = .156875193295039e+01;
+static double q4 = .361951990101499e+00;
+static double q5 = .325038868253937e-01;
+static double q6 = .667465618796164e-03;
+static double r0 = .422784335098467e+00;
+static double r1 = .848044614534529e+00;
+static double r2 = .565221050691933e+00;
+static double r3 = .156513060486551e+00;
+static double r4 = .170502484022650e-01;
+static double r5 = .497958207639485e-03;
+static double s1 = .124313399877507e+01;
+static double s2 = .548042109832463e+00;
+static double s3 = .101552187439830e+00;
+static double s4 = .713309612391000e-02;
+static double s5 = .116165475989616e-03;
+static double gamln1,w,x;
/*
..
.. Executable Statements ..
@@ -7848,25 +7848,25 @@ double Xgamm(double *a)
-----------------------------------------------------------------------
*/
{
-double d = .41893853320467274178e0;
-double pi = 3.1415926535898e0;
-double r1 = .820756370353826e-03;
-double r2 = -.595156336428591e-03;
-double r3 = .793650663183693e-03;
-double r4 = -.277777777770481e-02;
-double r5 = .833333333333333e-01;
-double p[7] = {
+static double d = .41893853320467274178e0;
+static double pi = 3.1415926535898e0;
+static double r1 = .820756370353826e-03;
+static double r2 = -.595156336428591e-03;
+static double r3 = .793650663183693e-03;
+static double r4 = -.277777777770481e-02;
+static double r5 = .833333333333333e-01;
+static double p[7] = {
.539637273585445e-03,.261939260042690e-02,.204493667594920e-01,
.730981088720487e-01,.279648642639792e+00,.553413866010467e+00,1.0e0
};
-double q[7] = {
+static double q[7] = {
-.832979206704073e-03,.470059485860584e-02,.225211131035340e-01,
-.170458969313360e+00,-.567902761974940e-01,.113062953091122e+01,1.0e0
};
-int K2 = 3;
-int K3 = 0;
-double Xgamm,bot,g,lnx,s,t,top,w,x,z;
-int i,j,m,n,T1;
+static int K2 = 3;
+static int K3 = 0;
+static double Xgamm,bot,g,lnx,s,t,top,w,x,z;
+static int i,j,m,n,T1;
/*
..
.. Executable Statements ..
@@ -7982,8 +7982,8 @@ S120:
void grat1(double *a,double *x,double *r,double *p,double *q,
double *eps)
{
-int K2 = 0;
-double a2n,a2nm1,am0,an,an0,b2n,b2nm1,c,cma,g,h,j,l,sum,t,tol,w,z,T1,T3;
+static int K2 = 0;
+static double a2n,a2nm1,am0,an,an0,b2n,b2nm1,c,cma,g,h,j,l,sum,t,tol,w,z,T1,T3;
/*
..
.. Executable Statements ..
@@ -8113,72 +8113,72 @@ void gratio(double *a,double *x,double *ans,double *qans,int *ind)
--------------------
*/
{
-double alog10 = 2.30258509299405e0;
-double d10 = -.185185185185185e-02;
-double d20 = .413359788359788e-02;
-double d30 = .649434156378601e-03;
-double d40 = -.861888290916712e-03;
-double d50 = -.336798553366358e-03;
-double d60 = .531307936463992e-03;
-double d70 = .344367606892378e-03;
-double rt2pin = .398942280401433e0;
-double rtpi = 1.77245385090552e0;
-double third = .333333333333333e0;
-double acc0[3] = {
+static double alog10 = 2.30258509299405e0;
+static double d10 = -.185185185185185e-02;
+static double d20 = .413359788359788e-02;
+static double d30 = .649434156378601e-03;
+static double d40 = -.861888290916712e-03;
+static double d50 = -.336798553366358e-03;
+static double d60 = .531307936463992e-03;
+static double d70 = .344367606892378e-03;
+static double rt2pin = .398942280401433e0;
+static double rtpi = 1.77245385090552e0;
+static double third = .333333333333333e0;
+static double acc0[3] = {
5.e-15,5.e-7,5.e-4
};
-double big[3] = {
+static double big[3] = {
20.0e0,14.0e0,10.0e0
};
-double d0[13] = {
+static double d0[13] = {
.833333333333333e-01,-.148148148148148e-01,.115740740740741e-02,
.352733686067019e-03,-.178755144032922e-03,.391926317852244e-04,
-.218544851067999e-05,-.185406221071516e-05,.829671134095309e-06,
-.176659527368261e-06,.670785354340150e-08,.102618097842403e-07,
-.438203601845335e-08
};
-double d1[12] = {
+static double d1[12] = {
-.347222222222222e-02,.264550264550265e-02,-.990226337448560e-03,
.205761316872428e-03,-.401877572016461e-06,-.180985503344900e-04,
.764916091608111e-05,-.161209008945634e-05,.464712780280743e-08,
.137863344691572e-06,-.575254560351770e-07,.119516285997781e-07
};
-double d2[10] = {
+static double d2[10] = {
-.268132716049383e-02,.771604938271605e-03,.200938786008230e-05,
-.107366532263652e-03,.529234488291201e-04,-.127606351886187e-04,
.342357873409614e-07,.137219573090629e-05,-.629899213838006e-06,
.142806142060642e-06
};
-double d3[8] = {
+static double d3[8] = {
.229472093621399e-03,-.469189494395256e-03,.267720632062839e-03,
-.756180167188398e-04,-.239650511386730e-06,.110826541153473e-04,
-.567495282699160e-05,.142309007324359e-05
};
-double d4[6] = {
+static double d4[6] = {
.784039221720067e-03,-.299072480303190e-03,-.146384525788434e-05,
.664149821546512e-04,-.396836504717943e-04,.113757269706784e-04
};
-double d5[4] = {
+static double d5[4] = {
-.697281375836586e-04,.277275324495939e-03,-.199325705161888e-03,
.679778047793721e-04
};
-double d6[2] = {
+static double d6[2] = {
-.592166437353694e-03,.270878209671804e-03
};
-double e00[3] = {
+static double e00[3] = {
.25e-3,.25e-1,.14e0
};
-double x00[3] = {
+static double x00[3] = {
31.0e0,17.0e0,9.7e0
};
-int K1 = 1;
-int K2 = 0;
-double a2n,a2nm1,acc,am0,amn,an,an0,apn,b2n,b2nm1,c,c0,c1,c2,c3,c4,c5,c6,
+static int K1 = 1;
+static int K2 = 0;
+static double a2n,a2nm1,acc,am0,amn,an,an0,apn,b2n,b2nm1,c,c0,c1,c2,c3,c4,c5,c6,
cma,e,e0,g,h,j,l,r,rta,rtx,s,sum,t,t1,tol,twoa,u,w,x0,y,z;
-int i,iop,m,max,n;
-double wk[20],T3;
-int T4,T5;
-double T6,T7;
+static int i,iop,m,max,n;
+static double wk[20],T3;
+static int T4,T5;
+static double T6,T7;
/*
..
.. Executable Statements ..
@@ -8507,7 +8507,7 @@ double gsumln(double *a,double *b)
-----------------------------------------------------------------------
*/
{
-double gsumln,x,T1,T2;
+static double gsumln,x,T1,T2;
/*
..
.. Executable Statements ..
@@ -8548,29 +8548,29 @@ double psi(double *xx)
---------------------------------------------------------------------
*/
{
-double dx0 = 1.461632144968362341262659542325721325e0;
-double piov4 = .785398163397448e0;
-double p1[7] = {
+static double dx0 = 1.461632144968362341262659542325721325e0;
+static double piov4 = .785398163397448e0;
+static double p1[7] = {
.895385022981970e-02,.477762828042627e+01,.142441585084029e+03,
.118645200713425e+04,.363351846806499e+04,.413810161269013e+04,
.130560269827897e+04
};
-double p2[4] = {
+static double p2[4] = {
-.212940445131011e+01,-.701677227766759e+01,-.448616543918019e+01,
-.648157123766197e+00
};
-double q1[6] = {
+static double q1[6] = {
.448452573429826e+02,.520752771467162e+03,.221000799247830e+04,
.364127349079381e+04,.190831076596300e+04,.691091682714533e-05
};
-double q2[4] = {
+static double q2[4] = {
.322703493791143e+02,.892920700481861e+02,.546117738103215e+02,
.777788548522962e+01
};
-int K1 = 3;
-int K2 = 1;
-double psi,aug,den,sgn,upper,w,x,xmax1,xmx0,xsmall,z;
-int i,m,n,nq;
+static int K1 = 3;
+static int K2 = 1;
+static double psi,aug,den,sgn,upper,w,x,xmax1,xmx0,xsmall,z;
+static int i,m,n,nq;
/*
..
.. Executable Statements ..
@@ -8728,8 +8728,8 @@ double rcomp(double *a,double *x)
-------------------
*/
{
-double rt2pin = .398942280401433e0;
-double rcomp,t,t1,u;
+static double rt2pin = .398942280401433e0;
+static double rcomp,t,t1,u;
/*
..
.. Executable Statements ..
@@ -8759,13 +8759,13 @@ double rexp(double *x)
-----------------------------------------------------------------------
*/
{
-double p1 = .914041914819518e-09;
-double p2 = .238082361044469e-01;
-double q1 = -.499999999085958e+00;
-double q2 = .107141568980644e+00;
-double q3 = -.119041179760821e-01;
-double q4 = .595130811860248e-03;
-double rexp,w;
+static double p1 = .914041914819518e-09;
+static double p2 = .238082361044469e-01;
+static double q1 = -.499999999085958e+00;
+static double q2 = .107141568980644e+00;
+static double q3 = -.119041179760821e-01;
+static double q4 = .595130811860248e-03;
+static double rexp,w;
/*
..
.. Executable Statements ..
@@ -8789,14 +8789,14 @@ double rlog(double *x)
-------------------
*/
{
-double a = .566749439387324e-01;
-double b = .456512608815524e-01;
-double p0 = .333333333333333e+00;
-double p1 = -.224696413112536e+00;
-double p2 = .620886815375787e-02;
-double q1 = -.127408923933623e+01;
-double q2 = .354508718369557e+00;
-double rlog,r,t,u,w,w1;
+static double a = .566749439387324e-01;
+static double b = .456512608815524e-01;
+static double p0 = .333333333333333e+00;
+static double p1 = -.224696413112536e+00;
+static double p2 = .620886815375787e-02;
+static double q1 = -.127408923933623e+01;
+static double q2 = .354508718369557e+00;
+static double rlog,r,t,u,w,w1;
/*
..
.. Executable Statements ..
@@ -8839,14 +8839,14 @@ double rlog1(double *x)
-----------------------------------------------------------------------
*/
{
-double a = .566749439387324e-01;
-double b = .456512608815524e-01;
-double p0 = .333333333333333e+00;
-double p1 = -.224696413112536e+00;
-double p2 = .620886815375787e-02;
-double q1 = -.127408923933623e+01;
-double q2 = .354508718369557e+00;
-double rlog1,h,r,t,w,w1;
+static double a = .566749439387324e-01;
+static double b = .456512608815524e-01;
+static double p0 = .333333333333333e+00;
+static double p1 = -.224696413112536e+00;
+static double p2 = .620886815375787e-02;
+static double q1 = -.127408923933623e+01;
+static double q2 = .354508718369557e+00;
+static double rlog1,h,r,t,w,w1;
/*
..
.. Executable Statements ..
@@ -8911,12 +8911,12 @@ double spmpar(int *i)
-----------------------------------------------------------------------
*/
{
-int K1 = 4;
-int K2 = 8;
-int K3 = 9;
-int K4 = 10;
-double spmpar,b,binv,bm1,one,w,z;
-int emax,emin,ibeta,m;
+static int K1 = 4;
+static int K2 = 8;
+static int K3 = 9;
+static int K4 = 10;
+static double spmpar,b,binv,bm1,one,w,z;
+static int emax,emin,ibeta,m;
/*
..
.. Executable Statements ..
@@ -8980,16 +8980,16 @@ double stvaln(double *p)
**********************************************************************
*/
{
-double xden[5] = {
+static double xden[5] = {
0.993484626060e-1,0.588581570495e0,0.531103462366e0,0.103537752850e0,
0.38560700634e-2
};
-double xnum[5] = {
+static double xnum[5] = {
-0.322232431088e0,-1.000000000000e0,-0.342242088547e0,-0.204231210245e-1,
-0.453642210148e-4
};
-int K1 = 5;
-double stvaln,sign,y,z;
+static int K1 = 5;
+static double stvaln,sign,y,z;
/*
..
.. Executable Statements ..
diff --git a/ipmpar.c b/ipmpar.c
index c0b2926..1762ed6 100644
--- a/ipmpar.c
+++ b/ipmpar.c
@@ -67,8 +67,8 @@ int ipmpar(int*);
*/
int ipmpar(int *i)
{
-int imach[11];
-int ipmpar;
+static int imach[11];
+static int ipmpar;
/* MACHINE CONSTANTS FOR AMDAHL MACHINES. */
/*
imach[1] = 2;
|
Copyright © 2001-2024 The PHP Group