Question

我试图将两个矩阵乘以MPI，但我的两个矩阵都是矢量形式。

vector <vector <int> > Mat1;以MxN顺序

vector <vector <int> > Mat2;这是NxL订单

，结果也是矢量形式

vector <vector <int> > Mat3;

这是MxL订单。

我正在使用MPI_BCAST将值广播为

MPI_Bcast(Mat2.data(), Mat2.size(), MPI_INT, 0, MPI_COMM_WORLD);

并使用分散和收集功能来接收数据

    MPI_Scatterv(Mat1.data(), CNTS, SNTS, MPI_INT, MatPartS.data(), MatPartS.size(), MPI_INT, 0, MPI_COMM_WORLD);


    for (int i = 0; i < myRowsSize; ++i)
    {
            for (int j = 0; j < L; ++j)
            {
                    ResultMatPart[i][j]=0;
                    for (int k = 0; k < N; ++k)
                    {
                             ResultMatPart[i][j] += MatPartS[i][k] * Mat2[k][j];
                    }
            }
    }               


    MPI_Gatherv(ResultMatPart.data(), ResultMatPart.size(), MPI_INT, ResultMat.data(), RCNTS, SCNTS, MPI_INT, 0, MPI_COMM_WORLD);

块引用

我可以编译程序但它崩溃并且不会产生任何输出。

它打印的信息是：

Fatal error in PMPI_Scatterv: Message truncated, error stack:
PMPI_Scatterv(376).....: MPI_Scatterv(sbuf=0x108fe00, scnts=0x108fef0, displs=0x108ffc0, MPI_INT, rbuf=0x1090270, rcount=4, MPI_INT, root=0, MPI_COMM_WORLD) failed
MPIR_Scatterv_impl(187): 
MPIR_Scatterv(106).....: 
MPIR_Localcopy(340)....: Message truncated; 64 bytes received but buffer size is 16

如何在MPI中传递多维STL向量的任何帮助？

非常感谢

取值

Answer 1

那些MPI函数期待一维数组。您需要将2D矢量收集到单个矢量中。这只是使用向量矢量来表示矩阵的原因之一。

看看我给别人的矢量矢量困境的答案。它使用单个向量提供简单的矩阵类：

Get the first column of a matrix represented by a vector of vectors

Answer 2

您始终可以使用行主顺序（或列主顺序）来存储矩阵（http://en.wikipedia.org/wiki/Row-major_order）。由于这是矩阵的一维表示，您可以使用散点图将其发送出去。

如果你真的需要一个双维数组，你应该看MPI_PACK MPI_UNPACK或MPI派生类型。

Answer 3

vector <vector <int> > Mat1;以MxN顺序

vector <vector <int> > Mat2;这是NxL订单

这甚至不能确保行的长度确实相等。

包含在具有重载[] -op的类中的平坦M * N大小的字段（例如，简单real_type*或std::valarray<real_type>）是一个更加健壮的模型，也非常适合对于MPI。

如何在MPI中使用多维STL向量

3 个答案: