1 a84cbb2a 2004-04-19 devnull #include <u.h>
2 a84cbb2a 2004-04-19 devnull #include <libc.h>
3 a84cbb2a 2004-04-19 devnull #include <bio.h>
4 a84cbb2a 2004-04-19 devnull #include <mach.h>
7 a84cbb2a 2004-04-19 devnull * These routines assume that if the number is representable
8 a84cbb2a 2004-04-19 devnull * in IEEE floating point, it will be representable in the native
9 a84cbb2a 2004-04-19 devnull * double format. Naive but workable, probably.
12 a84cbb2a 2004-04-19 devnull ieeeftoa64(char *buf, uint n, u32int h, u32int l)
14 a84cbb2a 2004-04-19 devnull double fr;
17 a84cbb2a 2004-04-19 devnull if (n <= 0)
18 a84cbb2a 2004-04-19 devnull return 0;
21 929ba9b9 2004-04-20 devnull if(h & (1UL<<31)){
22 a84cbb2a 2004-04-19 devnull *buf++ = '-';
23 929ba9b9 2004-04-20 devnull h &= ~(1UL<<31);
25 a84cbb2a 2004-04-19 devnull *buf++ = ' ';
27 a84cbb2a 2004-04-19 devnull if(l == 0 && h == 0)
28 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "0.");
29 a84cbb2a 2004-04-19 devnull exp = (h>>20) & ((1L<<11)-1L);
30 a84cbb2a 2004-04-19 devnull if(exp == 0)
31 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "DeN(%.8lux%.8lux)", h, l);
32 a84cbb2a 2004-04-19 devnull if(exp == ((1L<<11)-1L)){
33 a84cbb2a 2004-04-19 devnull if(l==0 && (h&((1L<<20)-1L)) == 0)
34 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "Inf");
36 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "NaN(%.8lux%.8lux)", h&((1L<<20)-1L), l);
38 a84cbb2a 2004-04-19 devnull exp -= (1L<<10) - 2L;
39 a84cbb2a 2004-04-19 devnull fr = l & ((1L<<16)-1L);
40 a84cbb2a 2004-04-19 devnull fr /= 1L<<16;
41 a84cbb2a 2004-04-19 devnull fr += (l>>16) & ((1L<<16)-1L);
42 a84cbb2a 2004-04-19 devnull fr /= 1L<<16;
43 a84cbb2a 2004-04-19 devnull fr += (h & (1L<<20)-1L) | (1L<<20);
44 a84cbb2a 2004-04-19 devnull fr /= 1L<<21;
45 a84cbb2a 2004-04-19 devnull fr = ldexp(fr, exp);
46 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "%.18g", fr);
50 a84cbb2a 2004-04-19 devnull ieeeftoa32(char *buf, uint n, u32int h)
52 a84cbb2a 2004-04-19 devnull double fr;
55 a84cbb2a 2004-04-19 devnull if (n <= 0)
56 a84cbb2a 2004-04-19 devnull return 0;
58 929ba9b9 2004-04-20 devnull if(h & (1UL<<31)){
59 a84cbb2a 2004-04-19 devnull *buf++ = '-';
60 929ba9b9 2004-04-20 devnull h &= ~(1UL<<31);
62 a84cbb2a 2004-04-19 devnull *buf++ = ' ';
64 a84cbb2a 2004-04-19 devnull if(h == 0)
65 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "0.");
66 a84cbb2a 2004-04-19 devnull exp = (h>>23) & ((1L<<8)-1L);
67 a84cbb2a 2004-04-19 devnull if(exp == 0)
68 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "DeN(%.8lux)", h);
69 a84cbb2a 2004-04-19 devnull if(exp == ((1L<<8)-1L)){
70 a84cbb2a 2004-04-19 devnull if((h&((1L<<23)-1L)) == 0)
71 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "Inf");
73 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "NaN(%.8lux)", h&((1L<<23)-1L));
75 a84cbb2a 2004-04-19 devnull exp -= (1L<<7) - 2L;
76 a84cbb2a 2004-04-19 devnull fr = (h & ((1L<<23)-1L)) | (1L<<23);
77 a84cbb2a 2004-04-19 devnull fr /= 1L<<24;
78 a84cbb2a 2004-04-19 devnull fr = ldexp(fr, exp);
79 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "%.9g", fr);
83 a84cbb2a 2004-04-19 devnull beieeeftoa32(char *buf, uint n, void *s)
85 a84cbb2a 2004-04-19 devnull return ieeeftoa32(buf, n, beswap4(*(u32int*)s));
89 a84cbb2a 2004-04-19 devnull beieeeftoa64(char *buf, uint n, void *s)
91 a84cbb2a 2004-04-19 devnull return ieeeftoa64(buf, n, beswap4(*(u32int*)s), beswap4(((u32int*)(s))[1]));
95 a84cbb2a 2004-04-19 devnull leieeeftoa32(char *buf, uint n, void *s)
97 a84cbb2a 2004-04-19 devnull return ieeeftoa32(buf, n, leswap4(*(u32int*)s));
101 a84cbb2a 2004-04-19 devnull leieeeftoa64(char *buf, uint n, void *s)
103 a84cbb2a 2004-04-19 devnull return ieeeftoa64(buf, n, leswap4(((u32int*)(s))[1]), leswap4(*(u32int*)s));
106 a84cbb2a 2004-04-19 devnull /* packed in 12 bytes, with s[2]==s[3]==0; mantissa starts at s[4]*/
108 a84cbb2a 2004-04-19 devnull beieeeftoa80(char *buf, uint n, void *s)
110 a84cbb2a 2004-04-19 devnull uchar *reg = (uchar*)s;
112 a84cbb2a 2004-04-19 devnull ulong x;
113 a84cbb2a 2004-04-19 devnull uchar ieee[8+8]; /* room for slop */
114 a84cbb2a 2004-04-19 devnull uchar *p, *q;
116 a84cbb2a 2004-04-19 devnull memset(ieee, 0, sizeof(ieee));
117 a84cbb2a 2004-04-19 devnull /* sign */
118 a84cbb2a 2004-04-19 devnull if(reg[0] & 0x80)
119 a84cbb2a 2004-04-19 devnull ieee[0] |= 0x80;
121 a84cbb2a 2004-04-19 devnull /* exponent */
122 a84cbb2a 2004-04-19 devnull x = ((reg[0]&0x7F)<<8) | reg[1];
123 a84cbb2a 2004-04-19 devnull if(x == 0) /* number is ±0 */
124 a84cbb2a 2004-04-19 devnull goto done;
125 a84cbb2a 2004-04-19 devnull if(x == 0x7FFF){
126 a84cbb2a 2004-04-19 devnull if(memcmp(reg+4, ieee+1, 8) == 0){ /* infinity */
127 a84cbb2a 2004-04-19 devnull x = 2047;
128 a84cbb2a 2004-04-19 devnull }else{ /* NaN */
129 a84cbb2a 2004-04-19 devnull x = 2047;
130 a84cbb2a 2004-04-19 devnull ieee[7] = 0x1; /* make sure */
132 a84cbb2a 2004-04-19 devnull ieee[0] |= x>>4;
133 a84cbb2a 2004-04-19 devnull ieee[1] |= (x&0xF)<<4;
134 a84cbb2a 2004-04-19 devnull goto done;
136 a84cbb2a 2004-04-19 devnull x -= 0x3FFF; /* exponent bias */
137 a84cbb2a 2004-04-19 devnull x += 1023;
138 a84cbb2a 2004-04-19 devnull if(x >= (1<<11) || ((reg[4]&0x80)==0 && x!=0))
139 a84cbb2a 2004-04-19 devnull return snprint(buf, n, "not in range");
140 a84cbb2a 2004-04-19 devnull ieee[0] |= x>>4;
141 a84cbb2a 2004-04-19 devnull ieee[1] |= (x&0xF)<<4;
143 a84cbb2a 2004-04-19 devnull /* mantissa */
144 a84cbb2a 2004-04-19 devnull p = reg+4;
145 a84cbb2a 2004-04-19 devnull q = ieee+1;
146 a84cbb2a 2004-04-19 devnull for(i=0; i<56; i+=8, p++, q++){ /* move one byte */
147 a84cbb2a 2004-04-19 devnull x = (p[0]&0x7F) << 1;
148 a84cbb2a 2004-04-19 devnull if(p[1] & 0x80)
150 a84cbb2a 2004-04-19 devnull q[0] |= x>>4;
151 a84cbb2a 2004-04-19 devnull q[1] |= (x&0xF)<<4;
154 a84cbb2a 2004-04-19 devnull return beieeeftoa64(buf, n, (void*)ieee);
159 a84cbb2a 2004-04-19 devnull leieeeftoa80(char *buf, uint n, void *s)
162 a84cbb2a 2004-04-19 devnull char *cp;
163 a84cbb2a 2004-04-19 devnull char b[12];
165 a84cbb2a 2004-04-19 devnull cp = (char*) s;
166 a84cbb2a 2004-04-19 devnull for(i=0; i<12; i++)
167 a84cbb2a 2004-04-19 devnull b[11-i] = *cp++;
168 a84cbb2a 2004-04-19 devnull return beieeeftoa80(buf, n, b);