pintos-os.org Git - pspp/blob - src/data/data-out.c

   1 /* PSPP - a program for statistical analysis.
   2    Copyright (C) 1997-9, 2000, 2006, 2009, 2011 Free Software Foundation, Inc.
   3
   4    This program is free software: you can redistribute it and/or modify
   5    it under the terms of the GNU General Public License as published by
   6    the Free Software Foundation, either version 3 of the License, or
   7    (at your option) any later version.
   8
   9    This program is distributed in the hope that it will be useful,
  10    but WITHOUT ANY WARRANTY; without even the implied warranty of
  11    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12    GNU General Public License for more details.
  13
  14    You should have received a copy of the GNU General Public License
  15    along with this program.  If not, see <http://www.gnu.org/licenses/>. */
  16
  17 #include <config.h>
  18
  19 #include "data/data-out.h"
  20
  21 #include <ctype.h>
  22 #include <float.h>
  23 #include <math.h>
  24 #include <stdint.h>
  25 #include <stdlib.h>
  26 #include <time.h>
  27 #include <unistr.h>
  28
  29 #include "data/calendar.h"
  30 #include "data/format.h"
  31 #include "data/settings.h"
  32 #include "data/value.h"
  33 #include "libpspp/assertion.h"
  34 #include "libpspp/cast.h"
  35 #include "libpspp/float-format.h"
  36 #include "libpspp/i18n.h"
  37 #include "libpspp/integer-format.h"
  38 #include "libpspp/message.h"
  39 #include "libpspp/misc.h"
  40 #include "libpspp/pool.h"
  41 #include "libpspp/str.h"
  42
  43 #include "gl/minmax.h"
  44
  45 #include "gettext.h"
  46 #define _(msgid) gettext (msgid)
  47 \f
  48 /* A representation of a number that can be quickly rounded to
  49    any desired number of decimal places (up to a specified
  50    maximum). */
  51 struct rounder
  52   {
  53     char string[64];    /* Magnitude of number with excess precision. */
  54     int integer_digits; /* Number of digits before decimal point. */
  55     int leading_nines;  /* Number of `9's or `.'s at start of string. */
  56     int leading_zeros;  /* Number of `0's or `.'s at start of string. */
  57     bool negative;      /* Is the number negative? */
  58   };
  59
  60 static void rounder_init (struct rounder *, double number, int max_decimals);
  61 static int rounder_width (const struct rounder *, int decimals,
  62                           int *integer_digits, bool *negative);
  63 static void rounder_format (const struct rounder *, int decimals,
  64                             char *output);
  65 \f
  66 typedef void data_out_converter_func (const union value *,
  67                                       const struct fmt_spec *,
  68                                       char *);
  69 #define FMT(NAME, METHOD, IMIN, OMIN, IO, CATEGORY) \
  70         static data_out_converter_func output_##METHOD;
  71 #include "format.def"
  72
  73 static bool output_decimal (const struct rounder *, const struct fmt_spec *,
  74                             bool require_affixes, char *);
  75 static bool output_scientific (double, const struct fmt_spec *,
  76                                bool require_affixes, char *);
  77
  78 static double power10 (int) PURE_FUNCTION;
  79 static double power256 (int) PURE_FUNCTION;
  80
  81 static void output_infinite (double, const struct fmt_spec *, char *);
  82 static void output_missing (const struct fmt_spec *, char *);
  83 static void output_overflow (const struct fmt_spec *, char *);
  84 static bool output_bcd_integer (double, int digits, char *);
  85 static void output_binary_integer (uint64_t, int bytes, enum integer_format,
  86                                    char *);
  87 static void output_hex (const void *, size_t bytes, char *);
  88 \f
  89
  90 static data_out_converter_func *const converters[FMT_NUMBER_OF_FORMATS] =
  91     {
  92 #define FMT(NAME, METHOD, IMIN, OMIN, IO, CATEGORY) output_##METHOD,
  93 #include "format.def"
  94     };
  95
  96 /* Converts the INPUT value, encoded in INPUT_ENCODING, according to format
  97    specification FORMAT, appending the output to OUTPUT in OUTPUT_ENCODING.
  98    However, binary formats (FMT_P, FMT_PK, FMT_IB, FMT_PIB, FMT_RB) yield the
  99    binary results, which may not be properly encoded for OUTPUT_ENCODING.
 100
 101    VALUE must be the correct width for FORMAT, that is, its width must be
 102    fmt_var_width(FORMAT).
 103
 104    INPUT_ENCODING can normally be obtained by calling dict_get_encoding() on
 105    the dictionary with which INPUT is associated.  ENCODING is only important
 106    when FORMAT's type is FMT_A. */
 107 void
 108 data_out_recode (const union value *input, const char *input_encoding,
 109                  const struct fmt_spec *format,
 110                  struct string *output, const char *output_encoding)
 111 {
 112   assert (fmt_check_output (format));
 113   if (format->type == FMT_A)
 114     {
 115       char *in = CHAR_CAST (char *, value_str (input, format->w));
 116       char *out = recode_string (output_encoding, input_encoding,
 117                                  in, format->w);
 118       ds_put_cstr (output, out);
 119       free (out);
 120     }
 121   else if (fmt_get_category (format->type) == FMT_CAT_BINARY)
 122     converters[format->type] (input, format,
 123                               ds_put_uninit (output, format->w));
 124   else
 125     {
 126       char *utf8_encoded = data_out (input, input_encoding, format);
 127       char *output_encoded = recode_string (output_encoding, UTF8,
 128                                             utf8_encoded, -1);
 129       ds_put_cstr (output, output_encoded);
 130       free (output_encoded);
 131       free (utf8_encoded);
 132     }
 133 }
 134
 135 static char *
 136 binary_to_utf8 (const char *in, struct pool *pool)
 137 {
 138   uint8_t *out = pool_alloc_unaligned (pool, strlen (in) * 2 + 1);
 139   uint8_t *p = out;
 140
 141   while (*in != '\0')
 142     {
 143       uint8_t byte = *in++;
 144       int mblen = u8_uctomb (p, byte, 2);
 145       assert (mblen > 0);
 146       p += mblen;
 147     }
 148   *p = '\0';
 149
 150   return CHAR_CAST (char *, out);
 151 }
 152
 153 /* Converts the INPUT value into a UTF-8 encoded string, according to format
 154    specification FORMAT.
 155
 156    VALUE must be the correct width for FORMAT, that is, its width must be
 157    fmt_var_width(FORMAT).
 158
 159    ENCODING must be the encoding of INPUT.  Normally this can be obtained by
 160    calling dict_get_encoding() on the dictionary with which INPUT is
 161    associated.  ENCODING is only important when FORMAT's type is FMT_A.
 162
 163    The return value is dynamically allocated, and must be freed by the caller.
 164    If POOL is non-null, then the return value is allocated on that pool.  */
 165 char *
 166 data_out_pool (const union value *input, const char *encoding,
 167                const struct fmt_spec *format, struct pool *pool)
 168 {
 169   assert (fmt_check_output (format));
 170   if (format->type == FMT_A)
 171     {
 172       char *in = CHAR_CAST (char *, value_str (input, format->w));
 173       return recode_string_pool (UTF8, encoding, in, format->w, pool);
 174     }
 175   else if (fmt_get_category (format->type) == FMT_CAT_BINARY)
 176     {
 177       char tmp[16];
 178
 179       assert (format->w + 1 <= sizeof tmp);
 180       converters[format->type] (input, format, tmp);
 181       return binary_to_utf8 (tmp, pool);
 182     }
 183   else
 184     {
 185       const struct fmt_number_style *style = settings_get_style (format->type);
 186       size_t size = format->w + style->extra_bytes + 1;
 187       char *output;
 188
 189       output = pool_alloc_unaligned (pool, size);
 190       converters[format->type] (input, format, output);
 191       return output;
 192     }
 193 }
 194
 195 char *
 196 data_out (const union value *input, const char *encoding, const struct fmt_spec *format)
 197 {
 198   return data_out_pool (input, encoding, format, NULL);
 199 }
 200
 201 \f
 202 /* Main conversion functions. */
 203
 204 /* Outputs F, COMMA, DOT, DOLLAR, PCT, E, CCA, CCB, CCC, CCD, and
 205    CCE formats. */
 206 static void
 207 output_number (const union value *input, const struct fmt_spec *format,
 208                char *output)
 209 {
 210   double number = input->f;
 211
 212   if (number == SYSMIS)
 213     output_missing (format, output);
 214   else if (!isfinite (number))
 215     output_infinite (number, format, output);
 216   else
 217     {
 218       if (format->type != FMT_E && fabs (number) < 1.5 * power10 (format->w))
 219         {
 220           struct rounder r;
 221           rounder_init (&r, number, format->d);
 222
 223           if (output_decimal (&r, format, true, output)
 224               || output_scientific (number, format, true, output)
 225               || output_decimal (&r, format, false, output))
 226             return;
 227         }
 228
 229       if (!output_scientific (number, format, false, output))
 230         output_overflow (format, output);
 231     }
 232 }
 233
 234 /* Outputs N format. */
 235 static void
 236 output_N (const union value *input, const struct fmt_spec *format,
 237           char *output)
 238 {
 239   double number = input->f * power10 (format->d);
 240   if (input->f == SYSMIS || number < 0)
 241     output_missing (format, output);
 242   else
 243     {
 244       char buf[128];
 245       number = fabs (round (number));
 246       if (number < power10 (format->w)
 247           && sprintf (buf, "%0*.0f", format->w, number) == format->w)
 248         memcpy (output, buf, format->w);
 249       else
 250         output_overflow (format, output);
 251     }
 252
 253   output[format->w] = '\0';
 254 }
 255
 256 /* Outputs Z format. */
 257 static void
 258 output_Z (const union value *input, const struct fmt_spec *format,
 259           char *output)
 260 {
 261   double number = input->f * power10 (format->d);
 262   char buf[128];
 263   if (input->f == SYSMIS)
 264     output_missing (format, output);
 265   else if (fabs (number) < power10 (format->w)
 266            && sprintf (buf, "%0*.0f", format->w,
 267                        fabs (round (number))) == format->w)
 268     {
 269       if (number < 0 && strspn (buf, "0") < format->w)
 270         {
 271           char *p = &buf[format->w - 1];
 272           *p = "}JKLMNOPQR"[*p - '0'];
 273         }
 274       memcpy (output, buf, format->w);
 275       output[format->w] = '\0';
 276     }
 277   else
 278     output_overflow (format, output);
 279 }
 280
 281 /* Outputs P format. */
 282 static void
 283 output_P (const union value *input, const struct fmt_spec *format,
 284           char *output)
 285 {
 286   if (output_bcd_integer (fabs (input->f * power10 (format->d)),
 287                           format->w * 2 - 1, output)
 288       && input->f < 0.0)
 289     output[format->w - 1] |= 0xd;
 290   else
 291     output[format->w - 1] |= 0xf;
 292 }
 293
 294 /* Outputs PK format. */
 295 static void
 296 output_PK (const union value *input, const struct fmt_spec *format,
 297            char *output)
 298 {
 299   output_bcd_integer (input->f * power10 (format->d), format->w * 2, output);
 300 }
 301
 302 /* Outputs IB format. */
 303 static void
 304 output_IB (const union value *input, const struct fmt_spec *format,
 305            char *output)
 306 {
 307   double number = round (input->f * power10 (format->d));
 308   if (input->f == SYSMIS
 309       || number >= power256 (format->w) / 2 - 1
 310       || number < -power256 (format->w) / 2)
 311     memset (output, 0, format->w);
 312   else
 313     {
 314       uint64_t integer = fabs (number);
 315       if (number < 0)
 316         integer = -integer;
 317       output_binary_integer (integer, format->w,
 318                              settings_get_output_integer_format (),
 319                              output);
 320     }
 321
 322   output[format->w] = '\0';
 323 }
 324
 325 /* Outputs PIB format. */
 326 static void
 327 output_PIB (const union value *input, const struct fmt_spec *format,
 328             char *output)
 329 {
 330   double number = round (input->f * power10 (format->d));
 331   if (input->f == SYSMIS
 332       || number < 0 || number >= power256 (format->w))
 333     memset (output, 0, format->w);
 334   else
 335     output_binary_integer (number, format->w,
 336                            settings_get_output_integer_format (), output);
 337
 338   output[format->w] = '\0';
 339 }
 340
 341 /* Outputs PIBHEX format. */
 342 static void
 343 output_PIBHEX (const union value *input, const struct fmt_spec *format,
 344                char *output)
 345 {
 346   double number = round (input->f);
 347   if (input->f == SYSMIS)
 348     output_missing (format, output);
 349   else if (input->f < 0 || number >= power256 (format->w / 2))
 350     output_overflow (format, output);
 351   else
 352     {
 353       char tmp[8];
 354       output_binary_integer (number, format->w / 2, INTEGER_MSB_FIRST, tmp);
 355       output_hex (tmp, format->w / 2, output);
 356     }
 357
 358 }
 359
 360 /* Outputs RB format. */
 361 static void
 362 output_RB (const union value *input, const struct fmt_spec *format,
 363            char *output)
 364 {
 365   double d = input->f;
 366   memcpy (output, &d, format->w);
 367
 368   output[format->w] = '\0';
 369 }
 370
 371 /* Outputs RBHEX format. */
 372 static void
 373 output_RBHEX (const union value *input, const struct fmt_spec *format,
 374               char *output)
 375 {
 376   double d = input->f;
 377
 378   output_hex (&d, format->w / 2, output);
 379 }
 380
 381 /* Outputs DATE, ADATE, EDATE, JDATE, SDATE, QYR, MOYR, WKYR,
 382    DATETIME, TIME, and DTIME formats. */
 383 static void
 384 output_date (const union value *input, const struct fmt_spec *format,
 385              char *output)
 386 {
 387   double number = input->f;
 388   int year, month, day, yday;
 389
 390   const char *template = fmt_date_template (format->type);
 391   size_t template_width = strlen (template);
 392   int excess_width = format->w - template_width;
 393
 394   char tmp[64];
 395   char *p = tmp;
 396
 397   assert (format->w >= template_width);
 398   if (number == SYSMIS)
 399     goto missing;
 400
 401   if (fmt_get_category (format->type) == FMT_CAT_DATE)
 402     {
 403       if (number <= 0)
 404         goto missing;
 405       calendar_offset_to_gregorian (number / 60. / 60. / 24.,
 406                                     &year, &month, &day, &yday);
 407       number = fmod (number, 60. * 60. * 24.);
 408     }
 409   else
 410     year = month = day = yday = 0;
 411
 412   while (*template != '\0')
 413     {
 414       int ch = *template;
 415       int count = 1;
 416       while (template[count] == ch)
 417         count++;
 418       template += count;
 419
 420       switch (ch)
 421         {
 422         case 'd':
 423           if (count < 3)
 424             p += sprintf (p, "%02d", day);
 425           else
 426             p += sprintf (p, "%03d", yday);
 427           break;
 428         case 'm':
 429           if (count < 3)
 430             p += sprintf (p, "%02d", month);
 431           else
 432             {
 433               static const char *const months[12] =
 434                 {
 435                   "JAN", "FEB", "MAR", "APR", "MAY", "JUN",
 436                   "JUL", "AUG", "SEP", "OCT", "NOV", "DEC",
 437                 };
 438               p = stpcpy (p, months[month - 1]);
 439             }
 440           break;
 441         case 'y':
 442           if (count >= 4 || excess_width >= 2)
 443             {
 444               if (year <= 9999)
 445                 p += sprintf (p, "%04d", year);
 446               else if (format->type == FMT_DATETIME)
 447                 p = stpcpy (p, "****");
 448               else
 449                 goto overflow;
 450             }
 451           else
 452             {
 453               int epoch =  settings_get_epoch ();
 454               int offset = year - epoch;
 455               if (offset < 0 || offset > 99)
 456                 goto overflow;
 457               p += sprintf (p, "%02d", abs (year) % 100);
 458             }
 459           break;
 460         case 'q':
 461           p += sprintf (p, "%d", (month - 1) / 3 + 1);
 462           break;
 463         case 'w':
 464           p += sprintf (p, "%2d", (yday - 1) / 7 + 1);
 465           break;
 466         case 'D':
 467           if (number < 0)
 468             *p++ = '-';
 469           number = fabs (number);
 470           p += sprintf (p, "%*.0f", count, floor (number / 60. / 60. / 24.));
 471           number = fmod (number, 60. * 60. * 24.);
 472           break;
 473         case 'H':
 474           if (number < 0)
 475             *p++ = '-';
 476           number = fabs (number);
 477           p += sprintf (p, "%0*.0f", count, floor (number / 60. / 60.));
 478           number = fmod (number, 60. * 60.);
 479           break;
 480         case 'M':
 481           p += sprintf (p, "%02d", (int) floor (number / 60.));
 482           number = fmod (number, 60.);
 483           excess_width = format->w - (p - tmp);
 484           if (excess_width < 0)
 485             goto overflow;
 486           if (excess_width == 3 || excess_width == 4
 487               || (excess_width >= 5 && format->d == 0))
 488             p += sprintf (p, ":%02d", (int) number);
 489           else if (excess_width >= 5)
 490             {
 491               int d = MIN (format->d, excess_width - 4);
 492               int w = d + 3;
 493               sprintf (p, ":%0*.*f", w, d, number);
 494               if (settings_get_decimal_char (FMT_F) != '.')
 495                 {
 496                   char *cp = strchr (p, '.');
 497                   if (cp != NULL)
 498                     *cp = settings_get_decimal_char (FMT_F);
 499                 }
 500               p += strlen (p);
 501             }
 502           break;
 503         case 'X':
 504           *p++ = ' ';
 505           break;
 506         default:
 507           assert (count == 1);
 508           *p++ = ch;
 509           break;
 510         }
 511     }
 512
 513   buf_copy_lpad (output, format->w, tmp, p - tmp, ' ');
 514   output[format->w] = '\0';
 515   return;
 516
 517  overflow:
 518   output_overflow (format, output);
 519   return;
 520
 521  missing:
 522   output_missing (format, output);
 523   return;
 524 }
 525
 526 /* Outputs WKDAY format. */
 527 static void
 528 output_WKDAY (const union value *input, const struct fmt_spec *format,
 529               char *output)
 530 {
 531   static const char *const weekdays[7] =
 532     {
 533       "SUNDAY", "MONDAY", "TUESDAY", "WEDNESDAY",
 534       "THURSDAY", "FRIDAY", "SATURDAY",
 535     };
 536
 537   if (input->f >= 1 && input->f < 8)
 538     {
 539       buf_copy_str_rpad (output, format->w,
 540                          weekdays[(int) input->f - 1], ' ');
 541       output[format->w] = '\0';
 542     }
 543   else
 544     {
 545       if (input->f != SYSMIS)
 546         msg (ME, _("Weekday number %f is not between 1 and 7."), input->f);
 547       output_missing (format, output);
 548     }
 549
 550 }
 551
 552 /* Outputs MONTH format. */
 553 static void
 554 output_MONTH (const union value *input, const struct fmt_spec *format,
 555               char *output)
 556 {
 557   static const char *const months[12] =
 558     {
 559       "JANUARY", "FEBRUARY", "MARCH", "APRIL", "MAY", "JUNE",
 560       "JULY", "AUGUST", "SEPTEMBER", "OCTOBER", "NOVEMBER", "DECEMBER",
 561     };
 562
 563   if (input->f >= 1 && input->f < 13)
 564     {
 565       buf_copy_str_rpad (output, format->w, months[(int) input->f - 1], ' ');
 566       output[format->w] = '\0';
 567     }
 568   else
 569     {
 570       if (input->f != SYSMIS)
 571         msg (ME, _("Month number %f is not between 1 and 12."), input->f);
 572       output_missing (format, output);
 573     }
 574
 575 }
 576
 577 /* Outputs A format. */
 578 static void
 579 output_A (const union value *input UNUSED,
 580           const struct fmt_spec *format UNUSED, char *output UNUSED)
 581 {
 582   NOT_REACHED ();
 583 }
 584
 585 /* Outputs AHEX format. */
 586 static void
 587 output_AHEX (const union value *input, const struct fmt_spec *format,
 588              char *output)
 589 {
 590   output_hex (value_str (input, format->w), format->w / 2, output);
 591 }
 592 \f
 593 /* Decimal and scientific formatting. */
 594
 595 /* If REQUEST plus the current *WIDTH fits within MAX_WIDTH,
 596    increments *WIDTH by REQUEST and return true.
 597    Otherwise returns false without changing *WIDTH. */
 598 static bool
 599 allocate_space (int request, int max_width, int *width)
 600 {
 601   assert (*width <= max_width);
 602   if (request + *width <= max_width)
 603     {
 604       *width += request;
 605       return true;
 606     }
 607   else
 608     return false;
 609 }
 610
 611 /* Tries to compose the number represented by R, in the style of
 612    FORMAT, into OUTPUT.  Returns true if successful, false on
 613    failure, which occurs if FORMAT's width is too narrow.  If
 614    REQUIRE_AFFIXES is true, then the prefix and suffix specified
 615    by FORMAT's style must be included; otherwise, they may be
 616    omitted to make the number fit. */
 617 static bool
 618 output_decimal (const struct rounder *r, const struct fmt_spec *format,
 619                 bool require_affixes, char *output)
 620 {
 621   const struct fmt_number_style *style =
 622     settings_get_style (format->type);
 623
 624   int decimals;
 625
 626   for (decimals = format->d; decimals >= 0; decimals--)
 627     {
 628       /* Formatted version of magnitude of NUMBER. */
 629       char magnitude[64];
 630
 631       /* Number of digits in MAGNITUDE's integer and fractional parts. */
 632       int integer_digits;
 633
 634       /* Amount of space within the field width already claimed.
 635          Initially this is the width of MAGNITUDE, then it is reduced
 636          in stages as space is allocated to prefixes and suffixes and
 637          grouping characters. */
 638       int width;
 639
 640       /* Include various decorations? */
 641       bool add_neg_prefix;
 642       bool add_affixes;
 643       bool add_grouping;
 644
 645       /* Position in output. */
 646       char *p;
 647
 648       /* Make sure there's room for the number's magnitude, plus
 649          the negative suffix, plus (if negative) the negative
 650          prefix. */
 651       width = rounder_width (r, decimals, &integer_digits, &add_neg_prefix);
 652       width += style->neg_suffix.width;
 653       if (add_neg_prefix)
 654         width += style->neg_prefix.width;
 655       if (width > format->w)
 656         continue;
 657
 658       /* If there's room for the prefix and suffix, allocate
 659          space.  If the affixes are required, but there's no
 660          space, give up. */
 661       add_affixes = allocate_space (fmt_affix_width (style),
 662                                     format->w, &width);
 663       if (!add_affixes && require_affixes)
 664         continue;
 665
 666       /* Check whether we should include grouping characters.
 667          We need room for a complete set or we don't insert any at all.
 668          We don't include grouping characters if decimal places were
 669          requested but they were all dropped. */
 670       add_grouping = (style->grouping != 0
 671                       && integer_digits > 3
 672                       && (format->d == 0 || decimals > 0)
 673                       && allocate_space ((integer_digits - 1) / 3,
 674                                          format->w, &width));
 675
 676       /* Format the number's magnitude. */
 677       rounder_format (r, decimals, magnitude);
 678
 679       /* Assemble number. */
 680       p = output;
 681       if (format->w > width)
 682         p = mempset (p, ' ', format->w - width);
 683       if (add_neg_prefix)
 684         p = stpcpy (p, style->neg_prefix.s);
 685       if (add_affixes)
 686         p = stpcpy (p, style->prefix.s);
 687       if (!add_grouping)
 688         p = mempcpy (p, magnitude, integer_digits);
 689       else
 690         {
 691           int i;
 692           for (i = 0; i < integer_digits; i++)
 693             {
 694               if (i > 0 && (integer_digits - i) % 3 == 0)
 695                 *p++ = style->grouping;
 696               *p++ = magnitude[i];
 697             }
 698         }
 699       if (decimals > 0)
 700         {
 701           *p++ = style->decimal;
 702           p = mempcpy (p, &magnitude[integer_digits + 1], decimals);
 703         }
 704       if (add_affixes)
 705         p = stpcpy (p, style->suffix.s);
 706       if (add_neg_prefix)
 707         p = stpcpy (p, style->neg_suffix.s);
 708       else
 709         p = mempset (p, ' ', style->neg_suffix.width);
 710
 711       assert (p >= output + format->w);
 712       assert (p <= output + format->w + style->extra_bytes);
 713       *p = '\0';
 714
 715       return true;
 716     }
 717   return false;
 718 }
 719
 720 /* Formats NUMBER into OUTPUT in scientific notation according to
 721    the style of the format specified in FORMAT. */
 722 static bool
 723 output_scientific (double number, const struct fmt_spec *format,
 724                    bool require_affixes, char *output)
 725 {
 726   const struct fmt_number_style *style =
 727     settings_get_style (format->type);
 728   int width;
 729   int fraction_width;
 730   bool add_affixes;
 731   char *p;
 732
 733   /* Allocate minimum required space. */
 734   width = 6 + style->neg_suffix.width;
 735   if (number < 0)
 736     width += style->neg_prefix.width;
 737   if (width > format->w)
 738     return false;
 739
 740   /* Check for room for prefix and suffix. */
 741   add_affixes = allocate_space (fmt_affix_width (style), format->w, &width);
 742   if (require_affixes && !add_affixes)
 743     return false;
 744
 745   /* Figure out number of characters we can use for the fraction,
 746      if any.  (If that turns out to be 1, then we'll output a
 747      decimal point without any digits following; that's what the
 748      # flag does in the call to sprintf, below.) */
 749   fraction_width = MIN (MIN (format->d + 1, format->w - width), 16);
 750   if (format->type != FMT_E && fraction_width == 1)
 751     fraction_width = 0;
 752   width += fraction_width;
 753
 754   /* Format (except suffix). */
 755   p = output;
 756   if (width < format->w)
 757     p = mempset (p, ' ', format->w - width);
 758   if (number < 0)
 759     p = stpcpy (p, style->neg_prefix.s);
 760   if (add_affixes)
 761     p = stpcpy (p, style->prefix.s);
 762   if (fraction_width > 0)
 763     sprintf (p, "%#.*E", fraction_width - 1, fabs (number));
 764   else
 765     sprintf (p, "%.0E", fabs (number));
 766
 767   /* The C locale always uses a period `.' as a decimal point.
 768      Translate to comma if necessary. */
 769   if (style->decimal != '.')
 770     {
 771       char *cp = strchr (p, '.');
 772       if (cp != NULL)
 773         *cp = style->decimal;
 774     }
 775
 776   /* Make exponent have exactly three digits, plus sign. */
 777   {
 778     char *cp = strchr (p, 'E') + 1;
 779     long int exponent = strtol (cp, NULL, 10);
 780     if (abs (exponent) > 999)
 781       return false;
 782     sprintf (cp, "%+04ld", exponent);
 783   }
 784
 785   /* Add suffixes. */
 786   p = strchr (p, '\0');
 787   if (add_affixes)
 788     p = stpcpy (p, style->suffix.s);
 789   if (number < 0)
 790     p = stpcpy (p, style->neg_suffix.s);
 791   else
 792     p = mempset (p, ' ', style->neg_suffix.width);
 793
 794   assert (p >= output + format->w);
 795   assert (p <= output + format->w + style->extra_bytes);
 796   *p = '\0';
 797
 798   return true;
 799 }
 800 \f
 801 /* Returns true if the magnitude represented by R should be
 802    rounded up when chopped off at DECIMALS decimal places, false
 803    if it should be rounded down. */
 804 static bool
 805 should_round_up (const struct rounder *r, int decimals)
 806 {
 807   int digit = r->string[r->integer_digits + decimals + 1];
 808   assert (digit >= '0' && digit <= '9');
 809   return digit >= '5';
 810 }
 811
 812 /* Initializes R for formatting the magnitude of NUMBER to no
 813    more than MAX_DECIMAL decimal places. */
 814 static void
 815 rounder_init (struct rounder *r, double number, int max_decimals)
 816 {
 817   assert (fabs (number) < 1e41);
 818   assert (max_decimals >= 0 && max_decimals <= 16);
 819   if (max_decimals == 0)
 820     {
 821       /* Fast path.  No rounding needed.
 822
 823          We append ".00" to the integer representation because
 824          round_up assumes that fractional digits are present.  */
 825       sprintf (r->string, "%.0f.00", fabs (round (number)));
 826     }
 827   else
 828     {
 829       /* Slow path.
 830
 831          This is more difficult than it really should be because
 832          we have to make sure that numbers that are exactly
 833          halfway between two representations are always rounded
 834          away from zero.  This is not what sprintf normally does
 835          (usually it rounds to even), so we have to fake it as
 836          best we can, by formatting with extra precision and then
 837          doing the rounding ourselves.
 838
 839          We take up to two rounds to format numbers.  In the
 840          first round, we obtain 2 digits of precision beyond
 841          those requested by the user.  If those digits are
 842          exactly "50", then in a second round we format with as
 843          many digits as are significant in a "double".
 844
 845          It might be better to directly implement our own
 846          floating-point formatting routine instead of relying on
 847          the system's sprintf implementation.  But the classic
 848          Steele and White paper on printing floating-point
 849          numbers does not hint how to do what we want, and it's
 850          not obvious how to change their algorithms to do so.  It
 851          would also be a lot of work. */
 852       sprintf (r->string, "%.*f", max_decimals + 2, fabs (number));
 853       if (!strcmp (r->string + strlen (r->string) - 2, "50"))
 854         {
 855           int binary_exponent, decimal_exponent, format_decimals;
 856           frexp (number, &binary_exponent);
 857           decimal_exponent = binary_exponent * 3 / 10;
 858           format_decimals = (DBL_DIG + 1) - decimal_exponent;
 859           if (format_decimals > max_decimals + 2)
 860             sprintf (r->string, "%.*f", format_decimals, fabs (number));
 861         }
 862     }
 863
 864   if (r->string[0] == '0')
 865     memmove (r->string, &r->string[1], strlen (r->string));
 866
 867   r->leading_zeros = strspn (r->string, "0.");
 868   r->leading_nines = strspn (r->string, "9.");
 869   r->integer_digits = strchr (r->string, '.') - r->string;
 870   r->negative = number < 0;
 871 }
 872
 873 /* Returns the number of characters required to format the
 874    magnitude represented by R to DECIMALS decimal places.
 875    The return value includes integer digits and a decimal point
 876    and fractional digits, if any, but it does not include any
 877    negative prefix or suffix or other affixes.
 878
 879    *INTEGER_DIGITS is set to the number of digits before the
 880    decimal point in the output, between 0 and 40.
 881
 882    If R represents a negative number and its rounded
 883    representation would include at least one nonzero digit,
 884    *NEGATIVE is set to true; otherwise, it is set to false. */
 885 static int
 886 rounder_width (const struct rounder *r, int decimals,
 887                int *integer_digits, bool *negative)
 888 {
 889   /* Calculate base measures. */
 890   int width = r->integer_digits;
 891   if (decimals > 0)
 892     width += decimals + 1;
 893   *integer_digits = r->integer_digits;
 894   *negative = r->negative;
 895
 896   /* Rounding can cause adjustments. */
 897   if (should_round_up (r, decimals))
 898     {
 899       /* Rounding up leading 9s adds a new digit (a 1). */
 900       if (r->leading_nines >= width)
 901         {
 902           width++;
 903           ++*integer_digits;
 904         }
 905     }
 906   else
 907     {
 908       /* Rounding down. */
 909       if (r->leading_zeros >= width)
 910         {
 911           /* All digits that remain after rounding are zeros.
 912              Therefore we drop the negative sign. */
 913           *negative = false;
 914           if (r->integer_digits == 0 && decimals == 0)
 915             {
 916               /* No digits at all are left.  We need to display
 917                  at least a single digit (a zero). */
 918               assert (width == 0);
 919               width++;
 920               *integer_digits = 1;
 921             }
 922         }
 923     }
 924   return width;
 925 }
 926
 927 /* Formats the magnitude represented by R into OUTPUT, rounding
 928    to DECIMALS decimal places.  Exactly as many characters as
 929    indicated by rounder_width are written.  No terminating null
 930    is appended. */
 931 static void
 932 rounder_format (const struct rounder *r, int decimals, char *output)
 933 {
 934   int base_width = r->integer_digits + (decimals > 0 ? decimals + 1 : 0);
 935   if (should_round_up (r, decimals))
 936     {
 937       if (r->leading_nines < base_width)
 938         {
 939           /* Rounding up.  This is the common case where rounding
 940              up doesn't add an extra digit. */
 941           char *p;
 942           memcpy (output, r->string, base_width);
 943           for (p = output + base_width - 1; ; p--)
 944             {
 945               assert (p >= output);
 946               if (*p == '9')
 947                 *p = '0';
 948               else if (*p >= '0' && *p <= '8')
 949                 {
 950                   (*p)++;
 951                   break;
 952                 }
 953               else
 954                 assert (*p == '.');
 955             }
 956         }
 957       else
 958         {
 959           /* Rounding up leading 9s causes the result to be a 1
 960              followed by a number of 0s, plus a decimal point. */
 961           char *p = output;
 962           *p++ = '1';
 963           p = mempset (p, '0', r->integer_digits);
 964           if (decimals > 0)
 965             {
 966               *p++ = '.';
 967               p = mempset (p, '0', decimals);
 968             }
 969           assert (p == output + base_width + 1);
 970         }
 971     }
 972   else
 973     {
 974       /* Rounding down. */
 975       if (r->integer_digits != 0 || decimals != 0)
 976         {
 977           /* Common case: just copy the digits. */
 978           memcpy (output, r->string, base_width);
 979         }
 980       else
 981         {
 982           /* No digits remain.  The output is just a zero. */
 983           output[0] = '0';
 984         }
 985     }
 986 }
 987 \f
 988 /* Helper functions. */
 989
 990 /* Returns 10**X. */
 991 static double PURE_FUNCTION
 992 power10 (int x)
 993 {
 994   static const double p[] =
 995     {
 996       1e0, 1e1, 1e2, 1e3, 1e4, 1e5, 1e6, 1e7, 1e8, 1e9,
 997       1e10, 1e11, 1e12, 1e13, 1e14, 1e15, 1e16, 1e17, 1e18, 1e19,
 998       1e20, 1e21, 1e22, 1e23, 1e24, 1e25, 1e26, 1e27, 1e28, 1e29,
 999       1e30, 1e31, 1e32, 1e33, 1e34, 1e35, 1e36, 1e37, 1e38, 1e39,
1000       1e40,
1001     };
1002   return x >= 0 && x < sizeof p / sizeof *p ? p[x] : pow (10.0, x);
1003 }
1004
1005 /* Returns 256**X. */
1006 static double PURE_FUNCTION
1007 power256 (int x)
1008 {
1009   static const double p[] =
1010     {
1011       1.0,
1012       256.0,
1013       65536.0,
1014       16777216.0,
1015       4294967296.0,
1016       1099511627776.0,
1017       281474976710656.0,
1018       72057594037927936.0,
1019       18446744073709551616.0
1020     };
1021   return x >= 0 && x < sizeof p / sizeof *p ? p[x] : pow (256.0, x);
1022 }
1023
1024 /* Formats non-finite NUMBER into OUTPUT according to the width
1025    given in FORMAT. */
1026 static void
1027 output_infinite (double number, const struct fmt_spec *format, char *output)
1028 {
1029   assert (!isfinite (number));
1030
1031   if (format->w >= 3)
1032     {
1033       const char *s;
1034
1035       if (isnan (number))
1036         s = "NaN";
1037       else if (isinf (number))
1038         s = number > 0 ? "+Infinity" : "-Infinity";
1039       else
1040         s = "Unknown";
1041
1042       buf_copy_str_lpad (output, format->w, s, ' ');
1043     }
1044   else
1045     output_overflow (format, output);
1046
1047   output[format->w] = '\0';
1048 }
1049
1050 /* Formats OUTPUT as a missing value for the given FORMAT. */
1051 static void
1052 output_missing (const struct fmt_spec *format, char *output)
1053 {
1054   memset (output, ' ', format->w);
1055
1056   if (format->type != FMT_N)
1057     {
1058       int dot_ofs = (format->type == FMT_PCT ? 2
1059                      : format->type == FMT_E ? 5
1060                      : 1);
1061       output[MAX (0, format->w - format->d - dot_ofs)] = '.';
1062     }
1063   else
1064     output[format->w - 1] = '.';
1065
1066   output[format->w] = '\0';
1067 }
1068
1069 /* Formats OUTPUT for overflow given FORMAT. */
1070 static void
1071 output_overflow (const struct fmt_spec *format, char *output)
1072 {
1073   memset (output, '*', format->w);
1074   output[format->w] = '\0';
1075 }
1076
1077 /* Converts the integer part of NUMBER to a packed BCD number
1078    with the given number of DIGITS in OUTPUT.  If DIGITS is odd,
1079    the least significant nibble of the final byte in OUTPUT is
1080    set to 0.  Returns true if successful, false if NUMBER is not
1081    representable.  On failure, OUTPUT is cleared to all zero
1082    bytes. */
1083 static bool
1084 output_bcd_integer (double number, int digits, char *output)
1085 {
1086   char decimal[64];
1087
1088   assert (digits < sizeof decimal);
1089
1090   output[DIV_RND_UP (digits, 2)] = '\0';
1091   if (number != SYSMIS
1092       && number >= 0.
1093       && number < power10 (digits)
1094       && sprintf (decimal, "%0*.0f", digits, round (number)) == digits)
1095     {
1096       const char *src = decimal;
1097       int i;
1098
1099       for (i = 0; i < digits / 2; i++)
1100         {
1101           int d0 = *src++ - '0';
1102           int d1 = *src++ - '0';
1103           *output++ = (d0 << 4) + d1;
1104         }
1105       if (digits % 2)
1106         *output = (*src - '0') << 4;
1107
1108       return true;
1109     }
1110   else
1111     {
1112       memset (output, 0, DIV_RND_UP (digits, 2));
1113       return false;
1114     }
1115 }
1116
1117 /* Writes VALUE to OUTPUT as a BYTES-byte binary integer of the
1118    given INTEGER_FORMAT. */
1119 static void
1120 output_binary_integer (uint64_t value, int bytes,
1121                        enum integer_format integer_format, char *output)
1122 {
1123   integer_put (value, integer_format, output, bytes);
1124 }
1125
1126 /* Converts the BYTES bytes in DATA to twice as many hexadecimal
1127    digits in OUTPUT. */
1128 static void
1129 output_hex (const void *data_, size_t bytes, char *output)
1130 {
1131   const uint8_t *data = data_;
1132   size_t i;
1133
1134   for (i = 0; i < bytes; i++)
1135     {
1136       static const char hex_digits[] = "0123456789ABCDEF";
1137       *output++ = hex_digits[data[i] >> 4];
1138       *output++ = hex_digits[data[i] & 15];
1139     }
1140   *output = '\0';
1141 }