[Rcpp-devel] armadillo mailing list?

baptiste Auguié baptiste.auguie at googlemail.com
Thu Jul 29 15:00:22 CEST 2010


On Jul 29, 2010, at 1:23 PM, Romain Francois wrote:

> Le 29/07/10 09:30, baptiste auguie a écrit :
>> Dear all,
>> 
>> I'm porting some R code to C++ using the armadillo framework (I'm
>> using operations on complex matrices, etc.). The final goal is of
>> course to wrap the code using RcppArmadillo. Currently, I'm still
>> struggling with my little C++ knowledge. Does anyone know of an
>> Armadillo mailing list where I could ask C++ questions?
> 
> From armadillo's faq http://arma.sourceforge.net/faq.html
> you can find a link to a forum:
> http://sourceforge.net/apps/phpbb/arma/

I don't know how I'd missed it, thanks!

It turns out that my example would work, one only needs to make sure that the replacement has the correct type (matrix) and dimensions:

this works

 mat A = "0 0 0; 0 0 0; 0 0 0;";
 A.submat(1,1,1,1) = ones(1);

Thanks,

baptiste

> 
>> Just in case, I'll ask here too. I have a small 3x3 matrix that I
>> generate in a nested for loop, and I want to insert this matrix as a
>> block inside a larger matrix. In R language, I'd want to use a
>> construct like,
>> 
>> A[ seq(ii, ii+2), seq(jj, jj+2)]<- matrix(rnorm(9), 3, 3)
>> 
>> I tried the following as a test,
>> 
>> mat A = "0 0 0; 0 0 0; 0 0 0;";
>> A.submat(1,1,1,1) = 1;
>> 
>> but it seems that .submat is only to view elements, not assign them
>> (the `[<-` method in R). Any tip?
>> 
>> Sincerely,
>> 
>> baptiste
> 
> That is the kind of thing we should add in Rcpp sugar at some point, but currently I don't know how to.
> 
> Romain
> 
> -- 
> Romain Francois
> Professional R Enthusiast
> +33(0) 6 28 91 30 30
> http://romainfrancois.blog.free.fr
> |- http://bit.ly/aryfrk : useR! 2010
> |- http://bit.ly/bc8jNi : Rcpp 0.8.4
> `- http://bit.ly/dz0RlX : bibtex 0.2-1
> 



More information about the Rcpp-devel mailing list