bitonic_withoutparallel.cpp

#include<bits/stdc++.h>
#include <omp.h>
using namespace std;


void compAndSwap(int a[], int i, int j, int dir)
{
    if (dir==(a[i]>a[j]))
        swap(a[i],a[j]);
}

void bitonicMerge(int a[], int low, int cnt, int dir)
{
    if (cnt>1)
    {
        int k = cnt/2;
        for (int i=low; i<low+k; i++)
        compAndSwap(a, i, i+k, dir);
        bitonicMerge(a, low, k, dir);
        bitonicMerge(a, low+k, k, dir);
    }
}

void bitonicSort(int a[],int low, int cnt, int dir)
{
    if (cnt>1)
    {
        int k = cnt/2;

        // #pragma omp parallel
        // {
            bitonicSort(a, low, k, 1);
            bitonicSort(a, low+k, k, 0);
            bitonicMerge(a,low, cnt, dir);            
        // }
    }
}

void sort(int a[], int N, int up)
{
    bitonicSort(a,0, N, up);
}

int main()
{   
    int N;
    cin >> N;
    int a[N];
    for(int i=0;i<N;i++){
        cin >> a[i];
    }

    int up = 1;   
    sort(a, N, up);

    printf("Sorted array: \n");
    for (int i=0; i<N; i++)
        printf("%d ", a[i]);
    return 0;
}