在C++中使用稀疏矩陣的Rcpp

Question

從RCCP開始不堪重負。我如何能夠像代碼中定義的那樣，使用（索引，讀取和賦值）稀疏矩陣，就像我能用'標準'矩陣一樣？

library('Matrix') 
library(Rcpp) 
library(inline) 

r <- matrix(seq(1,9,1),ncol=3,nrow=3) 
i <- Matrix(0, nrow = nrow(r) * ncol(r), ncol = nrow(r)*ncol(r), sparse=TRUE) 
fx <- cxxfunction(signature(x_ = "matrix" ,y_="dsCMatrix"), ' 
    NumericMatrix x(x_) ; 
    int nr = x.nrow(), nc = x.ncol() ; 
    for (int i = 0; i < nr; i++) { 
      for (int j = 1; j < nc; j++) { 
      x(i,j) = 1; 
     } 
    } 
    return wrap(x) ; 
', plugin = "Rcpp") 

fx(r,i) 

Answer 1

你最好的選擇可能是

• 在Rcpp Gallery的有關稀疏矩陣的帖子，並
• 的RCPP-devel郵件列表。