코드 설명 (MPICH)

Question

#include "mpi.h" 
#include <stdio.h> 
#include <math.h> 


double f(double a) 
{ 
    return (4.0/(1.0 + a*a)); 
} 

void main(int argc, char *argv[]) 
{ 
    int done = 0, n, myid, numprocs,i; 
    double PI25DT = 3.141592653589793238462643; 
    double mypi, pi, h, sum, x; 
    double startwtime, endwtime; 
    int namelen; 
    char processor_name[MPI_MAX_PROCESSOR_NAME];  
    MPI_Init(&argc,&argv); 
    MPI_Comm_size(MPI_COMM_WORLD,&numprocs); 
    MPI_Comm_rank(MPI_COMM_WORLD,&myid); 
    MPI_Get_processor_name(processor_name,&namelen);  
    fprintf(stderr,"Process %d on %s\n", 
     myid, processor_name); 
    fflush(stderr); 
    n = 0; 
    while (!done) 
    { 
     if (myid == 0) 
     { 
     printf("Enter the number of intervals: (0 quits) ");fflush(stdout); 
     scanf("%d",&n); 

     startwtime = MPI_Wtime(); 
     } 
     MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD); 
     if (n == 0) 
      done = 1; 
     else 
     { 
      h = 1.0/(double) n; 
      sum = 0.0; 
      for (i = myid + 1; i <= n; i += numprocs) 
      { 
       x = h * ((double)i - 0.5); 
       sum += f(x); 
      } 
      mypi = h * sum; 

      MPI_Reduce(&mypi, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD); 

      if (myid == 0) 
     { 
       printf("pi is approximately %.16f, Error is %.16f\n", 
      pi, fabs(pi - PI25DT)); 
     endwtime = MPI_Wtime(); 
     printf("wall clock time = %f\n", endwtime-startwtime);   
     } 
     } 
    } 
    MPI_Finalize(); 
} 

누구든지 위 코드를 어떻게 설명 할 수 있습니까? 나는 실험실에 있는데, 내가 그리워서 설명해달라고 부탁했는데 이해가 안된다.

Answer 1

이것은 pi를 계산하기 위해 MPralH 라이브러리에서 parrallel 처리를 사용하는 예입니다. 이 예제는 mpich 설치에 포함되어 있습니다. 본질적으로 pi를 계산하는 작업을 여러 프로세서 또는 호스트 사이에 분수를 더하여 나눕니다.