ucx/dlist.c

Fri, 12 Oct 2012 12:12:59 +0200

author
Mike Becker <universe@uap-core.de>
date
Fri, 12 Oct 2012 12:12:59 +0200
changeset 73
f15c7d6aebb9
parent 69
fb59270b1de3
child 87
bd444539cced
permissions
-rw-r--r--

lists now sort on heap to prevent stack overflows

universe@4 1 #include "dlist.h"
universe@4 2
universe@67 3 UcxDlist *restrict ucx_dlist_clone(UcxDlist *restrict l,
universe@67 4 copy_func fnc, void *data) {
universe@18 5 UcxDlist *ret = NULL;
universe@18 6 while (l != NULL) {
universe@18 7 if (fnc != NULL) {
universe@18 8 ret = ucx_dlist_append(ret, fnc(l->data, data));
universe@18 9 } else {
universe@18 10 ret = ucx_dlist_append(ret, l->data);
universe@18 11 }
universe@18 12 l = l->next;
universe@18 13 }
universe@18 14 return ret;
universe@18 15 }
universe@18 16
universe@67 17 int ucx_dlist_equals(const UcxDlist *l1, const UcxDlist *l2,
universe@67 18 cmp_func fnc, void* data) {
universe@18 19 if (l1 == l2) return 1;
universe@18 20
universe@18 21 while (l1 != NULL && l2 != NULL) {
universe@18 22 if (fnc == NULL) {
universe@18 23 if (l1->data != l2->data) return 0;
universe@18 24 } else {
universe@18 25 if (fnc(l1->data, l2->data, data) != 0) return 0;
universe@18 26 }
universe@18 27 l1 = l1->next;
universe@18 28 l2 = l2->next;
universe@18 29 }
universe@18 30
universe@18 31 return (l1 == NULL && l2 == NULL);
universe@18 32 }
universe@18 33
universe@8 34 void ucx_dlist_free(UcxDlist *l) {
universe@8 35 UcxDlist *e = l, *f;
universe@8 36 while (e != NULL) {
universe@8 37 f = e;
universe@8 38 e = e->next;
universe@8 39 free(f);
universe@8 40 }
universe@8 41 }
universe@8 42
universe@7 43 UcxDlist *ucx_dlist_append(UcxDlist *l, void *data) {
universe@8 44 UcxDlist *nl = (UcxDlist*) malloc(sizeof(UcxDlist));
universe@8 45 if (nl == NULL) return NULL;
universe@7 46
universe@8 47 nl->data = data;
universe@8 48 nl->next = NULL;
universe@8 49 if (l == NULL) {
universe@22 50 nl->prev = NULL;
universe@8 51 return nl;
universe@8 52 } else {
universe@8 53 UcxDlist *t = ucx_dlist_last(l);
universe@8 54 t->next = nl;
universe@8 55 nl->prev = t;
universe@8 56 return l;
universe@8 57 }
universe@7 58 }
universe@7 59
universe@7 60 UcxDlist *ucx_dlist_prepend(UcxDlist *l, void *data) {
universe@8 61 UcxDlist *nl = ucx_dlist_append(NULL, data);
universe@8 62 if (nl == NULL) return NULL;
universe@7 63
universe@8 64 if (l != NULL) {
universe@8 65 nl->next = l;
universe@8 66 l->prev = nl;
universe@8 67 }
universe@8 68 return nl;
universe@7 69 }
universe@7 70
universe@67 71 UcxDlist *ucx_dlist_concat(UcxDlist *restrict l1, UcxDlist *restrict l2) {
universe@8 72 if (l1 == NULL) {
universe@8 73 return l2;
universe@8 74 } else {
universe@8 75 UcxDlist *last = ucx_dlist_last(l1);
universe@8 76 last->next = l2;
universe@8 77 l2->prev = last;
universe@8 78 return l1;
universe@8 79 }
universe@7 80 }
universe@7 81
universe@67 82 UcxDlist *ucx_dlist_last(const UcxDlist *l) {
universe@7 83 if (l == NULL) return NULL;
universe@7 84
universe@67 85 const UcxDlist *e = l;
universe@7 86 while (e->next != NULL) {
universe@7 87 e = e->next;
universe@7 88 }
universe@67 89 return (UcxDlist*)e;
universe@7 90 }
universe@7 91
universe@67 92 UcxDlist *ucx_dlist_get(const UcxDlist *l, int index) {
universe@8 93 if (l == NULL) return NULL;
universe@8 94
universe@67 95 const UcxDlist *e = l;
universe@8 96 while (e->next != NULL && index > 0) {
universe@8 97 e = e->next;
universe@8 98 index--;
universe@8 99 }
universe@7 100
universe@67 101 return (UcxDlist*)(index == 0 ? e : NULL);
universe@7 102 }
universe@7 103
universe@67 104 size_t ucx_dlist_size(const UcxDlist *l) {
universe@7 105 if (l == NULL) return 0;
universe@7 106
universe@67 107 const UcxDlist *e = l;
universe@7 108 size_t s = 1;
universe@7 109 while (e->next != NULL) {
universe@7 110 e = e->next;
universe@7 111 s++;
universe@7 112 }
universe@7 113
universe@7 114 return s;
universe@7 115 }
universe@7 116
universe@37 117 UcxDlist *ucx_dlist_sort_merge(int length,
universe@67 118 UcxDlist* restrict ls, UcxDlist* restrict le, UcxDlist* restrict re,
universe@37 119 cmp_func fnc, void* data) {
universe@73 120
universe@73 121 UcxDlist** sorted = (UcxDlist**) malloc(sizeof(UcxDlist*)*length);
universe@37 122 UcxDlist *rc, *lc;
universe@35 123
universe@67 124 lc = ls; rc = le;
universe@37 125 int n = 0;
universe@67 126 while (lc && lc != le && rc != re) {
universe@37 127 if (fnc(lc->data, rc->data, data) <= 0) {
universe@37 128 sorted[n] = lc;
universe@37 129 lc = lc->next;
universe@37 130 } else {
universe@37 131 sorted[n] = rc;
universe@37 132 rc = rc->next;
universe@35 133 }
universe@37 134 n++;
universe@37 135 }
universe@67 136 while (lc && lc != le) {
universe@37 137 sorted[n] = lc;
universe@37 138 lc = lc->next;
universe@37 139 n++;
universe@37 140 }
universe@67 141 while (rc && rc != re) {
universe@37 142 sorted[n] = rc;
universe@37 143 rc = rc->next;
universe@37 144 n++;
universe@35 145 }
universe@35 146
universe@37 147 // Update pointer
universe@37 148 sorted[0]->prev = NULL;
universe@37 149 for (int i = 0 ; i < length-1 ; i++) {
universe@37 150 sorted[i]->next = sorted[i+1];
universe@37 151 sorted[i+1]->prev = sorted[i];
universe@37 152 }
universe@37 153 sorted[length-1]->next = NULL;
universe@35 154
universe@69 155 UcxDlist *ret = sorted[0];
universe@73 156 free(sorted);
universe@69 157 return ret;
universe@35 158 }
universe@35 159
universe@36 160 UcxDlist *ucx_dlist_sort(UcxDlist *l, cmp_func fnc, void *data) {
universe@35 161 if (l == NULL) {
universe@35 162 return NULL;
universe@35 163 }
universe@37 164
universe@37 165 UcxDlist *lc;
universe@37 166 int ln = 1;
universe@37 167
universe@67 168 UcxDlist *restrict ls = l, *restrict le, *restrict re;
universe@37 169 lc = ls;
universe@37 170 while (lc->next != NULL && fnc(lc->next->data, lc->data, data) > 0) {
universe@37 171 lc = lc->next;
universe@37 172 ln++;
universe@37 173 }
universe@37 174 le = lc->next;
universe@37 175
universe@67 176 if (le == NULL) {
universe@37 177 return l; // this list is already sorted :)
universe@37 178 } else {
universe@37 179 UcxDlist *rc;
universe@37 180 int rn = 1;
universe@67 181 rc = le;
universe@37 182 while (rc->next != NULL && fnc(rc->next->data, rc->data, data) > 0) {
universe@37 183 rc = rc->next;
universe@37 184 rn++;
universe@37 185 }
universe@37 186 re = rc->next;
universe@37 187
universe@37 188 // Something left? Sort it!
universe@37 189 UcxDlist *remainder = re;
universe@37 190 size_t remainder_length = ucx_dlist_size(remainder);
universe@37 191 if (remainder != NULL) {
universe@37 192 remainder = ucx_dlist_sort(remainder, fnc, data);
universe@37 193 }
universe@37 194
universe@37 195 // {ls,...,le->prev} and {rs,...,re->prev} are sorted - merge them
universe@37 196 UcxDlist *sorted = ucx_dlist_sort_merge(ln+rn,
universe@67 197 ls, le, re,
universe@37 198 fnc, data);
universe@37 199
universe@37 200 // merge sorted list with (also sorted) remainder
universe@37 201 l = ucx_dlist_sort_merge(ln+rn+remainder_length,
universe@67 202 sorted, remainder, NULL, fnc, data);
universe@37 203
universe@35 204 return l;
universe@35 205 }
universe@35 206 }
universe@35 207
universe@7 208 /* dlist specific functions */
universe@67 209 UcxDlist *ucx_dlist_first(const UcxDlist *l) {
universe@8 210 if (l == NULL) return NULL;
universe@7 211
universe@67 212 const UcxDlist *e = l;
universe@8 213 while (e->prev != NULL) {
universe@8 214 e = e->prev;
universe@8 215 }
universe@67 216 return (UcxDlist *)e;
olaf@13 217 }
universe@22 218
universe@22 219 UcxDlist *ucx_dlist_remove(UcxDlist *l, UcxDlist *e) {
universe@22 220 if (e->prev == NULL) {
olaf@31 221 if(e->next != NULL) {
olaf@31 222 e->next->prev = NULL;
olaf@31 223 l = e->next;
olaf@31 224 } else {
olaf@31 225 l = NULL;
olaf@31 226 }
olaf@31 227
universe@22 228 } else {
universe@22 229 e->prev->next = e->next;
universe@22 230 e->next->prev = e->prev;
universe@22 231 }
universe@22 232 free(e);
universe@22 233 return l;
universe@22 234 }

mercurial