diff options
author | Rich Felker <dalias@aerifal.cx> | 2012-03-16 21:01:34 -0400 |
---|---|---|
committer | Rich Felker <dalias@aerifal.cx> | 2012-03-16 21:01:34 -0400 |
commit | 93a50a26cd0f9efc59cc83daae7b2d916b327ab1 (patch) | |
tree | ada41a23c4b01f27d1d2cf8cd88948b061aa8f61 /src/math/jnf.c | |
parent | 2cbb24bba39ad3529315098b5619b6fec078eb82 (diff) | |
parent | 40305f74bd70a575ce73260227ed3b64e0625b13 (diff) | |
download | musl-93a50a26cd0f9efc59cc83daae7b2d916b327ab1.tar.gz musl-93a50a26cd0f9efc59cc83daae7b2d916b327ab1.tar.xz musl-93a50a26cd0f9efc59cc83daae7b2d916b327ab1.zip |
Merge remote branch 'nsz/master'
Diffstat (limited to 'src/math/jnf.c')
-rw-r--r-- | src/math/jnf.c | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/src/math/jnf.c b/src/math/jnf.c index 7db93ae7..648db32b 100644 --- a/src/math/jnf.c +++ b/src/math/jnf.c @@ -64,7 +64,7 @@ float jnf(int n, float x) if (n > 33) /* underflow */ b = zero; else { - temp = x*(float)0.5; + temp = 0.5f * x; b = temp; for (a=one,i=2; i<=n; i++) { a *= (float)i; /* a = n! */ @@ -106,13 +106,13 @@ float jnf(int n, float x) float q0,q1,h,tmp; int32_t k,m; - w = (n+n)/(float)x; - h = (float)2.0/(float)x; + w = (n+n)/x; + h = 2.0f/x; z = w+h; q0 = w; - q1 = w*z - (float)1.0; + q1 = w*z - 1.0f; k = 1; - while (q1 < (float)1.0e9) { + while (q1 < 1.0e9f) { k += 1; z += h; tmp = z*q1 - q0; @@ -135,7 +135,7 @@ float jnf(int n, float x) tmp = n; v = two/x; tmp = tmp*logf(fabsf(v*tmp)); - if (tmp < (float)8.8721679688e+01) { + if (tmp < 88.721679688f) { for (i=n-1,di=(float)(i+i); i>0; i--) { temp = b; b *= di; @@ -151,7 +151,7 @@ float jnf(int n, float x) a = temp; di -= two; /* scale b to avoid spurious overflow */ - if (b > (float)1e10) { + if (b > 1e10f) { a /= b; t /= b; b = one; |