mirror of
https://github.com/VCMP-SqMod/SqMod.git
synced 2024-11-08 00:37:15 +01:00
49 lines
1.3 KiB
Plaintext
49 lines
1.3 KiB
Plaintext
|
/* Sort function. */
|
||
|
/* from https://github.com/bel2125/sort_r */
|
||
|
|
||
|
static void
|
||
|
mg_sort(void *data,
|
||
|
size_t elemcount,
|
||
|
size_t elemsize,
|
||
|
int (*compfunc)(const void *data1, const void *data2, void *userarg),
|
||
|
void *userarg)
|
||
|
{
|
||
|
/* We cannot use qsort_r here. For a detailed reason, see
|
||
|
* https://github.com/civetweb/civetweb/issues/1048#issuecomment-1047093014
|
||
|
* https://stackoverflow.com/questions/39560773/different-declarations-of-qsort-r-on-mac-and-linux
|
||
|
*/
|
||
|
|
||
|
/* We use ShellSort here with this gap sequence: https://oeis.org/A102549 */
|
||
|
int A102549[9] = {1, 4, 10, 23, 57, 132, 301, 701, 1750};
|
||
|
int Aidx, gap, i, j, k;
|
||
|
void *tmp = alloca(elemsize);
|
||
|
|
||
|
for (Aidx = 8; Aidx >= 0; Aidx--) {
|
||
|
gap = A102549[Aidx];
|
||
|
if (gap > ((int)elemcount / 2)) {
|
||
|
continue;
|
||
|
}
|
||
|
for (i = 0; i < gap; i++) {
|
||
|
for (j = i; j < (int)elemcount; j += gap) {
|
||
|
memcpy(tmp, (void *)((ptrdiff_t)data + elemsize * j), elemsize);
|
||
|
|
||
|
for (k = j; k >= gap; k -= gap) {
|
||
|
void *cmp =
|
||
|
(void *)((ptrdiff_t)data + elemsize * (k - gap));
|
||
|
int cmpres = compfunc(cmp, tmp, userarg);
|
||
|
if (cmpres > 0) {
|
||
|
memcpy((void *)((ptrdiff_t)data + elemsize * k),
|
||
|
cmp,
|
||
|
elemsize);
|
||
|
} else {
|
||
|
break;
|
||
|
}
|
||
|
}
|
||
|
memcpy((void *)((ptrdiff_t)data + elemsize * k), tmp, elemsize);
|
||
|
}
|
||
|
}
|
||
|
}
|
||
|
}
|
||
|
|
||
|
/* end if sort.inl */
|