C++ MPI_Send（109）：无效列组的值为-2，但必须为非负且小于1
有人能帮我找到代码中的以下错误吗： MPI_发送中的致命错误：无效列组，错误堆栈： MPI_发送（190）：MPI_发送（buf=0x1123060，计数=40，MPI_浮点，dest=MPI_任意源，标记=25，MPI_通信世界）失败 MPI_Send（109）：无效列组的值为-2，但必须为非负且小于1“C++ MPI_Send（109）：无效列组的值为-2，但必须为非负且小于1,c++,mpi,C++,Mpi,有人能帮我找到代码中的以下错误吗： MPI_发送中的致命错误：无效列组，错误堆栈： MPI_发送（190）：MPI_发送（buf=0x1123060，计数=40，MPI_浮点，dest=MPI_任意源，标记=25，MPI_通信世界）失败 MPI_Send（109）：无效列组的值为-2，但必须为非负且小于1“ intmain（intargc，char*argv[]） { name=“quick.dat”； n=40； 自动启动=标准：：时钟：：高分辨率时钟：：现在（）； int i； 双sp； 双
intmain（intargc，char*argv[]）
{
name=“quick.dat”；
n=40；
自动启动=标准：：时钟：：高分辨率时钟：：现在（）；
int i；
双sp；
双*ap、*ae、*aw、*aww、*cte；
ap=新的双精度[n]；
ae=新的双精度[n]；
aw=新的双精度[n]；
aww=新的双精度[n]；
cte=新双精度[n]；
全国公共广播电台（nproc）的int my_排名；
MPI_状态；
MPI_Init（&argc，&argv）；
MPI通信等级（MPI通信世界和我的通信等级）；
MPI通信大小（MPI通信世界和nproc）；
如果（我的排名==0）
{
MPI_发送（&ap[0]，n，MPI_浮点，MPI_任意源，25，MPI_通信世界）；
MPI_发送（&ae[0]，n，MPI_浮点，MPI_任意源，25，MPI_通信世界）；
MPI_发送（&aw[0]，n，MPI_浮点，MPI_任意源，25，MPI_通信世界）；
MPI_发送（&aww[0]，n，MPI_浮点，MPI_任意源，25，MPI_通信世界）；
MPI_发送（&cte[0]，n，MPI_浮点，MPI_任意源，25，MPI_通信世界）；
对于（i=0；i<（n-1）/nproc+1；i++）
{
如果（i==0）
{
ae[i]=a（（i+1）*1.0/n）*D（（i+1）*1.0/n）+
1./3.0*a（i*1.0/n）*D（i*1.0/n）-
3.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）；
cte[i]=-（8./3.*a（i*1.0/n）*D（i*1.0/n）+
2.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）+
a（i*1.0/n）*u（i*1.0/n））；
ap[i]=-（ae[i]-cte[i]+a（（i+1）*1.0/n）*u（（i+1）*1.0/n）-
a（i*1.0/n）*u（i*1.0/n））；
继续；
}
如果（i==1）
{
aw[i]=a（i*1.0/n）*D（i*1.0/n）+
7./8.*a（i*1.0/n）*u（i*1.0/n）+
1.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）；
热膨胀系数[i]=0.25*a（i*1.0/n）*u（i*1.0/n）；
ae[i]=a（（i+1）*1.0/n）*D（（i+1）*1.0/n）-
3.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）；
ap[i]=-（ae[i]+aw[i]-cte[i]+a（（i+1）*1.0/n）*u（（i+1）*1.0/n）-
a（i*1.0/n）*u（i*1.0/n））；
继续；
}
如果（i==n-1）
{
aw[i]=a（i*1.0/n）*D（i*1.0/n）+
1./3.0*a（（i+1）*1.0/n）*D（（i+1）*1.0/n）+
6.0/8.0*a（i*1.0/n）*u（i*1.0/n）；
aww[i]=-1.0/8.0*a（i*1.0/n）*u（i*1.0/n）；
cte[i]=0；
sp=8.0/3.0*a（（i+1）*1.0/n）*D（（i+1）*1.0/n）-
（i+1）*1.0/n）*u（（i+1）*1.0/n）；
ap[i]=-（aww[i]+aw[i]+sp+a（（i+1）*1.0/n）*u（（i+1）*1.0/n）-
a（i*1.0/n）*u（i*1.0/n））；
继续；
}
如果（i！=0&&i！=1&&i！=n-1）
{
ae[i]=a（（i+1）*1.0/n）*D（（i+1）*1.0/n）-
3.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）；
aw[i]=a（i*1.0/n）*D（i*1.0/n）+
6./8.*a（i*1.0/n）*u（i*1.0/n）+
1.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）；
aww[i]=-1.0/8.0*a（i*1.0/n）*u（i*1.0/n）；
cte[i]=0；
ap[i]=-（aww[i]+aw[i]+ae[i]-cte[i]+
a（（i+1）*1.0/n）*u（（i+1）*1.0/n）-
a（i*1.0/n）*u（i*1.0/n））；
}
}
MPI_Recv（&ap[0]，n，MPI_FLOAT，MPI_ANY_SOURCE，25，MPI_COMM_WORLD，&status）；
MPI_Recv（&ae[0]，n，MPI_FLOAT，MPI_ANY_SOURCE，25，MPI_COMM_WORLD，&status）；
MPI_记录（&aw[0]，n，MPI_浮动，MPI_任意源，25，MPI_通信世界和状态）；
MPI_Recv（&aww[0]，n，MPI_FLOAT，MPI_ANY_SOURCE，25，MPI_COMM_WORLD，&status）；
MPI_Recv（&cte[0]，n，MPI_FLOAT，MPI_ANY_SOURCE，25，MPI_COMM_WORLD，&status）；
}
否则如果（我的排名！=0&&my的排名！=nproc-1）
{
MPI_Recv（&ap[0]，n，MPI_FLOAT，0，25，MPI_COMM_WORLD，&status）；
MPI_Recv（&ae[0]，n，MPI_FLOAT，0，25，MPI_COMM_WORLD，&status）；
MPI_记录（&aw[0]，n，MPI_浮动，0，25，MPI_通信世界和状态）；
MPI_Recv（&aww[0]，n，MPI_FLOAT，0，25，MPI_COMM_WORLD，&status）；
MPI_Recv（和cte[0]，n，MPI_FLOAT，0，25，MPI_COMM_WORLD，和status）；
对于（i=（（n-1）/nproc）*我的排名+1；i<（（n-1）/nproc）*（我的排名+1）+1；i++）
{
如果（i==0）
{
ae[i]=a（（i+1）*1.0/n）*D（（i+1）*1.0/n）+
1./3.0*a（i*1.0/n）*D（i*1.0/n）-
3.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）；
cte[i]=-（8./3.*a（i*1.0/n）*D（i*1.0/n）+
2.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）+
a（i*1.0/n）*u（i*1.0/n））；
ap[i]=-（ae[i]-cte[i]+a（（i+1）*1.0/n）*u（（i+1）*1.0/n）-
a（i*1.0/n）*u（i*1.0/n））；
继续；
}
如果（i==1）
{
aw[i]=a（i*1.0/n）*D（i*1.0/n）+
7./8.*a（i*1.0/n）*u（i*1.0/n）+
1.0/8.0*a（（i+1）*1.0/n）*u（（i+1）*1.0/n）；
cte[i]=
int main(int argc, char* argv[])
{
name = "quick.dat";
n = 40;
auto start = std::chrono::high_resolution_clock::now();
int i;
double sp;
double *ap, *ae, *aw, *aww, *cte;
ap = new double[n];
ae = new double[n];
aw = new double[n];
aww = new double[n];
cte = new double[n];
int my_rank, nproc;
MPI_Status status;
MPI_Init(&argc, &argv);
MPI_Comm_rank(MPI_COMM_WORLD, &my_rank);
MPI_Comm_size(MPI_COMM_WORLD, &nproc);
if (my_rank == 0)
{
MPI_Send(&ap[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD);
MPI_Send(&ae[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD);
MPI_Send(&aw[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD);
MPI_Send(&aww[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD);
MPI_Send(&cte[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD);
for (i = 0; i < (n - 1) / nproc + 1; i++)
{
if (i == 0)
{
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) +
1. / 3.0 * a(i * 1.0 / n) * D(i * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
cte[i] = -(8. / 3. * a(i * 1.0 / n) * D(i * 1.0 / n) +
2.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) +
a(i * 1.0 / n) * u(i * 1.0 / n));
ap[i] = -(ae[i] - cte[i] + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i == 1)
{
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
7. / 8. * a(i * 1.0 / n) * u(i * 1.0 / n) +
1.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
cte[i] = 0.25 * a(i * 1.0 / n) * u(i * 1.0 / n);
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
ap[i] = -(ae[i] + aw[i] - cte[i] + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i == n - 1)
{
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
1. / 3.0 * a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) +
6.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
aww[i] = -1.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
cte[i] = 0;
sp = 8.0 / 3.0 * a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
ap[i] = -(aww[i] + aw[i] + sp + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i != 0 && i != 1 && i != n - 1)
{
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
6. / 8. * a(i * 1.0 / n) * u(i * 1.0 / n) +
1.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
aww[i] = -1.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
cte[i] = 0;
ap[i] = -(aww[i] + aw[i] + ae[i] - cte[i] +
a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
}
}
MPI_Recv(&ap[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&ae[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&aw[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&aww[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&cte[0], n, MPI_FLOAT, MPI_ANY_SOURCE, 25, MPI_COMM_WORLD, &status);
}
else if (my_rank != 0 && my_rank != nproc - 1)
{
MPI_Recv(&ap[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&ae[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&aw[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&aww[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&cte[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
for (i = ((n - 1) / nproc) * my_rank + 1; i < ((n - 1) / nproc) * (my_rank + 1) + 1; i++)
{
if (i == 0)
{
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) +
1. / 3.0 * a(i * 1.0 / n) * D(i * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
cte[i] = -(8. / 3. * a(i * 1.0 / n) * D(i * 1.0 / n) +
2.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) +
a(i * 1.0 / n) * u(i * 1.0 / n));
ap[i] = -(ae[i] - cte[i] + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i == 1)
{
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
7. / 8. * a(i * 1.0 / n) * u(i * 1.0 / n) +
1.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
cte[i] = 0.25 * a(i * 1.0 / n) * u(i * 1.0 / n);
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
ap[i] = -(ae[i] + aw[i] - cte[i] + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i == n - 1)
{
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
1. / 3.0 * a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) +
6.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
aww[i] = -1.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
cte[i] = 0;
sp = 8.0 / 3.0 * a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
ap[i] = -(aww[i] + aw[i] + sp + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i != 0 && i != 1 && i != n - 1)
{
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
6. / 8. * a(i * 1.0 / n) * u(i * 1.0 / n) +
1.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
aww[i] = -1.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
cte[i] = 0;
ap[i] = -(aww[i] + aw[i] + ae[i] - cte[i] +
a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
}
}
MPI_Send(&ap[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD);
MPI_Send(&ae[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD);
MPI_Send(&aw[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD);
MPI_Send(&aww[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD);
MPI_Send(&cte[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD);
}
else if (my_rank == nproc - 1 && nproc != 1)
{
MPI_Recv(&ap[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&ae[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&aw[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&aww[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
MPI_Recv(&cte[0], n, MPI_FLOAT, 0, 25, MPI_COMM_WORLD, &status);
for (i = ((n - 1) / nproc) * my_rank + 1; i < n; i++)
{
if (i == 0)
{
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) +
1. / 3.0 * a(i * 1.0 / n) * D(i * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
cte[i] = -(8. / 3. * a(i * 1.0 / n) * D(i * 1.0 / n) +
2.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) +
a(i * 1.0 / n) * u(i * 1.0 / n));
ap[i] = -(ae[i] - cte[i] + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i == 1)
{
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
7. / 8. * a(i * 1.0 / n) * u(i * 1.0 / n) +
1.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
cte[i] = 0.25 * a(i * 1.0 / n) * u(i * 1.0 / n);
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
ap[i] = -(ae[i] + aw[i] - cte[i] + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i == n - 1)
{
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
1. / 3.0 * a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) +
6.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
aww[i] = -1.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
cte[i] = 0;
sp = 8.0 / 3.0 * a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
ap[i] = -(aww[i] + aw[i] + sp + a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
continue;
}
if (i != 0 && i != 1 && i != n - 1)
{
ae[i] = a((i + 1) * 1.0 / n) * D((i + 1) * 1.0 / n) -
3.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
aw[i] = a(i * 1.0 / n) * D(i * 1.0 / n) +
6. / 8. * a(i * 1.0 / n) * u(i * 1.0 / n) +
1.0 / 8.0 * a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n);
aww[i] = -1.0 / 8.0 * a(i * 1.0 / n) * u(i * 1.0 / n);
cte[i] = 0;
ap[i] = -(aww[i] + aw[i] + ae[i] - cte[i] +
a((i + 1) * 1.0 / n) * u((i + 1) * 1.0 / n) -
a(i * 1.0 / n) * u(i * 1.0 / n));
}
}
MPI_Send(&ap[0], n, MPI_FLOAT, 0, MPI_ANY_TAG, MPI_COMM_WORLD);
MPI_Send(&ae[0], n, MPI_FLOAT, 0, MPI_ANY_TAG, MPI_COMM_WORLD);
MPI_Send(&aw[0], n, MPI_FLOAT, 0, MPI_ANY_TAG, MPI_COMM_WORLD);
MPI_Send(&aww[0], n, MPI_FLOAT, 0, MPI_ANY_TAG, MPI_COMM_WORLD);
MPI_Send(&cte[0], n, MPI_FLOAT, 0, MPI_ANY_TAG, MPI_COMM_WORLD);
}
// double t22 = MPI_Wtime();
MPI_Finalize();
// printf("[IJK] Compute time [s] : %6.3f \n", t22-t11 );
iteration(aww, aw, ap, ae, cte);
auto end = std::chrono::high_resolution_clock::now();
std::chrono::duration<double> diff = end - start;
std::cout << "Time to Execute the FVM code "
<< " is : " << diff.count() << " s\n";
return 0;
}
MPI_Send(&ap[0], n, MPI_FLOAT, 0, MPI_COMM_WORLD);