web-dev-qa-db-ja.com

Rcpp関数はパッケージ化されるとRStudioをクラッシュさせます

単一の関数でRcpp.package.skeleton(myPackage)を使用してダミーのRcppパッケージを作成しました。

_// [[Rcpp::export]]
double triple_balance(const double& balance) {

  if(R_IsNA(balance)) {
    stop("balance is NA, please provide a numeric value");
  } else {
     double result = balance*3;
    return result;
  }
}
_

単体テストでは、間違った入力が渡されたときに関数がエラーを返すことを確認したいと思います。

_Rcpp::sourceCpp('src/triple_balance.cpp')
triple_balance("10")
_

Triple_balance( "10")のエラー:要求されたタイプと互換性がありません:[type = character; target = double]。

ただし、パッケージをインストールして再起動すると、triple_balance("10")を実行するとRStudioがクラッシュします。クラッシュは毎回発生するわけではなく、ライブラリの初期状態に依存しているようです。以前にパッケージの以前のバージョンをインストールしたかどうか。たとえば、関数定義にconstを付けても付けなくても。

このような一貫性のない行動を説明できるものは何ですか?

私のsessionInfo

_> sessionInfo()
R version 3.6.0 (2019-04-26)
Platform: x86_64-redhat-linux-gnu (64-bit)
Running under: Red Hat Enterprise Linux Server 7.8 (Maipo)

Matrix products: default
BLAS/LAPACK: /usr/lib64/R/lib/libRblas.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] rcppDummy_1.0

loaded via a namespace (and not attached):
[1] compiler_3.6.0 tools_3.6.0    Rcpp_1.0.2     Packrat_0.5.0 
_
3
Kasia Kulma

ユースケースでは

Rcpp::sourceCpp('src/triple_balance.cpp')
triple_balance("10")

あなたのインターフェースがdoubleを使用しているので、エラーはあなたのものです

double triple_balance(const double& balance)

つまり、ではありません"10"-これは文字ベクトル(または「文字列」)です!

ベクトル引数を使用すると、RcppでautoMAGICallyにそれをキャッチできます。

// [[Rcpp::export]]
NumericVector triple_balance(const NumericVector & balance) {
    return balance * 3;
}

これは、Rcppデータ構造がいつインスタンス化されるかをチェックします。

R> triple_balance(10)
[1] 30
R> triple_balance(10L)
[1] 30
R> triple_balance("10")
Error in triple_balance("10") : 
  Not compatible with requested type: [type=character; target=double].
R> 

何よりも、ベクターに対して自動的に機能します。

R> triple_balance(c(c(2,3), seq(10,30,by=10), log10(1:3)))
[1]  6.00000  9.00000 30.00000 60.00000 90.00000  0.00000  0.90309  1.43136
R> 

また、非有限値も処理します。

R> triple_balance(c(10, NA, NaN, Inf, -Inf))
[1]   30   NA  NaN  Inf -Inf
R> 

編集ちなみに、私があなたのバージョンを_origを追加しました。私も問題なく動作します(Ubuntu 20.04;現在のRcpp)

R> triple_balance_orig("10")
Error in triple_balance_orig("10") : 
  Not compatible with requested type: [type=character; target=double].
R> packageVersion("Rcpp")
[1] ‘1.0.4.11’
R> 

その動作は新しいコードではないので、RHELインストールがクラッシュすることに少し驚いています。

編集2:議論のために、パッケージ化も行いましたが、RStudioもダウンしません。最新のテストバージョン1.4.390を実行しています。 (パッケージとエラーの再コンパイルでこのようなクラッシュが発生したことがありました。これは発生する可能性があります。ここでは発生しないはずです。)

https://rstudio.cloud でコードを試し、別の環境をすばやく試すことができます。

5