1 b2cfc4e2 2003-09-30 devnull /*
 2 b2cfc4e2 2003-09-30 devnull  * The authors of this software are Rob Pike and Ken Thompson.
 3 b2cfc4e2 2003-09-30 devnull  *              Copyright (c) 2002 by Lucent Technologies.
 4 b2cfc4e2 2003-09-30 devnull  * Permission to use, copy, modify, and distribute this software for any
 5 b2cfc4e2 2003-09-30 devnull  * purpose without fee is hereby granted, provided that this entire notice
 6 b2cfc4e2 2003-09-30 devnull  * is included in all copies of any software which is or includes a copy
 7 b2cfc4e2 2003-09-30 devnull  * or modification of this software and in all copies of the supporting
 8 b2cfc4e2 2003-09-30 devnull  * documentation for such software.
 9 b2cfc4e2 2003-09-30 devnull  * THIS SOFTWARE IS BEING PROVIDED "AS IS", WITHOUT ANY EXPRESS OR IMPLIED
10 b2cfc4e2 2003-09-30 devnull  * WARRANTY.  IN PARTICULAR, NEITHER THE AUTHORS NOR LUCENT TECHNOLOGIES MAKE ANY
11 b2cfc4e2 2003-09-30 devnull  * REPRESENTATION OR WARRANTY OF ANY KIND CONCERNING THE MERCHANTABILITY
12 b2cfc4e2 2003-09-30 devnull  * OF THIS SOFTWARE OR ITS FITNESS FOR ANY PARTICULAR PURPOSE.
13 b2cfc4e2 2003-09-30 devnull  */
14 b2cfc4e2 2003-09-30 devnull #include <stdarg.h>
15 b2cfc4e2 2003-09-30 devnull #include <string.h>
16 b2cfc4e2 2003-09-30 devnull #include "utf.h"
17 b2cfc4e2 2003-09-30 devnull #include "fmt.h"
18 b2cfc4e2 2003-09-30 devnull #include "fmtdef.h"
19 b2cfc4e2 2003-09-30 devnull 
20 b2cfc4e2 2003-09-30 devnull /*
21 b2cfc4e2 2003-09-30 devnull  * this table might overflow 127-bit exponent representations.
22 b2cfc4e2 2003-09-30 devnull  * in that case, truncate it after 1.0e38.
23 b2cfc4e2 2003-09-30 devnull  * it is important to get all one can from this
24 b2cfc4e2 2003-09-30 devnull  * routine since it is used in atof to scale numbers.
25 b2cfc4e2 2003-09-30 devnull  * the presumption is that C converts fp numbers better
26 b2cfc4e2 2003-09-30 devnull  * than multipication of lower powers of 10.
27 b2cfc4e2 2003-09-30 devnull  */
28 b2cfc4e2 2003-09-30 devnull 
29 b2cfc4e2 2003-09-30 devnull static
30 b2cfc4e2 2003-09-30 devnull double	tab[] =
31 b2cfc4e2 2003-09-30 devnull {
32 b2cfc4e2 2003-09-30 devnull 	1.0e0, 1.0e1, 1.0e2, 1.0e3, 1.0e4, 1.0e5, 1.0e6, 1.0e7, 1.0e8, 1.0e9,
33 b2cfc4e2 2003-09-30 devnull 	1.0e10,1.0e11,1.0e12,1.0e13,1.0e14,1.0e15,1.0e16,1.0e17,1.0e18,1.0e19,
34 b2cfc4e2 2003-09-30 devnull 	1.0e20,1.0e21,1.0e22,1.0e23,1.0e24,1.0e25,1.0e26,1.0e27,1.0e28,1.0e29,
35 b2cfc4e2 2003-09-30 devnull 	1.0e30,1.0e31,1.0e32,1.0e33,1.0e34,1.0e35,1.0e36,1.0e37,1.0e38,1.0e39,
36 b2cfc4e2 2003-09-30 devnull 	1.0e40,1.0e41,1.0e42,1.0e43,1.0e44,1.0e45,1.0e46,1.0e47,1.0e48,1.0e49,
37 b2cfc4e2 2003-09-30 devnull 	1.0e50,1.0e51,1.0e52,1.0e53,1.0e54,1.0e55,1.0e56,1.0e57,1.0e58,1.0e59,
38 b2cfc4e2 2003-09-30 devnull 	1.0e60,1.0e61,1.0e62,1.0e63,1.0e64,1.0e65,1.0e66,1.0e67,1.0e68,1.0e69,
39 b2cfc4e2 2003-09-30 devnull };
40 b2cfc4e2 2003-09-30 devnull 
41 b2cfc4e2 2003-09-30 devnull double
42 b2cfc4e2 2003-09-30 devnull __fmtpow10(int n)
43 b2cfc4e2 2003-09-30 devnull {
44 b2cfc4e2 2003-09-30 devnull 	int m;
45 b2cfc4e2 2003-09-30 devnull 
46 b2cfc4e2 2003-09-30 devnull 	if(n < 0) {
47 b2cfc4e2 2003-09-30 devnull 		n = -n;
48 b2cfc4e2 2003-09-30 devnull 		if(n < (int)(sizeof(tab)/sizeof(tab[0])))
49 b2cfc4e2 2003-09-30 devnull 			return 1/tab[n];
50 b2cfc4e2 2003-09-30 devnull 		m = n/2;
51 b2cfc4e2 2003-09-30 devnull 		return __fmtpow10(-m) * __fmtpow10(m-n);
52 b2cfc4e2 2003-09-30 devnull 	}
53 b2cfc4e2 2003-09-30 devnull 	if(n < (int)(sizeof(tab)/sizeof(tab[0])))
54 b2cfc4e2 2003-09-30 devnull 		return tab[n];
55 b2cfc4e2 2003-09-30 devnull 	m = n/2;
56 b2cfc4e2 2003-09-30 devnull 	return __fmtpow10(m) * __fmtpow10(n-m);
57 b2cfc4e2 2003-09-30 devnull }