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

Многопоточность при сортировке массивов - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Компиляторы. Какие преимущества между "кнопкой" и консолью? http://www.cyberforum.ru/cpp-beginners/thread1147477.html
В наличии два компилятора: от vs и intel с оболочками и консоли. Какие преимущества между "кнопкой" и консолью?
C++ Поможите реализовать эту программу поиск ближайших объектов class SPHERE_SET_TREE { public: class POSITION; class ITEM; class NODE; // Структура задает элемент множества сфер struct SPHERE { http://www.cyberforum.ru/cpp-beginners/thread1147473.html
Логическое выражение преобразовать в бинарное дерево C++
Преобразовать логическое выражение в бинарное дерево и вывести его.Выражение может содержать логическое и,или,отрицание(&,|,! соответственно).Например: 0|(0 | 1)&X(значение х пользователь вводит руками).
Выяснить, верно ли, что в последовательности имеются пять идущих подряд букв е C++
Условия: Даны натуральное число n, символы s1, …, sn. Выяснить, верно ли, что в последовательности s1, …, sn имеются пять идущих подряд букв е. Примерное решение(моё, может быть абсолютно неправильным т.к. делал по примеру): #include <conio.h> #include <clocale> #include <iostream> #include <cmath> using namespace std; int main () {setlocale(LC_ALL, "Russian_Russia.1251"); int n,i,k=0;
C++ Советы по программе, которая должна сделать перевод литров в четверти http://www.cyberforum.ru/cpp-beginners/thread1147446.html
Народ, нужны ваши советы, надо зделать програму на языке C++, перевод литров в чверти. В одной чверте 3.08 литров. Я пытался зделать сам и оно получилось (один раз, без понятия почему потом не хотело), но потом начелов выдавать такое, до этой програмы мне такого не писало: "suchera.exe": Загружено: "D:\suchera\Debug\suchera.exe", Символы загружены. "suchera.exe": Загружено:...
C++ По заданным значениям х и n найти элементы ряда и вывести на экран 2. По заданным значениям х и n найти элементы ряда и вывести на экран. x1/2, x2/3, x3/4, ...... , xn/(n+1). подробнее

Показать сообщение отдельно
NEvOl
19 / 18 / 0
Регистрация: 13.08.2012
Сообщений: 726
18.04.2014, 11:02  [ТС]     Многопоточность при сортировке массивов
хм... с сортировкой я вроде разобрался, действительно работает быстрее, но вот меня беспокоит функция создания kd дерева, в ней так же происходит сортировка массивов, если я использую 3 потока для сортировки то построение выполняется в разы медленнее чем когда одним потоком сортирую, код далек от совершенства, но рабочий, подскажите пожалуйста почему так ? при проверке показало что если использовать 3 потока для сортировки то времени уходит около 1000 секунд, а при использовании одного потока для сортировки всего 20 с, это у меня вообще в голове не укладывается как так ?
переделал немного код:
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
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
void nodeKdTree::CreateTree(photon *p, nodePhoton *newP, int num, int iStart, int rootIndex, int iBranch, bool thr)
{
    if(num == 1)
    {
        this->pht = p[0];
        newP[iStart].pos = p[0].pos;
        newP[iStart].energy = p[0].energy;
        newP[iStart].l = -1;
        newP[iStart].r = -1;
        if(iBranch == 1)
            newP[rootIndex].r = iStart;
        else
            newP[rootIndex].l = iStart;
        this->left = NULL;
        this->right = NULL;
        this->key = 0;
        this->lNum = 0;
        this->rNum = 0;
        this->numInArr = 0;
        return;
    }
    photon *x = new photon[num];
    photon *y = new photon[num];
    photon *z = new photon[num];
 
    std::memcpy(x, p,sizeof(photon)*num);
    std::memcpy(y, p,sizeof(photon)*num);
    std::memcpy(z, p,sizeof(photon)*num);
 
    if(thr)
    {
        std::thread th1(qSortX, x, num-1);
        std::thread th2(qSortY, y, num-1);
        std::thread th3(qSortZ, z, num-1);
 
        th1.join();
        th2.join();
        th3.join();
    }
    else
    {
        qSortX(x, num-1);
        qSortY(y, num-1);
        qSortZ(z, num-1);
    }
 
    int dx = abs(x[0].pos.x - x[num-1].pos.x);
    int dy = abs(y[0].pos.y - y[num-1].pos.y);
    int dz = abs(z[0].pos.z - z[num-1].pos.z);
    this->left = new nodeKdTree();
    this->right = new nodeKdTree();
 
 
    int newIndex(-1);
    if((dx > dy && dx > dz) || (dx == dy && dx > dz) || (dx == dz && dx > dy))
    {
        int m = num/2;
        this->numInArr = m;
        int lnum;
        int rnum;
        this->pht = x[m];
        newIndex = iStart+m;
        newP[newIndex].pos = x[m].pos;
        newP[newIndex].energy = x[m].energy;
        if(iBranch == 1)
            newP[rootIndex].r = newIndex;
        if(iBranch == -1)
            newP[rootIndex].l = newIndex;
        if(rootIndex == -1)
        {
            newP[newIndex].l = -1;
            newP[newIndex].r = -1;
        }
        newP[newIndex].k = 1;
        this->key = 1;
        lnum = m;
        photon *lx = new photon[lnum];
        std::memcpy(lx, x, sizeof(photon)*lnum);
        photon *rx;
        if(num%2)
        {
            rnum = m;
            rx = new photon[rnum];
            std::memcpy(rx, x+m+1, sizeof(photon)*rnum);
        }
        else
        {
            rnum = num-m-1;
            if(rnum > 0)
            {
                rx = new photon[rnum];
                std::memcpy(rx, x+m+1, sizeof(photon)*rnum);
            }
        }
        this->left->CreateTree(lx, newP, lnum, iStart, newIndex, -1, thr);
        delete []lx;
        this->lNum = lnum;
        this->rNum = rnum;
        if(rnum > 0)
        {
            this->right->CreateTree(rx, newP, rnum, newIndex+1, newIndex, 1, thr);
            delete []rx;
        }
    }
    else
        if((dy > dx && dy > dz) || (dy == dz && dy > dx))
        {
            int m = num/2;
            this->numInArr = m;
            int lnum;
            int rnum;
            this->pht = y[m];
            newIndex = iStart+m;
            newP[newIndex].pos = y[m].pos;
            newP[newIndex].energy = y[m].energy;
            if(iBranch == 1)
                newP[rootIndex].r = newIndex;
            if(iBranch == -1)
                newP[rootIndex].l = newIndex;
            if(rootIndex == -1)
            {
                newP[newIndex].l = -1;
                newP[newIndex].r = -1;
            }
            newP[newIndex].k = 2;
            this->key = 2;
            lnum = m;
            photon *ly = new photon[lnum];
            std::memcpy(ly, y, sizeof(photon)*lnum);
            photon *ry;
            if(num%2)
            {
                rnum = m;
                ry = new photon[rnum];
                std::memcpy(ry, y+m+1, sizeof(photon)*rnum);
            }
            else
            {
                rnum = num-m-1;
                if(rnum > 0)
                {
                    ry = new photon[rnum];
                    std::memcpy(ry, y+m+1, sizeof(photon)*rnum);
                }
            }
            this->left->CreateTree(ly, newP, lnum, iStart, newIndex, -1, thr);
            delete []ly;
            this->lNum = lnum;
            this->rNum = rnum;
            if(rnum > 0)
            {
                this->right->CreateTree(ry, newP, rnum, newIndex+1, newIndex, 1, thr);
                delete []ry;
            }
        }
        else
            {
                int m = num/2;
                this->numInArr = m;
                int lnum;
                int rnum;
                this->pht = z[m];
                newIndex = iStart+m;
                newP[newIndex].pos = z[m].pos;
                newP[newIndex].energy = z[m].energy;
                if(iBranch == 1)
                    newP[rootIndex].r = newIndex;
                if(iBranch == -1)
                    newP[rootIndex].l = newIndex;
                if(rootIndex == -1)
                {
                    newP[newIndex].l = -1;
                    newP[newIndex].r = -1;
                }
                newP[newIndex].k = 3;
                this->key = 3;
                lnum = m;
                photon *lz = new photon[lnum];
                std::memcpy(lz, z, sizeof(photon)*lnum);
                photon *rz;
                if(num%2)
                {
                        rnum = m;
                        rz = new photon[rnum];
                        std::memcpy(rz, z+m+1, sizeof(photon)*rnum);
                }
                else
                {
                    rnum = num-m-1;
                    if(rnum > 0)
                    {
                        rz = new photon[rnum];
                        std::memcpy(rz, z+m+1, sizeof(photon)*rnum);
                    }
                }
                this->left->CreateTree(lz, newP, lnum, iStart, newIndex, -1, thr);
                delete []lz;
                this->lNum = lnum;
                this->rNum = rnum;
                if(rnum > 0)
                {
                    this->right->CreateTree(rz, newP,  rnum, newIndex+1, newIndex, 1, thr);
                    delete []rz;
                }
            }
    return;
}
 
void qSortZ(photon *p, int e)
{
    std::qsort(p, e+1, sizeof(photon), funcZ);
}
 
void qSortY(photon *p, int e)
{
    std::qsort(p, e+1, sizeof(photon), funcY);
}
 
 
void qSortX(photon *p, int e)
{
    std::qsort(p, e+1, sizeof(photon), funcX);
}
 
int funcX(const void* a, const void* b)
{
    if((*(photon*)a).pos.x > (*(photon*)b).pos.x)
        return 1;
    if((*(photon*)a).pos.x < (*(photon*)b).pos.x)
        return -1;
    return 0;
}
 
int funcY(const void* a, const void* b)
{
    if((*(photon*)a).pos.y > (*(photon*)b).pos.y)
        return 1;
    if((*(photon*)a).pos.y < (*(photon*)b).pos.y)
        return -1;
    return 0;
}
 
int funcZ(const void* a, const void* b)
{
    if((*(photon*)a).pos.z > (*(photon*)b).pos.z)
        return 1;
    if((*(photon*)a).pos.z < (*(photon*)b).pos.z)
        return -1;
    return 0;
}
пояснения:
в функцию передается исходный массив с координатами точек в 3D(p), новый массив который постепенно заполняется и представляет само kd дерево(newP), количество элементов в переданном массиве p (num), индекс начала переданного подмассива в исходном массиве p (iStart), индекс узла дерева чьи ветви рассчитываются в текущий момент в массиве newp (rootIndex), индекс правой - 1 и левой - -1 ветви (iBranch) и последним параметром является метод сортировки массивов 3-мя потоками или одним (thr). Суть функции заключается в том что сортируются координаты точек по 3-м осям (x, y, z) затем находится по какой из осей максимальное расстояние между крайними точками, массив разбивается на 2 части, центральный элемент сохраняется после этого остальные части массива (правая и левая) отправляются рекурсивно в эту же функцию и так пока весь массив не переберется. помогите пожалуйста разобраться.
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru