MPI_Cart_create、MPICHでNULLコミュニケータを取得する

Question

私の宿題（並列プログラミング）用にこのプログラムを書いていますが、実行中にエラーが発生しています。

これはソースです：

/**************************************** 
* Program:     * 
* ALL_TO_ALL_BROADCAST_MESH  * 
*      * 
* Author:     * 
* -------------------   * 
* -------------------   * 
****************************************/ 

/* 
* Program's using MPI_Send, MPI_Recv and cartesian topology functions 
* 
* Compile using: mpicc -o all_to_all_bc all_to_all_bc.c 
*/ 

#include <stdlib.h> 
#include <stdio.h> 
#include <math.h> 
#include <string.h> 
#include "mpi.h" 

/* Number of dimensions */ 
#define NUMB_DIMS 2 

/* Max value for the rand() function */ 
#define MAX_NUMB 128 

#define TRUE 1 
#define FALSE 0 

/* Checks if square root of n can be calculated  * 
* Returns the square root on success or -1 on failure */ 
int is_sqrt(int n); 

int main(int argc, char **argv) { 

    int my_rank;   /* Rank of the current process   */ 
    int size;   /* Number of processes in the mesh  */ 
    int mesh_size;   /* Size of a single mesh row   */ 
    int my_numb;   /* Message     */ 
    int *all_numbs;   /* All numbers owned by the current process */ 
    int numb_numbs;   /* Number of numbers in the all_numbs array */ 
    int left, right, up, down; /* Neighbourghs of the process as rank numbers */ 

    /* Neighbourghs as cartesian coords */ 
    int *my_coords, *left_coords, *right_coords, *up_coords, *down_coords; 

    /* Allocate memory for cartesian coords */ 
    my_coords = (int *)malloc(NUMB_DIMS*sizeof(int)); 
    left_coords = (int *)malloc(NUMB_DIMS*sizeof(int)); 
    right_coords = (int *)malloc(NUMB_DIMS*sizeof(int)); 
    up_coords = (int *)malloc(NUMB_DIMS*sizeof(int)); 
    down_coords = (int *)malloc(NUMB_DIMS*sizeof(int)); 

    MPI_Comm grid_comm; /* Mesh Communicator */ 

    int *dim_size = (int *)malloc(NUMB_DIMS*sizeof(int)); /* Sizes of mesh dimensions  */ 
    int *peroids = (int *)malloc(NUMB_DIMS*sizeof(int)); /* True if the column wraps around */ 
    int reorder;      /* I have no idea what this does ;) */ 

    int i,j; /* counters */ 

    /* MPI initialization */ 
    MPI_Init(&argc,&argv); 
    MPI_Comm_size(MPI_COMM_WORLD,&size); 
    MPI_Comm_rank(MPI_COMM_WORLD,&my_rank); 

    /* Check if p^2=N */ 
    if (mesh_size=is_sqrt(size)<0) { 
     perror("Square root of the number of processes can't be calculated!"); 
     MPI_Finalize(); 
     return -1; 
    } 

    /* Configure parameters for MPI_Cart_Create() */ 
    for (i=0;i<NUMB_DIMS;i++) { 
     dim_size[i]=mesh_size; 
     peroids[i]=TRUE; 
    } 
    reorder=FALSE; 

    /* create new communicator */ 
    if (MPI_Cart_create(MPI_COMM_WORLD,NUMB_DIMS,dim_size,peroids,reorder,&grid_comm)<0) { 
     perror("Unable to create new communicator!"); 
     MPI_Finalize(); 
     return -1; 
    } 

    /* calculate the coords and ranks */ 
    MPI_Cart_coords(grid_comm,my_rank,NUMB_DIMS,my_coords); 

    /* !!this block is for 2D grid only!! */ 
     memcpy(left_coords,my_coords,NUMB_DIMS*sizeof(int)); 
     memcpy(right_coords,my_coords,NUMB_DIMS*sizeof(int)); 
     memcpy(up_coords,my_coords,NUMB_DIMS*sizeof(int)); 
     memcpy(down_coords,my_coords,NUMB_DIMS*sizeof(int)); 

     left_coords[0]=(left_coords[0]-1)%mesh_size; 
     right_coords[0]=(right_coords[0]+1)%mesh_size; 
     up_coords[1]=(up_coords[1]+1)%mesh_size; 
     down_coords[1]=(down_coords[1]-1)%mesh_size; 

     MPI_Cart_rank(grid_comm,left_coords,&left); 
     MPI_Cart_rank(grid_comm,right_coords,&right); 
     MPI_Cart_rank(grid_comm,up_coords,&up); 
     MPI_Cart_rank(grid_comm,down_coords,&down); 
    /* !!this block is for 2D grid only!! */ 

    /* get a random number for this process */ 
    my_numb = rand()%MAX_NUMB; 

    all_numbs = (int *)malloc(size*sizeof(int)); 

    MPI_Status status; 

    /**************** 
    * BROADCASTING * 
    ****************/ 

    printf("I'm process (%d,%d),\nMy number is %d\n\n",my_coords[0],my_coords[1],my_numb); 

    all_numbs[0]=my_numb; 
    numb_numbs=1; 

    /* Communication along rows */ 
    for (i=0;i<mesh_size-1;i++) { 

     printf("I'm process (%d,%d),\nI'm sending %d to process (%d,%d)\n\n",my_coords[0],my_coords[1],my_numb,right_coords[0],right_coords[1]); 
     MPI_Send(&my_numb,sizeof(int),MPI_INT,right,99,grid_comm); 

     printf("I'm process (%d,%d),\nI'm recieving data from process (%d,%d)\n\n",my_coords[0],my_coords[1],left_coords[0],left_coords[1]); 
     MPI_Recv(&my_numb,sizeof(int),MPI_INT,left,99,grid_comm,&status); 

     all_numbs[numb_numbs]=my_numb; 
     numb_numbs++; 

     printf("My current numbers are:\n"); 
     for (j=0;j<numb_numbs;j++) 
      printf("%d, ",all_numbs[j]); 
     printf("\n\n"); 
    } 

    /* Communication along columns */ 
    my_numb=all_numbs[0]; 

    for (i=0;i<mesh_size-1;i++) { 

     printf("I'm process (%d,%d),\nI'm sending %d to process (%d,%d)\n\n",my_coords[0],my_coords[1],my_numb,up_coords[0],up_coords[1]); 
     MPI_Send(&my_numb,sizeof(int),MPI_INT,up,99,grid_comm); 

     printf("I'm process (%d,%d),\nI'm recieving data from process (%d,%d)\n\n",my_coords[0],my_coords[1],down_coords[0],down_coords[1]); 
     MPI_Recv(&my_numb,sizeof(int),MPI_INT,down,99,grid_comm,&status); 

     all_numbs[numb_numbs]=my_numb; 
     numb_numbs++; 

     printf("My current numbers are:\n"); 
     for (j=0;j<numb_numbs;j++) 
      printf("%d, ",all_numbs[j]); 
     printf("\n\n"); 
    } 

    printf("I'm process %d\nBroacasting performed!\nMy numbers are:\n"); 
    for (i=0;i<numb_numbs;i++) 
     printf("%d, ",all_numbs[i]); 
    printf("\n\n"); 

    MPI_Finalize(); 
    return 0; 
} 

int is_sqrt(int n) { 

    double a=n; 
    double b=sqrt(a); 
    double c=b*b; 

    int result; 

    if (a==c) { 
     result=(int)b; 
     return result; 
    } 
    else 
     return -1; 
} 

そして、ここでは、私が取得していますエラーです：

% mpirun -np 4 all_to_all_bc 
[cli_1]: [cli_3]: aborting job: 
Fatal error in MPI_Cart_coords: Invalid communicator, error stack: 
MPI_Cart_coords(130): MPI_Cart_coords(MPI_COMM_NULL, rank=3, maxdims=2, coords=0x6a9010) failed 
MPI_Cart_coords(74).: Null communicator 
aborting job: 
Fatal error in MPI_Cart_coords: Invalid communicator, error stack: 
MPI_Cart_coords(130): MPI_Cart_coords(MPI_COMM_NULL, rank=1, maxdims=2, coords=0x6a9010) failed 
MPI_Cart_coords(74).: Null communicator 
[cli_0]: aborting job: 
Fatal error in MPI_Cart_coords: Invalid communicator, error stack: 
MPI_Cart_coords(130): MPI_Cart_coords(MPI_COMM_NULL, rank=0, maxdims=2, coords=0x6a9010) failed 
MPI_Cart_coords(74).: Null communicator 
[cli_2]: aborting job: 
Fatal error in MPI_Cart_coords: Invalid communicator, error stack: 
MPI_Cart_coords(130): MPI_Cart_coords(MPI_COMM_NULL, rank=2, maxdims=2, coords=0x6a9010) failed 
MPI_Cart_coords(74).: Null communicator 
rank 3 in job 18 ---host---_58157 caused collective abort of all ranks 
    exit status of rank 3: killed by signal 9 
rank 2 in job 18 ---host---_58157 caused collective abort of all ranks 
    exit status of rank 2: return code 1 
rank 1 in job 18 ---host---_58157 caused collective abort of all ranks 
    exit status of rank 1: killed by signal 9 
rank 0 in job 18 ---host---_58157 caused collective abort of all ranks 
    exit status of rank 0: return code 1 

Answer 1

は編集：

はあなたのエラーメッセージを見る：

Fatal error in MPI_Cart_coords: Invalid communicator, error stack: 
MPI_Cart_coords(130): MPI_Cart_coords(MPI_COMM_NULL, rank=3, maxdims=2, coords=0x6a9010) failed 
MPI_Cart_coords(74).: Null communicator 

MPI_Cart_coordsのコミュニケータがnullのようです。スキャンバックすると、これは上のMPI_Cart_create〜5行の呼び出しによって初期化されます。 MPI_Cart_create manページから

（私の強調）：

説明

MPI_Cart_createは デカルトトポロジ情報が付加された新しいコミュニケータへのハンドルを返します。 reorder = falseの場合、 新しいグループの各プロセスのランクは、古いグループのランク と同じです。さもなければ、機能はプロセス を並べ替えるかもしれません（物理マシンに仮想トポロジ のよい埋め込みを選ぶために）。 デカルト・グリッド の合計サイズが通信グループのサイズより小さい場合、MPI_Comm_splitと同様に、一部のプロセス がMPI_COMM_NULLで戻されます。グループ より大きいグリッドが指定されている場合は、 コールに誤りがあります。

これはあなたの問題かもしれません。