r - rcpparmadlio の自動微分関数?

Question

関数 f の導関数を Rcpp で計算したい。https://cran.r-project.org/web/packages/StanHeaders/vignettes/stanmath.htmlで、スタン ヘッダーと rcppEigen を使用するリソースを見つけまし た。私のプログラムはすべてrcpparmadillioによってコーディングされているので、rcpparmadillioとstanヘッダーによる自動派生関数にどのようにアクセスできるのか疑問に思っています(他のADパッケージかもしれません)。

 #include <RcppArmadillo.h>

// [[Rcpp::plugins("cpp11")]]
// [[Rcpp::depends(RcppArmadillo)]]
#include <RcppArmadilloExtensions/sample.h>
using namespace Rcpp;
using namespace arma;
using namespace std;


/* function */
// [[Rcpp::export]]
doule f(arma::rowvec x){
 arma::mat yy=x.t()*x;
 double res=arma::sun(x);
 rerutn(res);
}