MPI_File_write C - PullRequest
Новая
       27

MPI_File_write C

0 голосов
/ 14 ноября 2018

Я должен использовать MPI для вычисления значения концентрации из трехмерной матрицы.

Я делю свою матрицу на кусочки в зависимости от значения числа процедур, которые я буду использовать.

Я получаю трехмерную матрицу меньшего размера, которая работает независимо, но я хочу, чтобы вычисляемые ими значения концентрации записывались в один и тот же файл, например строки от 0 до n-1, записывались proc 0, n - 2n-1 proc 1, ..

Я знаю, что могу использовать смещение в функции mpi_file_write_at, чтобы выполнить такую ​​вещь, но я не знаю, как это реализовать.

Я знаю, что простым решением было бы попросить proc 0 открыть файл, написать и закрыть его, а затем proc 1. открывается.

Не обращайте внимания на сам код, так как он еще не закончен 

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

int main(int argc, char *argv[])
{
    char *name = argv[1];

    double h, m, L, T_max, v_x, v_y, v_z, D, r_thresh;
    unsigned int S;

    FILE *file = fopen(name, "rb");


    if(!file)
    {
        fprintf(stderr, "Not able to open the file %s\n", name);
    }
    else
        printf("Able to open the file %s\n", name);

    fscanf(file, "%lf", &h);
    fscanf(file, "%lf", &m);
    fscanf(file, "%lf", &L);
    fscanf(file, "%lf", &T_max);
    fscanf(file, "%lf", &v_x);
    fscanf(file, "%lf", &v_y);
    fscanf(file, "%lf", &v_z);
    fscanf(file, "%lf", &D);
    fscanf(file, "%u", &S);
    fscanf(file, "%lf", &r_thresh);

    fclose(file);

    if(((L/h) % 2) != 0){
        printf("Uncorrect parameters L and h ! \n");
        return -1;
    }
    else{
        unsigned int edge = (L/h);
        unsigned int middle = edge/2;
        unsigned int N = edge+1;
    }
    if(((T_max/m) % 1) != 0){
        printf("Uncorrect parameters T_max and m ! \n");
        return -1;
    }
    else{
        unsigned int t_end = T_max/m;
    }

    unsigned int myrank, nbproc;
    MPI_Init(&argc, &argv) ;
    MPI_Comm_rank( MPI_COMM_WORLD , &myrank );
    MPI_Comm_size( MPI_COMM_WORLD , &nbproc );

    double w = N/nbproc;
    unsigned int width = round(w);
    unsigned int k_min, k_max;
    k_min = myrank*width;
    k_max = (myrank+1)*width-1;
    width = width+2; 
    if(nbproc==1){
        width = N;
    }
    else if(myrank==(nbproc-1)){
        k_max = N-1; 
        width = k_max-k_min+1; 
        width = width+1; 
    }
    else if(myrank==0){
        width = width-1; 
    }



    double ***C ;
    C = malloc(width * sizeof(*C)); 
    for(k=0 ; k < width ; k++){
        C[k] = malloc(N * sizeof(**C));
    }
    for(k=0 ; k < width ; k++){
        for(j=0 ; j < N ; j++){
            C[k][j] = malloc(N * sizeof(***C));
        }
    }


    unsigned int i,j,k,n;
    for(k=0 ; k<width ; k++){
        for(j=0 ; j<N ; j++){
            for(i=0 ; i<N ; i++){
                C[k][j][i]=0;
            }
        }
    }
    if(middle>=k_min && middle=<k_max){
        C[middle-k_min][middle][middle] = 1;
    }


    if(myrank==0 || nbproc==1){
        k_min = k_min +1 ;
    }
    if(myrank==(nbproc-1) || nbproc==1){
        k_max = k_max -1 ;
    }





    for(n=1 ; n<T_max ; n++){

        char *c = "c_";
        char buffer[10];
        itoa(n,buffer,2);
        strcat(c,buffer);

        // I want to write some value stored in a buffer now
        // computed by different process in the same file 
        // MPI_File_write_at ???


    }

    for(k=0 ; k < width ; k++){
        for(j=0 ; j < N ; j++){
            free(C[k][j]);
            }
        free(C[k]);
    }
    free(C);

    return 0;
}

1 Ответ

0 голосов
/ 30 ноября 2018

Каждый процесс в программе MPI имеет локальное адресное пространство.Файловый ввод / вывод изменяется следующим образом: файл является глобально доступным линейным потоком байтов.Вам нужно подумать о том, где в глобальном пространстве пишет каждый процесс.Пара способов сделать это:

  • иногда вы можете просто записать каждый процесс в соответствующий фрагмент файла.ранг ноль записывает в первый мегабайт, ранг 1 во 2 и т. д.Возможно, это слишком упрощенно для вашего трехмерного массива.

  • Для более сложной декомпозиции вы устанавливаете «представление файла»: вы создаете тип данных MPI, описывающий расположение данных на диске.в вашем случае MPI_TYPE_CREATE_SUBARRAY - это, вероятно, то, что вы хотите: глобальная трехмерная матрица, и каждый процесс вносит некоторый вклад в общую картину.

...