Question

这是我想要将2D数组与矢量数组相乘的代码：

#include<iostream>
#include<mpi.h>
using namespace std;
int v_array[10] ;
int ch_size, start, close;
int res ;
int rows, cols;

int main(int argc, char *argv[])
{
    int pro_id, tot_pros;
    MPI_Init(&argc, &argv); 
    MPI_Comm_rank(MPI_COMM_WORLD, &pro_id); 
    MPI_Comm_size(MPI_COMM_WORLD, &tot_pros); 

    if (pro_id == 0) {
        cout << "Enter rows and columns: ";
        cin >> rows >> cols;

        int **array = new int*[rows];
        int size1 = rows * cols;
        array[0] = new int[size1];
        for (int j = 1; j < rows; j++) {
            array[j] = &array[0][j*cols];
        }
        for (int i = 0; i < rows; i++) {
            v_array[i] = 1;
            for (int j = 0; j < cols; j++) {
                array[i][j] = 1;
            }
        }
        for (int i = 1; i < tot_pros; i++) {
            ch_size = (rows / (tot_pros - 1));
            start = (i - 1) * ch_size;
            if (((i + 1) == tot_pros) && ((rows % (tot_pros - 1)) != 0)) {
                close = rows;
            }
            else {
                close = start + ch_size;
            }
            MPI_Send(&start, 1, MPI_INT, i, 1, MPI_COMM_WORLD);
            MPI_Send(&close, 1, MPI_INT, i, 2, MPI_COMM_WORLD);
            MPI_Send(&cols, 1, MPI_INT, i, 4, MPI_COMM_WORLD);
            MPI_Send(&array[start][0], ch_size *cols, MPI_INT, i, 3, MPI_COMM_WORLD);
        }
    }
    else
     {
        int cols;
        MPI_Recv(&start, 1, MPI_INT, 0, 1, MPI_COMM_WORLD, MPI_STATUSES_IGNORE);
        MPI_Recv(&close, 1, MPI_INT, 0, 2, MPI_COMM_WORLD, MPI_STATUSES_IGNORE);
        MPI_Recv(&cols, 1, MPI_INT, 0, 4, MPI_COMM_WORLD, MPI_STATUSES_IGNORE);
        int **array = new int*[(close - start)*cols];
        MPI_Recv(array, (close - start) *cols , MPI_INT, 0, 3, MPI_COMM_WORLD, MPI_STATUSES_IGNORE);
        for (int i = start; i < close; i++) {
            for (int j = 0; j < cols; j++) {
                cout << array[i]<<array[j];
                    res += array[i][j] * v_array[i];
                    cout << res;
                }
            }
         }
    MPI_Finalize(); 
    return 0;
}

当我有静态数组时，同样的程序工作正常，但是动态我得到了这个错误。

E：\ MS（CS）\ 2nd Semester \ Parallel   编程\ program \ arr_multi \ Debug＆gt; mpiexec -n 4 arr_multi.exe输入   行和列：3 2

工作中止：[排名]消息

[0-1]终止

[2]进程退出而没有调用finalize

[3]终止了

----错误分析-----
     RAMISHA-PC上的
[2] arr_multi.exe过早结束，可能有   坠毁。退出代码0xc0000005

----错误分析-----

我声明了一个具有连续位置的数组，我的行在进程之间正确划分。我认为我的数据结构存在问题并尝试了许多解决方案但是徒劳无功。

Answer 1

首先，是调试MPI应用程序的方法，这应该是您的首要任务。多进程应用程序的一般方法是在应用程序开始时暂停，例如使用getchar()然后使用调试程序附加到每个进程here：

编译，链接并开始运行您的MPI程序（您可能希望在执行后续步骤时尽早发布read语句来保存程序。

附加到当前正在运行的MPI进程之一：调试 - 附加到进程会弹出一个列出可用进程的对话框。您应该会看到可执行文件的NUM个实例（N来自mpiexec -n NUM）。选择所有这些并单击 Attach 。您现在可以通过添加断点等进行调试。要在MPI进程之间移动，请使用代码清单上方的Process下拉菜单。

话虽如此，至少有一个问题出在这个部分：int **array = new int*[(close - start)*cols];（在应用程序的接收部分）。您分配第一个维度而不是第二个维度，因此第一个维度中的所有指针都是未初始化的。

将其更改为：

        int *array = new int[(close - start) * cols];
        MPI_Recv(array, (close - start) *cols, MPI_INT, 0, 3, MPI_COMM_WORLD, MPI_STATUSES_IGNORE);
        for (int i = start; i < close; i++) {
            for (int j = 0; j < cols; j++) {
                cout << array[(i - start) * cols];
                res += array[(i - start) * cols] * v_array[i];
                cout << res;
            }
        }
        delete[] array;

或者，如果您真的想要使用2D数组，请从发送部分复制初始化代码：

        int rows = close - start;
        int **array = new int*[rows];
        int size1 = rows * cols;
        array[0] = new int[size1];
        for (int j = 1; j < rows; j++) {
            array[j] = &array[0][j*cols];
        }

第二个问题是v_array作为全局变量未在接收器进程中初始化。请记住，在MPI中，每个过程都是一个独立的程序。因此，您应该始终初始化v_array，即不管pro_id。

        for (int i = 0; i < rows; i++) {
            v_array[i] = 1;
        }

程序提前结束，可能已经崩溃。退出代码0xc0000005

1 个答案: