Форум программистов, компьютерный форум, киберфорум
OpenMP
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.60/5: Рейтинг темы: голосов - 5, средняя оценка - 4.60
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105

Как реализовать параллельные вычисления рекурсивной функции с помощью openmp?

15.11.2023, 20:02. Показов 1467. Ответов 16

Студворк — интернет-сервис помощи студентам
Пытаюсь реализовать алгоритм с параллельной быстрой сортировкой с помощью OpenMP. Пробовал через секции делать, получал предупреждение
c6993: при анализе кода структуры openmp игнорируются; анализируется код одного потока
Поменял на задачи или task'и , всё равно выводит то же самое предупреждение. Может есть какой-то более простой способ распараллелить сортировку? код следующий:
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
#include<iostream>
#include<omp.h>
 
using std::cout;
using std::endl;
#define count_of_count 2000000
class ParallelQuickSort {
    //keep count of threads
    int k = 0;
 
private:
    //partitioning procedure
    int partition(int *arr, int l, int r) {
        int i = l + 1;
        int j = r;
        int key = arr[l];
        int temp;
        while (true) {
            while (i < r && key >= arr[i])
                i++;
            while (key < arr[j])
                j--;
            if (i < j) {
                temp = arr[i];
                arr[i] = arr[j];
                arr[j] = temp;
            }
            else {
                temp = arr[l];
                arr[l] = arr[j];
                arr[j] = temp;
                return j;
            }
        }
    }
 
public:
    void quickSort(int* arr, int l, int r) {
        if (l < r) {
            int p = partition(arr, l, r);
 
 
#pragma omp task
                {
                    k = k + 1;
                    quickSort(arr, l, p - 1);
                }
#pragma omp task
                {
                    k = k + 1;
                    quickSort(arr, p + 1, r);
                }
#pragma omp taskwait
 
        }
    }
    //prints array
    void printArr(int arr[], int n) {
        for (int i = 0; i < n; i++)
            cout << arr[i] << " ";
        cout << endl;
    }
    //run the whole procedure
    void run() {
        double endTime = 0, startTime = 0, totalTime = 0;
        int* dan;
        dan = new int[count_of_count];
        for (int i = 0; i < count_of_count; i++) {
            dan[i] = rand() % 100000;
        }
        startTime = time(NULL);
        quickSort(dan, 0, count_of_count - 1);
        endTime = time(NULL);
        cout << endTime - startTime << "\n";
    }
};
 
int main() {
    ParallelQuickSort pqs;
    pqs.run();
    return 0;
}
0
Лучшие ответы (1)
Programming
Эксперт
39485 / 9562 / 3019
Регистрация: 12.04.2006
Сообщений: 41,671
Блог
15.11.2023, 20:02
Ответы с готовыми решениями:

Clojure Реализовать цикл while с помощью рекурсивной функции
Здравствуйте, товарищи форумчане! Дошел до изучения макросов, необходимо реализовать одну из разновидностей цикла, я остановился на цикле...

Параллельные вычисления: произведение двух векторов (OpenMP)
Помогите с написанием программы вычисления произведения двух векторов с использованием параллельного алгоритма OpenMP) Добавлено через...

Параллельные вычисления с использованием параллельных секций в OpenMP
Доброго времени суток, ребята! Ребят, нужно вычислить интеграл на отрезке для функции. Нужно разбить вручную отрезок на несколько...

16
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
15.11.2023, 20:03  [ТС]
Вот такое предупреждение
Миниатюры
Как реализовать параллельные вычисления рекурсивной функции с помощью openmp?  
0
15.11.2023, 20:10

Не по теме:

Есть раздел OpenMP.

0
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
15.11.2023, 20:13  [ТС]
Цитата Сообщение от nmcf Посмотреть сообщение
Не по теме:
Есть раздел OpenMP.
Понял, перенесу
0
Asm/C++/Delphi/Py/PHP/VBA
 Аватар для Jin X
6808 / 2048 / 238
Регистрация: 14.12.2014
Сообщений: 4,299
Записей в блоге: 12
16.11.2023, 11:47
Могу только дать статью на эту тему (как раз с примером на OpenMP):
https://iq.opengenus.org/parallel-quicksort/
0
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
16.11.2023, 18:28  [ТС]
Цитата Сообщение от Jin X Посмотреть сообщение
Могу только дать статью на эту тему (как раз с примером на OpenMP):
https://iq.opengenus.org/parallel-quicksort/
Спасибо, но я как раз оттуда брал ранее код с секциями) и он тоже не работал) походу дело в настройках проекта
0
Asm/C++/Delphi/Py/PHP/VBA
 Аватар для Jin X
6808 / 2048 / 238
Регистрация: 14.12.2014
Сообщений: 4,299
Записей в блоге: 12
17.11.2023, 11:36
У меня код по ссылке работает, ваш валится.
0
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
17.11.2023, 17:21  [ТС]
А у вас нет предупреждения под omp ? У меня тоже работает кол по ссылке, но если я уберу все конструкции omp, программа будет работать точно также, как и с ними, т е ускорения нет никакого
0
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
18.11.2023, 14:21  [ТС]
Поменял код, решил распараллелить цикл, теперь он работает примерно в 4 раза медленнее, чем с 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
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
107
108
109
#include <iostream>
#include <omp.h>
#ifndef _OPENMP
static_assert(false, "openmp support required");
#endif
using namespace std;
#define count_of_count 20000000
// Function to swap two numbers a and b
void swap(int* a, int* b)
{
    int t = *a;
    *a = *b;
    *b = t;
}
 
// Function to perform the partitioning
// of array arr[]
int partition(int *arr, int start, int end)
{
    // Declaration
    int pivot = arr[end];
    int i = (start - 1);
 
#pragma omp parallel for num_threads(8)
    for (int j = start; j <= end - 1; j++) {
        if (arr[j] < pivot) {
#pragma omp critical
            {
                i++;
                swap(&arr[i], &arr[j]);
            }
            
        }
    }
    swap(&arr[i + 1], &arr[end]);
 
    // Returning the respective index
    return (i + 1);
}
 
// Function to perform QuickSort Algorithm
// using openmp
void quicksort(int *arr, int start, int end)
{
    // Declaration
    int index;
 
    if (start < end) {
 
        // Getting the index of pivot
        // by partitioning
        index = partition(arr, start, end);
        
        // Parallel sections
 
        
            
 
            
 
 
                
                    //cout << omp_get_num_threads() << "\n";
                    quicksort(arr, start, index - 1);
                
 
            
                    //cout << omp_get_num_threads() << "             222222\n";
                    quicksort(arr, index + 1, end);
                
            
        
    }
}
 
// Driver Code
int main()
{
    // Declaration
    int *dan;
    dan = new int[count_of_count];
    double endTime = 0, startTime = 0, totalTime = 0;
 
 
    cout << "Enter the array: \n";
 
    // Taking input that array
    for (int i = 0; i < count_of_count; i++) {
        dan[i] = rand() % 100000;
    }
 
    // Calling quicksort having parallel
    // code implementation
    startTime = time(NULL);
    quicksort(dan, 0, count_of_count - 1);
    endTime = time(NULL);
    cout << endTime - startTime << "\n";
    // Printing the sorted array
    cout << "Array after Sorting is: \n";
    /*for (int i = 0; i < count_of_count; i++) {
        if (i % 1000 == 0)
        {
            cout << dan[i] << " ";
        }
    }*/
    
 
    return 0;
}
0
7804 / 6568 / 2988
Регистрация: 14.04.2014
Сообщений: 28,705
18.11.2023, 21:10
Ты точно понимаешь, что делаешь? Насколько я знаю, директива "parallel for" разделяет цикл на части; используется, если итерации независимы.
0
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
19.11.2023, 13:11  [ТС]
Цитата Сообщение от nmcf Посмотреть сообщение
Ты точно понимаешь, что делаешь? Насколько я знаю, директива "parallel for" разделяет цикл на части; используется, если итерации независимы.
Выше используются секции, но это всё равно не помогает конкретно в данном случае
0
Asm/C++/Delphi/Py/PHP/VBA
 Аватар для Jin X
6808 / 2048 / 238
Регистрация: 14.12.2014
Сообщений: 4,299
Записей в блоге: 12
19.11.2023, 16:39
pashamem7894, я бы добавил условие, чтобы слишком короткие циклы не разбивались на потоки.
Плюс, зачем тут critical? Ведь мы параллельно обрабатываем левую и правую сторону, которые не пересекаются.
0
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
19.11.2023, 16:42  [ТС]
Цитата Сообщение от Jin X Посмотреть сообщение
pashamem7894, я бы добавил условие, чтобы слишком короткие циклы не разбивались на потоки.
Плюс, зачем тут critical? Ведь мы параллельно обрабатываем левую и правую сторону, которые не пересекаются.
__________________
На счёт critical, не знаю, с примера брал. А по поводу условия про слишком короткие циклы, цикл ж происходит внутри потока уже, так как же отменить разбиение ещё до того, как известно время...
0
Asm/C++/Delphi/Py/PHP/VBA
 Аватар для Jin X
6808 / 2048 / 238
Регистрация: 14.12.2014
Сообщений: 4,299
Записей в блоге: 12
19.11.2023, 16:46
Лучший ответ Сообщение было отмечено pashamem7894 как решение

Решение

#pragma omp task if p-1-l > 1000
#pragma omp task if r-(p+1) > 1000
1
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
19.11.2023, 17:19  [ТС]
Цитата Сообщение от Jin X Посмотреть сообщение
#pragma omp task if p-1-l > 1000
#pragma omp task if r-(p+1) > 1000
А это действительно сработало, код ускорился многократно и сейчас работает быстрее последовательного, спасибо вам большое) Буду теперь пробовать как-то установить количество потоков в тасках(хотя это вроде нельзя делать). Может у вас есть информация о том, как установить определённое количество потоков на таски?)
0
Asm/C++/Delphi/Py/PHP/VBA
 Аватар для Jin X
6808 / 2048 / 238
Регистрация: 14.12.2014
Сообщений: 4,299
Записей в блоге: 12
20.11.2023, 10:45
Цитата Сообщение от pashamem7894 Посмотреть сообщение
Может у вас есть информация о том, как установить определённое количество потоков на таски?)
Вы же делали #pragma omp parallel num_threads(8) — это не срабатывает?
Ещё есть omp_set_num_threads.
0
2 / 2 / 1
Регистрация: 06.10.2021
Сообщений: 105
22.11.2023, 18:44  [ТС]
Цитата Сообщение от Jin X Посмотреть сообщение
Вы же делали #pragma omp parallel num_threads(8) — это не срабатывает?
Ещё есть omp_set_num_threads.
Пробовал в тасках задавать. я так думаю, что 1 task = 1 потоку, никак иначе
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
inter-admin
Эксперт
29715 / 6470 / 2152
Регистрация: 06.03.2009
Сообщений: 28,500
Блог
22.11.2023, 18:44
Помогаю со студенческими работами здесь

Написать программу для вычисления выражения с помощью рекурсивной функции
Помогите написать программу C++ для вычисления выражения an. Вычисление организовать в виде рекурсивной функции. {a}_{n} =...

Перевести программу с паскаля на фортран, используя Параллельные вычисления и OpenMP
Есть несложная программа на паскале, надо ее перевести в фортран с помощью параллельных вычислений и библиотеки OpenMP и MPI. Готовых...

Параллельные вычисления: найти число пи с помощью ряда ряд Лейбница
пишу программу для вычисления числа пи с помощью ряда ряд Лейбница. (pi/4=1-1/3+1/5-1/7+... ) получаются всегда разные значения и не...

С помощью рекурсивной функции вывести значение функции sin(x) от А до B включая с шагом step
Дан прототип функции void print_tab (float A, float B, float step) Как с помощью рекурсии вывести значение функции sin(x) от А до B...

OPENmp параллельные потоки QT
#include &lt;iostream&gt; #include &quot;omp.h&quot; using namespace std; int main() { #pragma omp parallel { ...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
17
Ответ Создать тему
Новые блоги и статьи
PhpStorm 2025.3: WSL Terminal всегда стартует в ~
and_y87 14.12.2025
PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .
Как объединить две одинаковые БД Access с разными данными
VikBal 11.12.2025
Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.
Новый ноутбук
volvo 07.12.2025
Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro
Музыка, написанная Искусственным Интеллектом
volvo 04.12.2025
Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .
От async/await к виртуальным потокам в Python
IndentationError 23.11.2025
Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов
Argus19 22.11.2025
Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .
Сколько Государство потратило денег на меня, обеспечивая инсулином.
Programma_Boinc 20.11.2025
Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .
Ломающие изменения в C#.NStar Alpha
Etyuhibosecyu 20.11.2025
Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .
Мысли в слух
kumehtar 18.11.2025
Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .
Создание Single Page Application на фреймах
krapotkin 16.11.2025
Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2025, CyberForum.ru