c - MPI_Reduce() 程序出现段错误

Question

我开始研究使用 C 和 OpenMPI 库进行并行编程。所以现在一切看起来都太复杂了。

我正在尝试做一个单程序多数据:掌握:- 初始化一个数组- segmentation 它- 将相同大小的位发送到不同的进程(奴隶)奴隶:- 进程更改数组值并计算所有新元素的总和- 将更改后的数组发送回 master掌握:- 进行集体交流以收集和总结新值(value)的总和- 打印接收到的每个新数组的前五个元素和新值的全局总和。- 打印花费的时间。

这是我写的:

#include <stdio.h>
#include <stdlib.h>
#include "mpi.h"
#define ARRAYSIZE   16000000

int main (int argc, char *argv[]) {
    MPI_Init(&argc, &argv);
    int myrank;
    char name[100];
    int result;
    int size = 0;
    int number;
    MPI_Comm_rank(MPI_COMM_WORLD, &myrank);
    MPI_Request request;
    int buffer;
    int count;
    double t1, t2;
    MPI_Comm_size (MPI_COMM_WORLD,&size);
    int cells = ARRAYSIZE/size;
    float array[cells];

    if (myrank == 0) {
        t1 = MPI_Wtime();
        t2 = MPI_Wtime();
        MPI_Get_processor_name(name, &result);

        // ********************** INICIALIZANDO ARRAY **********************
        int     i;          /* loop variable */
        float   data[ARRAYSIZE];    /* the intial array */

        printf("Starting serial array example...\n");
        printf("Using array of %d floats. Requires %ld bytes\n",ARRAYSIZE,sizeof(data));

        /* Initialize the array */
        printf("Initializing array...\n");
        for(i=0; i<ARRAYSIZE; i++)
          data[i] =  i * 1.0;

        /* Print a few sample results */
        printf("Sample results\n");
        printf("   data[1]=%e\n",  data[1]);
        printf("   data[100]=%e\n",  data[100]);
        printf("   data[1000]=%e\n",  data[1000]);
        printf("   data[10000]=%e\n",  data[10000]);
        printf("   data[100000]=%e\n",  data[100000]);
        printf("   data[1000000]=%e\n",  data[1000000]);
        printf("\nAll Done!\n");
        // ********************** ARRAY INICIALIZADO **********************

        MPI_Comm_size (MPI_COMM_WORLD,&size);
        printf("Total of tasks: %d", size);
        printf("Each task process %d cells", ARRAYSIZE/size);
        int cells = ARRAYSIZE/size;
        int id_task;
        for(id_task = 0; id_task < size; id_task++) {
            //float array[cells];
            int i=0;
            for(i=0; i<cells; i++)
                array[i] =  i * (id_task+1.0);
            MPI_Send(&array[id_task*cells], cells, MPI_FLOAT, id_task, 0, MPI_COMM_WORLD);
        }

        printf("master: %d at processor: %s\n",myrank, name);

    }
    MPI_Recv(array, cells, MPI_FLOAT, 0, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
    int i;
    float resultado;
    for(i=0; i<cells; i++)
        array[i] =  i * (myrank+1.0);

    if(myrank!=0){
        MPI_Send(array, cells, MPI_FLOAT, 0, 0, MPI_COMM_WORLD);
    }

    MPI_Reduce(&array, &resultado, 1, MPI_FLOAT, MPI_SUM, 0, MPI_COMM_WORLD);

    if (myrank == 0) {
        int j;
        for(j=0; j<cells; j++){
            for(i=0; i<5; i++){
                printf("| %lf ",array[i*j]);
            }
        }
        //t1 = MPI_Wtime();
        t2 = MPI_Wtime();
        MPI_Get_processor_name(name, &result);
        printf("master: %d at processor: %s\ntime: %lf\n",myrank, name,t2-t1);
    }
    MPI_Finalize();

}

但是当我运行它时，我得到的只是一个“Segmentation Fault Error”。我已阅读段错误问题，但无法诊断为什么我的代码会出现这种情况。

提前致谢。(我很抱歉我写得不好，英语不是我的母语)

更新:我包含了一个 malloc 和一个 free，但在运行时我仍然有“mpirun 注意到节点 Brunos-Air 上 PID 0 的进程等级 0 在信号 11 上退出(段错误:11)。"

如果我尝试“$HOME/opt/usr/local/bin/mpicc -o -Wall -Wextra -g programa3-r1 ./programa3-r1.c”，它给出:“ld:无法与架构 x86_64 的主可执行文件“programa3-r1”链接 clang:错误:链接器命令失败，退出代码为 1(使用 -v 查看调用) "

Answer 1

for(id_task = 0; id_task < size; id_task++) {
    //float array[cells];
    int i=0;
    for(i=0; i<cells; i++)
        array[i] =  i * (id_task+1.0);
    MPI_Send(&array[id_task*cells], cells, MPI_FLOAT, id_task, 0, MPI_COMM_WORLD);
}

你应该MPI_Send() array (aka &array[0]) 而不是&array[id_task*cells]。这可能是崩溃的根本原因。

请注意，从 MPI 的角度来看，您的程序是不正确的，因为任务 0 MPI_Send() 到它自己，然后 MPI_Recv() 从它自己。对于小 消息，这可能会很好地工作，而对于长 消息，则可能会出现死锁。 小 和长 取决于(除其他事项外)您的 MPI 库，所使用的互连，因此您真的必须避免这种情况。

if(myrank!=0){
    MPI_Send(array, cells, MPI_FLOAT, 0, 0, MPI_COMM_WORLD);
}

在等级 0 上没有匹配的 MPI_Recv()，所以这也可能是死锁。

MPI_Reduce(&array, &resultado, 1, MPI_FLOAT, MPI_SUM, 0, MPI_COMM_WORLD);

这对所有等级的 array[0] 求和，并将结果存储在等级 0 的 resultado 中。我不确定这是否真的是您的想法(但我不知道你想要达到什么目的，所以这可能没问题。