|  | 
|  | 1 | +#include <stdio.h> | 
|  | 2 | +#include <stdint.h> | 
|  | 3 | +#include <stdlib.h> | 
|  | 4 | +#include <stdbool.h> | 
|  | 5 | + | 
|  | 6 | +// declare function here, linker will find this when linked to | 
|  | 7 | +// libx86simdsortcpp.so | 
|  | 8 | +void keyvalue_qsort_float_sizet(float*, size_t*, size_t); | 
|  | 9 | +void keyvalue_qsort_float_uint32(float*, uint32_t*, uint32_t); | 
|  | 10 | +void keyvalue_qsort_sizet_sizet(size_t*, size_t*, size_t); | 
|  | 11 | +void keyvalue_qsort_sizet_uint32(size_t*, uint32_t*, uint32_t); | 
|  | 12 | +void keyvalue_qsort_uint32_sizet(uint32_t*, size_t*, size_t); | 
|  | 13 | +void keyvalue_qsort_uint32_uint32(uint32_t*, uint32_t*, uint32_t); | 
|  | 14 | +void keyvalue_qsort_int32_sizet(int32_t*, size_t*, size_t); | 
|  | 15 | +void keyvalue_qsort_int32_uint32(int32_t*, uint32_t*, uint32_t); | 
|  | 16 | + | 
|  | 17 | +// struct definition, we will sort an array of these: | 
|  | 18 | +struct Point { | 
|  | 19 | +    int x; | 
|  | 20 | +    int y; | 
|  | 21 | +    float distance; | 
|  | 22 | +    size_t metric; | 
|  | 23 | +}; | 
|  | 24 | + | 
|  | 25 | +#define SWAP(a, b, type) {type temp = a; a = b; b = temp;} | 
|  | 26 | + | 
|  | 27 | +// Function to sort an array of objects: | 
|  | 28 | +void object_qsort(struct Point* arr, size_t size) | 
|  | 29 | +{ | 
|  | 30 | +    /* (1) Create and initialize arrays of key and value  */ | 
|  | 31 | +    size_t* key = malloc(size * sizeof(size_t)); | 
|  | 32 | +    size_t* arg = malloc(size * sizeof(size_t)); | 
|  | 33 | +    bool* done = malloc(size * sizeof(bool)); | 
|  | 34 | +    for (size_t ii = 0; ii < size; ++ii) { | 
|  | 35 | +        key[ii] = arr[ii].metric; | 
|  | 36 | +        arg[ii] = ii; | 
|  | 37 | +        done[ii] = false; | 
|  | 38 | +    } | 
|  | 39 | + | 
|  | 40 | +    /* (2) IndexSort using the keyvalue_qsort */ | 
|  | 41 | +    keyvalue_qsort_sizet_sizet(key, arg, size); | 
|  | 42 | + | 
|  | 43 | +    /* (3) Permute obj array in-place */ | 
|  | 44 | +    for (size_t ii = 0; ii < size; ++ii) { | 
|  | 45 | +        if (done[ii]) { continue; } | 
|  | 46 | +        done[ii] = true; | 
|  | 47 | +        size_t prev_j = ii; | 
|  | 48 | +        size_t jj = arg[ii]; | 
|  | 49 | +        while (ii != jj) { | 
|  | 50 | +            SWAP(arr[prev_j], arr[jj], struct Point); | 
|  | 51 | +            done[jj] = true; | 
|  | 52 | +            prev_j = jj; | 
|  | 53 | +            jj = arg[jj]; | 
|  | 54 | +        } | 
|  | 55 | +    } | 
|  | 56 | +    free(key); | 
|  | 57 | +    free(arg); | 
|  | 58 | +    free(done); | 
|  | 59 | +} | 
|  | 60 | + | 
|  | 61 | +int main() { | 
|  | 62 | +    const size_t size = 10; | 
|  | 63 | +    struct Point arr[size]; | 
|  | 64 | + | 
|  | 65 | +    // Initialize: | 
|  | 66 | +    for (size_t ii = 0; ii < size; ++ii) { | 
|  | 67 | +        arr[ii].distance = (float) rand() / RAND_MAX; | 
|  | 68 | +        arr[ii].metric =  rand() % 100; | 
|  | 69 | +    } | 
|  | 70 | + | 
|  | 71 | +    // sort: | 
|  | 72 | +    object_qsort(arr, size); | 
|  | 73 | + | 
|  | 74 | +    // check if it is sorted: | 
|  | 75 | +    printf("arr = "); | 
|  | 76 | +    for (size_t ii = 0; ii < size; ++ii) { | 
|  | 77 | +        printf("%ld, ", arr[ii].metric); | 
|  | 78 | +    } | 
|  | 79 | +    printf("\n"); | 
|  | 80 | +    return 0; | 
|  | 81 | +} | 
0 commit comments