Форум программистов, компьютерный форум CyberForum.ru
Наши страницы

OpenMP Для сортровки слиянием - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Визначте номер мінімального додатного елемента масиву. Перший елемент масиву дорівнює 1000. Розмір масиву 6 http://www.cyberforum.ru/cpp-beginners/thread1632118.html
Розробіть програму, в якій визначається масив цілих чисел заданого розміру. Перший елемент масиву задається, всі інші елемен¬ти вводяться із клавіатури, потім виконуються задані обчислення,...
C++ Визначте мінімальний елемент масиву, що кратний 5. Перший елемент масиву дорівнює 500. Розмір масиву 8 Розробіть програму, в якій визначається масив цілих чисел заданого розміру. Перший елемент масиву задається, всі інші елемен¬ти вводяться із клавіатури, потім виконуються задані обчислення,... http://www.cyberforum.ru/cpp-beginners/thread1632116.html
Обчисліть суму цілих чисел, которие вводяться з клавіатури, доки сума не стане більше 100 C++
Обчисліть суму цілих чисел, що вводяться з клавіатури, доки сума не стане більше 100. Язык форума русский. Тема закрыта.
Вывод содержимого каталога C++
Здравствуйте. Возникла следующая проблема: Нужно вывести в консоль содержимое указанного каталога. Проблема в том, что у меня получилось вывести на экран только содержимое каталога из которого...
C++ Ссылка на неразрешенный внешний символ http://www.cyberforum.ru/cpp-beginners/thread1632089.html
Здравствуйте. Пробую сделать что-то вроде map. Ошибка при создании объекта. #pragma once #include <iostream> template <typename Type, class List> class Map { public: Map(); ~Map(); Type...
C++ Метод простой итерации Найти приближенное решение уравнения x = f(x), где f(x) = α*exp(x), при α = 0, 2 методом простой итерации: xn+1 = f(xn). В качестве начального приближения взять x0 = 0. Итерации продолжать до тех... подробнее

Показать сообщение отдельно
Мыша
4 / 4 / 0
Регистрация: 03.03.2011
Сообщений: 192

OpenMP Для сортровки слиянием - C++

05.01.2016, 16:22. Просмотров 183. Ответов 0
Метки (Все метки)

хочу расспарлелисть работу сортировки слиянием но почему-то время с openMP становится больше. Вот тут два разных алгоритма сортировки: один рекурсивный другой без рекурсии. Не могу понять почему!

1:
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
template <typename Item>
void Merge(Item Mas[], int left, int right, int medium)
{
    int j = left;
    int k = medium + 1;
    int count = right - left + 1;
 
    if (count <= 1) return;
 
    Item *TmpMas = new Item[count];
 
    for (int i = 0; i < count; ++i) {
        if (j <= medium && k <= right) {
            if (Mas[j] < Mas[k])
                TmpMas[i] = Mas[j++];
            else
                TmpMas[i] = Mas[k++];
        } else {
            if (j <= medium)
                TmpMas[i] = Mas[j++];
            else
                TmpMas[i] = Mas[k++];
        }
    }
 
    j = 0;
    for (int i = left; i <= right; ++i) {
        Mas[i] = TmpMas[j++];
    }
    delete[] TmpMas;
}
 
template <typename Item>
void MergeSort(Item a[], int l, int r)
{
    int m;
 
    // Условие выхода из рекурсии
    if(l >= r) return;
 
    m = (l + r) / 2;
 
    // Рекурсивная сортировка полученных массивов
    #pragma omp parallel
    {
    #pragma omp sections
    {
        #pragma omp section
        {
            MergeSort(a, l, m);
        }
        #pragma omp section
        {
            MergeSort(a, m + 1, r);
        }
    }
    }
    Merge(a, l, r, m);
}
int main()
{
 
    double start_time, end_time, start_time1, end_time1; 
    
    int a[1000];
 
    for (int i=0; i<1000; i++)
        a[i]=rand()%100+1;
    cout<<endl;
 
    start_time = omp_get_wtime();
 
    MergeSort(a,0,1000);
 
    end_time = omp_get_wtime();
    cout<<"time: "<<end_time-start_time<<endl;
 
 
    return 0;
}

2:

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
template <class T>
T* merge(T *m1, T* m2, int l1, int l2)
{
    T* ret = new T[l1+l2];
    int n = 0;
    // Сливаем массивы, пока один не закончится
    while (l1 && l2) 
    {   
        if (*m1 < *m2) 
        {
            ret[n] = *m1;
            m1++;
            --l1;
        } else 
        {
            ret[n] = *m2;
            ++m2;
            --l2;
        }
        ++n;    
    }
    // Если закончился первый массив
    if (l1 == 0) 
    {
        #pragma omp for
        for (int i = 0; i < l2; ++i) 
        {
            ret[n++] = *m2++;
        }
    } else { // Если закончился второй массив
        #pragma omp for
        for (int i = 0; i < l1; ++i) 
        {
            ret[n++] = *m1++;
        }
    }
    return ret;
}
 
// Функция восходящего слияния
template <class T>
void mergeSort(T * mas, int len)
{
    int n = 1, l, ost;
    T * mas1;
    while (n < len) 
    {
        l = 0;
        while (l < len) 
        {
            if (l + n >= len) break;
            ost = (l + n * 2 > len) ? (len - (l + n)) : n;
            mas1 = merge(mas + l, mas + l + n, n, ost);
            #pragma omp for
            for (int i = 0; i < n + ost; ++i)
                mas[l+i] = mas1[i];
            delete [] mas1;
            l += n * 2;
        }
        n *= 2;
    }
}
 
 
int main()
{
 
    double start_time, end_time, start_time1, end_time1; 
    
    int a[1000];
 
    for (int i=0; i<1000; i++)
        a[i]=rand()%100+1;
    cout<<endl;
 
    start_time = omp_get_wtime();
 
    mergeSort(a,0,1000);
 
    end_time = omp_get_wtime();
    cout<<"time: "<<end_time-start_time<<endl;
 
 
    return 0;
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru