7cdd7d9c5ecb93a826c4946a85cf1a1165240e6b
[pspp-builds.git] / utilities / pspp-dump-sav.c
1 /* PSPP - a program for statistical analysis.
2    Copyright (C) 2007, 2008, 2009, 2010, 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 <ctype.h>
20 #include <errno.h>
21 #include <getopt.h>
22 #include <inttypes.h>
23 #include <limits.h>
24 #include <stdio.h>
25 #include <stdlib.h>
26
27 #include "data/val-type.h"
28 #include "libpspp/cast.h"
29 #include "libpspp/compiler.h"
30 #include "libpspp/float-format.h"
31 #include "libpspp/integer-format.h"
32 #include "libpspp/misc.h"
33
34 #include "gl/error.h"
35 #include "gl/minmax.h"
36 #include "gl/progname.h"
37 #include "gl/version-etc.h"
38 #include "gl/xalloc.h"
39
40 #define ID_MAX_LEN 64
41
42 struct sfm_reader
43   {
44     const char *file_name;
45     FILE *file;
46
47     int n_variable_records, n_variables;
48
49     int *var_widths;
50     size_t n_var_widths, allocated_var_widths;
51
52     enum integer_format integer_format;
53     enum float_format float_format;
54
55     bool compressed;
56     double bias;
57   };
58
59 static void read_header (struct sfm_reader *);
60 static void read_variable_record (struct sfm_reader *);
61 static void read_value_label_record (struct sfm_reader *);
62 static void read_document_record (struct sfm_reader *);
63 static void read_extension_record (struct sfm_reader *);
64 static void read_machine_integer_info (struct sfm_reader *,
65                                        size_t size, size_t count);
66 static void read_machine_float_info (struct sfm_reader *,
67                                      size_t size, size_t count);
68 static void read_mrsets (struct sfm_reader *, size_t size, size_t count);
69 static void read_display_parameters (struct sfm_reader *,
70                                      size_t size, size_t count);
71 static void read_long_var_name_map (struct sfm_reader *r,
72                                     size_t size, size_t count);
73 static void read_long_string_map (struct sfm_reader *r,
74                                   size_t size, size_t count);
75 static void read_datafile_attributes (struct sfm_reader *r,
76                                       size_t size, size_t count);
77 static void read_variable_attributes (struct sfm_reader *r,
78                                       size_t size, size_t count);
79 static void read_ncases64 (struct sfm_reader *, size_t size, size_t count);
80 static void read_character_encoding (struct sfm_reader *r,
81                                        size_t size, size_t count);
82 static void read_long_string_value_labels (struct sfm_reader *r,
83                                            size_t size, size_t count);
84 static void read_unknown_extension (struct sfm_reader *,
85                                     size_t size, size_t count);
86 static void read_compressed_data (struct sfm_reader *, int max_cases);
87
88 static struct text_record *open_text_record (
89   struct sfm_reader *, size_t size);
90 static void close_text_record (struct text_record *);
91 static bool read_variable_to_value_pair (struct text_record *,
92                                          char **key, char **value);
93 static char *text_tokenize (struct text_record *, int delimiter);
94 static bool text_match (struct text_record *text, int c);
95 static const char *text_parse_counted_string (struct text_record *);
96 static size_t text_pos (const struct text_record *);
97
98 static void usage (void);
99 static void sys_warn (struct sfm_reader *, const char *, ...)
100      PRINTF_FORMAT (2, 3);
101 static void sys_error (struct sfm_reader *, const char *, ...)
102      PRINTF_FORMAT (2, 3)
103      NO_RETURN;
104
105 static void read_bytes (struct sfm_reader *, void *, size_t);
106 static bool try_read_bytes (struct sfm_reader *, void *, size_t);
107 static int read_int (struct sfm_reader *);
108 static int64_t read_int64 (struct sfm_reader *);
109 static double read_float (struct sfm_reader *);
110 static void read_string (struct sfm_reader *, char *, size_t);
111 static void skip_bytes (struct sfm_reader *, size_t);
112 static void trim_spaces (char *);
113
114 int
115 main (int argc, char *argv[])
116 {
117   int max_cases = 0;
118   struct sfm_reader r;
119   int i;
120
121   set_program_name (argv[0]);
122
123   for (;;)
124     {
125       static const struct option long_options[] =
126         {
127           { "data",    optional_argument, NULL, 'd' },
128           { "help",    no_argument,       NULL, 'h' },
129           { "version", no_argument,       NULL, 'v' },
130           { NULL,      0,                 NULL, 0 },
131         };
132
133       int c;
134
135       c = getopt_long (argc, argv, "d::hv", long_options, NULL);
136       if (c == -1)
137         break;
138
139       switch (c)
140         {
141         case 'd':
142           max_cases = optarg ? atoi (optarg) : INT_MAX;
143           break;
144
145         case 'v':
146           version_etc (stdout, "pspp-dump-sav", PACKAGE_NAME, PACKAGE_VERSION,
147                        "Ben Pfaff", "John Darrington", NULL_SENTINEL);
148           exit (EXIT_SUCCESS);
149
150         case 'h':
151           usage ();
152           exit (EXIT_SUCCESS);
153
154         default:
155           exit (EXIT_FAILURE);
156         }
157     }
158
159   if (optind == argc)
160     error (1, 0, "at least one non-option argument is required; "
161            "use --help for help");
162
163   for (i = optind; i < argc; i++)
164     {
165       int rec_type;
166
167       r.file_name = argv[i];
168       r.file = fopen (r.file_name, "rb");
169       if (r.file == NULL)
170         error (EXIT_FAILURE, errno, "error opening `%s'", r.file_name);
171       r.n_variable_records = 0;
172       r.n_variables = 0;
173       r.n_var_widths = 0;
174       r.allocated_var_widths = 0;
175       r.var_widths = 0;
176       r.compressed = false;
177
178       if (argc - optind > 1)
179         printf ("Reading \"%s\":\n", r.file_name);
180       
181       read_header (&r);
182       while ((rec_type = read_int (&r)) != 999)
183         {
184           switch (rec_type)
185             {
186             case 2:
187               read_variable_record (&r);
188               break;
189
190             case 3:
191               read_value_label_record (&r);
192               break;
193
194             case 4:
195               sys_error (&r, "Misplaced type 4 record.");
196
197             case 6:
198               read_document_record (&r);
199               break;
200
201             case 7:
202               read_extension_record (&r);
203               break;
204
205             default:
206               sys_error (&r, "Unrecognized record type %d.", rec_type);
207             }
208         }
209       printf ("%08llx: end-of-dictionary record "
210               "(first byte of data at %08llx)\n",
211               (long long int) ftello (r.file),
212               (long long int) ftello (r.file) + 4);
213
214       if (r.compressed && max_cases > 0)
215         read_compressed_data (&r, max_cases);
216
217       fclose (r.file);
218     }
219   
220   return 0;
221 }
222
223 static void
224 read_header (struct sfm_reader *r)
225 {
226   char rec_type[5];
227   char eye_catcher[61];
228   uint8_t raw_layout_code[4];
229   int32_t layout_code;
230   int32_t nominal_case_size;
231   int32_t compressed;
232   int32_t weight_index;
233   int32_t ncases;
234   uint8_t raw_bias[8];
235   char creation_date[10];
236   char creation_time[9];
237   char file_label[65];
238
239   read_string (r, rec_type, sizeof rec_type);
240   read_string (r, eye_catcher, sizeof eye_catcher);
241
242   if (strcmp ("$FL2", rec_type) != 0)
243     sys_error (r, "This is not an SPSS system file.");
244
245   /* Identify integer format. */
246   read_bytes (r, raw_layout_code, sizeof raw_layout_code);
247   if ((!integer_identify (2, raw_layout_code, sizeof raw_layout_code,
248                           &r->integer_format)
249        && !integer_identify (3, raw_layout_code, sizeof raw_layout_code,
250                              &r->integer_format))
251       || (r->integer_format != INTEGER_MSB_FIRST
252           && r->integer_format != INTEGER_LSB_FIRST))
253     sys_error (r, "This is not an SPSS system file.");
254   layout_code = integer_get (r->integer_format,
255                              raw_layout_code, sizeof raw_layout_code);
256
257   nominal_case_size = read_int (r);
258   compressed = read_int (r);
259   weight_index = read_int (r);
260   ncases = read_int (r);
261
262   r->compressed = compressed != 0;
263
264   /* Identify floating-point format and obtain compression bias. */
265   read_bytes (r, raw_bias, sizeof raw_bias);
266   if (float_identify (100.0, raw_bias, sizeof raw_bias, &r->float_format) == 0)
267     {
268       sys_warn (r, "Compression bias is not the usual value of 100, or system "
269                 "file uses unrecognized floating-point format.");
270       if (r->integer_format == INTEGER_MSB_FIRST)
271         r->float_format = FLOAT_IEEE_DOUBLE_BE;
272       else
273         r->float_format = FLOAT_IEEE_DOUBLE_LE;
274     }
275   r->bias = float_get_double (r->float_format, raw_bias);
276
277   read_string (r, creation_date, sizeof creation_date);
278   read_string (r, creation_time, sizeof creation_time);
279   read_string (r, file_label, sizeof file_label);
280   trim_spaces (file_label);
281   skip_bytes (r, 3);
282
283   printf ("File header record:\n");
284   printf ("\t%17s: %s\n", "Product name", eye_catcher);
285   printf ("\t%17s: %"PRId32"\n", "Layout code", layout_code);
286   printf ("\t%17s: %"PRId32"\n", "Compressed", compressed);
287   printf ("\t%17s: %"PRId32"\n", "Weight index", weight_index);
288   printf ("\t%17s: %"PRId32"\n", "Number of cases", ncases);
289   printf ("\t%17s: %g\n", "Compression bias", r->bias);
290   printf ("\t%17s: %s\n", "Creation date", creation_date);
291   printf ("\t%17s: %s\n", "Creation time", creation_time);
292   printf ("\t%17s: \"%s\"\n", "File label", file_label);
293 }
294
295 static const char *
296 format_name (int format)
297 {
298   switch ((format >> 16) & 0xff)
299     {
300     case 1: return "A";
301     case 2: return "AHEX";
302     case 3: return "COMMA";
303     case 4: return "DOLLAR";
304     case 5: return "F";
305     case 6: return "IB";
306     case 7: return "PIBHEX";
307     case 8: return "P";
308     case 9: return "PIB";
309     case 10: return "PK";
310     case 11: return "RB";
311     case 12: return "RBHEX";
312     case 15: return "Z";
313     case 16: return "N";
314     case 17: return "E";
315     case 20: return "DATE";
316     case 21: return "TIME";
317     case 22: return "DATETIME";
318     case 23: return "ADATE";
319     case 24: return "JDATE";
320     case 25: return "DTIME";
321     case 26: return "WKDAY";
322     case 27: return "MONTH";
323     case 28: return "MOYR";
324     case 29: return "QYR";
325     case 30: return "WKYR";
326     case 31: return "PCT";
327     case 32: return "DOT";
328     case 33: return "CCA";
329     case 34: return "CCB";
330     case 35: return "CCC";
331     case 36: return "CCD";
332     case 37: return "CCE";
333     case 38: return "EDATE";
334     case 39: return "SDATE";
335     default: return "invalid";
336     }
337 }
338
339 /* Reads a variable (type 2) record from R and adds the
340    corresponding variable to DICT.
341    Also skips past additional variable records for long string
342    variables. */
343 static void
344 read_variable_record (struct sfm_reader *r)
345 {
346   int width;
347   int has_variable_label;
348   int missing_value_code;
349   int print_format;
350   int write_format;
351   char name[9];
352
353   printf ("%08llx: variable record #%d\n",
354           (long long int) ftello (r->file), r->n_variable_records++);
355
356   width = read_int (r);
357   has_variable_label = read_int (r);
358   missing_value_code = read_int (r);
359   print_format = read_int (r);
360   write_format = read_int (r);
361   read_string (r, name, sizeof name);
362   name[strcspn (name, " ")] = '\0';
363
364   if (width >= 0)
365     r->n_variables++;
366
367   if (r->n_var_widths >= r->allocated_var_widths)
368     r->var_widths = x2nrealloc (r->var_widths, &r->allocated_var_widths,
369                                 sizeof *r->var_widths);
370   r->var_widths[r->n_var_widths++] = width;
371
372   printf ("\tWidth: %d (%s)\n",
373           width,
374           width > 0 ? "string"
375           : width == 0 ? "numeric"
376           : "long string continuation record");
377   printf ("\tVariable label: %d\n", has_variable_label);
378   printf ("\tMissing values code: %d (%s)\n", missing_value_code,
379           (missing_value_code == 0 ? "no missing values"
380            : missing_value_code == 1 ? "one missing value"
381            : missing_value_code == 2 ? "two missing values"
382            : missing_value_code == 3 ? "three missing values"
383            : missing_value_code == -2 ? "one missing value range"
384            : missing_value_code == -3 ? "one missing value, one range"
385            : "bad value"));
386   printf ("\tPrint format: %06x (%s%d.%d)\n",
387           print_format, format_name (print_format),
388           (print_format >> 8) & 0xff, print_format & 0xff);
389   printf ("\tWrite format: %06x (%s%d.%d)\n",
390           write_format, format_name (write_format),
391           (write_format >> 8) & 0xff, write_format & 0xff);
392   printf ("\tName: %s\n", name);
393
394   /* Get variable label, if any. */
395   if (has_variable_label != 0 && has_variable_label != 1)
396     sys_error (r, "Variable label indicator field is not 0 or 1.");
397   if (has_variable_label == 1)
398     {
399       long long int offset = ftello (r->file);
400       size_t len, read_len;
401       char label[255 + 1];
402
403       len = read_int (r);
404
405       /* Read up to 255 bytes of label. */
406       read_len = MIN (sizeof label - 1, len);
407       read_string (r, label, read_len + 1);
408       printf("\t%08llx Variable label: \"%s\"\n", offset, label);
409
410       /* Skip unread label bytes. */
411       skip_bytes (r, len - read_len);
412
413       /* Skip label padding up to multiple of 4 bytes. */
414       skip_bytes (r, ROUND_UP (len, 4) - len);
415     }
416
417   /* Set missing values. */
418   if (missing_value_code != 0)
419     {
420       int i;
421
422       printf ("\t%08llx Missing values:", (long long int) ftello (r->file));
423       if (!width)
424         {
425           if (missing_value_code < -3 || missing_value_code > 3
426               || missing_value_code == -1)
427             sys_error (r, "Numeric missing value indicator field is not "
428                        "-3, -2, 0, 1, 2, or 3.");
429           if (missing_value_code < 0)
430             {
431               double low = read_float (r);
432               double high = read_float (r);
433               printf (" %g...%g", low, high);
434               missing_value_code = -missing_value_code - 2;
435             }
436           for (i = 0; i < missing_value_code; i++)
437             printf (" %g", read_float (r));
438         }
439       else if (width > 0)
440         {
441           if (missing_value_code < 1 || missing_value_code > 3)
442             sys_error (r, "String missing value indicator field is not "
443                        "0, 1, 2, or 3.");
444           for (i = 0; i < missing_value_code; i++)
445             {
446               char string[9];
447               read_string (r, string, sizeof string);
448               printf (" \"%s\"", string);
449             }
450         }
451       putchar ('\n');
452     }
453 }
454
455 static void
456 print_untyped_value (struct sfm_reader *r, char raw_value[8])
457 {
458   int n_printable;
459   double value;
460
461   value = float_get_double (r->float_format, raw_value);
462   for (n_printable = 0; n_printable < sizeof raw_value; n_printable++)
463     if (!isprint (raw_value[n_printable]))
464       break;
465
466   printf ("%g/\"%.*s\"", value, n_printable, raw_value);
467 }
468
469 /* Reads value labels from sysfile R and inserts them into the
470    associated dictionary. */
471 static void
472 read_value_label_record (struct sfm_reader *r)
473 {
474   int label_cnt, var_cnt;
475   int i;
476
477   printf ("%08llx: value labels record\n", (long long int) ftello (r->file));
478
479   /* Read number of labels. */
480   label_cnt = read_int (r);
481   for (i = 0; i < label_cnt; i++)
482     {
483       char raw_value[8];
484       unsigned char label_len;
485       size_t padded_len;
486       char label[256];
487
488       read_bytes (r, raw_value, sizeof raw_value);
489
490       /* Read label length. */
491       read_bytes (r, &label_len, sizeof label_len);
492       padded_len = ROUND_UP (label_len + 1, 8);
493
494       /* Read label, padding. */
495       read_bytes (r, label, padded_len - 1);
496       label[label_len] = 0;
497
498       printf ("\t");
499       print_untyped_value (r, raw_value);
500       printf (": \"%s\"\n", label);
501     }
502
503   /* Now, read the type 4 record that has the list of variables
504      to which the value labels are to be applied. */
505
506   /* Read record type of type 4 record. */
507   if (read_int (r) != 4)
508     sys_error (r, "Variable index record (type 4) does not immediately "
509                "follow value label record (type 3) as it should.");
510
511   /* Read number of variables associated with value label from type 4
512      record. */
513   printf ("\t%08llx: apply to variables", (long long int) ftello (r->file));
514   var_cnt = read_int (r);
515   for (i = 0; i < var_cnt; i++)
516     printf (" #%d", read_int (r));
517   putchar ('\n');
518 }
519
520 static void
521 read_document_record (struct sfm_reader *r)
522 {
523   int n_lines;
524   int i;
525
526   printf ("%08llx: document record\n", (long long int) ftello (r->file));
527   n_lines = read_int (r);
528   printf ("\t%d lines of documents\n", n_lines);
529
530   for (i = 0; i < n_lines; i++)
531     {
532       char line[81];
533       printf ("\t%08llx: ", (long long int) ftello (r->file));
534       read_string (r, line, sizeof line);
535       trim_spaces (line);
536       printf ("line %d: \"%s\"\n", i, line);
537     }
538 }
539
540 static void
541 read_extension_record (struct sfm_reader *r)
542 {
543   long long int offset = ftello (r->file);
544   int subtype = read_int (r);
545   size_t size = read_int (r);
546   size_t count = read_int (r);
547   size_t bytes = size * count;
548
549   printf ("%08llx: Record 7, subtype %d, size=%zu, count=%zu\n",
550           offset, subtype, size, count);
551
552   switch (subtype)
553     {
554     case 3:
555       read_machine_integer_info (r, size, count);
556       return;
557
558     case 4:
559       read_machine_float_info (r, size, count);
560       return;
561
562     case 5:
563       /* Variable sets information.  We don't use these yet.
564          They only apply to GUIs; see VARSETS on the APPLY
565          DICTIONARY command in SPSS documentation. */
566       break;
567
568     case 6:
569       /* DATE variable information.  We don't use it yet, but we
570          should. */
571       break;
572
573     case 7:
574     case 19:
575       read_mrsets (r, size, count);
576       return;
577
578     case 11:
579       read_display_parameters (r, size, count);
580       return;
581
582     case 13:
583       read_long_var_name_map (r, size, count);
584       return;
585
586     case 14:
587       read_long_string_map (r, size, count);
588       return;
589
590     case 16:
591       read_ncases64 (r, size, count);
592       return;
593
594     case 17:
595       read_datafile_attributes (r, size, count);
596       return;
597
598     case 18:
599       read_variable_attributes (r, size, count);
600       return;
601
602     case 20:
603       read_character_encoding (r, size, count);
604       return;
605
606     case 21:
607       read_long_string_value_labels (r, size, count);
608       return;
609
610     default:
611       sys_warn (r, "Unrecognized record type 7, subtype %d.", subtype);
612       read_unknown_extension (r, size, count);
613       return;
614     }
615
616   skip_bytes (r, bytes);
617 }
618
619 static void
620 read_machine_integer_info (struct sfm_reader *r, size_t size, size_t count)
621 {
622   long long int offset = ftello (r->file);
623   int version_major = read_int (r);
624   int version_minor = read_int (r);
625   int version_revision = read_int (r);
626   int machine_code = read_int (r);
627   int float_representation = read_int (r);
628   int compression_code = read_int (r);
629   int integer_representation = read_int (r);
630   int character_code = read_int (r);
631
632   printf ("%08llx: machine integer info\n", offset);
633   if (size != 4 || count != 8)
634     sys_error (r, "Bad size (%zu) or count (%zu) field on record type 7, "
635                "subtype 3.", size, count);
636
637   printf ("\tVersion: %d.%d.%d\n",
638           version_major, version_minor, version_revision);
639   printf ("\tMachine code: %d\n", machine_code);
640   printf ("\tFloating point representation: %d (%s)\n",
641           float_representation,
642           float_representation == 1 ? "IEEE 754"
643           : float_representation == 2 ? "IBM 370"
644           : float_representation == 3 ? "DEC VAX"
645           : "unknown");
646   printf ("\tCompression code: %d\n", compression_code);
647   printf ("\tEndianness: %d (%s)\n", integer_representation,
648           integer_representation == 1 ? "big"
649           : integer_representation == 2 ? "little" : "unknown");
650   printf ("\tCharacter code: %d\n", character_code);
651 }
652
653 /* Read record type 7, subtype 4. */
654 static void
655 read_machine_float_info (struct sfm_reader *r, size_t size, size_t count)
656 {
657   long long int offset = ftello (r->file);
658   double sysmis = read_float (r);
659   double highest = read_float (r);
660   double lowest = read_float (r);
661
662   printf ("%08llx: machine float info\n", offset);
663   if (size != 8 || count != 3)
664     sys_error (r, "Bad size (%zu) or count (%zu) on extension 4.",
665                size, count);
666
667   printf ("\tsysmis: %g\n", sysmis);
668   if (sysmis != SYSMIS)
669     sys_warn (r, "File specifies unexpected value %g as %s.",
670               sysmis, "SYSMIS");
671
672   printf ("\thighest: %g\n", highest);
673   if (highest != HIGHEST)
674     sys_warn (r, "File specifies unexpected value %g as %s.",
675               highest, "HIGHEST");
676
677   printf ("\tlowest: %g\n", lowest);
678   if (lowest != LOWEST)
679     sys_warn (r, "File specifies unexpected value %g as %s.",
680               lowest, "LOWEST");
681 }
682
683 /* Read record type 7, subtype 7. */
684 static void
685 read_mrsets (struct sfm_reader *r, size_t size, size_t count)
686 {
687   struct text_record *text;
688
689   printf ("%08llx: multiple response sets\n",
690           (long long int) ftello (r->file));
691   text = open_text_record (r, size * count);
692   for (;;)
693     {
694       const char *name;
695       enum { MRSET_MC, MRSET_MD } type;
696       bool cat_label_from_counted_values = false;
697       bool label_from_var_label = false;
698       const char *counted;
699       const char *label;
700       const char *variables;
701
702       name = text_tokenize (text, '=');
703       if (name == NULL)
704         break;
705
706       if (text_match (text, 'C'))
707         {
708           type = MRSET_MC;
709           counted = NULL;
710           if (!text_match (text, ' '))
711             {
712               sys_warn (r, "missing space following 'C' at offset %zu "
713                         "in mrsets record", text_pos (text));
714               break;
715             }
716         }
717       else if (text_match (text, 'D'))
718         {
719           type = MRSET_MD;
720         }
721       else if (text_match (text, 'E'))
722         {
723           char *number;
724
725           type = MRSET_MD;
726           cat_label_from_counted_values = true;
727
728           if (!text_match (text, ' '))
729             {
730               sys_warn (r, "Missing space following `%c' at offset %zu "
731                         "in MRSETS record", 'E', text_pos (text));
732               break;
733             }
734
735           number = text_tokenize (text, ' ');
736           if (!strcmp (number, "11"))
737             label_from_var_label = true;
738           else if (strcmp (number, "1"))
739             sys_warn (r, "Unexpected label source value `%s' "
740                       "following `E' at offset %zu in MRSETS record",
741                       number, text_pos (text));
742
743         }
744       else
745         {
746           sys_warn (r, "missing `C', `D', or `E' at offset %zu "
747                     "in mrsets record", text_pos (text));
748           break;
749         }
750
751       if (type == MRSET_MD)
752         {
753           counted = text_parse_counted_string (text);
754           if (counted == NULL)
755             break;
756         }
757
758       label = text_parse_counted_string (text);
759       if (label == NULL)
760         break;
761
762       variables = text_tokenize (text, '\n');
763       if (variables == NULL)
764         {
765           sys_warn (r, "missing variable names following label "
766                     "at offset %zu in mrsets record", text_pos (text));
767           break;
768         }
769
770       printf ("\t\"%s\": multiple %s set",
771               name, type == MRSET_MC ? "category" : "dichotomy");
772       if (counted != NULL)
773         printf (", counted value \"%s\"", counted);
774       if (cat_label_from_counted_values)
775         printf (", category labels from counted values");
776       if (label[0] != '\0')
777         printf (", label \"%s\"", label);
778       if (label_from_var_label)
779         printf (", label from variable label");
780       printf(", variables \"%s\"\n", variables);
781     }
782   close_text_record (text);
783 }
784
785 /* Read record type 7, subtype 11. */
786 static void
787 read_display_parameters (struct sfm_reader *r, size_t size, size_t count)
788 {
789   size_t n_vars;
790   bool includes_width;
791   size_t i;
792
793   printf ("%08llx: variable display parameters\n",
794           (long long int) ftello (r->file));
795   if (size != 4)
796     {
797       sys_warn (r, "Bad size %zu on extension 11.", size);
798       skip_bytes (r, size * count);
799       return;
800     }
801
802   n_vars = r->n_variables;
803   if (count == 3 * n_vars)
804     includes_width = true;
805   else if (count == 2 * n_vars)
806     includes_width = false;
807   else
808     {
809       sys_warn (r, "Extension 11 has bad count %zu (for %zu variables.",
810                 count, n_vars);
811       skip_bytes (r, size * count);
812       return;
813     }
814
815   for (i = 0; i < n_vars; ++i)
816     {
817       int measure = read_int (r);
818       int width = includes_width ? read_int (r) : 0;
819       int align = read_int (r);
820
821       printf ("\tVar #%zu: measure=%d (%s)", i, measure,
822               (measure == 1 ? "nominal"
823                : measure == 2 ? "ordinal"
824                : measure == 3 ? "scale"
825                : "invalid"));
826       if (includes_width)
827         printf (", width=%d", width);
828       printf (", align=%d (%s)\n", align,
829               (align == 0 ? "left"
830                : align == 1 ? "right"
831                : align == 2 ? "centre"
832                : "invalid"));
833     }
834 }
835
836 /* Reads record type 7, subtype 13, which gives the long name
837    that corresponds to each short name.  */
838 static void
839 read_long_var_name_map (struct sfm_reader *r, size_t size, size_t count)
840 {
841   struct text_record *text;
842   char *var;
843   char *long_name;
844
845   printf ("%08llx: long variable names (short => long)\n",
846           (long long int) ftello (r->file));
847   text = open_text_record (r, size * count);
848   while (read_variable_to_value_pair (text, &var, &long_name))
849     printf ("\t%s => %s\n", var, long_name);
850   close_text_record (text);
851 }
852
853 /* Reads record type 7, subtype 14, which gives the real length
854    of each very long string.  Rearranges DICT accordingly. */
855 static void
856 read_long_string_map (struct sfm_reader *r, size_t size, size_t count)
857 {
858   struct text_record *text;
859   char *var;
860   char *length_s;
861
862   printf ("%08llx: very long strings (variable => length)\n",
863           (long long int) ftello (r->file));
864   text = open_text_record (r, size * count);
865   while (read_variable_to_value_pair (text, &var, &length_s))
866     printf ("\t%s => %d\n", var, atoi (length_s));
867   close_text_record (text);
868 }
869
870 static bool
871 read_attributes (struct sfm_reader *r, struct text_record *text,
872                  const char *variable)
873 {
874   const char *key;
875   int index;
876
877   for (;;) 
878     {
879       key = text_tokenize (text, '(');
880       if (key == NULL)
881         return true;
882   
883       for (index = 1; ; index++)
884         {
885           /* Parse the value. */
886           const char *value = text_tokenize (text, '\n');
887           if (value == NULL) 
888             {
889               sys_warn (r, "%s: Error parsing attribute value %s[%d]",
890                         variable, key, index);
891               return false;
892             }
893           if (strlen (value) < 2
894               || value[0] != '\'' || value[strlen (value) - 1] != '\'')
895             sys_warn (r, "%s: Attribute value %s[%d] is not quoted: %s",
896                       variable, key, index, value);
897           else
898             printf ("\t%s: %s[%d] = \"%.*s\"\n",
899                     variable, key, index, (int) strlen (value) - 2, value + 1);
900
901           /* Was this the last value for this attribute? */
902           if (text_match (text, ')'))
903             break;
904         }
905
906       if (text_match (text, '/'))
907         return true; 
908     }
909 }
910
911 /* Read extended number of cases record. */
912 static void
913 read_ncases64 (struct sfm_reader *r, size_t size, size_t count)
914 {
915   int64_t unknown, ncases64;
916
917   if (size != 8)
918     {
919       sys_warn (r, "Bad size %zu for extended number of cases.", size);
920       skip_bytes (r, size * count);
921       return;
922     }
923   if (count != 2)
924     {
925       sys_warn (r, "Bad count %zu for extended number of cases.", size);
926       skip_bytes (r, size * count);
927       return;
928     }
929   unknown = read_int64 (r);
930   ncases64 = read_int64 (r);
931   printf ("%08llx: extended number of cases: "
932           "unknown=%"PRId64", ncases64=%"PRId64"\n",
933           (long long int) ftello (r->file), unknown, ncases64);
934 }
935
936 static void
937 read_datafile_attributes (struct sfm_reader *r, size_t size, size_t count) 
938 {
939   struct text_record *text;
940   
941   printf ("%08llx: datafile attributes\n", (long long int) ftello (r->file));
942   text = open_text_record (r, size * count);
943   read_attributes (r, text, "datafile");
944   close_text_record (text);
945 }
946
947 static void
948 read_character_encoding (struct sfm_reader *r, size_t size, size_t count)
949 {
950   long long int posn =  ftello (r->file);
951   char *encoding = xcalloc (size, count + 1);
952   read_string (r, encoding, count + 1);
953
954   printf ("%08llx: Character Encoding: %s\n", posn, encoding);
955 }
956
957 static void
958 read_long_string_value_labels (struct sfm_reader *r, size_t size, size_t count)
959 {
960   long long int start = ftello (r->file);
961
962   printf ("%08llx: long string value labels\n", start);
963   while (ftello (r->file) - start < size * count)
964     {
965       long long posn = ftello (r->file);
966       char var_name[ID_MAX_LEN + 1];
967       int var_name_len;
968       int n_values;
969       int width;
970       int i;
971
972       /* Read variable name. */
973       var_name_len = read_int (r);
974       if (var_name_len > ID_MAX_LEN)
975         sys_error (r, "Variable name length in long string value label "
976                    "record (%d) exceeds %d-byte limit.",
977                    var_name_len, ID_MAX_LEN);
978       read_string (r, var_name, var_name_len + 1);
979
980       /* Read width, number of values. */
981       width = read_int (r);
982       n_values = read_int (r);
983
984       printf ("\t%08llx: %s, width %d, %d values\n",
985               posn, var_name, width, n_values);
986
987       /* Read values. */
988       for (i = 0; i < n_values; i++)
989         {
990           char *value;
991           int value_length;
992
993           char *label;
994           int label_length;
995
996           posn = ftello (r->file);
997
998           /* Read value. */
999           value_length = read_int (r);
1000           value = xmalloc (value_length + 1);
1001           read_string (r, value, value_length + 1);
1002
1003           /* Read label. */
1004           label_length = read_int (r);
1005           label = xmalloc (label_length + 1);
1006           read_string (r, label, label_length + 1);
1007
1008           printf ("\t\t%08llx: \"%s\" (%d bytes) => \"%s\" (%d bytes)\n",
1009                   posn, value, value_length, label, label_length);
1010
1011           free (value);
1012           free (label);
1013         }
1014     }
1015 }
1016
1017 static void
1018 hex_dump (size_t offset, const void *buffer_, size_t buffer_size)
1019 {
1020   const uint8_t *buffer = buffer_;
1021
1022   while (buffer_size > 0)
1023     {
1024       size_t n = MIN (buffer_size, 16);
1025       size_t i;
1026
1027       printf ("%04zx", offset);
1028       for (i = 0; i < 16; i++)
1029         {
1030           if (i < n)
1031             printf ("%c%02x", i == 8 ? '-' : ' ', buffer[i]);
1032           else
1033             printf ("   ");
1034         }
1035
1036       printf (" |");
1037       for (i = 0; i < 16; i++)
1038         {
1039           unsigned char c = i < n ? buffer[i] : ' ';
1040           putchar (isprint (c) ? c : '.');
1041         }
1042       printf ("|\n");
1043
1044       offset += n;
1045       buffer += n;
1046       buffer_size -= n;
1047     }
1048 }
1049
1050 /* Reads and prints any type 7 record that we don't understand. */
1051 static void
1052 read_unknown_extension (struct sfm_reader *r, size_t size, size_t count)
1053 {
1054   unsigned char *buffer;
1055   size_t i;
1056
1057   if (size == 0 || count > 65536 / size)
1058     skip_bytes (r, size * count);
1059   else if (size != 1)
1060     {
1061       buffer = xmalloc (size);
1062       for (i = 0; i < count; i++)
1063         {
1064           read_bytes (r, buffer, size);
1065           hex_dump (i * size, buffer, size);
1066         }
1067       free (buffer);
1068     }
1069   else
1070     {
1071       buffer = xmalloc (count);
1072       read_bytes (r, buffer, count);
1073       if (memchr (buffer, 0, count) == 0)
1074         for (i = 0; i < count; i++)
1075           {
1076             unsigned char c = buffer[i];
1077
1078             if (c == '\\')
1079               printf ("\\\\");
1080             else if (c == '\n' || isprint (c))
1081               putchar (c);
1082             else
1083               printf ("\\%02x", c);
1084           }
1085       else
1086         hex_dump (0, buffer, count);
1087       free (buffer);
1088     }
1089 }
1090
1091 static void
1092 read_variable_attributes (struct sfm_reader *r, size_t size, size_t count) 
1093 {
1094   struct text_record *text;
1095   
1096   printf ("%08llx: variable attributes\n", (long long int) ftello (r->file));
1097   text = open_text_record (r, size * count);
1098   for (;;) 
1099     {
1100       const char *variable = text_tokenize (text, ':');
1101       if (variable == NULL || !read_attributes (r, text, variable))
1102         break; 
1103     }
1104   close_text_record (text);
1105 }
1106
1107 static void
1108 read_compressed_data (struct sfm_reader *r, int max_cases)
1109 {
1110   enum { N_OPCODES = 8 };
1111   uint8_t opcodes[N_OPCODES];
1112   long long int opcode_ofs;
1113   int opcode_idx;
1114   int case_num;
1115   int i;
1116
1117   read_int (r);
1118   printf ("\n%08llx: compressed data:\n", (long long int) ftello (r->file));
1119
1120   opcode_idx = N_OPCODES;
1121   opcode_ofs = 0;
1122   case_num = 0;
1123   for (case_num = 0; case_num < max_cases; case_num++)
1124     {
1125       printf ("%08llx: case %d's uncompressible data begins\n",
1126               (long long int) ftello (r->file), case_num);
1127       for (i = 0; i < r->n_var_widths; )
1128         {
1129           int width = r->var_widths[i];
1130           char raw_value[8];
1131           int opcode;
1132
1133           if (opcode_idx >= N_OPCODES)
1134             {
1135               opcode_ofs = ftello (r->file);
1136               if (i == 0)
1137                 {
1138                   if (!try_read_bytes (r, opcodes, 8))
1139                     return;
1140                 }
1141               else
1142                 read_bytes (r, opcodes, 8);
1143               opcode_idx = 0;
1144             }
1145           opcode = opcodes[opcode_idx];
1146           printf ("%08llx: variable %d: opcode %d: ",
1147                   opcode_ofs + opcode_idx, i, opcode);
1148
1149           switch (opcode)
1150             {
1151             default:
1152               printf ("%g", opcode - r->bias);
1153               if (width != 0)
1154                 printf (", but this is a string variable (width=%d)", width);
1155               printf ("\n");
1156               i++;
1157               break;
1158
1159             case 0:
1160               printf ("ignored padding\n");
1161               break;
1162
1163             case 252:
1164               printf ("end of data\n");
1165               return;
1166
1167             case 253:
1168               read_bytes (r, raw_value, 8);
1169               printf ("uncompressible data: ");
1170               print_untyped_value (r, raw_value);
1171               printf ("\n");
1172               i++;
1173               break;
1174
1175             case 254:
1176               printf ("spaces");
1177               if (width == 0)
1178                 printf (", but this is a numeric variable");
1179               printf ("\n");
1180               i++;
1181               break;
1182
1183             case 255:
1184               printf ("SYSMIS");
1185               if (width != 0)
1186                 printf (", but this is a string variable (width=%d)", width);
1187               printf ("\n");
1188               i++;
1189               break;
1190             }
1191
1192           opcode_idx++;
1193         }
1194     }
1195 }
1196 \f
1197 /* Helpers for reading records that consist of structured text
1198    strings. */
1199
1200 /* State. */
1201 struct text_record
1202   {
1203     struct sfm_reader *reader;  /* Reader. */
1204     char *buffer;               /* Record contents. */
1205     size_t size;                /* Size of buffer. */
1206     size_t pos;                 /* Current position in buffer. */
1207   };
1208
1209 /* Reads SIZE bytes into a text record for R,
1210    and returns the new text record. */
1211 static struct text_record *
1212 open_text_record (struct sfm_reader *r, size_t size)
1213 {
1214   struct text_record *text = xmalloc (sizeof *text);
1215   char *buffer = xmalloc (size + 1);
1216   read_bytes (r, buffer, size);
1217   buffer[size] = '\0';
1218   text->reader = r;
1219   text->buffer = buffer;
1220   text->size = size;
1221   text->pos = 0;
1222   return text;
1223 }
1224
1225 /* Closes TEXT and frees its storage.
1226    Not really needed, because the pool will free the text record anyway,
1227    but can be used to free it earlier. */
1228 static void
1229 close_text_record (struct text_record *text)
1230 {
1231   free (text->buffer);
1232   free (text);
1233 }
1234
1235 static char *
1236 text_tokenize (struct text_record *text, int delimiter)
1237 {
1238   size_t start = text->pos;
1239   while (text->pos < text->size
1240          && text->buffer[text->pos] != delimiter
1241          && text->buffer[text->pos] != '\0')
1242     text->pos++;
1243   if (start == text->pos)
1244     return NULL;
1245   text->buffer[text->pos++] = '\0';
1246   return &text->buffer[start];
1247 }
1248
1249 static bool
1250 text_match (struct text_record *text, int c) 
1251 {
1252   if (text->pos < text->size && text->buffer[text->pos] == c) 
1253     {
1254       text->pos++;
1255       return true;
1256     }
1257   else
1258     return false;
1259 }
1260
1261 /* Reads a integer value expressed in decimal, then a space, then a string that
1262    consists of exactly as many bytes as specified by the integer, then a space,
1263    from TEXT.  Returns the string, null-terminated, as a subset of TEXT's
1264    buffer (so the caller should not free the string). */
1265 static const char *
1266 text_parse_counted_string (struct text_record *text)
1267 {
1268   size_t start;
1269   size_t n;
1270   char *s;
1271
1272   start = text->pos;
1273   n = 0;
1274   while (isdigit ((unsigned char) text->buffer[text->pos]))
1275     n = (n * 10) + (text->buffer[text->pos++] - '0');
1276   if (start == text->pos)
1277     {
1278       sys_error (text->reader, "expecting digit at offset %zu in record",
1279                  text->pos);
1280       return NULL;
1281     }
1282
1283   if (!text_match (text, ' '))
1284     {
1285       sys_error (text->reader, "expecting space at offset %zu in record",
1286                  text->pos);
1287       return NULL;
1288     }
1289
1290   if (text->pos + n > text->size)
1291     {
1292       sys_error (text->reader, "%zu-byte string starting at offset %zu "
1293                  "exceeds record length %zu", n, text->pos, text->size);
1294       return NULL;
1295     }
1296
1297   s = &text->buffer[text->pos];
1298   if (s[n] != ' ')
1299     {
1300       sys_error (text->reader, "expecting space at offset %zu following "
1301                  "%zu-byte string", text->pos + n, n);
1302       return NULL;
1303     }
1304   s[n] = '\0';
1305   text->pos += n + 1;
1306   return s;
1307 }
1308
1309 /* Reads a variable=value pair from TEXT.
1310    Looks up the variable in DICT and stores it into *VAR.
1311    Stores a null-terminated value into *VALUE. */
1312 static bool
1313 read_variable_to_value_pair (struct text_record *text,
1314                              char **key, char **value)
1315 {
1316   *key = text_tokenize (text, '=');
1317   *value = text_tokenize (text, '\t');
1318   if (!*key || !*value)
1319     return false;
1320
1321   while (text->pos < text->size
1322          && (text->buffer[text->pos] == '\t'
1323              || text->buffer[text->pos] == '\0'))
1324     text->pos++;
1325   return true;
1326 }
1327
1328 /* Returns the current byte offset inside the TEXT's string. */
1329 static size_t
1330 text_pos (const struct text_record *text)
1331 {
1332   return text->pos;
1333 }
1334 \f
1335 static void
1336 usage (void)
1337 {
1338   printf ("\
1339 %s, a utility for dissecting system files.\n\
1340 Usage: %s [OPTION]... SYSFILE...\n\
1341 where each SYSFILE is the name of a system file.\n\
1342 \n\
1343 Options:\n\
1344   --data[=MAXCASES]   print (up to MAXCASES cases of) compressed data\n\
1345   --help              display this help and exit\n\
1346   --version           output version information and exit\n",
1347           program_name, program_name);
1348 }
1349
1350 /* Displays a corruption message. */
1351 static void
1352 sys_msg (struct sfm_reader *r, const char *format, va_list args)
1353 {
1354   printf ("\"%s\" near offset 0x%llx: ",
1355           r->file_name, (long long int) ftello (r->file));
1356   vprintf (format, args);
1357   putchar ('\n');
1358 }
1359
1360 /* Displays a warning for the current file position. */
1361 static void
1362 sys_warn (struct sfm_reader *r, const char *format, ...)
1363 {
1364   va_list args;
1365
1366   va_start (args, format);
1367   sys_msg (r, format, args);
1368   va_end (args);
1369 }
1370
1371 /* Displays an error for the current file position,
1372    marks it as in an error state,
1373    and aborts reading it using longjmp. */
1374 static void
1375 sys_error (struct sfm_reader *r, const char *format, ...)
1376 {
1377   va_list args;
1378
1379   va_start (args, format);
1380   sys_msg (r, format, args);
1381   va_end (args);
1382
1383   exit (EXIT_FAILURE);
1384 }
1385 \f
1386 /* Reads BYTE_CNT bytes into BUF.
1387    Returns true if exactly BYTE_CNT bytes are successfully read.
1388    Aborts if an I/O error or a partial read occurs.
1389    If EOF_IS_OK, then an immediate end-of-file causes false to be
1390    returned; otherwise, immediate end-of-file causes an abort
1391    too. */
1392 static inline bool
1393 read_bytes_internal (struct sfm_reader *r, bool eof_is_ok,
1394                      void *buf, size_t byte_cnt)
1395 {
1396   size_t bytes_read = fread (buf, 1, byte_cnt, r->file);
1397   if (bytes_read == byte_cnt)
1398     return true;
1399   else if (ferror (r->file))
1400     sys_error (r, "System error: %s.", strerror (errno));
1401   else if (!eof_is_ok || bytes_read != 0)
1402     sys_error (r, "Unexpected end of file.");
1403   else
1404     return false;
1405 }
1406
1407 /* Reads BYTE_CNT into BUF.
1408    Aborts upon I/O error or if end-of-file is encountered. */
1409 static void
1410 read_bytes (struct sfm_reader *r, void *buf, size_t byte_cnt)
1411 {
1412   read_bytes_internal (r, false, buf, byte_cnt);
1413 }
1414
1415 /* Reads BYTE_CNT bytes into BUF.
1416    Returns true if exactly BYTE_CNT bytes are successfully read.
1417    Returns false if an immediate end-of-file is encountered.
1418    Aborts if an I/O error or a partial read occurs. */
1419 static bool
1420 try_read_bytes (struct sfm_reader *r, void *buf, size_t byte_cnt)
1421 {
1422   return read_bytes_internal (r, true, buf, byte_cnt);
1423 }
1424
1425 /* Reads a 32-bit signed integer from R and returns its value in
1426    host format. */
1427 static int
1428 read_int (struct sfm_reader *r)
1429 {
1430   uint8_t integer[4];
1431   read_bytes (r, integer, sizeof integer);
1432   return integer_get (r->integer_format, integer, sizeof integer);
1433 }
1434
1435 /* Reads a 64-bit signed integer from R and returns its value in
1436    host format. */
1437 static int64_t
1438 read_int64 (struct sfm_reader *r)
1439 {
1440   uint8_t integer[8];
1441   read_bytes (r, integer, sizeof integer);
1442   return integer_get (r->integer_format, integer, sizeof integer);
1443 }
1444
1445 /* Reads a 64-bit floating-point number from R and returns its
1446    value in host format. */
1447 static double
1448 read_float (struct sfm_reader *r)
1449 {
1450   uint8_t number[8];
1451   read_bytes (r, number, sizeof number);
1452   return float_get_double (r->float_format, number);
1453 }
1454
1455 /* Reads exactly SIZE - 1 bytes into BUFFER
1456    and stores a null byte into BUFFER[SIZE - 1]. */
1457 static void
1458 read_string (struct sfm_reader *r, char *buffer, size_t size)
1459 {
1460   assert (size > 0);
1461   read_bytes (r, buffer, size - 1);
1462   buffer[size - 1] = '\0';
1463 }
1464
1465 /* Skips BYTES bytes forward in R. */
1466 static void
1467 skip_bytes (struct sfm_reader *r, size_t bytes)
1468 {
1469   while (bytes > 0)
1470     {
1471       char buffer[1024];
1472       size_t chunk = MIN (sizeof buffer, bytes);
1473       read_bytes (r, buffer, chunk);
1474       bytes -= chunk;
1475     }
1476 }
1477
1478 static void
1479 trim_spaces (char *s)
1480 {
1481   char *end = strchr (s, '\0');
1482   while (end > s && end[-1] == ' ')
1483     end--;
1484   *end = '\0';
1485 }