| 查看: 862 | 回复: 0 | |||
[交流]
MPI并行求解一维泊松方程，sendrecv函数有问题，代码如下，希望大家帮帮忙
|
|
//一维泊松方程 u(x)=x^3,cpu为四核 #include "stdio.h" #include "math.h" #define MPICH_SKIP_MPICXX //没有这句就会出现overriding virtual function differs from 'MPI::的编译错误 #include "mpi.h"//将函数库包含进来 #include <stdio.h> #pragma comment (lib, "mpi.lib"  //没有这句就会出现unresolved external symbol _MPI_Finalize的连接错误。double f(double x) {double y; y=6*x; return y; } int main( int argc, char *argv[]) { double l=1;//长度 double w=0,e=1;//左右边界条件 int num=4;//进程数 int n=5;//网格数 double dd=l/n/num;//网格大小 double u[7],u0[7]; double eps=0; int i,j; int myid, numprocs, namelen;//定义进程号、进程数、进程名 double startwtime, endwtime;//定义起止时间，为了计算效率 char processor_name[MPI_MAX_PROCESSOR_NAME];//用进程名最长的作为长度 MPI_Init(&argc,&argv);//初始化语句，每个程序都一模一样 MPI_Status status; MPI_Datatype onerow; startwtime = MPI_Wtime();//得到开始时间 MPI_Comm_size(MPI_COMM_WORLD,&numprocs);//得到进程数 MPI_Comm_rank(MPI_COMM_WORLD,&myid);//得到进程号 MPI_Get_processor_name(processor_name,&namelen);//得到进程名，其实进程名不是必须的 MPI_Type_contiguous(1,MPI_DOUBLE,&onerow); MPI_Type_commit( &onerow ); for(j=0;j<n+2;j++){ u[j]=0; u[0]=0; u[n+1]=1;}//初始化 te: eps=0; for(j=1;j<n+1;j++){ u0[j]=u[j];} //for(j=0;j<n+2;j++){ //u0[j]=pow(((j-0.5)*dd+myid*n*dd),3); //u0[j]=u[j];} int up=myid-1; if (up<0) {up=MPI_PROC_NULL;} int down=myid+1; if (down>1){ down=MPI_PROC_NULL;} //if(myid<3&&myid>0) { //MPI_Send(&u[1],1,MPI_DOUBLE,myid+1,1000,MPI_COMM_WORLD); ///MPI_Recv(&u[n+1],1,MPI_DOUBLE,myid,1000,MPI_COMM_WORLD,&status); //MPI_Send(&u[n],1,MPI_DOUBLE,myid-1,1000,MPI_COMM_WORLD); //MPI_Recv(&u[0],1,MPI_DOUBLE,myid,1000,MPI_COMM_WORLD,&status);} //if (myid>0) {} if(myid>0&&myid<3){ MPI_Sendrecv(&u[1],1,onerow,down,1000,&u[n+1],1,onerow,myid,1000,MPI_COMM_WORLD,&status);//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! MPI_Sendrecv(&u[n],1,onerow,up,1000,&u[0],1,myid,up,1000,MPI_COMM_WORLD,&status);}//!!!!!!!!!!!!!!!!!!!!!!!!!!!!! for(i=1;i<=n;i++) { double t=(i-0.5)*dd+myid*n*dd; if(i==1){ u=(u[i+1]+2*u[0]-f(t)*dd*dd)/3.0; continue;} if(i==(n)){ u=(2*u[n+1]+u[i-1]-f(t)*dd*dd)/3.0; continue;} u=(u[i-1]+u[i+1]-f(t)*dd*dd)/2; } //if(myid<num-1){u[myid][n+1]=/*u[myid+1][0];}*/(u[myid+1][1]+u[myid][n]-f((myid+1)*n*dd)*dd*dd)/2; //u[myid+1][0]=u[myid][n+1];} //if(myid>0){u[myid][0]=/*u[myid-1][n+1];}*/(u[myid][1]+u[myid-1][n+1]-f(myid*n*dd)*dd*dd)/2; //u[myid-1][n+1]=u[myid][0];} for(j=1;j<n+1;j++){ eps=eps+fabs(u0[j]-u[j]);} if(eps>0.1){goto te;} endwtime = MPI_Wtime();//得到终止时间 printf("wall clock time = %f\n", endwtime-startwtime);//得到墙上时间 MPI_Type_free( &onerow ); MPI_Finalize();//关闭MPI for(j=1;j<n+1;j++){ printf("u[%d]=%f\n",j,u[j]);} return 0; } |
回复此楼» 猜你喜欢
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有80人回复
-
物理学I论文润色/翻译怎么收费?
已经有70人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
40