1 /* e_jnf.c -- float version of e_jn.c.
2  * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
23 static const volatile float vone = 1, vzero = 0;
25 static const float
29 static const float zero  =  0.0000000000e+00;
31 float
32 jnf(int n, float x)  in jnf()
35 	float a, b, temp, di;  in jnf()
36 	float z, w;  in jnf()
56 	else if((float)n<=x) {  in jnf()
62 		b = b*((float)(i+i)/x) - a; /* avoid underflow */  in jnf()
73 		    temp = x*(float)0.5; b = temp;  in jnf()
75 			a *= (float)i;		/* a = n! */  in jnf()
110 		float t,v;  in jnf()
111 		float q0,q1,h,tmp; int32_t k,m;  in jnf()
112 		w  = (n+n)/(float)x; h = (float)2.0/(float)x;  in jnf()
113 		q0 = w;  z = w+h; q1 = w*z - (float)1.0; k=1;  in jnf()
114 		while(q1<(float)1.0e9) {  in jnf()
135 		if(tmp<(float)8.8721679688e+01) {  in jnf()
136 	    	    for(i=n-1,di=(float)(i+i);i>0;i--){  in jnf()
144 	    	    for(i=n-1,di=(float)(i+i);i>0;i--){  in jnf()
151 			if(b>(float)1e10) {  in jnf()
169 float
170 ynf(int n, float x)  in ynf()
174 	float a, b, temp;  in ynf()
196 	    b = ((float)(i+i)/x)*b - a;  in ynf()