MPI_Type_create_subarray и MPI_Send - PullRequest
Новая
       40

MPI_Type_create_subarray и MPI_Send

3 голосов
/ 07 апреля 2011

это мой первый вопрос здесь в stackoverflow.У меня есть два процесса, корень 0 и ведомый 1. Подчиненный выделяет 2d массив (CHUNK_ROWS + 2, CHUNK_COLUMNS + 2) и хочет отправить подмассив [CHUNK_ROWS] [CHUNK_COLUMNS].Root выделяет 2d массив (ROWS, COLUMNS) и получает подмассив, который хочет сохранить, из ddd [0] [0] и печатает его.В результате я получаю это неправильно .. Почему?Я знаю, что это бессмысленный код, но это всего лишь простая программа, которая должна помочь мне в более сложном проекте.Вот код: 

#include <mpi.h>
#include <iostream>
using namespace std;

#define ROWS 10
#define COLUMNS 10
#define CHUNK_ROWS 5
#define CHUNK_COLUMNS 5
#define TAG 0

int** alloca_matrice(int righe, int colonne)
{
int** matrice=NULL;
int i;

// per allocare la matrice devo fare in questo modo al fine di avere le righe contigue in memoria
// e poter così utilizzare il tipo 'colonna' che definisco con MPI_Type_Vector()
matrice = (int **)malloc(righe * sizeof(int*));

if(matrice != NULL){
    matrice[0] = (int *)malloc(righe*colonne*sizeof(int));
    if(matrice[0]!=NULL)
        for(i=1; i<righe; i++)
            matrice[i] = matrice[0]+i*colonne;
    else{
        free(matrice);
        matrice = NULL;
    }
}
else{
    matrice = NULL;
}
return matrice;
}


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

int my_id, numprocs,length,i,j;
int ndims, sizes[2],subsizes[2],starts[2];
int** DEBUG_CH;
int** ddd;
char name[BUFSIZ];
MPI_Datatype subarray;
MPI_Status status;
MPI_Init(&argc, &argv) ;                    // Chiamata obbligatoria di inizializzazione    
MPI_Comm_rank(MPI_COMM_WORLD, &my_id) ;     // Ottiene l'identificativo del processo
MPI_Comm_size(MPI_COMM_WORLD, &numprocs) ;  // Ottiene quanti processi sono attivi
MPI_Get_processor_name(name, &length);      // Il nome del nodo dove il processo ? in esecuzione



if(my_id==1){
    //creo una sottomatrice ripulita dalle ghost cells
    ndims=2;
    sizes[0] = CHUNK_ROWS+2;
    sizes[1] = CHUNK_COLUMNS+2;
    subsizes[0] = CHUNK_ROWS;
    subsizes[1] = CHUNK_COLUMNS;
    starts[0] = 1;
    starts[1] = 1;
    MPI_Type_create_subarray(ndims,sizes,subsizes,starts,MPI_ORDER_C,MPI_INT,&subarray);
    MPI_Type_commit(&subarray);

    DEBUG_CH = alloca_matrice(CHUNK_ROWS+2,CHUNK_COLUMNS+2);
    for(i=0; i<CHUNK_ROWS+2; i++){
        for(j=0; j<CHUNK_COLUMNS+2; j++){
            if(i==0 || i==CHUNK_ROWS+1 || j==0 || j==CHUNK_COLUMNS+1)
                DEBUG_CH[i][j] = 5;
            else
                DEBUG_CH[i][j] = 1;
        }
    }
    MPI_Send(DEBUG_CH,1,subarray,0,TAG,MPI_COMM_WORLD);
    MPI_Type_free(&subarray);
}
if(my_id==0){
    //creo una sottomatrice ripulita dalle ghost cells
    ndims=2;
    sizes[0] = ROWS;
    sizes[1] = COLUMNS;
    subsizes[0] = CHUNK_ROWS;
    subsizes[1] = CHUNK_COLUMNS;
    starts[0] = 0;
    starts[1] = 0;
    MPI_Type_create_subarray(ndims,sizes,subsizes,starts,MPI_ORDER_C,MPI_INT,&subarray);
    MPI_Type_commit(&subarray);

    ddd = alloca_matrice(ROWS,COLUMNS);
    MPI_Recv(ddd[0],1,subarray,1,TAG,MPI_COMM_WORLD,&status);
    MPI_Type_free(&subarray);
    for(i=0; i<CHUNK_ROWS; i++){
        for(j=0; j<CHUNK_COLUMNS; j++){
            printf("%d ",ddd[i][j]);
        }
        printf("\n");
    }
}
MPI_Finalize();                             // Chiusura di MPI.
return 0;
}

Заранее спасибо.

1 Ответ

2 голосов
/ 07 апреля 2011

Поздравляем, ваш код почти идеален, в MPI_Send есть только одна глупая ошибка, которую вы получили прямо в MPI_Recv.

Для отправки у вас есть 

MPI_Send(DEBUG_CH,1,subarray,0,TAG,MPI_COMM_WORLD);

, тогда как для Recv у вас есть 

MPI_Recv(ddd[0],1,subarray,1,TAG,MPI_COMM_WORLD,&status);

Второе верно. Вам необходимо отправить указатель на то место, откуда начинаются данные; это не DEBUG_CH, который является указателем на указатель на целые, а & (DEBUG_CH [0] [0]) или, что то же самое, DEBUG_CH [0], как вы это сделали с ddd.

...