Added result_class parameter to tab_double and updated all callers. Removed tab_fixed
[pspp] / src / language / stats / means.c
index 7ff286d7f3d1feae8d5680b79f21e6bb65ea8a59..f1726b9af1aef0d02773ea7aa32726aafd4569f4 100644 (file)
@@ -1,5 +1,5 @@
 /* PSPP - a program for statistical analysis.
-   Copyright (C) 2011, 2012 Free Software Foundation, Inc.
+   Copyright (C) 2011, 2012, 2013 Free Software Foundation, Inc.
 
    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
@@ -384,7 +384,7 @@ first_update (void *stat, double w UNUSED, double x)
   double *f = stat;
 
   if (*f == SYSMIS)
-     *f = x;
+    *f = x;
 }
 
 static double
@@ -502,11 +502,12 @@ parse_means_table_syntax (struct lexer *lexer, const struct means *cmd, struct m
   table->ii = 0;
   table->n_layers = 0;
   table->layers = NULL;
+  table->interactions = NULL;
 
   /* Dependent variable (s) */
-  if (!parse_variables_const (lexer, cmd->dict,
-                             &table->dep_vars, &table->n_dep_vars,
-                             PV_NO_DUPLICATE | PV_NUMERIC))
+  if (!parse_variables_const_pool (lexer, cmd->pool, cmd->dict,
+                                  &table->dep_vars, &table->n_dep_vars,
+                                  PV_NO_DUPLICATE | PV_NUMERIC))
     return false;
 
   /* Factor variable (s) */
@@ -517,10 +518,10 @@ parse_means_table_syntax (struct lexer *lexer, const struct means *cmd, struct m
          table->n_layers++;
           table->layers = 
            pool_realloc (cmd->pool, table->layers, 
-                     sizeof (*table->layers) * table->n_layers);
+                         sizeof (*table->layers) * table->n_layers);
 
-         if (!parse_variables_const 
-              (lexer, cmd->dict,
+         if (!parse_variables_const_pool 
+              (lexer, cmd->pool, cmd->dict,
                &table->layers[table->n_layers - 1].factor_vars,
                &table->layers[table->n_layers - 1].n_factor_vars,
                PV_NO_DUPLICATE))
@@ -538,7 +539,7 @@ parse_means_table_syntax (struct lexer *lexer, const struct means *cmd, struct m
   table->n_layers++;
   table->layers = 
     pool_realloc (cmd->pool, table->layers, 
-              sizeof (*table->layers) * table->n_layers);
+                 sizeof (*table->layers) * table->n_layers);
   table->layers[table->n_layers - 1].factor_vars = NULL;
   table->layers[table->n_layers - 1].n_factor_vars = 0;
 
@@ -633,10 +634,10 @@ cmd_means (struct lexer *lexer, struct dataset *ds)
       if (lex_match_id (lexer, "MISSING"))
        {
          /*
-            If no MISSING subcommand is specified, each combination of
-            a dependent variable and categorical variables is handled
-            separately.
-          */
+           If no MISSING subcommand is specified, each combination of
+           a dependent variable and categorical variables is handled
+           separately.
+         */
          lex_match (lexer, T_EQUALS);
          if (lex_match_id (lexer, "INCLUDE"))
            {
@@ -661,17 +662,17 @@ cmd_means (struct lexer *lexer, struct dataset *ds)
            }
          else if (lex_match_id (lexer, "DEPENDENT"))
            /*
-            Use the command "/MISSING=DEPENDENT" to
-            include user-missing values for the categorical variables, 
-            while excluding them for the dependent variables.
+             Use the command "/MISSING=DEPENDENT" to
+             include user-missing values for the categorical variables, 
+             while excluding them for the dependent variables.
 
-            Cases are dropped only when user-missing values
-            appear in dependent  variables.  User-missing
-            values for categorical variables are treated according to
-            their face value.
+             Cases are dropped only when user-missing values
+             appear in dependent  variables.  User-missing
+             values for categorical variables are treated according to
+             their face value.
 
-            Cases are ALWAYS dropped when System Missing values appear 
-            in the categorical variables.
+             Cases are ALWAYS dropped when System Missing values appear 
+             in the categorical variables.
            */
            {
              means.dep_exclude = MV_ANY;
@@ -698,7 +699,7 @@ cmd_means (struct lexer *lexer, struct dataset *ds)
                  int x;
                  means.cells =
                    pool_realloc (means.pool, means.cells,
-                             (means.n_cells += n_C) * sizeof (*means.cells));
+                                 (means.n_cells += n_C) * sizeof (*means.cells));
 
                  for (x = 0; x < n_C; ++x)
                    means.cells[means.n_cells - (n_C - 1 - x) - 1] = x;
@@ -711,7 +712,7 @@ cmd_means (struct lexer *lexer, struct dataset *ds)
                {
                  means.cells =
                    pool_realloc (means.pool, means.cells,
-                             (means.n_cells += 3) * sizeof (*means.cells));
+                                 (means.n_cells += 3) * sizeof (*means.cells));
                  
                  means.cells[means.n_cells - 2 - 1] = MEANS_MEAN;
                  means.cells[means.n_cells - 1 - 1] = MEANS_N;
@@ -725,7 +726,7 @@ cmd_means (struct lexer *lexer, struct dataset *ds)
                        {
                          means.cells =
                            pool_realloc (means.pool, means.cells,
-                                     ++means.n_cells * sizeof (*means.cells));
+                                         ++means.n_cells * sizeof (*means.cells));
                        
                          means.cells[means.n_cells - 1] = k;
                          break;
@@ -749,28 +750,28 @@ cmd_means (struct lexer *lexer, struct dataset *ds)
 
 
   for (t = 0; t < means.n_tables; ++t)
-  {
-    struct mtable *table = &means.table[t];
+    {
+      struct mtable *table = &means.table[t];
 
-    table->interactions =
-      pool_calloc (means.pool, table->n_layers, sizeof (*table->interactions));
+      table->interactions =
+       pool_calloc (means.pool, table->n_layers, sizeof (*table->interactions));
 
-    table->summary =
-      pool_calloc (means.pool, table->n_dep_vars * table->n_layers, sizeof (*table->summary));
+      table->summary =
+       pool_calloc (means.pool, table->n_dep_vars * table->n_layers, sizeof (*table->summary));
 
-    for (l = 0; l < table->n_layers; ++l)
-      {
-        int v;
-        const struct layer *lyr = &table->layers[l];
-       const int n_vars = lyr->n_factor_vars;
-        table->interactions[l] = interaction_create (NULL);
-        for (v = 0 ; v < n_vars ; ++v)
-          {
-            interaction_add_variable (table->interactions[l],
-                                      lyr->factor_vars[v]);
-          }
-      }
-  }
+      for (l = 0; l < table->n_layers; ++l)
+       {
+         int v;
+         const struct layer *lyr = &table->layers[l];
+         const int n_vars = lyr->n_factor_vars;
+         table->interactions[l] = interaction_create (NULL);
+         for (v = 0; v < n_vars ; ++v)
+           {
+             interaction_add_variable (table->interactions[l],
+                                       lyr->factor_vars[v]);
+           }
+       }
+    }
 
   {
     struct casegrouper *grouper;
@@ -786,12 +787,33 @@ cmd_means (struct lexer *lexer, struct dataset *ds)
     ok = proc_commit (ds) && ok;
   }
 
+  for (t = 0; t < means.n_tables; ++t)
+    {
+      int l;
+      struct mtable *table = &means.table[t];
+      if (table->interactions)
+       for (l = 0; l < table->n_layers; ++l)
+         {
+           interaction_destroy (table->interactions[l]);
+         }
+    }
 
   pool_destroy (means.pool);
   return CMD_SUCCESS;
 
-error:
+ error:
 
+  for (t = 0; t < means.n_tables; ++t)
+    {
+      int l;
+      struct mtable *table = &means.table[t];
+      if (table->interactions)
+       for (l = 0; l < table->n_layers; ++l)
+         {
+           interaction_destroy (table->interactions[l]);
+         }  
+    }
+  
   pool_destroy (means.pool);
   return CMD_FAILURE;
 }
@@ -827,6 +849,22 @@ struct per_cat_data
   bool warn;
 };
 
+
+static void 
+destroy_n (const void *aux1 UNUSED, void *aux2, void *user_data)
+{
+  struct mtable *table = aux2;
+  int v;
+  struct per_cat_data *per_cat_data = user_data;
+  struct per_var_data *pvd = per_cat_data->pvd;
+
+  for (v = 0; v < table->n_dep_vars; ++v)
+    {
+      struct per_var_data *pp = &pvd[v];
+      moments1_destroy (pp->mom);
+    }
+}
+
 static void *
 create_n (const void *aux1, void *aux2)
 {
@@ -938,19 +976,20 @@ run_means (struct means *cmd, struct casereader *input,
   struct payload payload;
   payload.create = create_n;
   payload.update = update_n;
-  payload.destroy = calculate_n;
+  payload.calculate = calculate_n;
+  payload.destroy = destroy_n;
   
   for (t = 0; t < cmd->n_tables; ++t)
-  {
-    struct mtable *table = &cmd->table[t];
-    table->cats
-      = categoricals_create (table->interactions,
-                            table->n_layers, wv, cmd->exclude);
+    {
+      struct mtable *table = &cmd->table[t];
+      table->cats
+       = categoricals_create (table->interactions,
+                              table->n_layers, wv, cmd->dep_exclude, cmd->exclude);
 
-    categoricals_set_payload (table->cats, &payload, cmd, table);
-  }
+      categoricals_set_payload (table->cats, &payload, cmd, table);
+    }
 
-  for (reader = casereader_clone (input);
+  for (reader = input;
        (c = casereader_read (reader)) != NULL; case_unref (c))
     {
       for (t = 0; t < cmd->n_tables; ++t)
@@ -1075,7 +1114,7 @@ output_case_processing_summary (const struct mtable *table)
            table->summary[row].non_missing;
 
          tab_double (t, 1, row + heading_rows,
-                     0, table->summary[row].non_missing, &F_8_0);
+                     0, table->summary[row].non_missing, NULL, RC_INTEGER);
 
          tab_text_format (t, 2, row + heading_rows,
                           0, _("%g%%"), 
@@ -1083,7 +1122,7 @@ output_case_processing_summary (const struct mtable *table)
 
 
          tab_double (t, 3, row + heading_rows,
-                     0, table->summary[row].missing, &F_8_0);
+                     0, table->summary[row].missing, NULL, RC_INTEGER);
 
 
          tab_text_format (t, 4, row + heading_rows,
@@ -1093,7 +1132,7 @@ output_case_processing_summary (const struct mtable *table)
 
          tab_double (t, 5, row + heading_rows,
                      0, table->summary[row].missing + 
-                     table->summary[row].non_missing, &F_8_0);
+                     table->summary[row].non_missing, NULL, RC_INTEGER);
 
          tab_text_format (t, 6, row + heading_rows,
                           0, _("%g%%"), 
@@ -1162,7 +1201,7 @@ output_report (const struct means *cmd,  int iact_idx,
          tab_text (t, 0,
                    heading_rows + dv * n_cats,
                    TAB_RIGHT | TAT_TITLE,
-                   var_get_name (table->dep_vars[dv])
+                   var_to_string (table->dep_vars[dv])
                    );
 
          if ( dv > 0)
@@ -1202,7 +1241,7 @@ output_report (const struct means *cmd,  int iact_idx,
 
              tab_double (t, heading_columns + i,
                          heading_rows + grp + dv * n_cats,
-                         0, result, 0);
+                         0, result, NULL, RC_OTHER);
            }
        }
     }