1 /****************************************************************
2
3 The author of this software is David M. Gay.
4
5 Copyright (C) 2009 by David M. Gay
6 All Rights Reserved
7
8 Permission to use, copy, modify, and distribute this software and
9 its documentation for any purpose and without fee is hereby
10 granted, provided that the above copyright notice appear in all
11 source-code copies and that both that the copyright notice and this
12 permission notice and warranty disclaimer appear in supporting
13 documentation.
14
15 THE AUTHOR DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE,
16 INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS.
17 IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY SPECIAL, INDIRECT OR
18 CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF
19 USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR
20 OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR
21 PERFORMANCE OF THIS SOFTWARE.
22
23 ****************************************************************/
24 #include "stdio1.h"
25 #include "gdtoa.h"
26 #include <string.h>
27
28 #undef allow_Quad
29 #undef want_Quad
30 #undef want_Ux
31 #define want_LD
32 typedef union Ud {double x; unsigned int u[2]; } Ud;
33 #ifdef __x86_64 /*{{*/
34 #define want_Ux
35 #ifndef NO_GDTOA_i386_Quad /*{*/
36 typedef union UQ {__float128 x; unsigned int u[4]; } UQ;
37 #define allow_Quad(x) x
38 #define want_Quad
39 #endif /*}*/
40 #else /*}{*/
41 #ifdef __i386 /*{{*/
42 #define want_Ux
43 #else /*}{*/
44 #ifdef __sparc /*{{*/
45 typedef union UQ {long double x; unsigned int u[4]; } Ux;
46 #else /*}{*/
47 #ifdef __INTEL_COMPILER /*{*/
48 #undef want_Quad
49 #undef want_Ux
50 #undef want_LD
51 #endif /*}*/
52 #endif /*}}*/
53 #endif /*}}*/
54 #endif /*}}*/
55
56 #ifndef allow_Quad
57 #define allow_Quad(x) /*nothing*/
58 #endif
59
60 #ifdef want_Ux /*{{*/
61 typedef union Ux {long double x; unsigned short u[5]; } Ux;
62 #else /*}{*/
63 #ifdef __sparc
64 #define want_Ux
65 #endif
66 #endif /*}}*/
67
68 int
main(void)69 main(void)
70 {
71 Ud d;
72 allow_Quad(UQ q;)
73 char *b, buf[256], fmt[32], *s;
74 #ifdef want_Ux
75 Ux x;
76 x.x = 0.;
77 #endif
78 int k;
79
80 k = 0;
81 strcpy(fmt, "%.g");
82 d.x = 0.;
83 allow_Quad(q.x = 0.;)
84 while(fgets(buf, sizeof(buf), stdin)) {
85 for(b = buf; *b && *b != '\n'; ++b);
86 *b = 0;
87 if (b == buf)
88 continue;
89 b = buf;
90 if (*b == '%') {
91 for(k = 0; *b > ' '; ++b)
92 #ifdef want_LD /*{{*/
93 switch(*b) {
94 case 'L':
95 k = 1;
96 #ifdef want_Quad
97 break;
98 case 'q':
99 if (k >= 1)
100 k = 2;
101 #endif
102 }
103 #else /*}{*/
104 ;
105 #endif /*}}*/
106 if (*b)
107 *b++ = 0;
108 if (b - buf < sizeof(fmt)) {
109 strcpy(fmt, buf);
110 }
111 }
112 if (*b) {
113 switch(k) {
114 case 0:
115 d.x = strtod(b,&s);
116 break;
117 case 1:
118 #ifdef want_Ux
119 #ifdef __sparc
120 strtopQ(b,&s,&x.x);
121 #else
122 strtopx(b,&s,&x.x);
123 #endif
124 #else
125 strtopQ(b,&s,&q.x);
126 #endif
127 break;
128 allow_Quad(case 2: strtopQ(b,&s,&q.x);)
129 }
130 if (*s)
131 printf("Ignoring \"%s\"\n", s);
132 }
133 switch(k) {
134 case 0:
135 printf("d.x = %.g = #%x %x; %s ==> ", d.x, d.u[1], d.u[0], fmt);
136 printf(fmt, d.x);
137 break;
138 case 1:
139 #ifdef __sparc
140 printf("x.x = %.Lg = #%x %x %x %x; %s ==> ", x.x,
141 x.u[0], x.u[1], x.u[2], x.u[3], fmt);
142 #else
143 printf("x.x = %.Lg = #%x %x %x %x %x; %s ==> ", x.x,
144 x.u[4], x.u[3], x.u[2], x.u[1], x.u[0], fmt);
145 #endif
146 printf(fmt, x.x);
147 #ifdef want_Quad
148 break;
149 case 2:
150 printf("q.x = %.Lqg = #%x %x %x %x; %s ==> ", q.x,
151 q.u[3], q.u[2], q.u[1], q.u[0], fmt);
152 printf(fmt, q.x);
153 #endif
154 }
155 putchar('\n');
156 }
157 return 0;
158 }
159