Как это распараллелить?

[nemato]

Начал учить опен мп протокол.
Как можно найти расстояние между N точек?

Как это распараллелить?
int main() {
for (i=0; i<NP-1; i++ )
for (j=i+1;j<NP;j++)
d=distance(x1,y1,x2,y2);
}

[rlib]

Так, например:

#include <stdio.h>
#include <omp.h>

#define NP 21221

int main() {
	int i, j;
	unsigned long k, sumk, expk;
	sumk = 0;

#pragma	omp parallel private(j,k) 
	{
		k = 0;
#pragma omp for
		for (i=0; i<NP-1; ++i) {
			for (j=i+1;j<NP; ++j) {
				++k;
			}
		}
#pragma omp atomic
		sumk += k;
		
	}
	if ( NP % 2 == 0 ) 
		expk = (NP/2)*(NP-1);
	else
		expk = ((NP-1)/2)*NP;
	printf("sumk=%lu\nExpected k=%lu\n", sumk, expk);
	system("PAUSE");
	return 0;
}

[nemato]

спасибо!
можно откаментировать?

[rlib]

У вас два цикла. Верхний (i) автоматом разбивается OMP на кол-во доступных потоков (равен кол-ву ваших cpu cores по дефолту). Внутренний цикл пробегается каждым из потоков. Поэтому делаем переменную j частной для каждого потока, чтобы избежать самой распространенной oшибки паралельного програминга - racing condition между потоками. В k подсчитываем кол-во итераций каждого потока (поэтому k тоже частный per поток). expk (expected k) равен кол-ву дистанций между точками, следовательно это число, которому должнa быть равна сумма всех k из всех потоков - сделано исключительно для проверки правильности работы алгоритма.
#pragma omp atomic необходима для исключения racing condition опять, но уже после выхода из цыкла каждого из потоков. Можно пользовать critical section (#pragma omp critical), но тогда будет куда тормознутее.