xingzhe

//
//进程之间相互问候

// void Hello(void);
//
// int main(int argc, char* argv[])
// {
//  int me, option, namelen, size;
//
//  char processor_name[MPI_MAX_PROCESSOR_NAME];
//  MPI_Init(&argc, &argv);
//  MPI_Comm_rank(MPI_COMM_WORLD, &me);
//  MPI_Comm_size(MPI_COMM_WORLD, &size);
//
//  if (size < 2){
//   //如果总的进程数小于2， 则出错推出
//   fprintf(stderr, "system require at least 2 processes.");
//   MPI_Abort(MPI_COMM_WORLD, 1);
//  }
//
//  //get current name
//     MPI_Get_processor_name(processor_name, &namelen);
//
//  fprintf(stderr, "Process &d is alive on &s\n", me, processor_name);
//
//  //sync
//  MPI_Barrier(MPI_COMM_WORLD);
//
//  //greet each other
//  Hello();
//
//  //
//  MPI_Finalize();
//
//  return 0;
// }
//
// void Hello(void){
//  
//  int nproc, me;
//  int type =1;
//  int buffer[2], node;
//  MPI_Status status;
//  MPI_Comm_rank(MPI_COMM_WORLD, &me);
//  MPI_Comm_size(MPI_COMM_WORLD, &nproc);
//  
//  if (me == 0){
//   //进程0负责打印提示信息
//   printf("\nHello test from all to all\n");
//   fflush(stdout);
//  }
//  
//  //循环对每个进程进行问候
//  for (node = 0; node < nproc; node++){ 
//   if (node != me){
//    buffer[0]= me;
//    buffer[1]= node;//greet node
//    
//    //send information to node
//    MPI_Send(buffer, 2, MPI_INT, node, type, MPI_COMM_WORLD);
//    
//    //receive information from node
//    //question  1. blocked or not blocked
//             //question  2. when MPI_Recv returns, the information must be sent from node！
//    MPI_Recv(buffer, 2, MPI_INT, node, type, MPI_COMM_WORLD, &status);
//    
//    //error if receiving information is not sent to me or
//    //????
//    if ((buffer[0] != node )||(buffer[1] != me)){
//     (void)fprintf(stderr, "Hello: %d!=%d or %d!=d%\n", buffer[0], node, buffer[1], me);
//    }else{
//     //receiving informaton
//     printf("Hello from %d  to %d\n", me, node);
//    }
//    
//    
//    
//    fflush(stdout);
//   }
//  }
//  
// }

//接收其他进程的数据
// int main(int argc, char** argv){
//  int rank;
//  int size;
//  int i;
//  int buf[1];
//
//  MPI_Status status;
//  MPI_Init(&argc, &argv);
//  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
//  MPI_Comm_size(MPI_COMM_WORLD, &size);
//
//
//  if (rank == 0){
//   for (i = 0; i < 100*(size-1); i++){
//    MPI_Recv(buf, 1, MPI_INT, MPI_ANY_SOURCE, MPI_ANY_TAG, MPI_COMM_WORLD, &status);
//    printf("Msg=%d from%d with tag%d\n", buf[0],status.MPI_SOURCE, status.MPI_TAG);
//   }
//  }else{
//   for (i=0; i < 100; i++){
//    buf[0]=rank+i;
//    MPI_Send(buf, 1, MPI_INT, 0, i, MPI_COMM_WORLD);
//   
//   }
//  }
//
//  MPI_Finalize();
//}

//对等模式的MPI程序设计
// const int totalsize = 16;
// const int mysize = totalsize/4;
// const int steps = 10;
// int aa[2][3]={{1,1,1},{3,3,3}};
// int a[totalsize][mysize + 2];
// int b[totalsize][mysize + 2];
//
// int n, myid, numprocs, i, j, rc;
//
// int main(int argc, char** argv){
//  int begin_col;
//  int end_col;
//  int ierr;
//  
//  MPI_Status status;
//  
//  MPI_Init(&argc, &argv);
//  MPI_Comm_rank(MPI_COMM_WORLD, &myid);
//  MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
//  
//  //数组初始化
//  for (j = 1; j < mysize+2; j++){
//   for (i=1; i < totalsize; i++){
//    a[i][j] = 0.0;
//   }
//  }
//  
//  //
//  if (myid == 0){
//   for (i = 1; i < totalsize; i++){
//    a[i][2] = 8.0;
//   }
//  }
//  
//  //
//     if (myid == 3){
//   for (i = 1; i < totalsize; i++){
//    a[i][mysize + 1] = 8.0;
//   }
//  }
//  
//  //
//  for (i = 1; i < mysize; i++){
//   a[1][i] = 8.0;
//   a[totalsize][i] = 8.0;
//  }
//
//
//  //jacobi迭代部分
//  for (n = 1; n < steps; n++){
//   //从右侧的邻居得到数据
//   if (myid < 3){
//    MPI_Recv(&a[1][mysize + 2],totalsize, MPI_INT, myid + 1, 10, MPI_COMM_WORLD, &status);
//   }
//
//   //向左侧的邻居发送数据
//   if (myid > 0){
//    MPI_Send(&a[1][2], totalsize, MPI_INT, myid, 10, MPI_COMM_WORLD);
//   }
//
//   //向右侧邻居发送数据
//   if (myid < 3){
//    MPI_Send(&a[1][mysize+1], totalsize, MPI_INT, myid + 1, 10, MPI_COMM_WORLD);
//   }
//
//   //向左侧邻居接收数据
//   if (myid > 0)
//   {
//    MPI_Recv(&a[1][1], totalsize, MPI_INT, myid - 1, 10, MPI_COMM_WORLD, &status);
//   }
//
//   //
//   begin_col = 2;
//   end_col = mysize + 1;
//   //
//   if (myid == 0){
//    begin_col = 3;
//   }
//         //
//   if (myid == 3){
//    end_col = mysize;
//   }
//   
//   //
//         for(j = begin_col; j < end_col; j++){
//    for (i = 2; i < totalsize - 1; i++){
//     b[i][j] = (a[i][j+1] + a[i][j-1] + a[i+1][j] + a[i-1][j])*0.25;
//    }
//   }
//
//   for(j = begin_col; j < end_col; j++){
//    for (i = 2; totalsize-1; i++){
//     a[i][j]= b[i][j];
//
//    }
//   }
//  }
//
//  for (i = 2; i < totalsize; i++){
//   //printf()
//  }
//
//  MPI_Finalize();
//  return 0;
//}

主从模式的MPI程序设计
void master_io();
void slave_io();
#define MSG_EXIT 1
#define MSG_PRINT_ORDERED 2
#define MSG_PRINT_UNORDERED 3

int main(int argc, char** argv){
 int rank, size;
 MPI_Init(&argc, &argv);
 MPI_Comm_rank(MPI_COMM_WORLD, &rank);
 if (rank == 0)
 {
  master_io();
 }else{
  slave_io();
 }

 MPI_Finalize();
}

void master_io(){
 int i, j, size, nslave, firstmsg;
 char  buf[256], buf2[256];
 MPI_Status status;
 MPI_Comm_size(MPI_COMM_WORLD, &size);
 nslave = size -1;
 while (nslave > 0){
  //从任意进程接收任意标识的消息
  MPI_Recv(buf, 256, MPI_CHAR, MPI_ANY_SOURCE, MPI_ANY_TAG,MPI_COMM_WORLD, &status);
  switch(status.MPI_TAG)
  {
  case MSG_EXIT: nslave--; break;
  case MSG_PRINT_UNORDERED:fputs(buf, stdout);break;
  case MSG_PRINT_ORDERED:
   firstmsg = status.MPI_SOURCE;
   for (i=1; i< size; i++){
    if (i == firstmsg){
     fputs(buf, stdout);
    }else{
     //指定源进程和消息标识
     MPI_Recv(buf2, 256, MPI_CHAR, i, MSG_PRINT_ORDERED, MPI_COMM_WORLD, &status);
     fputs(buf2, stdout);
    }     
   }
   break;
  }
 }
}

void slave_io(){
 char buf[256];
 int rank;
 MPI_Comm_rank(MPI_COMM_WORLD, & rank);
 sprintf(buf, "Hello from slave&d, ordered print\n", rank);

 //send ordered msg
 MPI_Send(buf, strlen(buf)+1, MPI_CHAR, 0, MSG_PRINT_ORDERED, MPI_COMM_WORLD);

 sprintf(buf, "Goodbye from slave %d, ordered print\n", rank);
 //send unordered msg
 MPI_Send(buf, strlen(buf)+1, MPI_CHAR, 0, MSG_PRINT_UNORDERED, MPI_COMM_WORLD);

 //send exit msg
 MPI_Send(buf, 0, MPI_CHAR, 0, MSG_EXIT, MPI_COMM_WORLD);

}