X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Fmath%2Fcategoricals.c;h=25de2cc6ea18a0d5240edb39d239a8959264686a;hb=53d725a9a75cf1644a4cbc9a7d08855ac02b65e4;hp=576e19b1d322be30f2ee878aeda65e33083347d2;hpb=3025fc73778a733bfd2aeefd57da1b1c1955a675;p=pspp-builds.git diff --git a/src/math/categoricals.c b/src/math/categoricals.c index 576e19b1..25de2cc6 100644 --- a/src/math/categoricals.c +++ b/src/math/categoricals.c @@ -1,5 +1,5 @@ /* PSPP - a program for statistical analysis. - Copyright (C) 2009 Free Software Foundation, Inc. + Copyright (C) 2009, 2010 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 @@ -26,6 +26,7 @@ #include #include #include +#include #include @@ -41,7 +42,6 @@ struct value_node Can be used as an index into an array */ }; - struct var_params { /* A map indexed by a union values */ @@ -62,6 +62,18 @@ struct var_params }; +/* Comparison function to sort the reverse_value_map in ascending order */ +static int +compare_value_node (const void *vn1_, const void *vn2_, const void *aux) +{ + const struct value_node * const *vn1 = vn1_; + const struct value_node * const *vn2 = vn2_; + const struct var_params *vp = aux; + + return value_compare_3way (&(*vn1)->value, &(*vn2)->value, var_get_width (vp->var)); +} + + struct categoricals { /* The weight variable */ @@ -115,7 +127,7 @@ categoricals_destroy ( struct categoricals *cat) } -#if 1 +#if 0 void categoricals_dump (const struct categoricals *cat) { @@ -299,7 +311,7 @@ categoricals_total (const struct categoricals *cat) /* This function must be called *before* any call to categoricals_get_*_by subscript an *after* all calls to categoricals_update */ void -categoricals_done (struct categoricals *cat) +categoricals_done (const struct categoricals *cat_) { /* Implementation Note: Whilst this function is O(n) in cat->n_cats_total, in most uses it will be more efficient that using a tree based structure, since it @@ -307,6 +319,7 @@ categoricals_done (struct categoricals *cat) 1 call of O(n) + 10^9 calls of O(1) is better than 10^9 calls of O(log n). */ + struct categoricals *cat = CONST_CAST (struct categoricals *, cat_); int v; int idx = 0; cat->reverse_variable_map = pool_calloc (cat->pool, @@ -330,6 +343,10 @@ categoricals_done (struct categoricals *cat) vp->reverse_value_map[vn->subscript] = vn; } + /* For some purposes (eg CONTRASTS in ONEWAY) the values need to be sorted */ + sort (vp->reverse_value_map, vp->n_cats, sizeof (const struct value_node *), + compare_value_node, vp); + /* Populate the reverse variable map. */ for (i = 0; i < vp->n_cats; ++i) @@ -338,7 +355,6 @@ categoricals_done (struct categoricals *cat) assert (cat->n_vars <= cat->n_vp); - // categoricals_dump (cat); }