如何从Fortran调用R函数？

Question

根据 http://gallery.rcpp.org/articles/r-function-from-c++,Rcpp 允许用户从C++调用R函数.Fortran中有类似的东西,以便人们可以在Fortran代码中调用R函数吗？

Answer 1

我相信你应该看看RFortran？AFAIK 它是 R 对象的唯一 Fortran 绑定，而且它是开源的。

编辑

根据下面的评论，我不知道 RFortran 仅适用于 Windows。因此，为了获得更便携的答案，让我们创建一个示例，在其中我将利用 RInside ，这使得在 C++ 代码中嵌入 R 函数变得更简单。我也会使用iso_c_binding与 C 的接口。

测试C.cpp

#include <iostream>
#include <RInside.h>                    

void helloR_(int argc, char *argv[], const char *msg);

extern "C" void helloR(int argc, char *argv[], const char *msg) {

    // create an embedded R instance
    RInside R(argc, argv);              

    // convert to string for RInside assignment
    std::string txt = std::string(msg);

    // C++ Notice
    std::cout << "This is C++, " << txt << std::endl;

    // Assign string to R object
    R.assign(txt, "txt");

    // eval the string, give R notice
    R.parseEvalQ("cat('This is R, ', txt, '\n')");
}

测试F.f

  PROGRAM MAIN
    USE iso_c_binding
    IMPLICIT NONE
    INTEGER :: argc
    CHARACTER(len=32) :: arg
    CHARACTER(len=32) :: msg

    INTERFACE
      SUBROUTINE R_FUN(argc, arg, msg) bind(C, name="helloR")
        USE iso_c_binding
        INTEGER(kind=c_int), INTENT(IN) :: argc
        CHARACTER(kind=c_char), INTENT(IN) :: arg(*)
        CHARACTER(kind = C_CHAR), INTENT(IN) :: msg(*)
      END SUBROUTINE R_FUN
    END INTERFACE

    print *, "Fortran Calling RInside"
    CALL R_FUN (argc, arg, "Hello World"//C_NULL_CHAR)

  END PROGRAM

Fortran 编译很简单：

gfortran -c testF.f

C++ 编译有点棘手，因为您必须知道 R、Rcpp 和 RInside 的包含目录所在的位置。

g++ testC.cpp -c -I/path/to/RInside/include -I/path/to/Rcpp/include -I/usr/share/R/include

然后，您需要提供正确的库和标志的链接lgfortran。

g++ -o fcr testF.o testC.o -L/usr/lib/R/lib -lR -L/path/to/RInside/lib/ -lRInside -L/path/to/Rcpp/libs/ -Wl,-rpath,/home/path/to/RInside/lib/ -lRInside -Wl,-rpath,/path/to/Rcpp/libs/ -lgfortran

现在我有一个小程序演示如何从 Fortran 访问 R 函数

./fcr
 Fortran Calling RInside
This is C++, Hello World
This is R,  Hello World