Пирамидальная сортировка работает очень медленно

#include <string>
#include <sstream>
#include <iostream>
#include <assert.h>
#include <exception>
#include <stdexcept>
#include <algorithm>
#include <time.h>
#include <iomanip>  
#include <math.h>
#include <stdio.h>
 
const int N = 25000000;
 
double a[N];
 
 
void fill(double * arr){
    for (int i = 0; i < N; ++i)
    {
        arr[i] = sin((double)i);
    }
}
 
 
 
int getIndexRightChild(int index) {
    return 2 * index + 1;
}
 
int getIndexLeftChild(int index) {
    return 2 * index + 2;
}
void swapTwoElements(double *arr, int first, int second)
{
    double temp = arr[first];
    arr[first] = arr[second];
    arr[second] = temp;
}
/*
  O(lg(n))
*/
void maxHeapify(double * arr, int index, int heapSize) {
    int indexLeft = getIndexLeftChild(index);  // O(1)
    int indexRight = getIndexRightChild(index); // O(1)
    int largest; // O (1)
    if (heapSize > indexLeft &&  arr[indexLeft] > arr[index]) // O(1)
        largest = indexLeft; // O (1)
    else //O(1)
        largest = index; // O(1)
    if (heapSize > indexRight &&  arr[indexRight] > arr[largest]) // O(1)
        largest = indexRight; // O(1)
    if (largest != index) //O1)
    { 
        swapTwoElements(arr, index, largest); // O(1)
        maxHeapify(arr, largest, heapSize); 
    }
}
 
/*
 
 O(n)
*/
void buildMaxHeap(double * arr, int heapSize){
    int side = floor(heapSize / 2);   // O(1)
    for (int i = side; i > -1; --i)
    {
        maxHeapify(arr, i, heapSize); // O(lg(n))
    }
}
 
 
/*
  O (n * lg(n))
*/
void heapSort(double *arr, int Num){
    int heapSize = Num;  // O(1)
    int size = Num;  /// O(1)
    buildMaxHeap(arr, heapSize);  //O(n)
 
    for (int i = (size - 1); i >= 1; --i) { // O(n)
        swapTwoElements(arr, 0, i);
        heapSize--;
        maxHeapify(arr, 0, heapSize); // O(lg(n))
    }
     
 
}
 
 
 
int _tmain(int argc, _TCHAR* argv[])
{
    fill(a);
    float fTimeStart = clock() / (float)CLOCKS_PER_SEC;
    //std::sort(a, a + N);
    heapSort(a, N);
    //mergeSort(a, 0, N);
    float fTimeStop = clock() / (float)CLOCKS_PER_SEC;
    printf("time == %f secund\n", fTimeStop - fTimeStart);
}

//#include <string>
//#include <sstream>
//#include <iostream>
//#include <assert.h>
//#include <exception>
//#include <stdexcept>
//#include <algorithm>
//#include <iomanip>  
#include <time.h>
#include <math.h>
#include <stdio.h>