/****************************************************************

The author of this software is David M. Gay.

Copyright (C) 2009 by David M. Gay
All Rights Reserved

Permission to use, copy, modify, and distribute this software and
its documentation for any purpose and without fee is hereby
granted, provided that the above copyright notice appear in all
source-code copies and that both that the copyright notice and this
permission notice and warranty disclaimer appear in supporting
documentation.

THE AUTHOR DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE,
INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS.
IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY SPECIAL, INDIRECT OR
CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF
USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR
OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR
PERFORMANCE OF THIS SOFTWARE.

****************************************************************/
#include "stdio1.h"
#include "gdtoa.h"
#include <string.h>

#undef allow_Quad
#undef want_Quad
#undef want_Ux
#define want_LD
typedef union Ud {double x; unsigned int u[2]; } Ud;
#ifdef __x86_64 /*{{*/
#define want_Ux
#ifndef NO_GDTOA_i386_Quad /*{*/
typedef union UQ {__float128 x; unsigned int u[4]; } UQ;
#define allow_Quad(x) x
#define want_Quad
#endif /*}*/
#else /*}{*/
#ifdef __i386 /*{{*/
#define want_Ux
#else /*}{*/
#ifdef __sparc /*{{*/
typedef union UQ {long double x; unsigned int u[4]; } Ux;
#else /*}{*/
#ifdef __INTEL_COMPILER /*{*/
#undef want_Quad
#undef want_Ux
#undef want_LD
#endif /*}*/
#endif /*}}*/
#endif /*}}*/
#endif /*}}*/

#ifndef allow_Quad
#define allow_Quad(x) /*nothing*/
#endif

#ifdef want_Ux /*{{*/
typedef union Ux {long double x; unsigned short u[5]; } Ux;
#else /*}{*/
#ifdef __sparc
#define want_Ux
#endif
#endif /*}}*/

 int
main(void)
{
        Ud d;
        allow_Quad(UQ q;)
        char *b, buf[256], fmt[32], *s;
#ifdef want_Ux
        Ux x;
        x.x = 0.;
#endif
        int k;

        k = 0;
        strcpy(fmt, "%.g");
        d.x = 0.;
        allow_Quad(q.x = 0.;)
        while(fgets(buf, sizeof(buf), stdin)) {
                for(b = buf; *b && *b != '\n'; ++b);
                *b = 0;
                if (b == buf)
                        continue;
                b = buf;
                if (*b == '%') {
                        for(k = 0; *b > ' '; ++b)
#ifdef want_LD /*{{*/
                                switch(*b) {
                                  case 'L':
                                        k = 1;
#ifdef want_Quad
                                        break;
                                  case 'q':
                                        if (k >= 1)
                                                k = 2;
#endif
                                  }
#else /*}{*/
                                ;
#endif /*}}*/
                        if (*b)
                                *b++ = 0;
                        if (b - buf < sizeof(fmt)) {
                                strcpy(fmt, buf);
                                }
                        }
                if (*b) {
                        switch(k) {
                          case 0:
                                d.x = strtod(b,&s);
                                break;
                          case 1:
#ifdef want_Ux
#ifdef __sparc
                                strtopQ(b,&s,&x.x);
#else
                                strtopx(b,&s,&x.x);
#endif
#else
                                strtopQ(b,&s,&q.x);
#endif
                                break;
                          allow_Quad(case 2: strtopQ(b,&s,&q.x);)
                          }
                        if (*s)
                                printf("Ignoring \"%s\"\n", s);
                        }
                switch(k) {
                        case 0:
                                printf("d.x = %.g = #%x %x; %s ==> ", d.x, d.u[1], d.u[0], fmt);
                                printf(fmt, d.x);
                                break;
                        case 1:
#ifdef __sparc
                                printf("x.x = %.Lg = #%x %x %x %x; %s ==> ", x.x,
                                        x.u[0], x.u[1], x.u[2], x.u[3], fmt);
#else
                                printf("x.x = %.Lg = #%x %x %x %x %x; %s ==> ", x.x,
                                        x.u[4], x.u[3], x.u[2], x.u[1], x.u[0], fmt);
#endif
                                printf(fmt, x.x);
#ifdef want_Quad
                                break;
                        case 2:
                                printf("q.x = %.Lqg = #%x %x %x %x; %s ==> ", q.x,
                                        q.u[3], q.u[2], q.u[1], q.u[0], fmt);
                                printf(fmt, q.x);
#endif
                        }
                putchar('\n');
                }
        return 0;
        }