X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Fmath%2Fsort.c;h=8719d877e096fcc9fc0d0f911d462d7bd0566725;hb=5c3291dc396b795696e94f47780308fd7ace6fc4;hp=aa2afe80951b7a31608f6e4d1bf4d03c8807de0a;hpb=43b1296aafe7582e7dbe6c2b6a8b478d7d9b0fcf;p=pspp-builds.git diff --git a/src/math/sort.c b/src/math/sort.c index aa2afe80..8719d877 100644 --- a/src/math/sort.c +++ b/src/math/sort.c @@ -1,5 +1,5 @@ /* PSPP - a program for statistical analysis. - Copyright (C) 1997-9, 2000, 2006 Free Software Foundation, Inc. + Copyright (C) 1997-9, 2000, 2006, 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 @@ -20,17 +20,18 @@ #include -#include #include #include #include #include #include -#include +#include #include #include #include +#include "xalloc.h" + #include "gettext.h" #define _(msgid) gettext (msgid) @@ -40,42 +41,44 @@ int max_buffers = INT_MAX; struct sort_writer { - struct case_ordering *ordering; + struct caseproto *proto; + struct subcase ordering; struct merge *merge; struct pqueue *pqueue; struct casewriter *run; casenumber run_id; - struct ccase run_end; + struct ccase *run_end; }; static struct casewriter_class sort_casewriter_class; -static struct pqueue *pqueue_create (const struct case_ordering *); +static struct pqueue *pqueue_create (const struct subcase *, + const struct caseproto *); static void pqueue_destroy (struct pqueue *); static bool pqueue_is_full (const struct pqueue *); static bool pqueue_is_empty (const struct pqueue *); static void pqueue_push (struct pqueue *, struct ccase *, casenumber); -static void pqueue_pop (struct pqueue *, struct ccase *, casenumber *); +static struct ccase *pqueue_pop (struct pqueue *, casenumber *); static void output_record (struct sort_writer *); struct casewriter * -sort_create_writer (struct case_ordering *ordering) +sort_create_writer (const struct subcase *ordering, + const struct caseproto *proto) { struct sort_writer *sort; sort = xmalloc (sizeof *sort); - sort->ordering = case_ordering_clone (ordering); - sort->merge = merge_create (ordering); - sort->pqueue = pqueue_create (ordering); + sort->proto = caseproto_ref (proto); + subcase_clone (&sort->ordering, ordering); + sort->merge = merge_create (ordering, proto); + sort->pqueue = pqueue_create (ordering, proto); sort->run = NULL; sort->run_id = 0; - case_nullify (&sort->run_end); - - case_ordering_destroy (ordering); + sort->run_end = NULL; - return casewriter_create (&sort_casewriter_class, sort); + return casewriter_create (proto, &sort_casewriter_class, sort); } static void @@ -88,9 +91,9 @@ sort_casewriter_write (struct casewriter *writer UNUSED, void *sort_, if (pqueue_is_full (sort->pqueue)) output_record (sort); - next_run = (case_is_null (&sort->run_end) - || case_ordering_compare_cases (c, &sort->run_end, - sort->ordering) < 0); + next_run = (sort->run_end == NULL + || subcase_compare_3way (&sort->ordering, c, + &sort->ordering, sort->run_end) < 0); pqueue_push (sort->pqueue, c, sort->run_id + (next_run ? 1 : 0)); } @@ -99,11 +102,12 @@ sort_casewriter_destroy (struct casewriter *writer UNUSED, void *sort_) { struct sort_writer *sort = sort_; - case_ordering_destroy (sort->ordering); + subcase_destroy (&sort->ordering); merge_destroy (sort->merge); pqueue_destroy (sort->pqueue); casewriter_destroy (sort->run); - case_destroy (&sort->run_end); + case_unref (sort->run_end); + caseproto_unref (sort->proto); free (sort); } @@ -116,8 +120,7 @@ sort_casewriter_convert_to_reader (struct casewriter *writer, void *sort_) if (sort->run == NULL && sort->run_id == 0) { /* In-core sort. */ - sort->run = mem_writer_create (case_ordering_get_value_cnt ( - sort->ordering)); + sort->run = mem_writer_create (sort->proto); sort->run_id = 1; } while (!pqueue_is_empty (sort->pqueue)) @@ -134,12 +137,12 @@ sort_casewriter_convert_to_reader (struct casewriter *writer, void *sort_) static void output_record (struct sort_writer *sort) { - struct ccase min_case; + struct ccase *min_case; casenumber min_run_id; - pqueue_pop (sort->pqueue, &min_case, &min_run_id); + min_case = pqueue_pop (sort->pqueue, &min_run_id); #if 0 - printf ("\toutput: %f to run %d\n", case_num_idx (&min_case, 0), min_run_id); + printf ("\toutput: %f to run %d\n", case_num_idx (min_case, 0), min_run_id); #endif if (sort->run_id != min_run_id && sort->run != NULL) @@ -149,15 +152,13 @@ output_record (struct sort_writer *sort) } if (sort->run == NULL) { - sort->run = tmpfile_writer_create (case_ordering_get_value_cnt ( - sort->ordering)); + sort->run = tmpfile_writer_create (sort->proto); sort->run_id = min_run_id; } - case_destroy (&sort->run_end); - case_clone (&sort->run_end, &min_case); - - casewriter_write (sort->run, &min_case); + case_unref (sort->run_end); + sort->run_end = case_ref (min_case); + casewriter_write (sort->run, min_case); } static struct casewriter_class sort_casewriter_class = @@ -168,20 +169,34 @@ static struct casewriter_class sort_casewriter_class = }; /* Reads all the cases from INPUT. Sorts the cases according to - ORDERING. Returns the sorted cases in a new casereader, or a - null pointer if an I/O error occurs. Both INPUT and ORDERING - are destroyed upon return, regardless of success. */ + ORDERING. Returns the sorted cases in a new casereader. */ struct casereader * -sort_execute (struct casereader *input, struct case_ordering *ordering) +sort_execute (struct casereader *input, const struct subcase *ordering) { - struct casewriter *output = sort_create_writer (ordering); + struct casewriter *output = + sort_create_writer (ordering, casereader_get_proto (input)); casereader_transfer (input, output); return casewriter_make_reader (output); } + +/* Reads all the cases from INPUT. Sorts the cases in ascending + order according to VARIABLE. Returns the sorted cases in a + new casereader. */ +struct casereader * +sort_execute_1var (struct casereader *input, const struct variable *var) +{ + struct subcase sc; + struct casereader *reader; + + subcase_init_var (&sc, var, SC_ASCEND); + reader = sort_execute (input, &sc); + subcase_destroy (&sc); + return reader; +} struct pqueue { - struct case_ordering *ordering; + struct subcase ordering; struct pqueue_record *records; size_t record_cnt; size_t record_cap; @@ -191,7 +206,7 @@ struct pqueue struct pqueue_record { casenumber id; - struct ccase c; + struct ccase *c; casenumber idx; }; @@ -199,14 +214,13 @@ static int compare_pqueue_records_minheap (const void *a, const void *b, const void *pq_); static struct pqueue * -pqueue_create (const struct case_ordering *ordering) +pqueue_create (const struct subcase *ordering, const struct caseproto *proto) { struct pqueue *pq; pq = xmalloc (sizeof *pq); - pq->ordering = case_ordering_clone (ordering); - pq->record_cap - = get_workspace_cases (case_ordering_get_value_cnt (ordering)); + subcase_clone (&pq->ordering, ordering); + pq->record_cap = settings_get_workspace_cases (proto); if (pq->record_cap > max_buffers) pq->record_cap = max_buffers; else if (pq->record_cap < min_buffers) @@ -225,12 +239,11 @@ pqueue_destroy (struct pqueue *pq) { while (!pqueue_is_empty (pq)) { - struct ccase c; casenumber id; - pqueue_pop (pq, &c, &id); - case_destroy (&c); + struct ccase *c = pqueue_pop (pq, &id); + case_unref (c); } - case_ordering_destroy (pq->ordering); + subcase_destroy (&pq->ordering); free (pq->records); free (pq); } @@ -257,15 +270,15 @@ pqueue_push (struct pqueue *pq, struct ccase *c, casenumber id) r = &pq->records[pq->record_cnt++]; r->id = id; - case_move (&r->c, c); + r->c = c; r->idx = pq->idx++; push_heap (pq->records, pq->record_cnt, sizeof *pq->records, compare_pqueue_records_minheap, pq); } -static void -pqueue_pop (struct pqueue *pq, struct ccase *c, casenumber *id) +static struct ccase * +pqueue_pop (struct pqueue *pq, casenumber *id) { struct pqueue_record *r; @@ -276,7 +289,7 @@ pqueue_pop (struct pqueue *pq, struct ccase *c, casenumber *id) r = &pq->records[pq->record_cnt]; *id = r->id; - case_move (c, &r->c); + return r->c; } /* Compares record-run tuples A and B on id, then on case data, @@ -290,7 +303,7 @@ compare_pqueue_records_minheap (const void *a_, const void *b_, const struct pqueue *pq = pq_; int result = a->id < b->id ? -1 : a->id > b->id; if (result == 0) - result = case_ordering_compare_cases (&a->c, &b->c, pq->ordering); + result = subcase_compare_3way (&pq->ordering, a->c, &pq->ordering, b->c); if (result == 0) result = a->idx < b->idx ? -1 : a->idx > b->idx; return -result;