X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Fmath%2Fnp.c;h=e61bf58e117aaf915d596bc26b8875152f0f0deb;hb=edf539fc797fedcdf7a71894f707d3c9728335d6;hp=e189b47091faee8304ee1a557042cd0f2cff873c;hpb=015e221b0f8578afee769528572c76387f26c629;p=pspp-builds.git diff --git a/src/math/np.c b/src/math/np.c index e189b470..e61bf58e 100644 --- a/src/math/np.c +++ b/src/math/np.c @@ -1,5 +1,5 @@ /* PSPP - a program for statistical analysis. - Copyright (C) 2008 Free Software Foundation, Inc. + Copyright (C) 2008, 2009 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 @@ -16,14 +16,18 @@ #include #include "np.h" -#include -#include -#include -#include + #include -#include +#include +#include + #include #include +#include +#include +#include + +#include "xalloc.h" static void destroy (struct statistic *stat) @@ -37,7 +41,7 @@ static void acc (struct statistic *s, const struct ccase *cx UNUSED, double c, double cc, double y) { - struct ccase cp; + struct ccase *cp; struct np *np = (struct np *) s; double rank = np->prev_cc + (c + 1) / 2.0; @@ -56,13 +60,11 @@ acc (struct statistic *s, const struct ccase *cx UNUSED, maximize (&np->y_max, y); minimize (&np->y_min, y); - case_create (&cp, n_NP_IDX); - - case_data_rw_idx (&cp, NP_IDX_Y)->f = y; - case_data_rw_idx (&cp, NP_IDX_NS)->f = ns; - case_data_rw_idx (&cp, NP_IDX_DNS)->f = dns; - - casewriter_write (np->writer, &cp); + cp = case_create (casewriter_get_proto (np->writer)); + case_data_rw_idx (cp, NP_IDX_Y)->f = y; + case_data_rw_idx (cp, NP_IDX_NS)->f = ns; + case_data_rw_idx (cp, NP_IDX_DNS)->f = dns; + casewriter_write (np->writer, cp); np->prev_cc = cc; } @@ -74,6 +76,8 @@ np_create (const struct moments1 *m) struct np *np = xzalloc (sizeof (*np)); struct statistic *stat = (struct statistic *) np; struct order_stats *os = (struct order_stats *) np; + struct caseproto *proto; + int i; np->prev_cc = 0; @@ -84,7 +88,11 @@ np_create (const struct moments1 *m) np->y_min = np->ns_min = np->dns_min = DBL_MAX; np->y_max = np->ns_max = np->dns_max = -DBL_MAX; - np->writer = autopaging_writer_create (n_NP_IDX); + proto = caseproto_create (); + for (i = 0; i < n_NP_IDX; i++) + proto = caseproto_add_width (proto, 0); + np->writer = autopaging_writer_create (proto); + caseproto_unref (proto); os->k = 0; stat->destroy = destroy;