Параллельная программа на С++ с использованием MPI не работает на нескольких потоках

@raizor4 · Регистрация: 03.03.2020

Студворк — интернет-сервис помощи студентам

разрабатываю для курсовой работы программу в параллельном коде. за основу взято готовое решение на паскале и переведено на с++. суть такая, что запараллелен самым простым способом на данном этапе процесс заполнения массива, фрагмент кода следующий:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
if (RankProcessa==0) lbound = RankProcessa*piece+1;
else lbound = RankProcessa*piece+1;
rbound=lbound+piece-1;
if (RankProcessa==ChisloProc-1 && Ny%2==0) rbound —;
for (int i=1; i<=Nx;i++)
for (int j=lbound;j<=rbound;j++)
T[i][j]=To;
if (RankProcessa>0)
{
for (int i=1; i<=Ny;j++)
MPI_Send(&T[i][lbound], piece, MPI_DOUBLE, 0, Tag, MPI_COMM_WORLD);
}
if (RankProcessa==0)
{
for (int i=1; i<ChisloProc;i++)
{
for (int j=1; j<Ny;j++)
MPI_Recv(&T[j][i*piece+1], 1, MPI_DOUBLE, i, Tag, MPI_COMM_WORLD, &status);
}
}

в mpich2 на нескольких процессах программа зацикливается и не работает. я подозреваю, что где-то на этапе пересылки сообщений процессы не вызывают send или recv. на одном процессе в mpich и вижуал студио при этом работает. Просьба подсказать, где ошибка и подсказать вариант распараллеливания обратного цикла в данной программе, к примеру for (i=Nx-1;i>=1;i--). весь код соответственно ниже:

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
#include <mpi.h>
#include <iostream>
#include <cmath>
#include <iomanip>
#include <stdio.h>
#include <fstream>
 
using namespace std;
 
int ChisloProc=0, RankProcessa=0;
const int mf = 102;
 
int main(int argc, char *argv[])
{
const int Tag=0;
MPI_Init (&argc, &argv);
MPI_Comm_size (MPI_COMM_WORLD, &ChisloProc);
MPI_Comm_rank (MPI_COMM_WORLD, &RankProcessa);
MPI_Status status;
ofstream fout("output.txt");
double alpha[mf+1], beta[mf+1], T[mf+1][mf+1];
double startwtime=0.0, endwtime;
int i,j,Nx,Ny;
double ai,bi,ci,fi,a,lamda,ro,c,hx,hy,tau,t_end,time,To,L,H,Th,Tc;
if (RankProcessa==0)
cout«"vvedite Nx,Ny"«endl;
scanf("%d%d",&Nx,&Ny);
//fflush(stdout);
cout«"vvedite t_end,L,H,lamda,ro,c,Th,Tc,To"«endl;
cin»t_end»L»H»lamda»ro»c»Th»Tc»To;
//scanf("%lf%lf%lf%lf%lf%lf%lf%lf%lf",&t_end,&L,&H,&lamda,&ro,&c,&Th,&Tc,&To);
fflush(stdout);
hx=L/(Nx-1);
hy=H/(Ny-1);
a=lamda/(ro*c);
tau=t_end/100.0;
int piece=Ny/ChisloProc;
int lbound, rbound;
if (RankProcessa==0)
startwtime=MPI_Wtime();
if (RankProcessa==0) lbound = RankProcessa*piece+1;
else lbound = RankProcessa*piece+1;
rbound=lbound+piece-1;
if (RankProcessa==ChisloProc-1 && Ny%2==0) rbound —;
for (int i=1; i<=Nx;i++)
for (int j=lbound;j<=rbound;j++)
T[i][j]=To;
if (RankProcessa>0)
{
for (int i=1; i<=Ny;j++)
MPI_Send(&T[i][lbound], piece, MPI_DOUBLE, 0, Tag, MPI_COMM_WORLD);
}
if (RankProcessa==0)
{
for (int i=1; i<ChisloProc;i++)
{
for (int j=1; j<Ny;j++)
MPI_Recv(&T[j][i*piece+1], 1, MPI_DOUBLE, i, Tag, MPI_COMM_WORLD, &status);
}
}
time=0;
if (RankProcessa==0)
{
while (time<t_end)
{
time+=tau;
for (int j=1;j<=Ny;j++)
{
alpha[1]=0.0;
beta[1]=Th;
 
for (int i=2; i<=Nx-1;i++)
{
ai=lamda/pow(hx,2);
bi=2.0*lamda/pow(hx,2)+ro*c/tau;
ci=lamda/pow(hx,2);
fi=(-1)*ro*c*T[i][j]/tau;
alpha[i]=ai/(bi-ci*alpha[i-1]);
beta[i]=(ci*beta[i-1]-fi)/(bi-ci*alpha[i-1]);
}
 
T[Nx][j]=Tc;
 
for (i=Nx-1;i>=1;i--)
T[i][j]=alpha[i]*T[i+1][j]+beta[i];
}
for (int i=2;i<=Nx-1;i++)
{
alpha[1]=2.0*a*tau/(2.0*a*tau+pow(hy,2));
beta[1]=pow(hy,2)*T[i][1]/(2.0*a*tau+pow(hy,2));
 
for(int j=2; j<=Ny-1;j++)
{
ai=lamda/pow(hy,2);
bi=2.0*lamda/pow(hy,2)+ro*c/tau;
ci=lamda/pow(hy,2);
fi=(-1)*ro*c*T[i][j]/tau;
alpha[j]=ai/(bi-ci*alpha[j-1]);
beta[j]=(ci*beta[j-1]-fi)/(bi-ci*alpha[j-1]);
}
T[i][Ny]=(2.0*a*tau*beta[Ny-1]+pow(hy,2)*T[i][Ny])/(2.0*a*tau*(1.0-alpha[Ny-1])+pow(hy,2));
 
for (int j=Ny-1;j>=1;j--)
T[i][j]=alpha[j]*T[i][j+1]+beta[j];
}
}
}
endwtime=MPI_Wtime();
if (RankProcessa==0)
{
fout«"Длина пластины L= "«L«endl;
fout«"Толщина пластины Н= "«H«endl;
fout«"Число узлов по координате х в пластине Nx= "«Nx«endl;
fout«"Число узлов по координате у в пластине Nу= "«Ny«endl;
fout«"Коэффициент теплопроводности материала пластины lamda= "«lamda«endl;
fout«"Плотность материала пластины ro= "«ro«endl;
fout«"Теплоемкость материала пластины с= "«c«endl;
fout«"Начальная температура То= "«To«endl;
fout«"Температура на границе х = 0 области решения Th= "«Th«endl;
fout«"Температура на границе х = L области решения Tc= "«Tc«endl;
fout«"Результат получен с шагом по координате х, hx= "«hx«endl;
fout«"Результат получен с шагом по координате y, hy= "«hy«endl;
fout«"Результат получен с шагом по времени tau= "«tau«endl;
fout«"Температурное поле в момент времени t= "«t_end«endl;
for (int i=1; i<=Nx;i++)
{
for (int j=1; j<=Ny;j++)
fout«hx*(i-1)«" | "«hy*(j-1)«" | "«T[i][j]«" | ";
fout«endl;
}
fout«endwtime-startwtime«" секунд."«endl;
}
 
MPI_Finalize();
fout.close();
return 0;
}

Новые блоги и статьи Все статьи Все блоги /
Вывод диалогового окна перед закрытием, если документ не проведён Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .	Программный контроль заполнения реквизита табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать контроль заполнения реквизита "ПричинаСписания". . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .
Вывод данных в справочнике через динамический список Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .	Программное заполнения текстового поля в реквизите формы документа Maks 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .	К слову об оптимизации kumehtar 01.04.2026 Вспоминаю начало 2000-х, университет, когда я писал на Delphi. Тогда среди программистов на форумах активно обсуждали аккуратную работу с памятью: нужно было следить за переменными, вовремя. . .	Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .