lib/libc/gdtoa/test/pftest.c

   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         }