X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Fmath%2Fnp.c;h=ccaa51e22d24a0e96a91996f7f597df547d3c6ce;hb=81579d9e9f994fb2908f50af41c3eb033d216e58;hp=e61bf58e117aaf915d596bc26b8875152f0f0deb;hpb=bd17d2af982332ee1791998361b1ac6731fe14fa;p=pspp-builds.git diff --git a/src/math/np.c b/src/math/np.c index e61bf58e..ccaa51e2 100644 --- a/src/math/np.c +++ b/src/math/np.c @@ -1,5 +1,5 @@ /* PSPP - a program for statistical analysis. - Copyright (C) 2008, 2009 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,25 +15,27 @@ along with this program. If not, see . */ #include -#include "np.h" + +#include "math/np.h" #include #include #include -#include -#include -#include -#include -#include +#include "data/case.h" +#include "data/casewriter.h" +#include "libpspp/cast.h" +#include "libpspp/compiler.h" +#include "libpspp/misc.h" +#include "math/moments.h" -#include "xalloc.h" +#include "gl/xalloc.h" static void destroy (struct statistic *stat) { - struct order_stats *os = (struct order_stats *) stat; - free (os); + struct np *np = UP_CAST (stat, struct np, parent.parent); + free (np); } @@ -42,7 +44,7 @@ acc (struct statistic *s, const struct ccase *cx UNUSED, double c, double cc, double y) { struct ccase *cp; - struct np *np = (struct np *) s; + struct np *np = UP_CAST (s, struct np, parent.parent); double rank = np->prev_cc + (c + 1) / 2.0; double ns = gsl_cdf_ugaussian_Pinv (rank / ( np->n + 1 )); @@ -69,13 +71,13 @@ acc (struct statistic *s, const struct ccase *cx UNUSED, np->prev_cc = cc; } -struct order_stats * +struct np * np_create (const struct moments1 *m) { double variance; struct np *np = xzalloc (sizeof (*np)); - struct statistic *stat = (struct statistic *) np; - struct order_stats *os = (struct order_stats *) np; + struct order_stats *os = &np->parent; + struct statistic *stat = &os->parent; struct caseproto *proto; int i; @@ -98,5 +100,5 @@ np_create (const struct moments1 *m) stat->destroy = destroy; stat->accumulate = acc; - return os; + return np; }