Loading...
--- Libc/Libc-320/stdio/hexfloat.c
+++ /dev/null
@@ -1,192 +0,0 @@
-/*
- * Copyright (c) 1999 Apple Computer, Inc. All rights reserved.
- *
- * @APPLE_LICENSE_HEADER_START@
- *
- * Copyright (c) 1999-2003 Apple Computer, Inc. All Rights Reserved.
- *
- * This file contains Original Code and/or Modifications of Original Code
- * as defined in and that are subject to the Apple Public Source License
- * Version 2.0 (the 'License'). You may not use this file except in
- * compliance with the License. Please obtain a copy of the License at
- * http://www.opensource.apple.com/apsl/ and read it before using this
- * file.
- *
- * The Original Code and all software distributed under the License are
- * distributed on an 'AS IS' basis, WITHOUT WARRANTY OF ANY KIND, EITHER
- * EXPRESS OR IMPLIED, AND APPLE HEREBY DISCLAIMS ALL SUCH WARRANTIES,
- * INCLUDING WITHOUT LIMITATION, ANY WARRANTIES OF MERCHANTABILITY,
- * FITNESS FOR A PARTICULAR PURPOSE, QUIET ENJOYMENT OR NON-INFRINGEMENT.
- * Please see the License for the specific language governing rights and
- * limitations under the License.
- *
- * @APPLE_LICENSE_HEADER_END@
- */
-
-#ifdef HEXFLOAT
-
-/*
- * hexfloat.c provides routines that vfprintf and vfwprintf can call to convert
- * floating point numbers to hexidecimal, as used by the %a and %A conversions.
- * This is necessarily dependent not only on the floating point data format,
- * but also byte order and existence of long double type.
- *
- * The union hexdouble represents the IEEE-754 double precision format, while
- * union hexlongdouble represents the IEEE-754 extended double precision
- * format.
- */
-
-#include <stdio.h>
-#include <math.h>
-#include <limits.h>
-
-#define EXPBIAS 1023
-#define EXPMIN -1022
-#define EXPSPECIAL 2047
-#define FRACTHEX 13
-
-union hexdouble {
- double d;
- struct {
-#if defined(__ppc__)
- unsigned int sign:1;
- unsigned int exp:11;
- unsigned long long fract:52;
-#elif defined(__i386__)
- unsigned long long fract:52;
- unsigned int exp:11;
- unsigned int sign:1;
-#else
-#error Unsupported architecture
-#endif
- } s;
-};
-
-#if !__TYPE_LONGDOUBLE_IS_DOUBLE
-#ifdef __i386__
-
-#define LEXPBIAS 16383
-#define LEXPMIN -16382
-#define LEXPSPECIAL 32767
-#define LFRACTHEX 16
-
-union hexlongdouble {
- long double d;
- struct {
- unsigned long long fract:63;
- unsigned int i:1;
- unsigned int exp:15;
- unsigned int sign:1;
- } s;
-};
-#endif /* __i386__ */
-#endif /* !__TYPE_LONGDOUBLE_IS_DOUBLE */
-
-int
-__hdtoa(double d, const char *xdigs, int prec, char *cp,
- int *expt, int *signflag, char **dtoaend)
-{
- union hexdouble u;
- char buf[FRACTHEX];
- char *hp;
- int i;
- long long fract;
- //char *start = cp; //DEBUG
-
- u.d = d;
- //printf("\nsign=%d exp=%x fract=%llx\n", u.s.sign, u.s.exp, u.s.fract); //DEBUG
- *signflag = u.s.sign;
- fract = u.s.fract;
- switch (u.s.exp) {
- case EXPSPECIAL: /* NaN or Inf */
- *expt = INT_MAX;
- *cp = (fract ? 'N' : 'I');
- return 0;
- case 0: /* Zero or denormalized */
- *cp++ = '0';
- *expt = (fract ? EXPMIN : 0);
- break;
- default: /* Normal numbers */
- *cp++ = '1';
- *expt = u.s.exp - EXPBIAS;
- break;
- }
- if (prec < 0)
- prec = FRACTHEX;
- //printf("prec=%d expt=%d\n", prec, *expt); //DEBUG
- if (prec > 0) {
- int dig = (prec > FRACTHEX ? FRACTHEX : prec);
- int zero = prec - dig;
- int shift = FRACTHEX - dig;
- if (shift > 0)
- fract >>= (shift << 2);
- for (hp = buf + dig, i = dig; i > 0; i--) {
- *--hp = xdigs[fract & 0xf];
- fract >>= 4;
- }
- strncpy(cp, hp, dig);
- cp += dig;
- while(zero-- > 0)
- *cp++ = '0';
- }
- *dtoaend = cp;
- //while (start < cp) putchar(*start++); //DEBUG
- //putchar('\n'); //DEBUG
- return prec;
-}
-
-#if !__TYPE_LONGDOUBLE_IS_DOUBLE
-#ifdef __i386__
-int
-__hldtoa(long double d, const char *xdigs, int prec, char *cp,
- int *expt, int *signflag, char **dtoaend)
-{
- union hexlongdouble u;
- char buf[LFRACTHEX];
- char *hp;
- int i;
- unsigned long long fract;
- //char *start = cp; //DEBUG
-
- u.d = d;
- //printf("d=%Lg u.d=%Lg\n", d, u.d); //DEBUG
- //printf("\nsign=%d exp=%x fract=%llx\n", u.s.sign, u.s.exp, u.s.fract); //DEBUG
- *signflag = u.s.sign;
- fract = (u.s.fract << 1);
- switch (u.s.exp) {
- case LEXPSPECIAL: /* NaN or Inf */
- *expt = INT_MAX;
- *cp = (fract ? 'N' : 'I');
- return 0;
- default: /* Normal or denormalized */
- *cp++ = u.s.i ? '1' : '0';
- *expt = u.s.exp - LEXPBIAS;
- break;
- }
- if (prec < 0)
- prec = LFRACTHEX;
- //printf("prec=%d expt=%d\n", prec, *expt); //DEBUG
- if (prec > 0) {
- int dig = (prec > LFRACTHEX ? LFRACTHEX : prec);
- int zero = prec - dig;
- int shift = LFRACTHEX - dig;
- if (shift > 0)
- fract >>= (shift << 2);
- for (hp = buf + dig, i = dig; i > 0; i--) {
- *--hp = xdigs[fract & 0xf];
- fract >>= 4;
- }
- strncpy(cp, hp, dig);
- cp += dig;
- while(zero-- > 0)
- *cp++ = '0';
- }
- *dtoaend = cp;
- //while (start < cp) putchar(*start++); //DEBUG
- //putchar('\n'); //DEBUG
- return prec;
-}
-#endif /* __i386__ */
-#endif /* !__TYPE_LONGDOUBLE_IS_DOUBLE */
-
-#endif /* HEXFLOAT */