Я запустил программу OpenMP для выполнения метода Якоби, и она работала очень хорошо, 2 потока выполнялись немного быстрее, чем в 2 раза по сравнению с 1 потоком, а 4 потока — в 2 раза быстрее, чем 1 поток. Я чувствовал, что все работает отлично... пока не достиг ровно 20, 22 и 24 потоков. Я продолжал ломать его, пока у меня не появилась эта простая программа
#include <stdio.h>
#include <omp.h>
int main(int argc, char *argv[]) {
int i, n, maxiter, threads, nsquared, execs = 0;
double begin, end;
if (argc != 4) {
printf("4 args\n");
return 1;
}
n = atoi(argv[1]);
threads = atoi(argv[2]);
maxiter = atoi(argv[3]);
omp_set_num_threads(threads);
nsquared = n * n;
begin = omp_get_wtime();
while (execs < maxiter) {
#pragma omp parallel for
for (i = 0; i < nsquared; i++) {
//do nothing
}
execs++;
}
end = omp_get_wtime();
printf("%f seconds\n", end - begin);
return 0;
}
И вот некоторый вывод для разных номеров потоков:
./a.out 500 1 1000
0.6765799 seconds
./a.out 500 8 1000
0.0851808 seconds
./a.out 500 20 1000
19.5467 seconds
./a.out 500 22 1000
21.2296 seconds
./a.out 500 24 1000
20.1268 seconds
./a.out 500 26 1000
0.1363 seconds
Я бы понял большое замедление, если бы оно продолжалось для всех потоков, следующих за 20, потому что я полагал, что это будет накладным расходом потока (хотя я чувствовал, что это было немного экстремально). Но даже изменение n оставляет времена 20, 22 и 24 неизменными. Изменение maxiter на 100 уменьшает его примерно до 1,9 секунды, 2,2 секунды, ..., что означает, что замедление вызывает только создание потока, а не внутренняя итерация.
Это как-то связано с попыткой ОС создать потоки, которых у нее нет? Если это что-то значит, omp_get_num_procs()
возвращает 24, и это на процессорах Intel Xeon (значит, 24 включает гиперпоточность?)
Спасибо за помощь.