Hi Dirk,
I am sorry for using the wrong address, it will not happen again.
Furthermore, I think that back in the day when I programmed my code the
interface I used was the recommended one (at least this is what I thought it
was). It was right around the time when the paper in JStatSoft from Doug Bates
and you came out and I tried to strongly follow the recommendations and used
the RcppEigen package skeleton. If the interface has changed meanwhile, I will
do my very best to update it and will look into the vignettes.
Thanks for the pointer,
Henrik
Am 15.10.2014 um 23:56 schrieb Dirk Eddelbuettel:
(Resending, corrected headers)
Henrik,
The list address is rcpp-de...@r-forge.wu-wien.ac.at Would you mind
using that? The gmane redirector messes up automated filtering. Thank you!
Moreover, it also seems to have swallowed my reply. Please do not create
extra work for us. Thanks.
On 15 October 2014 at 23:15, Henrik Singmann wrote:
|=20
|=20
| Hi Dirk and list,
|=20
| I have now managed to isolate the problem. The cpp file posted below =
can on both Windows and Linux crash R if (a) compiled by winbuilder (wi=
ndows only) or if build via devtools::load_all() or from within Rstudio=
("Build & Reload") and (b) when using the appropriate input.
|=20
| The following input crashes R:
| S <- 3
| Ineq <- structure(0, .Dim =3D c(1L, 1L))
| =20
| But this input is okay:
| S <- 3
| Ineq <- structure(c(-1, 1, 0), .Dim =3D c(1L, 3L))
|=20
| =20
| The c++ function producing the problem is (where SS =3D S and ineq =3D=
Ineq):
|=20
| #### code follows ######
|=20
| SEXP determinant2( SEXP SS, SEXP ineq){
| BEGIN_RCPP
I *strongly* recommend you look into the Rcpp Attributes vignette, and
the hundreds of examples in the list archive and on, say, the Rcpp Gallery,
and use an an interface such as
// [[Rcpp::export]]
Rcpp::NumericMatrix foo(Rcpp::NumericMatrix a, Rcpp::NumericMatrix n) { ...
(and I used base Rcpp for simplicity, RcppEigen works the same way).
We had several cases where manual use of RNGScope gets in the way. That may
not be the use, but it makes it easier for everybody to play along --- if you
use the same format we all use.
Otherwise -- thanks, and I'll try to take a look later.
Dirk
PS Your example is not self-contained. What arguments do you supply to
determinant2 ?
--
Dr. Henrik Singmann
PostDoc
Universität Zürich, Schweiz
http://singmann.org
_______________________________________________
Rcpp-devel mailing list
Rcpp-devel@lists.r-forge.r-project.org
https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel
