Оптимальная сортировка без дублирования данных

#include <stdio.h>
#include <stdlib.h>
#include <stddef.h>
#include <time.h>
 
// Правильная сортировка qsort для положительных и отрицательных значений float
int compare (const void * a, const void * b)
{
  float fa = *(const float*) a;
  float fb = *(const float*) b;
  return (fa < fb) - (fa > fb);
}
 
int main(){
 
 // Динамический массив с реальными данными является координатами точек на
 // графике: номер ячейки в массиве - число по оси X; записанные в ячейку
 // данные - число по оси Y.
 // В действительности, реальный массив содержит около 200 миллионов элементов float.
 ptrdiff_t elements_quantity = 1000;
 
 // Выделяем память под динамический массив
 float *x = (float*) malloc(elements_quantity*sizeof(float));
 // Проверка результата выделения памяти
 if ( !x ) {
        printf("Error: can't allocate memory!\n");
        getchar();  //задержка консоли
 }
 
 // Для примера заполним массив случайными данными
 srand(time(0));
 for(ptrdiff_t i=0; i<elements_quantity; i++) {
        x[i]=( rand() % 101 - 50 ) / 10.0;
//        printf("%f\n",x[i]);
 }
 
 // Создадим динамический массив для записи результатов вычислений
 ptrdiff_t *extracts;
 // Размер динамического массива
 ptrdiff_t size_of_extracts = 0;
 // Резервирование первой ячейки под массив
 extracts = (ptrdiff_t*)malloc(sizeof(ptrdiff_t));
 if ( !extracts ) {
    printf("Error: can't allocate memory!\n");
    getchar();  //задержка консоли
 }
 
 // Выполнить некие вычисления над всеми элементами реального массива x[]
 for(ptrdiff_t j=0; j<elements_quantity; j++) {
     // Результат вычислений поместить в виде номеров ячеек массива x[] в другой
     // динамический массив extracts
     if( x[j] >= -2.01 && x[j] <= 2.01 ){
 
         extracts[size_of_extracts] = j;
 
         // Переопределение памяти бОльшего размера
         size_of_extracts++;
         extracts = (ptrdiff_t*) realloc(extracts, sizeof(ptrdiff_t)*size_of_extracts+1);
         if ( !extracts ) {
             printf("Error: can't reallocate memory!\n");
             getchar();  //задержка консоли
         }
     }
 }
 
 // Вывод содержимого массива x[] в соответствие с результатами вычислений, помещённых в extracts[]
 // При работе с реальными данными размер массива extracts[] составляет около 30% от массива x[]
 for(ptrdiff_t k=0; k<size_of_extracts; k++) {
    printf("Номер ячейки в массиве и значение x[%td]=%.2f\n", extracts[k], x[extracts[k]]);
 }
 
 // Тут блок, сортирующий содержимое массива extracts в порядке убывания значений из массива x[]
 // Сортировка через qsort
 //
 // qsort(a, size_of_extracts, sizeof(extracts[0]), compare);
 
 // Вывод после сортировки
 for(ptrdiff_t k=0; k<size_of_extracts; k++) {
    printf("Номер ячейки в массиве и значение x[%td]=%.2f\n", extracts[k], x[extracts[k]]);
 }
 
 // Освобождение выделенной памяти
 free(x);
 free(extracts);
 
 return 0;
 
}

#pragma omp parallel for
        for(size_t i=max_elements_quantity-1; i >= 0; --i)