13 #define SIZE_THRESHOLD 32
15 static int icmp(
const void *aa,
const void *bb)
23 static int fcmp(
const void *aa,
const void *bb)
35 static int dcmp(
const void *aa,
const void *bb)
47 int f_median(
int argc,
const int *argt,
void **args)
51 bool use_heap =
false;
56 for (i = 1; i <= argc; i++)
57 if (argt[i] != argt[0])
63 CELL *a = stack_array;
72 CELL *a1 = &a[(argc - 1) / 2];
73 CELL *a2 = &a[argc / 2];
78 for (j = 0; j < argc && !nv; j++) {
88 qsort(a, argc,
sizeof(
CELL), icmp);
89 res[i] = (*a1 + *a2) / 2;
100 FCELL *a = stack_array;
107 FCELL *res = args[0];
109 FCELL *a1 = &a[(argc - 1) / 2];
110 FCELL *a2 = &a[argc / 2];
112 for (i = 0; i <
columns; i++) {
115 for (j = 0; j < argc && !nv; j++) {
125 qsort(a, argc,
sizeof(
FCELL), fcmp);
126 res[i] = (*a1 + *a2) / 2;
137 DCELL *a = stack_array;
144 DCELL *res = args[0];
146 DCELL *a1 = &a[(argc - 1) / 2];
147 DCELL *a2 = &a[argc / 2];
149 for (i = 0; i <
columns; i++) {
152 for (j = 0; j < argc && !nv; j++) {
162 qsort(a, argc,
sizeof(
DCELL), dcmp);
163 res[i] = (*a1 + *a2) / 2;
void G_free(void *)
Free allocated memory.
size_t Rast_cell_size(RASTER_MAP_TYPE)
Returns size of a raster cell in bytes.