Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread Conrad Sand
Hi Baptiste, Thanks for the bug report. I'll take a look at the underlying issues. On 31 May 2011 09:48, baptiste auguie wrote: > Sorry, I had too many spaces in the vectors. The conclusion is the same > though. > > library(inline) > require( RcppArmadillo ) > > ## let's calculate this produc

Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread Conrad Sand
On 31 May 2011 05:02, Savitsky, Terrance wrote: > Dr. Sanderson, I've been able to verify that my issue resides in the > inv() function (that I typically apply to small p x p matrices, where p > = 3 - 10). In particular, the new version/algorithm induces numerical > instability. This is a bit st

Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread baptiste auguie
Sorry, I had too many spaces in the vectors. The conclusion is the same though. library(inline) require( RcppArmadillo ) ## let's calculate this product c(-1-1i, 1-1i) %*% c(1i, -1i) ## 0-2i ## trans() with v0.2.19 fx <- cxxfunction( signature() , ' arma::cx_colvec A = "(-1,-1) (+1,-1)

Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread baptiste auguie
After a couple hours switching back and forth between versions, I found a problem with the scalar product of complex vectors. Here is a minimal example, library(inline) require( RcppArmadillo ) ## let's calculate this product c(-1-1i, 1-1i) %*% c(1i, -1i) ## 0-2i ## trans() with v0.2.19

Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread baptiste auguie
I was a bit optimistic yesterday; I have not yet been able to produce a minimal example but it seems trans/strans was not the end of the story in my code. Something else was broken with the change to 0.2.21.Thankfully the end results are very visibly wrong. Best regards, Baptiste On 31 May 2011

Re: [Rcpp-devel] package .onLoad multiple modules

2011-05-30 Thread baptiste auguie
Hi, On 31 May 2011 07:55, Romain Francois wrote: > Le 30/05/11 01:56, baptiste auguie a écrit : >> >> Thanks, that is helpful. I think I have nailed down the problem. >> >> I defined a wrapper at R level that used a C++ function from a module >> named "cd", >> >> circular_dichroism_spectrum<- fun

Re: [Rcpp-devel] Error: .onLoad failed in loadNamespace()

2011-05-30 Thread Romain Francois
Le 28/05/11 04:04, Dirk Eddelbuettel a écrit : Hi Laurent, On 27 May 2011 at 17:17, Laurent Gatto wrote: | Dear all, | | A basic packages with Rcpp modules produces the warning described | below, that I would ideally, with your help, get rid of. | | Rscript -e "require(Rcpp); Rcpp.package.skele

Re: [Rcpp-devel] package .onLoad multiple modules

2011-05-30 Thread Romain Francois
Le 30/05/11 01:56, baptiste auguie a écrit : Thanks, that is helpful. I think I have nailed down the problem. I defined a wrapper at R level that used a C++ function from a module named "cd", circular_dichroism_spectrum<- function( ... ){ [...] # some argument processing res<- cd$circula

Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread Savitsky, Terrance
Dr. Sanderson, I've been able to verify that my issue resides in the inv() function (that I typically apply to small p x p matrices, where p = 3 - 10). In particular, the new version/algorithm induces numerical instability. I've not yet tested that the results of the inv() computation are general

Re: [Rcpp-devel] RcppArmadillo fails on FreeBSD

2011-05-30 Thread Dirk Eddelbuettel
On 30 May 2011 at 17:21, Rainer Hurling wrote: | > I do not understand what you are asking: "as meant in RcppArmadilloConfig" ? | > Neither Romain, Doug nor I use a *BSD variant (if we ignore OS X as a BSD-derivative). | > So nothing in RcppArmadillo explicitly enables or disables *BSD. | | Than

Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread Savitsky, Terrance
Hello Dr. Sanderson. Thank you, so much, for your work; along with the work of Dirk and Romain, your tools make it possible for me to have impact I otherwise would never achieve. Anyway, I separately tested my (typical) use of trans() in a simple code and received no problems. I have various Rcp

Re: [Rcpp-devel] RcppArmadillo fails on FreeBSD

2011-05-30 Thread Rainer Hurling
On 30.05.2011 16:24 (UTC+1), Dirk Eddelbuettel wrote: On 30 May 2011 at 15:43, Rainer Hurling wrote: | I tried RcppArmadillo_0.2.21.tar.gz with R-devel from today on FreeBSD | 9.0-CURRENT (amd64). The build fails with the following messages: | | -- | #R CM

Re: [Rcpp-devel] RcppArmadillo fails on FreeBSD

2011-05-30 Thread Dirk Eddelbuettel
On 30 May 2011 at 15:43, Rainer Hurling wrote: | I tried RcppArmadillo_0.2.21.tar.gz with R-devel from today on FreeBSD | 9.0-CURRENT (amd64). The build fails with the following messages: | | -- | #R CMD INSTALL RcppArmadillo_0.2.21.tar.gz | * installing t

[Rcpp-devel] RcppArmadillo fails on FreeBSD

2011-05-30 Thread Rainer Hurling
I tried RcppArmadillo_0.2.21.tar.gz with R-devel from today on FreeBSD 9.0-CURRENT (amd64). The build fails with the following messages: -- #R CMD INSTALL RcppArmadillo_0.2.21.tar.gz * installing to library '/usr/local/lib/R/library' * installing *source*

Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread Conrad Sand
On 30 May 2011, Terrance Savitsky wrote: > Hello, After upgrading to the 0.2.21 release of RcppArmadillo, my > previously working code (across many functions) ceased working (on a > Windows XP installation). I re-installed the previous version (0.2.20) > from CRAN (via a server location not yet up

Re: [Rcpp-devel] trans() changed in latest RcppArmadillo

2011-05-30 Thread Dirk Eddelbuettel
On 29 May 2011 at 21:26, Savitsky, Terrance wrote: | Hello, After upgrading to the 0.2.21 release of RcppArmadillo, my previously | working code (across many functions) ceased working (on a Windows XP | installation). I re-installed the previous version (0.2.20) from CRAN (via a | server location