Rcpp参照によって渡された更新行列がRで返されます。

Question

Rcppを使用して参照渡しの方法に関する多くの例を確認しました。私は例えばthis非常に大きい参照してください。しかし、私は1つの質問があります。 Rのオブジェクトとして行列があり、エントリ[1,1]に1を加えたいとします。行列がCppにある場合の例を私が見たが、私はリターンステートメントを使用せずにRで更新を返したい。

この

は私がリストで行った例であり、それは非常によく

//[[Rcpp::export]] 
void test(List& a){ 
a(0)=0; 
} 

私はマトリックスとsimilarely行う必要があります動作します。

//[[Rcpp::export]] 
void test(arma::mat& a){ 
    a(0,0)=0; 
} 

2番目の列は更新されずにリストが更新されます。

誰でも手伝ってもらえますか？

Answer 1

これはおそらく悪い習慣であることを再確認してみましょう。 voidを使用しないで、変更されたオブジェクトを返す - より一般的なアプローチです。

つまり、どちらの方法でも動作させることができます。 RcppArmadilloの場合は、（明示的な）参照を渡します。私は希望の動作この短い例のうち

> sourceCpp("/tmp/so.cpp") 

> M1 <- M2 <- matrix(0, 2, 2) 

> bar(M1) 

> M1 
    [,1] [,2] 
[1,] 42 0 
[2,] 0 0 

> foo(M2) 

> M2 
    [,1] [,2] 
[1,] 42 0 
[2,] 0 0 
> 

を得る：

#include <RcppArmadillo.h> 

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

// [[Rcpp::export]] 
void bar(Rcpp::NumericMatrix M) { 
    M(0,0) = 42; 
} 

// [[Rcpp::export]] 
void foo(arma::mat M) { 
    M(0,0) = 42; 
} 

/*** R 
M1 <- M2 <- matrix(0, 2, 2) 

bar(M1) 
M1 

foo(M2) 
M2 
*/