diff options
| author | Jeff Johnston <jjohnstn@redhat.com> | 2000-05-26 22:42:39 +0000 |
|---|---|---|
| committer | Jeff Johnston <jjohnstn@redhat.com> | 2000-05-26 22:42:39 +0000 |
| commit | 8887021868a19a9feb8d4cdb22148823ee42579b (patch) | |
| tree | c3ac1cbf18fa5bd246106ad5d916cb7c03c048b9 /newlib/libm/math/ef_j0.c | |
| parent | e8844e5d9f6a07f7846ad0766cc3ef6e08dce8f5 (diff) | |
| download | cygnal-8887021868a19a9feb8d4cdb22148823ee42579b.tar.gz cygnal-8887021868a19a9feb8d4cdb22148823ee42579b.tar.bz2 cygnal-8887021868a19a9feb8d4cdb22148823ee42579b.zip | |
2000-05-26 Marek Michalkiewicz <marekm@linux.org.pl>
* libm/common/s_expm1.c (expm1): Add curly braces, avoid warnings.
* libm/common/s_log1p.c (log1p): Likewise.
* libm/common/s_scalbn.c (scalbn): Likewise.
* libm/math/e_log.c: Likewise.
* libm/math/e_asin.c: Likewise.
* libm/math/ef_asin.c: Likewise.
* libm/math/e_j0.c (pzero, qzero): Remove redundant test.
* libm/math/e_j1.c (pone, qone): Likewise.
* libm/math/ef_j0.c (pzerof, qzerof): Likewise.
* libm/math/ef_j1.c (ponef, qonef): Likewise.
* libm/mathfp/e_j0.c (pzero, qzero): Likewise.
* libm/mathfp/e_j1.c (pone, qone): Likewise.
* libm/mathfp/ef_j0.c (pzerof, qzerof): Likewise.
* libm/mathfp/ef_j1.c (ponef, qonef): Likewise.
Diffstat (limited to 'newlib/libm/math/ef_j0.c')
| -rw-r--r-- | newlib/libm/math/ef_j0.c | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/newlib/libm/math/ef_j0.c b/newlib/libm/math/ef_j0.c index 5ae6f308f..e7ee6e748 100644 --- a/newlib/libm/math/ef_j0.c +++ b/newlib/libm/math/ef_j0.c @@ -295,7 +295,7 @@ static float pS2[5] = { if(ix>=0x41000000) {p = pR8; q= pS8;} else if(ix>=0x40f71c58){p = pR5; q= pS5;} else if(ix>=0x4036db68){p = pR3; q= pS3;} - else if(ix>=0x40000000){p = pR2; q= pS2;} + else {p = pR2; q= pS2;} z = one/(x*x); r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5])))); s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*q[4])))); @@ -305,7 +305,7 @@ static float pS2[5] = { /* For x >= 8, the asymptotic expansions of qzero is * -1/8 s + 75/1024 s^3 - ..., where s = 1/x. - * We approximate pzero by + * We approximate qzero by * qzero(x) = s*(-1.25 + (R/S)) * where R = qR0 + qR1*s^2 + qR2*s^4 + ... + qR5*s^10 * S = 1 + qS0*s^2 + ... + qS5*s^12 @@ -431,7 +431,7 @@ static float qS2[6] = { if(ix>=0x41000000) {p = qR8; q= qS8;} else if(ix>=0x40f71c58){p = qR5; q= qS5;} else if(ix>=0x4036db68){p = qR3; q= qS3;} - else if(ix>=0x40000000){p = qR2; q= qS2;} + else {p = qR2; q= qS2;} z = one/(x*x); r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5])))); s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*(q[4]+z*q[5]))))); |