Thanks so much, Conrad! Best regards,
Baptiste On 31 May 2011 23:37, Conrad Sand <conradsand.r...@gmail.com> wrote: > I've found the cause of the issue (involves handling of small matrices). > > The fix is easy (and already done in the SVN repo), but I currently > don't have the time to roll out another release. I'll aim to release > a fix on the weekend, which will also give me time to double-check if > there is a problem in inv(). > > In the meantime I recommend sticking to the previous version of > Armadillo and RcppArmadillo. > > btw, for dot products I recommend using the dot() and cdot() functions > -- they're generally faster than going through the multiplication > operator. > > > On 31 May 2011 11:39, Conrad Sand <conradsand.r...@gmail.com> wrote: >> Hi Baptiste, >> >> Thanks for the bug report. I'll take a look at the underlying issues. >> >> >> On 31 May 2011 09:48, baptiste auguie <baptiste.aug...@googlemail.com> wrote: >>> 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);", B = "(+0,1) (+0,-1);" ; >>> arma::cx_colvec res = trans(A) * B; >>> >>> return wrap( res ) ; >>> ', plugin = "RcppArmadillo" ) >>> >>> fx() >>> ## 0-2i >>> >>> ## strans() with v0.2.21 >>> fx <- cxxfunction( signature() , ' >>> arma::cx_colvec A = "(-1,-1) (+1,-1);", B = "(+0,1) (+0,-1);" ; >>> arma::cx_colvec res = strans(A) * B; >>> return wrap( res ) ; >>> ', plugin = "RcppArmadillo" ) >>> >>> fx() >>> ## 1-1i >> > _______________________________________________ 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