通过Rcpp实现R向Fortran传递多维数组及编译错误解决

1. 修改代码解决编译错误的具体方案 #

针对你遇到的函数调用不匹配、指针类型转换、可变长度数组禁用等问题，按以下步骤调整：

• 正确声明Fortran子程序的C接口
  Fortran编译器会对子程序名做修饰（比如GCC默认转小写加下划线mytest_，Intel编译器可能转大写MYTEST），所以在Rcpp代码里必须用extern "C"声明兼容接口，避免C++的名字修饰导致函数找不到：

  extern "C" {
      // 假设Fortran子程序为 subroutine mytest(arr, n)，参数是double数组和整数长度
      void mytest_(double* arr, int* n); 
  }
  

  如果编译仍提示找不到函数，可通过nm mytest.o查看Fortran目标文件的符号名，再调整接口中的函数名。

• 正确传递数组与类型转换
  Rcpp的NumericVector/NumericMatrix底层是连续内存，直接用.data()获取double*指针即可，长度参数要传地址（Fortran默认传址调用）：

  #include <Rcpp.h>
  using namespace Rcpp;
  
  extern "C" {
      void mytest_(double* arr, int* n);
  }
  
  // [[Rcpp::export]]
  NumericVector call_mytest(NumericVector input) {
      int n = input.size();
      // 直接传递指针和长度的地址
      mytest_(input.data(), &n);
      return input;
  }
  

  如果是多维数组，需把每个维度的长度都作为参数传给Fortran子程序，比如2D数组要传递行数和列数。

• 解决可变长度数组问题
  如果Fortran子程序使用了可变长度数组（比如real*8 :: arr(n)），首先确保子程序已将数组长度作为参数传入（比如subroutine mytest(arr, n)）；其次在包的Makevars文件中启用Fortran 90+标准，避免编译器禁用可变长度数组：

  FCFLAGS = -std=f95 -O2
  OBJECTS = mytest.o Rcppmytest.o
  

2. 必须关注列/行主序差异 #

R和Fortran的数组都是列主序（按列存储），而C/C++是行主序（按行存储），这个差异在处理多维数组时必须重视：

• 直接传递R的NumericMatrix给Fortran时，无需调整顺序，两者存储逻辑一致；
• 若在C代码中手动构造多维数组再传给Fortran，必须反转维度顺序（比如C里的[row][col]要转为Fortran里的[col][row]），否则数组内容会错乱；
• 传递多维数组的维度参数时也要对应调整：比如R矩阵是3行4列，传给Fortran时要先传列数4，再传行数3。

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