为什么在尝试实施OpenMPI时设置单元格数

Question

我正在尝试学习如何使用OpenMPI并遇到了这个示例代码

#include "mpi.h”
int main(int argc, char **argv)
{
  // Call MPI initialization
  MPI_Init(&argc, &argv);

  // Get my processor ID and number of processor
  int myProcID;
  int numProcs;
  MPI_Comm_rank(MPI_COMM_WORLD, &myProcID);
  MPI_Comm_size(MPI_COMM_WORLD, &numProcs);

  // Set number of cells from command line argument
  int numCells  = atoi(argv[1]);

 <etc…>

  // Find which piece of the integral this 
  // processor is responsible for
  int numCellsPerProc = numCells/numProcs;
  int myStart = myProcID * numCellsPerProc;
  int myEnd   = myStart + numCellsPerProc;

  // Account for unequal load by making sure 
  // last processor has correct end
  if (myProcID == numProcs - 1) myEnd = numCells;
  // Loop over cells and compute integral 
  // using trapezoidal rule
  double myResult = 0.0;
  for (int i = myStart; i < myEnd; ++i)
  {
    double xL = xMin + i*dx;
    double xR = xL + dx;
    myResult += 0.5 * (myFunction(xL)+myFunction(xR)) * dx;
  }

  // Sum result across processors
  double totalResult;
  MPI_Reduce(&myResult, &totalResult, 1, MPI_DOUBLE, MPI_SUM,
             0, MPI_COMM_WORLD);

  // Print the result, but only from root processor
  if (myProcID == 0)
  {
    std::cout << "result = ";
    std::cout << std::fixed << std::setprecision(10)
              << totalResult << std::endl;
  }

  // Call MPI_Finalize
  MPI_Finalize();

  return 0;
}

<etc>

当谈到处理器的实际架构时，请原谅我的无知。为什么示例代码设置单元格数？我认为每个处理器作为一个整体负责一次一个工作？ 我根本不明白这些台词......

  // Set number of cells from command line argument
  int numCells  = atoi(argv[1]);

 <etc…>

  // Find which piece of the integral this 
  // processor is responsible for
  int numCellsPerProc = numCells/numProcs;
  int myStart = myProcID * numCellsPerProc;
  int myEnd   = myStart + numCellsPerProc

Answer 1

它取决于命令行参数 - argv[1] - 每个节点将有多少个作业（例如，在OpenMPI中，您可以通过-N指定每个节点的作业数） 。此外，您可以生成线程以使用多核处理器。

实际上，您正在计算积分。
您将积分间隔拆分为numProcs个部分，因此每个作业计算他的部分，并且最后，所有这些都归结为减少。

（单词cell - 在此上下文中不是一个好的变量名称）