распараллеливание циклов c++

[Denis++]

Добрый день!
Подскажите, пожалуйста, как распараллелить циклы через mpi.

Вот небольшой код.

i=1;
			do{	//Левый край по координате ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~	
				SumL[1]=0;
				SumC[1]=0;
				k=1;
				do{	SumL[1]=SumL[1]+inversL[i][k]*(2*(ph[k][2]-ph[k][1])-2*H*XP);
					SumC[1]=SumC[1]+C[i][k]*V[k][1];
				}while(++k<=N);
				V1[i][1]=TP*(c1*SumL[1]-beta*V[i][1]-sin(ph[i][1])+I+Irand[i]);
				P1[i][1]=TP*SumC[1];
				//Внутренние точки по координате ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
				x=2;
				do{	SumL[x]=0;
					SumC[x]=0;
					k=1;
					do{	SumL[x]=SumL[x]
							+inversL[i][k]*(ph[k][x+1]-2*ph[k][x]+ph[k][x-1]);
						SumC[x]=SumC[x]+C[i][k]*V[k][x];
					}while(++k<=N);
					V1[i][x]=TP*(c1*SumL[x]-beta*V[i][x]-sin(ph[i][x])+I+Irand[i]);
					P1[i][x]=TP*SumC[x];
				}while(++x<=xN-1);
				//Правый край по координате ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
				SumL[xN]=0;
				SumC[xN]=0;
				k=1;
				do{	SumL[xN]=SumL[xN]+inversL[i][k]*(2*(ph[k][xN-1]-ph[k][xN])+2*H*XP);
					SumC[xN]=SumC[xN]+C[i][k]*V[k][xN];
				}while(++k<=N);
				V1[i][xN]=TP*(c1*SumL[xN]-beta*V[i][xN]-sin(ph[i][xN])+I+Irand[i]);
				P1[i][xN]=TP*SumC[xN];
			}while(++i<=N);

[Denis++]

Помочь по данному вопросу видимо никто не сможет?

[p51x]

Самый просто вариант - используйте директивы для mpi. Посложнее - надо знать алгоритм или вначале тупо разбить циклы по счетчику, если они независимы.

[Denis++]

Самый просто вариант - используйте директивы для mpi.

А можно по подробнее, спасибо.

Есть вариант использовать MPI_Send и MPI_Recv т.е. передача с блокировкой.

[Denis++]

Нету вариантов?

[p51x]

В смысле? Я ж уже вам говорил: у вас есть цикл от 1 до Н + есть м потоков, вот разбиваете [1;Н] на м промежутков, считаете все это в разных потоках, ждете, синхронизируете, получаете результат

ну или меняете алгоритм и там раскидываете на потоки