Home | History | Annotate | Line # | Download | only in test 
      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
     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