Rcppの複素ガンマ関数

Question

from the docを理解しているので、boost C++ライブラリのtgamma関数を複素数で評価することができます。

私はRcppで使用しようとしています。これは私のコードです：

// [[Rcpp::depends(BH)]] 
#include <Rcpp.h> 
#include <boost/math/special_functions/gamma.hpp> 

// [[Rcpp::export]] 
std::complex<double> gamma_boost(std::complex<double> z) { 
    std::complex<double> result = tgamma(z); 
    return result; 
} 

このコードはコンパイルされません。

cannot convert 'std::complex<double>' to 'double' for argument '1' to 'double tgamma(double)' 

Answer 1

あなたがここで間違ってtgamma()があり、boost::math::tgamma(...)を試してみてください。私はエラーを取得します。そして/またはstd::complexにテンプレートを作成する必要があるかもしれません。

私の通常のアプローチは、まずコマンドラインで何かを稼働させ、Rcpp経由でRにそのコードを添付することです。

Answer 2

基本的に間違った機能を呼び出しています。

名前空間を指定しなかったため、のため、ADLは標準ライブラリのstd::tgammaを見つけました。

パラメータとしてstd::complexを使用しないため、コンパイラエラーが発生します。代わりにboost::math::tgammaが必要です。

しかしブーストのtgammaもstd::complex種類をサポートしていませんので、あなたが別のライブラリを使用するか、それを自分で実装する必要があります。