[Rcpp-devel] tracking volatile bug

2014-11-18 Thread Serguei Sokol
Hi everybody, I am facing a volatile bug which appears and disappears in identical calls on identical data. My RcppArmadilla code is expm_cpp() (it is obtained by rex2arma tool which I have presented before, cf. attached file). It is compared and benchmarked vs its R counterpart expm.higham() (a

Re: [Rcpp-devel] tracking volatile bug

2014-11-18 Thread Dirk Eddelbuettel
Serguei, First, I think we mentioned to you after your first post on rex2arma that all serious Rcpp development is generally done in packages. You should probably embrace that mode too. Many of us prefer it. Second, 'expm' does not need to be reimplemented. See my RcppKalman package (on GitHub

Re: [Rcpp-devel] tracking volatile bug

2014-11-18 Thread Romain François
> Le 18 nov. 2014 à 14:42, Serguei Sokol a écrit : > > Hi everybody, > > I am facing a volatile bug which appears and disappears in > identical calls on identical data. > > My RcppArmadilla code is expm_cpp() (it is obtained by rex2arma tool > which I have presented before, cf. attached file).

Re: [Rcpp-devel] tracking volatile bug

2014-11-18 Thread Serguei Sokol
Le 18/11/2014 16:39, Romain François a écrit : Le 18 nov. 2014 à 14:42, Serguei Sokol a écrit : /* call this R command: source("FileName") */ int errorOccurred; SEXP e = lang2(install("source"), mkString("FileName")); /* mkChar instead of mkString would lead to this runtime er

Re: [Rcpp-devel] tracking volatile bug

2014-11-18 Thread William Dunlap
> I will try "gctorture(TRUE)" suggested by Martin. Try running it under valgrind as well. E.g., % R --debugger=valgrind ==15338== Memcheck, a memory error detector ==15338== Copyright (C) 2002-2011, and GNU GPL'd, by Julian Seward et al. ==15338== Using Valgrind-3.7.0 and LibVEX; rerun w

[Rcpp-devel] Bug in loadRcppClass/loadModule?

2014-11-18 Thread Aaron Polhamus
Dear list, In the process of writing a comprehensive unit testing application for Rcpp I may have come across a bug in the code. It seems to me that the following block should execute just fine, exporting the C++ class "norm" to the global environment: require(inline) require(Rcpp) inc <- ' usin