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
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
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)
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
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
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
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
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
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
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
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
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
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
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*
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
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
16 matches
Mail list logo