14 #define SIZE_THRESHOLD 32
16 static int icmp(
const void *aa,
const void *bb)
24 static int fcmp(
const void *aa,
const void *bb)
36 static int dcmp(
const void *aa,
const void *bb)
48 int f_nmedian(
int argc,
const int *argt,
void **args)
52 bool use_heap =
false;
57 for (i = 1; i <= argc; i++)
58 if (argt[i] != argt[0])
64 CELL *a = stack_array;
79 for (j = 0; j < argc; j++) {
90 qsort(a, n,
sizeof(
CELL), icmp);
95 *resc = (*resc + a1) / 2;
108 FCELL *a = stack_array;
115 FCELL *res = args[0];
120 for (i = 0; i <
columns; i++) {
123 for (j = 0; j < argc; j++) {
134 qsort(a, n,
sizeof(
FCELL), fcmp);
139 *resc = (*resc + a1) / 2;
152 DCELL *a = stack_array;
159 DCELL *res = args[0];
164 for (i = 0; i <
columns; i++) {
167 for (j = 0; j < argc; j++) {
178 qsort(a, n,
sizeof(
DCELL), dcmp);
183 *resc = (*resc + a1) / 2;
void G_free(void *)
Free allocated memory.
size_t Rast_cell_size(RASTER_MAP_TYPE)
Returns size of a raster cell in bytes.