pintos-os.org Git - pspp/blob - src/data/sys-file-writer.c

   1 /* PSPP - a program for statistical analysis.
   2    Copyright (C) 1997-9, 2000, 2006, 2007 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 "sys-file-writer.h"
  20 #include "sys-file-private.h"
  21
  22 #include <ctype.h>
  23 #include <errno.h>
  24 #include <fcntl.h>
  25 #include <stdint.h>
  26 #include <stdlib.h>
  27 #include <sys/stat.h>
  28 #include <time.h>
  29 #include <unistd.h>
  30
  31 #include <libpspp/alloc.h>
  32 #include <libpspp/float-format.h>
  33 #include <libpspp/integer-format.h>
  34 #include <libpspp/magic.h>
  35 #include <libpspp/message.h>
  36 #include <libpspp/misc.h>
  37 #include <libpspp/str.h>
  38 #include <libpspp/version.h>
  39
  40 #include <data/case.h>
  41 #include <data/casewriter-provider.h>
  42 #include <data/casewriter.h>
  43 #include <data/dictionary.h>
  44 #include <data/file-handle-def.h>
  45 #include <data/format.h>
  46 #include <data/missing-values.h>
  47 #include <data/settings.h>
  48 #include <data/short-names.h>
  49 #include <data/value-labels.h>
  50 #include <data/variable.h>
  51
  52 #include "minmax.h"
  53 #include "unlocked-io.h"
  54
  55 #include "gettext.h"
  56 #define _(msgid) gettext (msgid)
  57
  58 /* Compression bias used by PSPP.  Values between (1 -
  59    COMPRESSION_BIAS) and (251 - COMPRESSION_BIAS) inclusive can be
  60    compressed. */
  61 #define COMPRESSION_BIAS 100
  62
  63 /* System file writer. */
  64 struct sfm_writer
  65   {
  66     struct file_handle *fh;     /* File handle. */
  67     FILE *file;                 /* File stream. */
  68
  69     bool compress;              /* 1=compressed, 0=not compressed. */
  70     casenumber case_cnt;        /* Number of cases written so far. */
  71
  72     /* Compression buffering.
  73
  74        Compressed data is output as groups of 8 1-byte opcodes
  75        followed by up to 8 (depending on the opcodes) 8-byte data
  76        items.  Data items and opcodes arrive at the same time but
  77        must be reordered for writing to disk, thus a small amount
  78        of buffering here. */
  79     uint8_t opcodes[8];         /* Buffered opcodes. */
  80     int opcode_cnt;             /* Number of buffered opcodes. */
  81     uint8_t data[8][8];         /* Buffered data. */
  82     int data_cnt;               /* Number of buffered data items. */
  83
  84     /* Variables. */
  85     struct sfm_var *sfm_vars;   /* Variables. */
  86     size_t sfm_var_cnt;         /* Number of variables. */
  87     size_t segment_cnt;         /* Number of variables including extra segments
  88                                    for long string variables. */
  89   };
  90
  91 static struct casewriter_class sys_file_casewriter_class;
  92
  93 static void write_header (struct sfm_writer *, const struct dictionary *);
  94 static void write_variable (struct sfm_writer *, const struct variable *);
  95 static void write_value_labels (struct sfm_writer *,
  96                                 struct variable *, int idx);
  97 static void write_integer_info_record (struct sfm_writer *);
  98 static void write_float_info_record (struct sfm_writer *);
  99
 100 static void write_longvar_table (struct sfm_writer *w,
 101                                  const struct dictionary *dict);
 102
 103 static void write_vls_length_table (struct sfm_writer *w,
 104                               const struct dictionary *dict);
 105
 106
 107 static void write_variable_display_parameters (struct sfm_writer *w,
 108                                                const struct dictionary *dict);
 109
 110 static void write_documents (struct sfm_writer *, const struct dictionary *);
 111
 112 static void write_int (struct sfm_writer *, int32_t);
 113 static inline void convert_double_to_output_format (double, uint8_t[8]);
 114 static void write_float (struct sfm_writer *, double);
 115 static void write_string (struct sfm_writer *, const char *, size_t);
 116 static void write_bytes (struct sfm_writer *, const void *, size_t);
 117 static void write_zeros (struct sfm_writer *, size_t);
 118 static void write_spaces (struct sfm_writer *, size_t);
 119 static void write_value (struct sfm_writer *, const union value *, int width);
 120
 121 static void write_case_uncompressed (struct sfm_writer *, struct ccase *);
 122 static void write_case_compressed (struct sfm_writer *, struct ccase *);
 123 static void flush_compressed (struct sfm_writer *);
 124 static void put_cmp_opcode (struct sfm_writer *, uint8_t);
 125 static void put_cmp_number (struct sfm_writer *, double);
 126 static void put_cmp_string (struct sfm_writer *, const void *, size_t);
 127
 128 bool write_error (const struct sfm_writer *);
 129 bool close_writer (struct sfm_writer *);
 130
 131 /* Returns default options for writing a system file. */
 132 struct sfm_write_options
 133 sfm_writer_default_options (void)
 134 {
 135   struct sfm_write_options opts;
 136   opts.create_writeable = true;
 137   opts.compress = get_scompression ();
 138   opts.version = 3;
 139   return opts;
 140 }
 141
 142 /* Opens the system file designated by file handle FH for writing
 143    cases from dictionary D according to the given OPTS.  If
 144    COMPRESS is nonzero, the system file will be compressed.
 145
 146    No reference to D is retained, so it may be modified or
 147    destroyed at will after this function returns.  D is not
 148    modified by this function, except to assign short names. */
 149 struct casewriter *
 150 sfm_open_writer (struct file_handle *fh, struct dictionary *d,
 151                  struct sfm_write_options opts)
 152 {
 153   struct sfm_writer *w = NULL;
 154   mode_t mode;
 155   int fd;
 156   int idx;
 157   int i;
 158
 159   /* Check version. */
 160   if (opts.version != 2 && opts.version != 3)
 161     {
 162       msg (ME, _("Unknown system file version %d. Treating as version %d."),
 163            opts.version, 3);
 164       opts.version = 3;
 165     }
 166
 167   /* Create file. */
 168   mode = S_IRUSR | S_IRGRP | S_IROTH;
 169   if (opts.create_writeable)
 170     mode |= S_IWUSR | S_IWGRP | S_IWOTH;
 171   fd = open (fh_get_file_name (fh), O_WRONLY | O_CREAT | O_TRUNC, mode);
 172   if (fd < 0)
 173     goto open_error;
 174
 175   /* Open file handle. */
 176   if (!fh_open (fh, FH_REF_FILE, "system file", "we"))
 177     goto error;
 178
 179   /* Create and initialize writer. */
 180   w = xmalloc (sizeof *w);
 181   w->fh = fh;
 182   w->file = fdopen (fd, "w");
 183
 184   w->compress = opts.compress;
 185   w->case_cnt = 0;
 186
 187   w->opcode_cnt = w->data_cnt = 0;
 188
 189   /* Figure out how to map in-memory case data to on-disk case
 190      data.  Also count the number of segments.  Very long strings
 191      occupy multiple segments, otherwise each variable only takes
 192      one segment. */
 193   w->segment_cnt = sfm_dictionary_to_sfm_vars (d, &w->sfm_vars,
 194                                                &w->sfm_var_cnt);
 195
 196   /* Check that file create succeeded. */
 197   if (w->file == NULL)
 198     {
 199       close (fd);
 200       goto open_error;
 201     }
 202
 203   /* Write the file header. */
 204   write_header (w, d);
 205
 206   /* Write basic variable info. */
 207   short_names_assign (d);
 208   for (i = 0; i < dict_get_var_cnt (d); i++)
 209     write_variable (w, dict_get_var (d, i));
 210
 211   /* Write out value labels. */
 212   idx = 0;
 213   for (i = 0; i < dict_get_var_cnt (d); i++)
 214     {
 215       struct variable *v = dict_get_var (d, i);
 216
 217       write_value_labels (w, v, idx);
 218       idx += sfm_width_to_octs (var_get_width (v));
 219     }
 220
 221   if (dict_get_documents (d) != NULL)
 222     write_documents (w, d);
 223
 224   write_integer_info_record (w);
 225   write_float_info_record (w);
 226
 227   write_variable_display_parameters (w, d);
 228
 229   if (opts.version >= 3)
 230     write_longvar_table (w, d);
 231
 232   write_vls_length_table (w, d);
 233
 234   /* Write end-of-headers record. */
 235   write_int (w, 999);
 236   write_int (w, 0);
 237
 238   if (write_error (w))
 239     goto error;
 240
 241   return casewriter_create (dict_get_next_value_idx (d),
 242                             &sys_file_casewriter_class, w);
 243
 244  error:
 245   close_writer (w);
 246   return NULL;
 247
 248  open_error:
 249   msg (ME, _("Error opening \"%s\" for writing as a system file: %s."),
 250        fh_get_file_name (fh), strerror (errno));
 251   goto error;
 252 }
 253
 254 /* Returns value of X truncated to two least-significant digits. */
 255 static int
 256 rerange (int x)
 257 {
 258   if (x < 0)
 259     x = -x;
 260   if (x >= 100)
 261     x %= 100;
 262   return x;
 263 }
 264
 265 /* Calculates the offset of data for TARGET_VAR from the
 266    beginning of each case's data for dictionary D.  The return
 267    value is in "octs" (8-byte units). */
 268 static int
 269 calc_oct_idx (const struct dictionary *d, struct variable *target_var)
 270 {
 271   int oct_idx;
 272   int i;
 273
 274   oct_idx = 0;
 275   for (i = 0; i < dict_get_var_cnt (d); i++)
 276     {
 277       struct variable *var = dict_get_var (d, i);
 278       if (var == target_var)
 279         break;
 280       oct_idx += sfm_width_to_octs (var_get_width (var));
 281     }
 282   return oct_idx;
 283 }
 284
 285 /* Write the sysfile_header header to system file W. */
 286 static void
 287 write_header (struct sfm_writer *w, const struct dictionary *d)
 288 {
 289   char prod_name[61];
 290   char creation_date[10];
 291   char creation_time[9];
 292   const char *file_label;
 293   struct variable *weight;
 294
 295   time_t t;
 296
 297   /* Record-type code. */
 298   write_string (w, "$FL2", 4);
 299
 300   /* Product identification. */
 301   snprintf (prod_name, sizeof prod_name, "@(#) SPSS DATA FILE %s - %s",
 302             version, host_system);
 303   write_string (w, prod_name, 60);
 304
 305   /* Layout code. */
 306   write_int (w, 2);
 307
 308   /* Number of `union value's per case. */
 309   write_int (w, calc_oct_idx (d, NULL));
 310
 311   /* Compressed? */
 312   write_int (w, w->compress);
 313
 314   /* Weight variable. */
 315   weight = dict_get_weight (d);
 316   write_int (w, weight != NULL ? calc_oct_idx (d, weight) + 1 : 0);
 317
 318   /* Number of cases.  We don't know this in advance, so we write
 319      -1 to indicate an unknown number of cases.  Later we can
 320      come back and overwrite it with the true value. */
 321   write_int (w, -1);
 322
 323   /* Compression bias. */
 324   write_float (w, COMPRESSION_BIAS);
 325
 326   /* Creation date and time. */
 327   if (time (&t) == (time_t) -1)
 328     {
 329       strcpy (creation_date, "01 Jan 70");
 330       strcpy (creation_time, "00:00:00");
 331     }
 332   else
 333     {
 334       static const char *month_name[12] =
 335         {
 336           "Jan", "Feb", "Mar", "Apr", "May", "Jun",
 337           "Jul", "Aug", "Sep", "Oct", "Nov", "Dec",
 338         };
 339       struct tm *tmp = localtime (&t);
 340       int day = rerange (tmp->tm_mday);
 341       int mon = rerange (tmp->tm_mon + 1);
 342       int year = rerange (tmp->tm_year);
 343       int hour = rerange (tmp->tm_hour + 1);
 344       int min = rerange (tmp->tm_min + 1);
 345       int sec = rerange (tmp->tm_sec + 1);
 346
 347       snprintf (creation_date, sizeof creation_date,
 348                 "%02d %s %02d", day, month_name[mon - 1], year);
 349       snprintf (creation_time, sizeof creation_time,
 350                 "%02d:%02d:%02d", hour - 1, min - 1, sec - 1);
 351     }
 352   write_string (w, creation_date, 9);
 353   write_string (w, creation_time, 8);
 354
 355   /* File label. */
 356   file_label = dict_get_label (d);
 357   if (file_label == NULL)
 358     file_label = "";
 359   write_string (w, file_label, 64);
 360
 361   /* Padding. */
 362   write_zeros (w, 3);
 363 }
 364
 365 /* Write format spec FMT to W, after adjusting it to be
 366    compatible with the given WIDTH. */
 367 static void
 368 write_format (struct sfm_writer *w, struct fmt_spec fmt, int width)
 369 {
 370   assert (fmt_check_output (&fmt));
 371   assert (sfm_width_to_segments (width) == 1);
 372
 373   if (width > 0)
 374     fmt_resize (&fmt, width);
 375   write_int (w, (fmt_to_io (fmt.type) << 16) | (fmt.w << 8) | fmt.d);
 376 }
 377
 378 /* Write a string continuation variable record for each 8-byte
 379    section beyond the initial 8 bytes, for a variable of the
 380    given WIDTH. */
 381 static void
 382 write_variable_continuation_records (struct sfm_writer *w, int width)
 383 {
 384   int position;
 385
 386   assert (sfm_width_to_segments (width) == 1);
 387   for (position = 8; position < width; position += 8)
 388     {
 389       write_int (w, 2);   /* Record type. */
 390       write_int (w, -1);  /* Width. */
 391       write_int (w, 0);   /* No variable label. */
 392       write_int (w, 0);   /* No missing values. */
 393       write_int (w, 0);   /* Print format. */
 394       write_int (w, 0);   /* Write format. */
 395       write_zeros (w, 8);   /* Name. */
 396     }
 397 }
 398
 399 /* Write the variable record(s) for variable V to system file
 400    W. */
 401 static void
 402 write_variable (struct sfm_writer *w, const struct variable *v)
 403 {
 404   int width = var_get_width (v);
 405   int segment_cnt = sfm_width_to_segments (width);
 406   int seg0_width = sfm_segment_alloc_width (width, 0);
 407   const struct missing_values *mv = var_get_missing_values (v);
 408   int i;
 409
 410   /* Record type. */
 411   write_int (w, 2);
 412
 413   /* Width. */
 414   write_int (w, seg0_width);
 415
 416   /* Variable has a variable label? */
 417   write_int (w, var_has_label (v));
 418
 419   /* Number of missing values.  If there is a range, then the
 420      range counts as 2 missing values and causes the number to be
 421      negated. */
 422   write_int (w, mv_has_range (mv) ? 2 - mv_n_values (mv) : mv_n_values (mv));
 423
 424   /* Print and write formats. */
 425   write_format (w, *var_get_print_format (v), seg0_width);
 426   write_format (w, *var_get_write_format (v), seg0_width);
 427
 428   /* Short name.
 429      The full name is in a translation table written
 430      separately. */
 431   write_string (w, var_get_short_name (v, 0), 8);
 432
 433   /* Value label. */
 434   if (var_has_label (v))
 435     {
 436       const char *label = var_get_label (v);
 437       size_t padded_len = ROUND_UP (MIN (strlen (label), 255), 4);
 438       write_int (w, padded_len);
 439       write_string (w, label, padded_len);
 440     }
 441
 442   /* Write the missing values, if any, range first. */
 443   if (mv_has_range (mv))
 444     {
 445       double x, y;
 446       mv_peek_range (mv, &x, &y);
 447       write_float (w, x);
 448       write_float (w, y);
 449     }
 450   for (i = 0; i < mv_n_values (mv); i++)
 451     {
 452       union value value;
 453       mv_peek_value (mv, &value, i);
 454       write_value (w, &value, seg0_width);
 455     }
 456
 457   write_variable_continuation_records (w, seg0_width);
 458
 459   /* Write additional segments for very long string variables. */
 460   for (i = 1; i < segment_cnt; i++)
 461     {
 462       int seg_width = sfm_segment_alloc_width (width, i);
 463       struct fmt_spec fmt = fmt_for_output (FMT_A, MAX (seg_width, 1), 0);
 464
 465       write_int (w, 2);           /* Variable record. */
 466       write_int (w, seg_width);   /* Width. */
 467       write_int (w, 0);           /* No variable label. */
 468       write_int (w, 0);           /* No missing values. */
 469       write_format (w, fmt, seg_width); /* Print format. */
 470       write_format (w, fmt, seg_width); /* Write format. */
 471       write_string (w, var_get_short_name (v, i), 8);
 472
 473       write_variable_continuation_records (w, seg_width);
 474     }
 475 }
 476
 477 /* Writes the value labels for variable V having system file
 478    variable index IDX to system file W. */
 479 static void
 480 write_value_labels (struct sfm_writer *w, struct variable *v, int idx)
 481 {
 482   const struct val_labs *val_labs;
 483   struct val_labs_iterator *i;
 484   struct val_lab *vl;
 485
 486   val_labs = var_get_value_labels (v);
 487   if (val_labs == NULL)
 488     return;
 489
 490   /* Value label record. */
 491   write_int (w, 3);             /* Record type. */
 492   write_int (w, val_labs_count (val_labs));
 493   for (vl = val_labs_first_sorted (val_labs, &i); vl != NULL;
 494        vl = val_labs_next (val_labs, &i))
 495     {
 496       uint8_t len = MIN (strlen (vl->label), 255);
 497
 498       write_value (w, &vl->value, var_get_width (v));
 499       write_bytes (w, &len, 1);
 500       write_bytes (w, vl->label, len);
 501       write_zeros (w, REM_RND_UP (len + 1, 8));
 502     }
 503
 504   /* Value label variable record. */
 505   write_int (w, 4);             /* Record type. */
 506   write_int (w, 1);             /* Number of variables. */
 507   write_int (w, idx + 1);       /* Variable's dictionary index. */
 508 }
 509
 510 /* Writes record type 6, document record. */
 511 static void
 512 write_documents (struct sfm_writer *w, const struct dictionary *d)
 513 {
 514   size_t line_cnt = dict_get_document_line_cnt (d);
 515
 516   write_int (w, 6);             /* Record type. */
 517   write_int (w, line_cnt);
 518   write_bytes (w, dict_get_documents (d), line_cnt * DOC_LINE_LENGTH);
 519 }
 520
 521 /* Write the alignment, width and scale values. */
 522 static void
 523 write_variable_display_parameters (struct sfm_writer *w,
 524                                    const struct dictionary *dict)
 525 {
 526   int i;
 527
 528   write_int (w, 7);             /* Record type. */
 529   write_int (w, 11);            /* Record subtype. */
 530   write_int (w, 4);             /* Data item (int32) size. */
 531   write_int (w, w->segment_cnt * 3); /* Number of data items. */
 532
 533   for (i = 0; i < dict_get_var_cnt (dict); ++i)
 534     {
 535       struct variable *v = dict_get_var (dict, i);
 536       int width = var_get_width (v);
 537       int segment_cnt = sfm_width_to_segments (width);
 538       int measure = (var_get_measure (v) == MEASURE_NOMINAL ? 1
 539                      : var_get_measure (v) == MEASURE_ORDINAL ? 2
 540                      : 3);
 541       int alignment = (var_get_alignment (v) == ALIGN_LEFT ? 0
 542                        : var_get_alignment (v) == ALIGN_RIGHT ? 1
 543                        : 2);
 544       int i;
 545
 546       for (i = 0; i < segment_cnt; i++)
 547         {
 548           int width_left = width - sfm_segment_effective_offset (width, i);
 549           write_int (w, measure);
 550           write_int (w, (i == 0 ? var_get_display_width (v)
 551                          : var_default_display_width (width_left)));
 552           write_int (w, alignment);
 553         }
 554     }
 555 }
 556
 557 /* Writes the table of lengths for very long string variables. */
 558 static void
 559 write_vls_length_table (struct sfm_writer *w,
 560                         const struct dictionary *dict)
 561 {
 562   struct string map;
 563   int i;
 564
 565   ds_init_empty (&map);
 566   for (i = 0; i < dict_get_var_cnt (dict); ++i)
 567     {
 568       const struct variable *v = dict_get_var (dict, i);
 569       if (sfm_width_to_segments (var_get_width (v)) > 1)
 570         ds_put_format (&map, "%s=%05d%c\t",
 571                        var_get_short_name (v, 0), var_get_width (v), 0);
 572     }
 573   if (!ds_is_empty (&map))
 574     {
 575       write_int (w, 7);         /* Record type. */
 576       write_int (w, 14);        /* Record subtype. */
 577       write_int (w, 1);         /* Data item (char) size. */
 578       write_int (w, ds_length (&map)); /* Number of data items. */
 579       write_bytes (w, ds_data (&map), ds_length (&map));
 580     }
 581   ds_destroy (&map);
 582 }
 583
 584 /* Writes the long variable name table. */
 585 static void
 586 write_longvar_table (struct sfm_writer *w, const struct dictionary *dict)
 587 {
 588   struct string map;
 589   size_t i;
 590
 591   ds_init_empty (&map);
 592   for (i = 0; i < dict_get_var_cnt (dict); i++)
 593     {
 594       struct variable *v = dict_get_var (dict, i);
 595
 596       if (i)
 597         ds_put_char (&map, '\t');
 598       ds_put_format (&map, "%s=%s",
 599                      var_get_short_name (v, 0), var_get_name (v));
 600     }
 601
 602   write_int (w, 7);             /* Record type. */
 603   write_int (w, 13);            /* Record subtype. */
 604   write_int (w, 1);             /* Data item (char) size. */
 605   write_int (w, ds_length (&map)); /* Number of data items. */
 606   write_bytes (w, ds_data (&map), ds_length (&map));
 607
 608   ds_destroy (&map);
 609 }
 610
 611 /* Write integer information record. */
 612 static void
 613 write_integer_info_record (struct sfm_writer *w)
 614 {
 615   int version_component[3];
 616   int float_format;
 617
 618   /* Parse the version string. */
 619   memset (version_component, 0, sizeof version_component);
 620   sscanf (bare_version, "%d.%d.%d",
 621           &version_component[0], &version_component[1], &version_component[2]);
 622
 623   /* Figure out the floating-point format. */
 624   if (FLOAT_NATIVE_64_BIT == FLOAT_IEEE_DOUBLE_LE
 625       || FLOAT_NATIVE_64_BIT == FLOAT_IEEE_DOUBLE_BE)
 626     float_format = 1;
 627   else if (FLOAT_NATIVE_64_BIT == FLOAT_Z_LONG)
 628     float_format = 2;
 629   else if (FLOAT_NATIVE_64_BIT == FLOAT_VAX_D)
 630     float_format = 3;
 631   else
 632     abort ();
 633
 634   /* Write record. */
 635   write_int (w, 7);             /* Record type. */
 636   write_int (w, 3);             /* Record subtype. */
 637   write_int (w, 4);             /* Data item (int32) size. */
 638   write_int (w, 8);             /* Number of data items. */
 639   write_int (w, version_component[0]);
 640   write_int (w, version_component[1]);
 641   write_int (w, version_component[2]);
 642   write_int (w, -1);          /* Machine code. */
 643   write_int (w, float_format);
 644   write_int (w, 1);           /* Compression code. */
 645   write_int (w, INTEGER_NATIVE == INTEGER_MSB_FIRST ? 1 : 2);
 646   write_int (w, 2);           /* 7-bit ASCII. */
 647 }
 648
 649 /* Write floating-point information record. */
 650 static void
 651 write_float_info_record (struct sfm_writer *w)
 652 {
 653   write_int (w, 7);             /* Record type. */
 654   write_int (w, 4);             /* Record subtype. */
 655   write_int (w, 8);             /* Data item (flt64) size. */
 656   write_int (w, 3);             /* Number of data items. */
 657   write_float (w, SYSMIS);      /* System-missing value. */
 658   write_float (w, HIGHEST);     /* Value used for HIGHEST in missing values. */
 659   write_float (w, LOWEST);      /* Value used for LOWEST in missing values. */
 660 }
 661 \f
 662 /* Writes case C to system file W. */
 663 static void
 664 sys_file_casewriter_write (struct casewriter *writer, void *w_,
 665                            struct ccase *c)
 666 {
 667   struct sfm_writer *w = w_;
 668
 669   if (ferror (w->file))
 670     {
 671       casewriter_force_error (writer);
 672       case_destroy (c);
 673       return;
 674     }
 675
 676   w->case_cnt++;
 677
 678   if (!w->compress)
 679     write_case_uncompressed (w, c);
 680   else
 681     write_case_compressed (w, c);
 682
 683   case_destroy (c);
 684 }
 685
 686 /* Destroys system file writer W. */
 687 static void
 688 sys_file_casewriter_destroy (struct casewriter *writer, void *w_)
 689 {
 690   struct sfm_writer *w = w_;
 691   if (!close_writer (w))
 692     casewriter_force_error (writer);
 693 }
 694
 695 /* Returns true if an I/O error has occurred on WRITER, false otherwise. */
 696 bool
 697 write_error (const struct sfm_writer *writer)
 698 {
 699   return ferror (writer->file);
 700 }
 701
 702 /* Closes a system file after we're done with it.
 703    Returns true if successful, false if an I/O error occurred. */
 704 bool
 705 close_writer (struct sfm_writer *w)
 706 {
 707   bool ok;
 708
 709   if (w == NULL)
 710     return true;
 711
 712   ok = true;
 713   if (w->file != NULL)
 714     {
 715       /* Flush buffer. */
 716       if (w->opcode_cnt > 0)
 717         flush_compressed (w);
 718       fflush (w->file);
 719
 720       ok = !write_error (w);
 721
 722       /* Seek back to the beginning and update the number of cases.
 723          This is just a courtesy to later readers, so there's no need
 724          to check return values or report errors. */
 725       if (ok && w->case_cnt <= INT32_MAX && !fseek (w->file, 80, SEEK_SET))
 726         {
 727           write_int (w, w->case_cnt);
 728           clearerr (w->file);
 729         }
 730
 731       if (fclose (w->file) == EOF)
 732         ok = false;
 733
 734       if (!ok)
 735         msg (ME, _("An I/O error occurred writing system file \"%s\"."),
 736              fh_get_file_name (w->fh));
 737     }
 738
 739   fh_close (w->fh, "system file", "we");
 740
 741   free (w->sfm_vars);
 742   free (w);
 743
 744   return ok;
 745 }
 746
 747 /* System file writer casewriter class. */
 748 static struct casewriter_class sys_file_casewriter_class =
 749   {
 750     sys_file_casewriter_write,
 751     sys_file_casewriter_destroy,
 752     NULL,
 753   };
 754 \f
 755 /* Writes case C to system file W, without compressing it. */
 756 static void
 757 write_case_uncompressed (struct sfm_writer *w, struct ccase *c)
 758 {
 759   size_t i;
 760
 761   for (i = 0; i < w->sfm_var_cnt; i++)
 762     {
 763       struct sfm_var *v = &w->sfm_vars[i];
 764
 765       if (v->width == 0)
 766         write_float (w, case_num_idx (c, v->case_index));
 767       else
 768         {
 769           write_bytes (w, case_str_idx (c, v->case_index) + v->offset,
 770                        v->width);
 771           write_spaces (w, v->padding);
 772         }
 773     }
 774 }
 775
 776 /* Writes case C to system file W, with compression. */
 777 static void
 778 write_case_compressed (struct sfm_writer *w, struct ccase *c)
 779 {
 780   size_t i;
 781
 782   for (i = 0; i < w->sfm_var_cnt; i++)
 783     {
 784       struct sfm_var *v = &w->sfm_vars[i];
 785
 786       if (v->width == 0)
 787         {
 788           double d = case_num_idx (c, v->case_index);
 789           if (d == SYSMIS)
 790             put_cmp_opcode (w, 255);
 791           else if (d >= 1 - COMPRESSION_BIAS
 792                    && d <= 251 - COMPRESSION_BIAS
 793                    && d == (int) d)
 794             put_cmp_opcode (w, (int) d + COMPRESSION_BIAS);
 795           else
 796             {
 797               put_cmp_opcode (w, 253);
 798               put_cmp_number (w, d);
 799             }
 800         }
 801       else
 802         {
 803           int offset = v->offset;
 804           int width, padding;
 805
 806           /* This code properly deals with a width that is not a
 807              multiple of 8, by ensuring that the final partial
 808              oct (8 byte unit) is treated as padded with spaces
 809              on the right. */
 810           for (width = v->width; width > 0; width -= 8, offset += 8)
 811             {
 812               const void *data = case_str_idx (c, v->case_index) + offset;
 813               int chunk_size = MIN (width, 8);
 814               if (!memcmp (data, "        ", chunk_size))
 815                 put_cmp_opcode (w, 254);
 816               else
 817                 {
 818                   put_cmp_opcode (w, 253);
 819                   put_cmp_string (w, data, chunk_size);
 820                 }
 821             }
 822
 823           /* This code deals properly with padding that is not a
 824              multiple of 8 bytes, by discarding the remainder,
 825              which was already effectively padded with spaces in
 826              the previous loop.  (Note that v->width + v->padding
 827              is always a multiple of 8.) */
 828           for (padding = v->padding / 8; padding > 0; padding--)
 829             put_cmp_opcode (w, 254);
 830         }
 831     }
 832 }
 833
 834 /* Flushes buffered compressed opcodes and data to W.
 835    The compression buffer must not be empty. */
 836 static void
 837 flush_compressed (struct sfm_writer *w)
 838 {
 839   assert (w->opcode_cnt > 0 && w->opcode_cnt <= 8);
 840
 841   write_bytes (w, w->opcodes, w->opcode_cnt);
 842   write_zeros (w, 8 - w->opcode_cnt);
 843
 844   write_bytes (w, w->data, w->data_cnt * sizeof *w->data);
 845
 846   w->opcode_cnt = w->data_cnt = 0;
 847 }
 848
 849 /* Appends OPCODE to the buffered set of compression opcodes in
 850    W.  Flushes the compression buffer beforehand if necessary. */
 851 static void
 852 put_cmp_opcode (struct sfm_writer *w, uint8_t opcode)
 853 {
 854   if (w->opcode_cnt >= 8)
 855     flush_compressed (w);
 856
 857   w->opcodes[w->opcode_cnt++] = opcode;
 858 }
 859
 860 /* Appends NUMBER to the buffered compression data in W.  The
 861    buffer must not be full; the way to assure that is to call
 862    this function only just after a call to put_cmp_opcode, which
 863    will flush the buffer as necessary. */
 864 static void
 865 put_cmp_number (struct sfm_writer *w, double number)
 866 {
 867   assert (w->opcode_cnt > 0);
 868   assert (w->data_cnt < 8);
 869
 870   convert_double_to_output_format (number, w->data[w->data_cnt++]);
 871 }
 872
 873 /* Appends SIZE bytes of DATA to the buffered compression data in
 874    W, followed by enough spaces to pad the output data to exactly
 875    8 bytes (thus, SIZE must be no greater than 8).  The buffer
 876    must not be full; the way to assure that is to call this
 877    function only just after a call to put_cmp_opcode, which will
 878    flush the buffer as necessary. */
 879 static void
 880 put_cmp_string (struct sfm_writer *w, const void *data, size_t size)
 881 {
 882   assert (w->opcode_cnt > 0);
 883   assert (w->data_cnt < 8);
 884   assert (size <= 8);
 885
 886   memset (w->data[w->data_cnt], ' ', 8);
 887   memcpy (w->data[w->data_cnt], data, size);
 888   w->data_cnt++;
 889 }
 890 \f
 891 /* Writes 32-bit integer X to the output file for writer W. */
 892 static void
 893 write_int (struct sfm_writer *w, int32_t x)
 894 {
 895   write_bytes (w, &x, sizeof x);
 896 }
 897
 898 /* Converts NATIVE to the 64-bit format used in output files in
 899    OUTPUT. */
 900 static inline void
 901 convert_double_to_output_format (double native, uint8_t output[8])
 902 {
 903   /* If "double" is not a 64-bit type, then convert it to a
 904      64-bit type.  Otherwise just copy it. */
 905   if (FLOAT_NATIVE_DOUBLE != FLOAT_NATIVE_64_BIT)
 906     float_convert (FLOAT_NATIVE_DOUBLE, &native, FLOAT_NATIVE_64_BIT, output);
 907   else
 908     memcpy (output, &native, sizeof native);
 909 }
 910
 911 /* Writes floating-point number X to the output file for writer
 912    W. */
 913 static void
 914 write_float (struct sfm_writer *w, double x)
 915 {
 916   uint8_t output[8];
 917   convert_double_to_output_format (x, output);
 918   write_bytes (w, output, sizeof output);
 919 }
 920
 921 /* Writes contents of VALUE with the given WIDTH to W, padding
 922    with zeros to a multiple of 8 bytes.
 923    To avoid a branch, and because we don't actually need to
 924    support it, WIDTH must be no bigger than 8. */
 925 static void
 926 write_value (struct sfm_writer *w, const union value *value, int width)
 927 {
 928   assert (width <= 8);
 929   if (width == 0)
 930     write_float (w, value->f);
 931   else
 932     {
 933       write_bytes (w, value->s, width);
 934       write_zeros (w, 8 - width);
 935     }
 936 }
 937
 938 /* Writes null-terminated STRING in a field of the given WIDTH to
 939    W.  If WIDTH is longer than WIDTH, it is truncated; if WIDTH
 940    is narrowed, it is padded on the right with spaces. */
 941 static void
 942 write_string (struct sfm_writer *w, const char *string, size_t width)
 943 {
 944   size_t data_bytes = MIN (strlen (string), width);
 945   size_t pad_bytes = width - data_bytes;
 946   write_bytes (w, string, data_bytes);
 947   while (pad_bytes-- > 0)
 948     putc (' ', w->file);
 949 }
 950
 951 /* Writes SIZE bytes of DATA to W's output file. */
 952 static void
 953 write_bytes (struct sfm_writer *w, const void *data, size_t size)
 954 {
 955   fwrite (data, 1, size, w->file);
 956 }
 957
 958 /* Writes N zeros to W's output file. */
 959 static void
 960 write_zeros (struct sfm_writer *w, size_t n)
 961 {
 962   while (n-- > 0)
 963     putc (0, w->file);
 964 }
 965
 966 /* Writes N spaces to W's output file. */
 967 static void
 968 write_spaces (struct sfm_writer *w, size_t n)
 969 {
 970   while (n-- > 0)
 971     putc (' ', w->file);
 972 }