如何在R中调用基于Rcpp的MPI C++多处理器程序且无需Rmpi？

可以实现，但不能在单进程R会话中直接调用函数来启动多处理器任务，必须通过MPI启动器创建多进程R实例来执行你的函数。具体步骤如下：

1. 调整C++代码 #

修改原代码以适配多进程启动场景，避免重复初始化MPI和错误调用终止函数：

#include <mpi.h>
#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
int hello_world() {
    int provided;
    // 仅在MPI未初始化时启动环境
    if (!MPI_Initialized(&provided)) {
        MPI_Init_thread(NULL, NULL, MPI_THREAD_FUNNELED, &provided);
    }
  
    int world_size;
    MPI_Comm_size(MPI_COMM_WORLD, &world_size);
  
    int world_rank;
    MPI_Comm_rank(MPI_COMM_WORLD, &world_rank);
  
    // 同步进程输出，避免打印内容乱序
    MPI_Barrier(MPI_COMM_WORLD);
    Rcout << "Hello world from processor rank " << world_rank << " out of " << world_size <<  " processors\n";
    MPI_Barrier(MPI_COMM_WORLD);
  
    // 禁止手动调用MPI_Finalize，R进程退出时会自动清理MPI环境
    // MPI_Finalize();
  
    return 0;
}

2. 编译安装R包 #

通过常规Rcpp包编译流程完成安装，比如在终端执行：

R CMD INSTALL Rmpihelloworld

或在RStudio中通过包管理界面完成安装。

3. 用MPI启动器启动多进程任务 #

在终端中执行以下命令（将x替换为需要的处理器数量）：

mpirun -np x Rscript -e "library(Rmpihelloworld); hello_world()"

预期输出 #

运行后会看到所有进程的输出，例如使用4个处理器时：

Hello world from processor rank 0 out of 4 processors
Hello world from processor rank 1 out of 4 processors
Hello world from processor rank 2 out of 4 processors
Hello world from processor rank 3 out of 4 processors

核心逻辑说明 #

• 本质是通过mpirun/mpiexec创建多个独立的R进程，每个进程加载你的包并执行hello_world，MPI会在这些进程间建立通信域，从而实现多处理器并行。
• 不要在函数内调用MPI_Finalize，否则会导致R进程异常崩溃，MPI环境会在R进程退出时自动清理。
• 后续如需实现进程间数据通信，直接在C++代码中使用标准MPI通信函数（如MPI_Send/MPI_Recv）即可，无需额外依赖Rmpi。

内容的提问来源于stack exchange，提问作者Chris Cheung