Index: trunk/minix/lib/gnu/ieee_float/Makedepend-ack =================================================================== --- trunk/minix/lib/gnu/ieee_float/Makedepend-ack (revision 9) +++ (revision ) @@ -1,3 +1,0 @@ -depend-ack: - rm .depend-ack - touch .depend-ack Index: trunk/minix/lib/gnu/ieee_float/Makedepend-gnu =================================================================== --- trunk/minix/lib/gnu/ieee_float/Makedepend-gnu (revision 9) +++ (revision ) @@ -1,9 +1,0 @@ -depend-gnu: - rm .depend-gnu - touch .depend-gnu - mkdep 'gcc -O -I.. -E' frexp.c | sed -e '//d' -e '//d' -e 's:^\(.\):../../obj-gnu/./gnu/ieee_float/\1:' >> .depend-gnu - mkdep 'gcc -O -I.. -E' isinf.c | sed -e '//d' -e '//d' -e 's:^\(.\):../../obj-gnu/./gnu/ieee_float/\1:' >> .depend-gnu - mkdep 'gcc -O -I.. -E' isnan.c | sed -e '//d' -e '//d' -e 's:^\(.\):../../obj-gnu/./gnu/ieee_float/\1:' >> .depend-gnu - mkdep 'gcc -O -I.. -E' ldexp.c | sed -e '//d' -e '//d' -e 's:^\(.\):../../obj-gnu/./gnu/ieee_float/\1:' >> .depend-gnu - mkdep 'gcc -O -I.. -E' modf.c | sed -e '//d' -e '//d' -e 's:^\(.\):../../obj-gnu/./gnu/ieee_float/\1:' >> .depend-gnu - mkdep 'gcc -O -I.. -E' fphook.c | sed -e '//d' -e '//d' -e 's:^\(.\):../../obj-gnu/./gnu/ieee_float/\1:' >> .depend-gnu Index: trunk/minix/lib/gnu/ieee_float/Makefile =================================================================== --- trunk/minix/lib/gnu/ieee_float/Makefile (revision 9) +++ (revision ) @@ -1,60 +1,0 @@ -#Generated from ./gnu/ieee_float/Makefile.in -all: all-ack - -all-ack: -all-gnu: - -makefiles: Makefile -Makedepend-ack Makedepend-gnu: - sh ../../generate.sh ./gnu/ieee_float ../../obj-ack/ ../../obj-gnu - -Makefile: Makefile.in Makedepend-ack Makedepend-gnu - sh ../../generate.sh ./gnu/ieee_float ../../obj-ack/ ../../obj-gnu - @echo - @echo *Attention* - @echo Makefile is regenerated... rerun command to see changes - @echo *Attention* - @echo - -depend: depend-ack -all-gnu: ../../obj-gnu/libc.a - -../../obj-gnu/libc.a: ../../obj-gnu/./gnu/ieee_float/frexp.o -../../obj-gnu/libc.a: ../../obj-gnu/./gnu/ieee_float/isinf.o -../../obj-gnu/libc.a: ../../obj-gnu/./gnu/ieee_float/isnan.o -../../obj-gnu/libc.a: ../../obj-gnu/./gnu/ieee_float/ldexp.o -../../obj-gnu/libc.a: ../../obj-gnu/./gnu/ieee_float/modf.o -../../obj-gnu/libc.a: ../../obj-gnu/./gnu/ieee_float/fphook.o - -../../obj-gnu/libc.a: - gar cr ../../obj-gnu/libc.a $? - -../../obj-gnu/./gnu/ieee_float/frexp.o: frexp.c - gcc -O -I.. -c -o ../../obj-gnu/./gnu/ieee_float/frexp.o frexp.c - -../../obj-gnu/./gnu/ieee_float/isinf.o: isinf.c - gcc -O -I.. -c -o ../../obj-gnu/./gnu/ieee_float/isinf.o isinf.c - -../../obj-gnu/./gnu/ieee_float/isnan.o: isnan.c - gcc -O -I.. -c -o ../../obj-gnu/./gnu/ieee_float/isnan.o isnan.c - -../../obj-gnu/./gnu/ieee_float/ldexp.o: ldexp.c - gcc -O -I.. -c -o ../../obj-gnu/./gnu/ieee_float/ldexp.o ldexp.c - -../../obj-gnu/./gnu/ieee_float/modf.o: modf.c - gcc -O -I.. -c -o ../../obj-gnu/./gnu/ieee_float/modf.o modf.c - -../../obj-gnu/./gnu/ieee_float/fphook.o: fphook.c - gcc -O -I.. -c -o ../../obj-gnu/./gnu/ieee_float/fphook.o fphook.c - - - - -clean:: - rm -f ../../obj-gnu/./gnu/ieee_float/* - -include Makedepend-ack -include .depend-ack - -include Makedepend-gnu -include .depend-gnu Index: trunk/minix/lib/gnu/ieee_float/Makefile.in =================================================================== --- trunk/minix/lib/gnu/ieee_float/Makefile.in (revision 9) +++ (revision ) @@ -1,14 +1,0 @@ - -CFLAGS="-O -I.." - -LIBRARIES=libc - -libc_FILES=" \ - frexp.c \ - isinf.c \ - isnan.c \ - ldexp.c \ - modf.c \ - fphook.c" - -TYPE=gnu Index: trunk/minix/lib/gnu/ieee_float/fphook.c =================================================================== --- trunk/minix/lib/gnu/ieee_float/fphook.c (revision 9) +++ (revision ) @@ -1,195 +1,0 @@ -/* - * fltpr.c - print floating point numbers - */ -/* $Header: /cvsup/minix/src/lib/gnu/ieee_float/fphook.c,v 1.1 2005/10/10 15:28:15 beng Exp $ */ - -#ifndef NOFLOAT -#include -#include -#include "../stdio/loc_incl.h" -int _fp_hook = 1; - -static char * -_pfloat(long double r, register char *s, int n, int flags) -{ - register char *s1; - int sign, dp; - register int i; - - s1 = _fcvt(r, n, &dp, &sign); - if (sign) - *s++ = '-'; - else if (flags & FL_SIGN) - *s++ = '+'; - else if (flags & FL_SPACE) - *s++ = ' '; - - if (dp<=0) - *s++ = '0'; - for (i=dp; i>0; i--) - if (*s1) *s++ = *s1++; - else *s++ = '0'; - if (((i=n) > 0) || (flags & FL_ALT)) - *s++ = '.'; - while (++dp <= 0) { - if (--i<0) - break; - *s++ = '0'; - } - while (--i >= 0) - if (*s1) *s++ = *s1++; - else *s++ = '0'; - return s; -} - -static char * -_pscien(long double r, register char *s, int n, int flags) -{ - int sign, dp; - register char *s1; - - s1 = _ecvt(r, n + 1, &dp, &sign); - if (sign) - *s++ = '-'; - else if (flags & FL_SIGN) - *s++ = '+'; - else if (flags & FL_SPACE) - *s++ = ' '; - - *s++ = *s1++; - if ((n > 0) || (flags & FL_ALT)) - *s++ = '.'; - while (--n >= 0) - if (*s1) *s++ = *s1++; - else *s++ = '0'; - *s++ = 'e'; - if ( r != 0 ) --dp ; - if ( dp<0 ) { - *s++ = '-' ; dp= -dp ; - } else { - *s++ = '+' ; - } - if (dp >= 100) { - *s++ = '0' + (dp / 100); - dp %= 100; - } - *s++ = '0' + (dp/10); - *s++ = '0' + (dp%10); - return s; -} - -#define NDIGINEXP(exp) (((exp) >= 100 || (exp) <= -100) ? 3 : 2) -#define LOW_EXP -4 -#define USE_EXP(exp, ndigits) (((exp) < LOW_EXP + 1) || (exp >= ndigits + 1)) - -static char * -_gcvt(long double value, int ndigit, char *s, int flags) -{ - int sign, dp; - register char *s1, *s2; - register int i; - register int nndigit = ndigit; - - s1 = _ecvt(value, ndigit, &dp, &sign); - s2 = s; - if (sign) *s2++ = '-'; - else if (flags & FL_SIGN) - *s2++ = '+'; - else if (flags & FL_SPACE) - *s2++ = ' '; - - if (!(flags & FL_ALT)) - for (i = nndigit - 1; i > 0 && s1[i] == '0'; i--) - nndigit--; - - if (USE_EXP(dp,ndigit)) { - /* Use E format */ - dp--; - *s2++ = *s1++; - if ((nndigit > 1) || (flags & FL_ALT)) *s2++ = '.'; - while (--nndigit > 0) *s2++ = *s1++; - *s2++ = 'e'; - if (dp < 0) { - *s2++ = '-'; - dp = -dp; - } - else *s2++ = '+'; - s2 += NDIGINEXP(dp); - *s2 = 0; - for (i = NDIGINEXP(dp); i > 0; i--) { - *--s2 = dp % 10 + '0'; - dp /= 10; - } - return s; - } - /* Use f format */ - if (dp <= 0) { - if (*s1 != '0') { - /* otherwise the whole number is 0 */ - *s2++ = '0'; - *s2++ = '.'; - } - while (dp < 0) { - dp++; - *s2++ = '0'; - } - } - for (i = 1; i <= nndigit; i++) { - *s2++ = *s1++; - if (i == dp) *s2++ = '.'; - } - if (i <= dp) { - while (i++ <= dp) *s2++ = '0'; - *s2++ = '.'; - } - if ((s2[-1]=='.') && !(flags & FL_ALT)) s2--; - *s2 = '\0'; - return s; -} - -char * -_f_print(va_list *ap, int flags, char *s, char c, int precision) -{ - register char *old_s = s; - long double ld_val; - - if (flags & FL_LONGDOUBLE) ld_val = va_arg(*ap, long double); - else ld_val = (long double) va_arg(*ap, double); - - switch(c) { - case 'f': - s = _pfloat(ld_val, s, precision, flags); - break; - case 'e': - case 'E': - s = _pscien(ld_val, s, precision , flags); - break; - case 'g': - case 'G': - s = _gcvt(ld_val, precision, s, flags); - s += strlen(s); - break; - } - if ( c == 'E' || c == 'G') { - while (*old_s && *old_s != 'e') old_s++; - if (*old_s == 'e') *old_s = 'E'; - } - return s; -} -#endif /* NOFLOAT */ -/* $Header: /cvsup/minix/src/lib/gnu/ieee_float/fphook.c,v 1.1 2005/10/10 15:28:15 beng Exp $ */ - -#include -#include "../ansi/ext_fmt.h" - -void _str_ext_cvt(const char *s, char **ss, struct EXTEND *e); -double _ext_dbl_cvt(struct EXTEND *e); - -double -strtod(const char *p, char **pp) -{ - struct EXTEND e; - - _str_ext_cvt(p, pp, &e); - return _ext_dbl_cvt(&e); -} Index: trunk/minix/lib/gnu/ieee_float/frexp.c =================================================================== --- trunk/minix/lib/gnu/ieee_float/frexp.c (revision 9) +++ (revision ) @@ -1,60 +1,0 @@ -/* -libc/ieee_float/frexp.c - -Created: Oct 14, 1993 by Philip Homburg - -Implementation of frexp that directly manipulates the exponent bits in an -ieee float -*/ - -#include -#include - -#include "ieee_float.h" - -double frexp(value, eptr) -double value; -int *eptr; -{ - struct f64 *f64p; - int exp, exp_bias; - double factor; - - f64p= (struct f64 *)&value; - exp_bias= 0; - - exp= F64_GET_EXP(f64p); - if (exp == F64_EXP_MAX) - { /* Either infinity or Nan */ - *eptr= 0; - return value; - } - if (exp == 0) - { - /* Either 0 or denormal */ - if (F64_GET_MANT_LOW(f64p) == 0 && - F64_GET_MANT_HIGH(f64p) == 0) - { - *eptr= 0; - return value; - } - - /* Multiply by 2^64 */ - factor= 65536.0; /* 2^16 */ - factor *= factor; /* 2^32 */ - factor *= factor; /* 2^64 */ - value *= factor; - exp_bias= 64; - exp= F64_GET_EXP(f64p); - } - - exp= exp - F64_EXP_BIAS - exp_bias + 1; - *eptr= exp; - F64_SET_EXP(f64p, F64_EXP_BIAS-1); - - return value; -} - -/* - * $PchId: frexp.c,v 1.3 1996/02/22 21:01:39 philip Exp $ - */ Index: trunk/minix/lib/gnu/ieee_float/ieee_float.h =================================================================== --- trunk/minix/lib/gnu/ieee_float/ieee_float.h (revision 9) +++ (revision ) @@ -1,51 +1,0 @@ -/* -libc/ieee_float/ieee_float.h - -Created: Oct 14, 1993 by Philip Homburg - -Define structures and macros for manipulating IEEE floats -*/ - -#ifndef IEEE_FLOAT_H -#define IEEE_FLOAT_H - -#define isnan __IsNan - -struct f64 -{ - u32_t low_word; - u32_t high_word; -}; - -#define F64_SIGN_SHIFT 31 -#define F64_SIGN_MASK 1 - -#define F64_EXP_SHIFT 20 -#define F64_EXP_MASK 0x7ff -#define F64_EXP_BIAS 1023 -#define F64_EXP_MAX 2047 - -#define F64_MANT_SHIFT 0 -#define F64_MANT_MASK 0xfffff - -#define F64_GET_SIGN(fp) (((fp)->high_word >> F64_SIGN_SHIFT) & \ - F64_SIGN_MASK) -#define F64_GET_EXP(fp) (((fp)->high_word >> F64_EXP_SHIFT) & \ - F64_EXP_MASK) -#define F64_SET_EXP(fp, val) ((fp)->high_word= ((fp)->high_word & \ - ~(F64_EXP_MASK << F64_EXP_SHIFT)) | \ - (((val) & F64_EXP_MASK) << F64_EXP_SHIFT)) - -#define F64_GET_MANT_LOW(fp) ((fp)->low_word) -#define F64_SET_MANT_LOW(fp, val) ((fp)->low_word= (val)) -#define F64_GET_MANT_HIGH(fp) (((fp)->high_word >> F64_MANT_SHIFT) & \ - F64_MANT_MASK) -#define F64_SET_MANT_HIGH(fp, val) ((fp)->high_word= ((fp)->high_word & \ - ~(F64_MANT_MASK << F64_MANT_SHIFT)) | \ - (((val) & F64_MANT_MASK) << F64_MANT_SHIFT)) - -#endif /* IEEE_FLOAT_H */ - -/* - * $PchId: ieee_float.h,v 1.3 1996/02/22 21:01:39 philip Exp $ - */ Index: trunk/minix/lib/gnu/ieee_float/isinf.c =================================================================== --- trunk/minix/lib/gnu/ieee_float/isinf.c (revision 9) +++ (revision ) @@ -1,31 +1,0 @@ -/* -libc/ieee_float/isinf.c - -Created: Oct 14, 1993 by Philip Homburg - -Implementation of isinf that directly tests the bits in an ieee float -*/ - -#define _MINIX_SOURCE - -#include -#include - -#include "ieee_float.h" - -int isinf(value) -double value; -{ - struct f64 *f64p; - int exp; - - f64p= (struct f64 *)&value; - exp= F64_GET_EXP(f64p); - if (exp != F64_EXP_MAX) - return 0; - return F64_GET_MANT_LOW(f64p) == 0 && F64_GET_MANT_HIGH(f64p) == 0; -} - -/* - * $PchId: isinf.c,v 1.3 1996/02/22 21:01:39 philip Exp $ - */ Index: trunk/minix/lib/gnu/ieee_float/isnan.c =================================================================== --- trunk/minix/lib/gnu/ieee_float/isnan.c (revision 9) +++ (revision ) @@ -1,31 +1,0 @@ -/* -libc/ieee_float/isnan.c - -Created: Oct 14, 1993 by Philip Homburg - -Implementation of isnan that directly tests the bits in an ieee float -*/ - -#define _MINIX_SOURCE - -#include -#include - -#include "ieee_float.h" - -int isnan(value) -double value; -{ - struct f64 *f64p; - int exp; - - f64p= (struct f64 *)&value; - exp= F64_GET_EXP(f64p); - if (exp != F64_EXP_MAX) - return 0; - return F64_GET_MANT_LOW(f64p) != 0 || F64_GET_MANT_HIGH(f64p) != 0; -} - -/* - * $PchId: isnan.c,v 1.3 1996/02/22 21:01:39 philip Exp $ - */ Index: trunk/minix/lib/gnu/ieee_float/ldexp.c =================================================================== --- trunk/minix/lib/gnu/ieee_float/ldexp.c (revision 9) +++ (revision ) @@ -1,99 +1,0 @@ -/* -libc/ieee_float/ldexp.c - -Created: Oct 14, 1993 by Philip Homburg - -Implementation of ldexp that directly manipulates the exponent bits in an -ieee float -*/ - -#include -#include -#include - -#include "ieee_float.h" - -double ldexp(value, exp) -double value; -int exp; -{ - struct f64 *f64p; - int oldexp, exp_bias; - double factor; - - f64p= (struct f64 *)&value; - exp_bias= 0; - - oldexp= F64_GET_EXP(f64p); - if (oldexp == F64_EXP_MAX) - { /* Either infinity or Nan */ - return value; - } - if (oldexp == 0) - { - /* Either 0 or denormal */ - if (F64_GET_MANT_LOW(f64p) == 0 && - F64_GET_MANT_HIGH(f64p) == 0) - { - return value; - } - } - - /* If exp is too large (> 2*F64_EXP_MAX) or too small - * (< -2*F64_EXP_MAX) return HUGE_VAL or 0. This prevents overflows - * in exp if exp is really weird - */ - if (exp >= 2*F64_EXP_MAX) - { - errno= ERANGE; - return HUGE_VAL; - } - if (exp <= -2*F64_EXP_MAX) - { - errno= ERANGE; - return 0; - } - - /* Normalize a denormal */ - if (oldexp == 0) - { - /* Multiply by 2^64 */ - factor= 65536.0; /* 2^16 */ - factor *= factor; /* 2^32 */ - factor *= factor; /* 2^64 */ - value *= factor; - exp= -64; - oldexp= F64_GET_EXP(f64p); - } - - exp= oldexp + exp; - if (exp >= F64_EXP_MAX) - { /* Overflow */ - errno= ERANGE; - return HUGE_VAL; - } - if (exp > 0) - { - /* Normal */ - F64_SET_EXP(f64p, exp); - return value; - } - /* Denormal, or underflow. */ - exp += 64; - F64_SET_EXP(f64p, exp); - /* Divide by 2^64 */ - factor= 65536.0; /* 2^16 */ - factor *= factor; /* 2^32 */ - factor *= factor; /* 2^64 */ - value /= factor; - if (value == 0.0) - { - /* Underflow */ - errno= ERANGE; - } - return value; -} - -/* - * $PchId: ldexp.c,v 1.3 1996/02/22 21:01:39 philip Exp $ - */ Index: trunk/minix/lib/gnu/ieee_float/modf.c =================================================================== --- trunk/minix/lib/gnu/ieee_float/modf.c (revision 9) +++ (revision ) @@ -1,61 +1,0 @@ -/* -libc/ieee_float/modf.c - -Created: Oct 14, 1993 by Philip Homburg - -Implementation of modf that directly manipulates the exponent bits in an -ieee float -*/ - -#include -#include - -#include "ieee_float.h" - -double modf(value, iptr) -double value; -double *iptr; -{ - struct f64 *f64p; - double tmp; - int exp; - int mask_bits; - u32_t mant; - - f64p= (struct f64 *)&value; - - exp= F64_GET_EXP(f64p); - exp -= F64_EXP_BIAS; - if (exp < 0) - { - *iptr= 0; - return value; - } - mask_bits= 52-exp; - if (mask_bits <= 0) - { - *iptr= value; - return 0; - } - tmp= value; - if (mask_bits >= 32) - { - F64_SET_MANT_LOW(f64p, 0); - mask_bits -= 32; - mant= F64_GET_MANT_HIGH(f64p); - mant &= ~((1 << mask_bits)-1); - F64_SET_MANT_HIGH(f64p, mant); - } - else - { - mant= F64_GET_MANT_LOW(f64p); - mant &= ~((1 << mask_bits)-1); - F64_SET_MANT_LOW(f64p, mant); - } - *iptr= value; - return tmp-value; -} - -/* - * $PchId: modf.c,v 1.3 1996/02/22 21:01:39 philip Exp $ - */