php.net |  support |  documentation |  report a bug |  advanced search |  search howto |  statistics |  random bug |  login
Return to Bug #58394
Patch use-current-dcdflib revision 2016-01-20 23:53 UTC by y dot uchiyama dot 1015 at gmail dot com

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;
 
PHP Copyright © 2001-2024 The PHP Group
All rights reserved.
Last updated: Tue Jul 23 17:01:30 2024 UTC