矩陣向量乘法在MPI和C

我想通過使用MPI和C的向量乘以方陣。我必須使用MPI_Allgather將矩陣的所有部分發送到所有進程。這是我到目前爲止矩陣向量乘法在MPI和C

#include "mpi.h" 
#include <stdlib.h> 
#include <stdio.h> 
#include <math.h> 
#include <string.h> 
#include <time.h> 
#define DIM 500 

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

     int i, j, n; 
     int nlocal;  /* Number of locally stored rows of A */ 
     double *fb; 
     double a[DIM*DIM], b[DIM], x[DIM];  /* Will point to a buffer that stores the entire vector b */ 
     int npes, myrank; 
     MPI_Status status; 

     MPI_Init(&argc,&argv); 

    /* Get information about the communicator */ 
    MPI_Comm_rank(MPI_COMM_WORLD, &myrank); 
    MPI_Comm_size(MPI_COMM_WORLD, &npes); 

    /* Allocate the memory that will store the entire vector b */ 
    fb = (double*)malloc(n*sizeof(double)); 

    nlocal = n/npes; 

    /* Gather the entire vector b on each processor using MPI's ALLGATHER operation */ 
    MPI_Allgather(b, nlocal, MPI_DOUBLE, fb, nlocal, MPI_DOUBLE, MPI_COMM_WORLD); 

    /* Perform the matrix-vector multiplication involving the locally stored submatrix */ 
    for (i=0; i<nlocal; i++) { 
     x[i] = 0.0; 
     for (j=0; j<n; j++) 
     x[i] += a[i*n+j]*fb[j]; 
    } 


    free(fb); 

    MPI_Finalize(); 
}//end main

OK現在，它的工作原理！它編譯，但然後我得到一個mpi allgather內部錯誤！我也嘗試了其他解決方案。

Fatal error in MPI_Allgather: Internal MPI error!, error stack: 
MPI_Allgather(961).......: MPI_Allgather(sbuf=0xa1828, scount=407275437, MPI_DOU                            BLE, rbuf=0xf61d0008, rcount=407275437, MPI_DOUBLE, MPI_COMM_WORLD) failed 
MPIR_Allgather_impl(807).: 
MPIR_Allgather(766)......: 
MPIR_Allgather_intra(560): 
MPIR_Localcopy(357)......: memcpy arguments alias each other, dst=0xb8513d70 src                            =0xa1828 len=-1036763800 
Fatal error in MPI_Allgather: Internal MPI error!, error stack: 
MPI_Allgather(961).......: MPI_Allgather(sbuf=0xa1828, scount=407275437, MPI_DOU                            BLE, rbuf=0xf61d0008, rcount=407275437, MPI_DOUBLE, MPI_COMM_WORLD) failed 
MPIR_Allgather_impl(807).: 
MPIR_Allgather(766)......: 
MPIR_Allgather_intra(560): 
MPIR_Localcopy(357)......: memcpy arguments alias each other, dst=0x3cb9b840 src                            =0xa1828 len=-1036763800 
Fatal error in MPI_Allgather: Internal MPI error!, error stack: 
MPI_Allgather(961).......: MPI_Allgather(sbuf=0xa1828, scount=407275437, MPI_DOU                            BLE, rbuf=0xf61d0008, rcount=407275437, MPI_DOUBLE, MPI_COMM_WORLD) failed 
MPIR_Allgather_impl(807).: 
MPIR_Allgather(766)......: 
MPIR_Allgather_intra(560): 
MPIR_Localcopy(357)......: memcpy arguments alias each other, dst=0x7a857ad8 src                            =0xa1828 len=-1036763800

任何人都可以幫我嗎？

來源

2012-12-06 Mário Santana

'a'是一個二維數組，但您正在訪問它，因爲它只有一個維度。語法是'a [??] [??]'而不是'a [??]'，但是你放置什麼取決於你的數據佈局。 – user786653

我喜歡一維，因爲它在不同的過程中通過行分解。 –

我試過，但仍然給我那個錯誤。不知道如何打開這個圓環 –

任何方式要訪問使用一個維（行或列主要），可以分配二維陣列作爲a[DIM*DIM]代替a[DIM][DIM]和訪問此以線性方式的二維陣列。此解決方案正在爲我工作。

來源

2012-12-08 09:27:34 Srini

它工作！謝謝！但現在我得到一個內部的mpi錯誤 –

矩陣向量乘法在MPI和C

回答

相關問題