/* PSPP - computes sample statistics.
Copyright (C) 1997-9, 2000, 2006 Free Software Foundation, Inc.
- Written by Ben Pfaff <blp@gnu.org>.
This program is free software; you can redistribute it and/or
modify it under the terms of the GNU General Public License as
#include <config.h>
#include "sys-file-writer.h"
-#include "sfm-private.h"
+#include "sys-file-private.h"
#include <ctype.h>
#include <errno.h>
#include "case.h"
#include "dictionary.h"
#include "file-handle-def.h"
+#include "format.h"
+#include "missing-values.h"
#include "settings.h"
#include "value-labels.h"
#include "variable.h"
-#include "stat-macros.h"
#include "minmax.h"
#include "gettext.h"
#define _(msgid) gettext (msgid)
+/* Find 64-bit floating-point type. */
+#if SIZEOF_FLOAT == 8
+ #define flt64 float
+ #define FLT64_MAX FLT_MAX
+#elif SIZEOF_DOUBLE == 8
+ #define flt64 double
+ #define FLT64_MAX DBL_MAX
+#elif SIZEOF_LONG_DOUBLE == 8
+ #define flt64 long double
+ #define FLT64_MAX LDBL_MAX
+#else
+ #error Which one of your basic types is 64-bit floating point?
+#endif
+
+/* Figure out SYSMIS value for flt64. */
+#include <libpspp/magic.h>
+#if SIZEOF_DOUBLE == 8
+#define second_lowest_flt64 second_lowest_value
+#else
+#error Must define second_lowest_flt64 for your architecture.
+#endif
+
+/* Record Type 1: General Information. */
+struct sysfile_header
+ {
+ char rec_type[4] ; /* 00: Record-type code, "$FL2". */
+ char prod_name[60] ; /* 04: Product identification. */
+ int32_t layout_code ; /* 40: 2. */
+ int32_t nominal_case_size ; /* 44: Number of `value's per case.
+ Note: some systems set this to -1 */
+ int32_t compress ; /* 48: 1=compressed, 0=not compressed. */
+ int32_t weight_idx ; /* 4c: 1-based index of weighting var, or 0. */
+ int32_t case_cnt ; /* 50: Number of cases, -1 if unknown. */
+ flt64 bias ; /* 54: Compression bias (100.0). */
+ char creation_date[9] ; /* 5c: `dd mmm yy' creation date of file. */
+ char creation_time[8] ; /* 65: `hh:mm:ss' 24-hour creation time. */
+ char file_label[64] ; /* 6d: File label. */
+ char padding[3] ; /* ad: Ignored padding. */
+ } ATTRIBUTE((packed)) ;
+
+/* Record Type 2: Variable. */
+struct sysfile_variable
+ {
+ int32_t rec_type ; /* 2. */
+ int32_t type ; /* 0=numeric, 1-255=string width,
+ -1=continued string. */
+ int32_t has_var_label ; /* 1=has a variable label, 0=doesn't. */
+ int32_t n_missing_values ; /* Missing value code of -3,-2,0,1,2, or 3. */
+ int32_t print ; /* Print format. */
+ int32_t write ; /* Write format. */
+ char name[SHORT_NAME_LEN] ; /* Variable name. */
+ /* The rest of the structure varies. */
+ } ATTRIBUTE((packed)) ;
+
/* Compression bias used by PSPP. Values between (1 -
COMPRESSION_BIAS) and (251 - COMPRESSION_BIAS) inclusive can be
compressed. */
var_flt64_cnt (const struct variable *v)
{
assert(sizeof(flt64) == MAX_SHORT_STRING);
- return width_to_bytes(var_get_width (v)) / MAX_SHORT_STRING ;
+ return sfm_width_to_bytes(var_get_width (v)) / MAX_SHORT_STRING ;
}
static inline int
struct sfm_var *sv = &w->vars[i];
sv->width = var_get_width (dv);
/* spss compatibility nonsense */
- if ( var_is_very_long_string (dv) )
+ if ( var_get_width (dv) >= MIN_VERY_LONG_STRING )
w->has_vls = true;
- sv->fv = dv->fv;
+ sv->fv = var_get_case_index (dv);
sv->flt64_cnt = var_flt64_cnt (dv);
}
int wcount = var_get_width (v);
do {
- struct variable var_cont = *v;
+ struct variable *var_cont = var_clone (v);
+ var_set_short_name (var_cont, var_get_short_name (v));
if ( var_is_alpha (v))
{
if ( 0 != count )
{
- var_clear_missing_values (&var_cont);
- var_set_short_name (&var_cont,
+ var_clear_missing_values (var_cont);
+ var_set_short_name (var_cont,
cont_var_name (var_get_short_name (v),
count));
- var_clear_label (&var_cont);
+ var_clear_label (var_cont);
w->var_cnt_vls++;
}
count++;
- if ( wcount > MAX_LONG_STRING )
+ if ( wcount >= MIN_VERY_LONG_STRING )
{
- var_set_width (&var_cont, MAX_LONG_STRING);
+ var_set_width (var_cont, MIN_VERY_LONG_STRING - 1);
wcount -= EFFECTIVE_LONG_STRING_LENGTH;
}
else
{
- var_set_width (&var_cont, wcount);
- wcount -= var_get_width (&var_cont);
+ var_set_width (var_cont, wcount);
+ wcount -= var_get_width (var_cont);
}
}
- write_variable (w, &var_cont);
+ write_variable (w, var_cont);
+ var_destroy (var_cont);
} while(wcount > 0);
}
const char *label = var_get_label (v);
sv.rec_type = 2;
- sv.type = MIN(var_get_width (v), MAX_LONG_STRING);
+ sv.type = MIN (var_get_width (v), MIN_VERY_LONG_STRING - 1);
sv.has_var_label = label != NULL;
mv_copy (&mv, var_get_missing_values (v));
memset (&sv.write, 0, sizeof sv.write);
memset (&sv.name, 0, sizeof sv.name);
- pad_count = DIV_RND_UP (MIN(var_get_width (v), MAX_LONG_STRING),
+ pad_count = DIV_RND_UP (MIN(var_get_width (v), MIN_VERY_LONG_STRING - 1),
(int) sizeof (flt64)) - 1;
for (i = 0; i < pad_count; i++)
buf_write (w, &sv, sizeof sv);
int32_t vars[1] ;
} ATTRIBUTE((packed));
+ const struct val_labs *val_labs;
struct val_labs_iterator *i;
struct value_label_rec *vlr;
struct var_idx_rec vir;
size_t vlr_size;
flt64 *loc;
- if (!val_labs_count (v->val_labs))
+ val_labs = var_get_value_labels (v);
+ if (val_labs == NULL)
return;
/* Pass 1: Count bytes. */
vlr_size = (sizeof (struct value_label_rec)
- + sizeof (flt64) * (val_labs_count (v->val_labs) - 1));
- for (vl = val_labs_first (v->val_labs, &i); vl != NULL;
- vl = val_labs_next (v->val_labs, &i))
+ + sizeof (flt64) * (val_labs_count (val_labs) - 1));
+ for (vl = val_labs_first (val_labs, &i); vl != NULL;
+ vl = val_labs_next (val_labs, &i))
vlr_size += ROUND_UP (strlen (vl->label) + 1, sizeof (flt64));
/* Pass 2: Copy bytes. */
vlr = xmalloc (vlr_size);
vlr->rec_type = 3;
- vlr->n_labels = val_labs_count (v->val_labs);
+ vlr->n_labels = val_labs_count (val_labs);
loc = vlr->labels;
- for (vl = val_labs_first_sorted (v->val_labs, &i); vl != NULL;
- vl = val_labs_next (v->val_labs, &i))
+ for (vl = val_labs_first_sorted (val_labs, &i); vl != NULL;
+ vl = val_labs_next (val_labs, &i))
{
size_t len = strlen (vl->label);
int32_t n_lines ; /* Number of lines of documents. */
} ATTRIBUTE((packed)) rec_6;
- const char *documents;
- size_t n_lines;
+ const char * documents = dict_get_documents (d);
+ size_t doc_bytes = strlen (documents);
- documents = dict_get_documents (d);
- n_lines = strlen (documents) / 80;
+ assert (doc_bytes % 80 == 0);
rec_6.rec_type = 6;
- rec_6.n_lines = n_lines;
+ rec_6.n_lines = doc_bytes / 80;
buf_write (w, &rec_6, sizeof rec_6);
- buf_write (w, documents, 80 * n_lines);
+ buf_write (w, documents, 80 * rec_6.n_lines);
}
/* Write the alignment, width and scale values */
v = dict_get_var(dict, i);
- params.measure = var_get_measure (v);
+ params.measure = (var_get_measure (v) == MEASURE_NOMINAL ? 1
+ : var_get_measure (v) == MEASURE_ORDINAL ? 2
+ : 3);
params.width = var_get_display_width (v);
- params.align = var_get_alignment (v);
+ params.align = (var_get_alignment (v) == ALIGN_LEFT ? 0
+ : var_get_alignment (v) == ALIGN_RIGHT ? 1
+ : 2);
buf_write (w, ¶ms, sizeof(params));
while (wcount > 0)
{
- params.width = wcount > MAX_LONG_STRING ? 32 : wcount;
+ params.width = wcount >= MIN_VERY_LONG_STRING ? 32 : wcount;
buf_write (w, ¶ms, sizeof(params));
{
const struct variable *v = dict_get_var (dict, i);
- if ( var_get_width (v) <= MAX_LONG_STRING )
+ if ( var_get_width (v) < MIN_VERY_LONG_STRING )
continue;
ds_put_format (&vls_length_map, "%s=%05d",
if (v->width == 0)
{
- *bounce_cur = case_num (c, v->fv);
+ *bounce_cur = case_num_idx (c, v->fv);
bounce_cur += v->flt64_cnt;
}
else
{ int ofs = 0;
while (ofs < v->width)
{
- int chunk = MIN (MAX_LONG_STRING, v->width - ofs);
+ int chunk = MIN (MIN_VERY_LONG_STRING - 1, v->width - ofs);
int nv = DIV_RND_UP (chunk, sizeof (flt64));
buf_copy_rpad ((char *) bounce_cur, nv * sizeof (flt64),
- case_data (c, v->fv)->s + ofs, chunk);
+ case_data_idx (c, v->fv)->s + ofs, chunk);
bounce_cur += nv;
ofs += chunk;
}