简短的MPI程序（C代码），root中的声明变量不工作

Question

因此，我对我的小MPI程序有一个疑问。 （我是编程的新手，所以它可能只是一些初学者的错误）

我制作了一个程序，该程序通过使用公式来计算PI：（1/n） * 4/（1+（1+（i/n）^2）的总和）;

唯一的问题是，我必须定义root函数中的迭代次数，但是一旦我设置了围绕它的任何形式的牙套，程序就无法正常工作。有什么方法可以在根部外定义“ n”，但在根函数内给出一个值？还是只是一些牙套问题，如果我正确设置它们，它仍然可以正常工作。

谢谢预先的

代码：（ 问题开始于“ if（process_rank == root）” - 是的，现在没有括号）

#include <stdio.h>
#include <stdlib.h> 
#include <unistd.h>
#include <mpi.h>
#include <math.h>

#define ROOT 0

void print_time(double time)
{
    printf("TIME:%.10f\n", time * 1000.0); // Print execution time in ms    
}

int main(int argc, char *argv[])
{
    int communicator_size, process_rank; //n
    double pi_appr; //time_elapsed
    double PI_ORIGINAL = 3.141592653589793238462643; // Original pi value for comparison

    double i, n; //size, error, rank
    double result=0.0, sum=0.0, begin=0.0, end=0.0; //pi=0.0

    MPI_Init(&argc, &argv); //error=
    MPI_Comm_size(MPI_COMM_WORLD, &communicator_size); //
    MPI_Comm_rank(MPI_COMM_WORLD, &process_rank); //&process_rank

    //Synchronize all processes and get the begin time
    MPI_Barrier(MPI_COMM_WORLD);
    begin = MPI_Wtime();

    if(process_rank == ROOT)
        
        n = 1000000; // Root defines number of computation iterations 
    

    n = 1000000; //if i dont declare n here again it doesnt work

    //Each process will caculate a part of the sum
    for (i=process_rank; i<n; i+=communicator_size) 
    {
        result += 4/(1+((i/n)*(i/n))); // for some reason pow() didnt work
    }

    //now we some up the results with MPI_Reduce
    MPI_Reduce(&result, &sum, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);
    
    //Synchronize all processes and get the end time
    MPI_Barrier(MPI_COMM_WORLD);
    end = MPI_Wtime();

    if (process_rank == ROOT)
    {
        pi_appr = (1/n)*sum; //calculate 

        printf("%f\n", end-begin); //we print the time by substracting the begin- from the end-time
        printf("Computed pi: %.16f (Error = %.16f)\n", pi_appr, fabs(pi_appr - PI_ORIGINAL));
    }

    MPI_Finalize();

    return 0;
}

Answer 1

您正在误解MPI的工作原理。它具有独立的过程和分布式内存。因此，如果您仅在一个过程中初始化变量，则不会在其他过程中初始化它。您可以使用mpi_bcast调用来传达从根部到其他进程的值。