1

我试图在 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 向量？

非常感谢

小号

4

3 回答 3

1

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

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

获取由向量向量表示的矩阵的第一列

于 2013-05-22T23:39:13.833 回答
0

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

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

于 2013-05-23T06:05:10.307 回答
0

vector <vector <int> > Mat1;哪个按 MxN 顺序

vector <vector <int> > Mat2;哪个是 NxL 订单

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

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

于 2013-05-23T19:08:06.130 回答