Прагма библиотеки omp C++ для ⇐ C++

[Гость]

I have this code, but I need help to use OpenMP

Код: Выделить всё

#pragma

s to increase its speed. I want to parallelize the

Код: Выделить всё

for

loops over variables

Код: Выделить всё

i

and

Код: Выделить всё

j

. The values of

Код: Выделить всё

n

and

Код: Выделить всё

m

can be much larger, too. For example,

Код: Выделить всё

n = 1000

,

Код: Выделить всё

m = 500

.
My code:

Код: Выделить всё

for (int i = 0; i < n; ++i) {
for (int j = 0; j < m; ++j) {
matrix[i][j] = {static_cast(distribution(gen)),
static_cast(distribution(gen))};
}
}

#pragma omp parallel
{
n = 50;
m = 50;
ReS = 0;
ImS = 0;

for (int nu = 0; nu < 1; nu++) {
#pragma omp for num_threads(4)
for (int i = 0; i < n; i++) {
for (int j = 0; j < m; j++) {
double* matrixPtr = &matrix[i][j].first;
double* matrixPtrSecond = &matrix[i][j].second;

for (double theta = 0; theta < PI; theta += theta_plus) {
for (double phi = 0; phi < 2 * PI; phi += theta_plus) {
double angle = *matrixPtr * cos(phi) * sin(theta) +
*matrixPtrSecond * sin(phi) * sin(theta);
ReS += cos(angle);
ImS += sin(angle);
}
}
}
}
}
}

Источник: https://stackoverflow.com/questions/781 ... pragma-for