X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Fmath%2Fpercentiles.c;h=9055279b5afa1db5f54db6a0cb062eaf734d78bf;hb=refs%2Fheads%2Fpivot-table2;hp=bf99de163ffbaafe2af8711685e9a640a0256784;hpb=9a331fe64eb814ae5c1322e21717a04fb254bf65;p=pspp diff --git a/src/math/percentiles.c b/src/math/percentiles.c index bf99de163f..9055279b5a 100644 --- a/src/math/percentiles.c +++ b/src/math/percentiles.c @@ -1,5 +1,5 @@ /* PSPP - a program for statistical analysis. - Copyright (C) 2008 Free Software Foundation, Inc. + Copyright (C) 2008, 2009, 2011 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 @@ -15,21 +15,22 @@ along with this program. If not, see . */ #include -#include "percentiles.h" -#include +#include "math/percentiles.h" + +#include "data/casereader.h" +#include "data/val-type.h" +#include "data/variable.h" +#include "libpspp/assertion.h" +#include "libpspp/cast.h" +#include "math/order-stats.h" + +#include "gl/xalloc.h" #include "gettext.h" #define _(msgid) gettext (msgid) #define N_(msgid) msgid -#include -#include -#include -#include -#include - - const char *const ptile_alg_desc[] = { "", N_("HAverage"), @@ -44,11 +45,9 @@ const char *const ptile_alg_desc[] = { double percentile_calculate (const struct percentile *ptl, enum pc_alg alg) { - struct percentile *mutable = (struct percentile *) ptl; + struct percentile *mutable = CONST_CAST (struct percentile *, ptl); const struct order_stats *os = &ptl->parent; - assert (os->cc == ptl->w); - if ( ptl->g1 == SYSMIS) mutable->g1 = (os->k[0].tc - os->k[0].cc) / os->k[0].c_p1; @@ -154,18 +153,19 @@ percentile_calculate (const struct percentile *ptl, enum pc_alg alg) static void destroy (struct statistic *stat) { - struct order_stats *os = (struct order_stats *) stat; + struct percentile *ptl = UP_CAST (stat, struct percentile, parent.parent); + struct order_stats *os = &ptl->parent; free (os->k); - free (os); + free (ptl); } -struct order_stats * +struct percentile * percentile_create (double p, double W) { struct percentile *ptl = xzalloc (sizeof (*ptl)); - struct order_stats *os = (struct order_stats *) ptl; - struct statistic *stat = (struct statistic *) ptl; + struct order_stats *os = &ptl->parent; + struct statistic *stat = &os->parent; assert (p >= 0); assert (p <= 1.0); @@ -174,7 +174,7 @@ percentile_create (double p, double W) ptl->w = W; os->n_k = 2; - os->k = xcalloc (sizeof (*os->k), 2); + os->k = xcalloc (2, sizeof (*os->k)); os->k[0].tc = W * p; os->k[1].tc = (W + 1.0) * p; @@ -186,6 +186,6 @@ percentile_create (double p, double W) stat->destroy = destroy; - return os; + return ptl; }