[Rcpp-devel] RcppArmadillo and 'trunc' problem

Dirk Eddelbuettel edd at debian.org
Wed May 30 13:38:45 CEST 2012


Hi Alexios,

On 30 May 2012 at 12:00, alexios ghalanos wrote:
| Hi,
| 
| While writing some extra functions into the rmgarch package (on r-forge) 
| which uses RcppArmadillo, I encountered the error relating to 'trunc' 
| and 'ftrunc' already discussed here:
| 
| http://lists.r-forge.r-project.org/pipermail/rcpp-devel/2010-September/001123.html

That's a long time ago, a few things changed.  That old email was regarding
OS X, your sessionInfo() shows Windows.
 
| This occurred for the first time once I made use of Rf_exp and Rf_rchisq 
| in code mixed with RcppArmadillo.

Can you provide a self-contained example, maybe via inline as we often do here?
 
| While the workaround described in the previous discussion works (i.e. 
| undefine the two macros in sugar) for me, it will not work for my 
| package once I release the update on r-forge.
| 
| Could you please suggest some solution or some mistake I am making which 
| gives rise to this problem?

We will certainly try to help, but so far we have nothing reproducible -- the
following works, and combines RcppArmadillo with rexp and rchisq:

R> library(inline)
R> 
R> f <- cxxfunction(signature(), plugin="RcppArmadillo", body='
+   arma::colvec X(2);
+   RNGScope scope;
+   X(0) = Rf_rexp(1);
+   X(1) = Rf_rchisq(1);
+   return Rcpp::wrap(X);
+ ')
R> set.seed(42); f()
         [,1]
[1,] 0.198337
[2,] 0.755413
R> set.seed(42); f()
         [,1]
[1,] 0.198337
[2,] 0.755413
R> set.seed(42); f()
         [,1]
[1,] 0.198337
[2,] 0.755413
R> 

Dirk
 
| sessionInfo ()
| R version 2.15.0 (2012-03-30)
| Platform: x86_64-pc-mingw32/x64 (64-bit)
| 
| attached base packages:
| [1] stats     graphics  grDevices utils     datasets  methods   base
| 
| other attached packages:
| [1] RcppArmadillo_0.3.2.0 Rcpp_0.9.10
| 
| loaded via a namespace (and not attached):
| [1] tools_2.15.0
| 
| 
| Thanks and Regards,
| 
| Alexios
| _______________________________________________
| Rcpp-devel mailing list
| Rcpp-devel at lists.r-forge.r-project.org
| https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel

-- 
Dirk Eddelbuettel | edd at debian.org | http://dirk.eddelbuettel.com  


More information about the Rcpp-devel mailing list