Fortran、C++をRと統合する 質問する

テスト関数.cpp

ここで、この Fortran コードを C++ コードから呼び出す必要があります。これを行うときは、いくつかの点を考慮する必要があります。

1. Fortran 引数は値ではなく参照によって渡されます。

2. は別のファイルで定義されているためMULTIPLY、コンパイラが引数と戻り値の型を認識できるように、C++ ファイルで宣言する必要があります。

   a. C++ ファイルの Fortran 関数を宣言するときに、関数名の大文字と小文字を区別せず、アンダースコアを追加します。これは、Fortran コンパイラがデフォルトでこれを行うためです。

   b. 関数をextern "C"リンケージ仕様内で宣言する必要があります。C++ コンパイラーは、オーバーロードが許可されているため、通常、関数名を一意の識別子として使用することはできませんが、Fortran 関数を呼び出す場合は、リンケージ仕様で達成されるものとまったく同じことを行う必要がありますextern "C"(たとえば、このSOの答え）。

#include "RcppArmadillo.h"

// [[Rcpp::depends(RcppArmadillo)]]

// First we'll declare the MULTIPLY Fortran function
// as multiply_ in an extern "C" linkage specification
// making sure to have the arguments passed as pointers.
extern "C" {
    double multiply_(double *x, double *y);
}

// Now our C++ function
// [[Rcpp::export]]
Rcpp::NumericVector test_function(Rcpp::NumericVector x) {
    // Get the size of the vector
    int n = x.size();
    // Create a new vector for our result
    Rcpp::NumericVector result(n);
    for ( int i = 0; i < n; ++i ) {
        // And for each element of the vector,
        // store as doubles the element and the index
        double starting_value = x[i], multiplier = (double)i;
        // Now we can call the Fortran function,
        // being sure to pass the address of the variables
        result[i] = multiply_(&starting_value, &multiplier);
    }
    return result;
}

mixedlang::test_function(0:9)
# [1]  0  1  4  9 16 25 36 49 64 81