Fatal编程技术网
  • c/
  • 我的MPI jacobi迭代程序给出了错误的结果

我的MPI jacobi迭代程序给出了错误的结果

c mpi

我的MPI jacobi迭代程序给出了错误的结果,c,mpi,linear-algebra,C,Mpi,Linear Algebra,在C中运行以下程序时，我遇到了上述错误 该程序首先生成3x3数组，然后进行Jacobi迭代。它使用MPI库。我不知道代码的哪些部分是错误的 #include <stdio.h> #include <string.h> #include <mpi.h> #include <math.h> // l2-norm // #include <time.h> int main(int argc, char **argv) { int numpro

在C中运行以下程序时，我遇到了上述错误

该程序首先生成3x3数组，然后进行Jacobi迭代。它使用MPI库。我不知道代码的哪些部分是错误的

#include <stdio.h>
#include <string.h>
#include <mpi.h>
#include <math.h> // l2-norm //
#include <time.h>
int main(int argc, char **argv)
{
int numprocs, myid;
MPI_Init(&argc,&argv);
MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
MPI_Comm_rank(MPI_COMM_WORLD, &myid);

double a[3][3];
double b[3];
double x[3]={0};
double xa[3]={0};
double xnew[3]={0};
double y[3]={0};
float sigancha;
time_t startTime=0, endTime=0;
int n=3;
int i, j =0;
int k=0;
int o;
int hoessu=300;
int minhoessu=300;
double sum=1;
int numsent =0;
int ans;
int row;
MPI_Status status;
int sender;
int po;
double *buffer;
/* synchronization */
MPI_Barrier(MPI_COMM_WORLD);
for (i=0; i<n; i++){
        b[i]=i*100;
        for (j=0; j<n; j++) {
                a[i][j]=((i+j)%10);
                if (i==j) {a[i][j]+=5000;}
        }
        x[i]=b[i]/a[i][i];
}
/* run if sum is greater than 0.0002 */

for (k=0; k<hoessu&&sum>0.0002||k<minhoessu; k++) {

numsent = 0;
for (o=myid+1; o<n+1; o+=numprocs) {
    i=o-1;
    xa[i]=b[i]+a[i][i]*x[i];
    for (j=0; j<n; j++) {
    xa[i]-=a[i][j]*x[j];
    }
xnew[i]=xa[i]/a[i][i];
/*send xnew[i] to master*/
MPI_Send(&xnew[i],1,MPI_DOUBLE,0,i,MPI_COMM_WORLD);
}

if (myid == 0){
/*get xnew[i]*/
for (i=0; i<n; i++) {
MPI_Recv(&ans, 1, MPI_DOUBLE, MPI_ANY_SOURCE, MPI_ANY_TAG, MPI_COMM_WORLD, &status);
sender = status.MPI_SOURCE;
row = status.MPI_TAG;
xnew[row] = ans;
}

/*calculates sum at master*/

for (j=0; j<n; j++){
    sum=0.0;
    sum+=(xnew[j]-x[j])*(xnew[j]-x[j]);
    x[j]=xnew[j];
}
sum=pow(sum,0.5);
MPI_Bcast(&x[0], n, MPI_DOUBLE, 0, MPI_COMM_WORLD);
}
}

if (myid == 0){
        endTime=clock();
        sigancha=(float)(endTime-startTime)/(CLOCKS_PER_SEC);
        printf("finished\n");
        for (j=0; j<n; j++) {
                printf("x[%d]=%fl\n",j+1,xnew[j]);
        }

        printf("iteration; %d times itereation are done. \n l2-norm, error is %fl .\n %f seceonds are used. \n",k,sum,sigancha);
}
        MPI_Finalize();

}
结果

finished
x[1]=-1736884775.000000l
x[2]=-370936800.000000l
x[3]=2118301216.000000l
iteration; 300 times itereation are done. 
 l2-norm, error is 34332272.000000l .
 0.020000 seceonds are used.
然而，该结果并非预期结果。如果这个程序工作得很好，它应该给出和串行雅可比迭代相同的结果。 是的

x[1]=-0.000020l
x[2]=-0.019968l
x[3]=0.399956l
我不知道为什么这个程序会产生错误的结果。

\include

#include <stdio.h>
#include <string.h>
#include <mpi.h>
#include <math.h> // l2-norm //
#include <time.h>
int main(int argc, char **argv)
{
int numprocs, myid;
MPI_Init(&argc,&argv);
MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
MPI_Comm_rank(MPI_COMM_WORLD, &myid);

double a[700][700];
double b[700];
double x[700]={0};
double xa[700]={0};
double xnew[700]={0};
double y[700]={0};
float sigancha;
time_t startTime=0, endTime=0;
int n=700;
int i, j =0;
int k=0;
int o;
int hoessu=300;
int minhoessu=300;
double sum=1;
int numsent =0;
int ans;
int row;
MPI_Status status;
int sender;
int po;
double *buffer;
/* synchronization */
MPI_Barrier(MPI_COMM_WORLD);
for (i=0; i<n; i++){
        b[i]=i*100;
        for (j=0; j<n; j++) {
                a[i][j]=((i+j)%10);
                if (i==j) {a[i][j]+=10000;}
        }
        x[i]=b[i]/a[i][i];
}
/* run if sum is greater than 0.0002 */

for (k=0; k<hoessu&&sum>0.0002||k<minhoessu; k++) {

numsent = 0;
for (o=myid+1; o<n+1; o+=numprocs) {
    i=o-1;
    xa[i]=b[i]+a[i][i]*x[i];
    for (j=0; j<n; j++) {
    xa[i]-=a[i][j]*x[j];
    }
xnew[i]=xa[i]/a[i][i];
/*send xnew[i] to master*/
ans=xnew[i];
MPI_Allgather(&xnew[i],1,MPI_DOUBLE,&xnew[i],1,MPI_DOUBLE,MPI_COMM_WORLD);
}

if (myid == 0){

/*calculates sum at master*/

for (j=0; j<n; j++){
    sum=0.0;
    sum+=(xnew[j]-x[j])*(xnew[j]-x[j]);
    x[j]=xnew[j];
}
sum=pow(sum,0.5);
MPI_Bcast(&x[0], n, MPI_DOUBLE, 0, MPI_COMM_WORLD);
}
}

if (myid == 0){
        endTime=clock();
        sigancha=(float)(endTime-startTime)/(CLOCKS_PER_SEC);
        printf("finished\n");
        for (j=0; j<n; j++) {
                printf("x[%d]=%fl\n",j+1,xnew[j]);
        }

        printf("iteration; %d times itereation are done. \n l2-norm, error is %fl .\n %f seceonds are used. \n",k,sum,sigancha);
}
        MPI_Finalize();

}


#包括
#包括
#包含//l2范数//
#包括
int main（int argc，字符**argv）
{
int numprocs，myid；
MPI_Init（&argc，&argv）；
MPI通信大小（MPI通信世界和numprocs）；
MPI通信等级（MPI通信世界和myid）；
双a[700][700]；
双b[700]；
双x[700]={0}；
双xa[700]={0}；
双xnew[700]={0}；
双y[700]={0}；
浮子；
开始时间=0，结束时间=0；
int n=700；
int i，j=0；
int k=0；
INTO；
int-hoesu=300；
int minhoesu=300；
双和=1；
int numsent=0；
INTANS；
int行；
MPI_状态；
int发送器；
国际邮政；
双*缓冲器；
/*同步*/
MPI_屏障（MPI_通信世界）；

对于（i=0；该消息不言自明！您的应用程序
MPI\u Send（）
对
i+1
进行排序，并且该排序很可能大于或等于
num\u procs
，尤其是因为您硬编码了
n-3
并运行了单个MPI任务。哦，我应该将i+1更改为0。但这还不够。观察到错误的结果。请正确缩进代码，并编辑带有obs描述的问题注释和预期的结果。注释很难阅读，但是从我可以看到的，<代码>错误是0 000000。<代码>看起来很好。请考虑下面的任何新问题。总是包括一个正确的格式。一旦解决了其中的一部分，不要改变你的问题，并且遇到完全不同的错误。如果你自己找到了解决方案。，将其作为您自己问题的答案发布。

#include <stdio.h>
#include <string.h>
#include <mpi.h>
#include <math.h> // l2-norm //
#include <time.h>
int main(int argc, char **argv)
{
int numprocs, myid;
MPI_Init(&argc,&argv);
MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
MPI_Comm_rank(MPI_COMM_WORLD, &myid);

double a[700][700];
double b[700];
double x[700]={0};
double xa[700]={0};
double xnew[700]={0};
double y[700]={0};
float sigancha;
time_t startTime=0, endTime=0;
int n=700;
int i, j =0;
int k=0;
int o;
int hoessu=300;
int minhoessu=300;
double sum=1;
int numsent =0;
int ans;
int row;
MPI_Status status;
int sender;
int po;
double *buffer;
/* synchronization */
MPI_Barrier(MPI_COMM_WORLD);
for (i=0; i<n; i++){
        b[i]=i*100;
        for (j=0; j<n; j++) {
                a[i][j]=((i+j)%10);
                if (i==j) {a[i][j]+=10000;}
        }
        x[i]=b[i]/a[i][i];
}
/* run if sum is greater than 0.0002 */

for (k=0; k<hoessu&&sum>0.0002||k<minhoessu; k++) {

numsent = 0;
for (o=myid+1; o<n+1; o+=numprocs) {
    i=o-1;
    xa[i]=b[i]+a[i][i]*x[i];
    for (j=0; j<n; j++) {
    xa[i]-=a[i][j]*x[j];
    }
xnew[i]=xa[i]/a[i][i];
/*send xnew[i] to master*/
ans=xnew[i];
MPI_Allgather(&xnew[i],1,MPI_DOUBLE,&xnew[i],1,MPI_DOUBLE,MPI_COMM_WORLD);
}

if (myid == 0){

/*calculates sum at master*/

for (j=0; j<n; j++){
    sum=0.0;
    sum+=(xnew[j]-x[j])*(xnew[j]-x[j]);
    x[j]=xnew[j];
}
sum=pow(sum,0.5);
MPI_Bcast(&x[0], n, MPI_DOUBLE, 0, MPI_COMM_WORLD);
}
}

if (myid == 0){
        endTime=clock();
        sigancha=(float)(endTime-startTime)/(CLOCKS_PER_SEC);
        printf("finished\n");
        for (j=0; j<n; j++) {
                printf("x[%d]=%fl\n",j+1,xnew[j]);
        }

        printf("iteration; %d times itereation are done. \n l2-norm, error is %fl .\n %f seceonds are used. \n",k,sum,sigancha);
}
        MPI_Finalize();

}