d0p1/ack
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

several bug fixes after test-suite

Browse source
This commit is contained in:
eck 1990-03-29 09:05:21 +00:00
parent f27e9ca63e
commit 34d54a20b1
2 changed files with 13 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
lang/cem/libcc.ansi/math/exp.c
View file

@ -38,11 +38,11 @@ exp(double x)
if (x <= M_LN_MIN_D) { if (x <= M_LN_MIN_D) {
if (x < M_LN_MIN_D) errno = ERANGE; if (x < M_LN_MIN_D) errno = ERANGE;
return DBL_MIN; return 0.0;
} }
if (x >= M_LN_MAX_D) { if (x >= M_LN_MAX_D) {
if (x > M_LN_MAX_D) errno = ERANGE; if (x > M_LN_MAX_D) errno = ERANGE;
return DBL_MAX; return HUGE_VAL;
} }
if (negative) x = -x; if (negative) x = -x;

11
lang/cem/libcc.ansi/math/ldexp.c
View file

@ -5,6 +5,8 @@
/* $Header$ */ /* $Header$ */
#include <math.h> #include <math.h>
#include <float.h>
#include <errno.h>
double double
ldexp(double fl, int exp) ldexp(double fl, int exp)
@ -12,6 +14,7 @@ ldexp(double fl, int exp)
int sign = 1; int sign = 1;
int currexp; int currexp;
if (fl == 0.0) return 0.0;
if (fl<0) { if (fl<0) {
fl = -fl; fl = -fl;
sign = -1; sign = -1;
@ -19,6 +22,10 @@ ldexp(double fl, int exp)
fl = frexp(fl,&currexp); fl = frexp(fl,&currexp);
exp += currexp; exp += currexp;
if (exp > 0) { if (exp > 0) {
if (exp > DBL_MAX_EXP) {
errno = ERANGE;
return sign * HUGE_VAL;
}
while (exp>30) { while (exp>30) {
fl *= (double) (1L << 30); fl *= (double) (1L << 30);
exp -= 30; exp -= 30;
@ -26,6 +33,10 @@ ldexp(double fl, int exp)
fl *= (double) (1L << exp); fl *= (double) (1L << exp);
} }
else { else {
/* number need not be normalized */
if (exp < DBL_MIN_EXP - DBL_MANT_DIG) {
return 0.0;
}
while (exp<-30) { while (exp<-30) {
fl /= (double) (1L << 30); fl /= (double) (1L << 30);
exp += 30; exp += 30;