-
Notifications
You must be signed in to change notification settings - Fork 0
/
OpenMP_Quick_Sort.c
107 lines (80 loc) · 1.64 KB
/
OpenMP_Quick_Sort.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
#include <time.h>
#include <stdio.h>
#include <stdlib.h>
#include <omp.h>
int j;
void swap(int* a, int* b)
{
int t = *a;
*a = *b;
*b = t;
}
int partition (int arr[], int low, int high)
{
int pivot = arr[high];
int i = (low - 1);
#pragma omp schedule(static) parallel num_threads(4)
{
#pragma omp for
for (j = low; j <= high- 1; j++)
{
if (arr[j] <= pivot)
{
i++;
swap(&arr[i], &arr[j]);
}
}
swap(&arr[i + 1], &arr[high]);
}
return (i + 1);
}
void quickSort(int arr[], int low, int high)
{
if (low < high)
{
int pi = partition(arr, low, high);
quickSort(arr, low, pi - 1);
quickSort(arr, pi + 1, high);
}
}
void printArray(int arr[], int size)
{
int i;
#pragma omp schedule(static) parallel num_threads(2)
{
#pragma omp for
for (i=0; i < size; i++)
printf("%d ", arr[i]);
}
printf("\n");
}
int main()
{
clock_t begin, end;
double total_t;
begin = clock();
int num=0;
FILE * fp;
fp = fopen("1-25000.txt","r");
fscanf(fp, "%d", &num);
int *arr = malloc(num * sizeof(int));
int n=num;
int i=0;
#pragma omp schedule(static) parallel num_threads(2)
{
#pragma omp for
for ( i=0; i < num; i++){
fscanf(fp, "%d", &arr[i]);
}
}
quickSort(arr, 0, n-1);
end= clock();
total_t = (double)(end - begin );
printArray(arr,num);
printf("Total numbers sorted: %d \n",num);
printf("Initial time taken by CPU: %d \n", begin );
printf("End time taken by CPU: %d \n", end );
printf("Total time taken by CPU: %f seconds\n", total_t );
free(arr);
return 0;
}