[Rcpp-devel] RcppEigen: Getting and setting dimnames of sparse matrix
Dirk Eddelbuettel
edd at debian.org
Sat Jan 19 00:17:40 CET 2013
Hi Søren,
On 18 January 2013 at 23:07, Søren Højsgaard wrote:
| Dear all,
| I put a sparse "dgCMatrix" into RcppEigen code as a sparse matrix. In Rcpp I create a new matrix of same dimensions as the input matrix. I would like to copy the dimnames from the input matrix to the output matrix. Can this be done? (It is essential that I can stick to working with sparse matrices). I know that I can set the dimensions on the R-side, but that involves AFAIK copying the object which takes time.
"It should".
I have not worked all that much with sparse matrices, but dimnames etc are
merely attributes at the R level. And we read and write these just as easily
at the C++ level.
Have a look at few of the recent posts at http://gallery.rcpp.org
in particular
http://gallery.rcpp.org/articles/getting-attributes-for-xts-example
and
http://gallery.rcpp.org/articles/creating-xts-from-c%2B%2B
where I, respectively, "get" and "set" the attributes that make an xts
object: dimnames, an index (itself with class and time zone attributes) and
timezone and class info.
For sparse matrices the mechanics may be slighty different due to the S4-ness
of things but the rest should work. I put a feeble attempt at working with
sparse matrices here:
http://gallery.rcpp.org/articles/armadillo-sparse-matrix/
Maybe we can do something similar with richer 'matrix factory' functions for
Eigen?
Dirk
|
| Example:
|
| src<-'
| using namespace Rcpp;
| typedef Eigen::SparseMatrix<double> SpMat;
| SpMat X(as<SpMat>(XX_));
| return(wrap(X));
| '
|
| names_ <- cxxfunction(signature(XX_="matrix"), plugin="RcppEigen",
| body=src)
|
| ugM <- new("dgCMatrix"
| , i = c(1L, 0L, 2L, 1L)
| , p = c(0L, 1L, 3L, 4L)
| , Dim = c(3L, 3L)
| , Dimnames = list(c("a", "b", "c"), c("a", "b", "c"))
| , x = c(1, 1, 1, 1)
| , factors = list())
|
| > ugM
| 3 x 3 sparse Matrix of class "dgCMatrix"
| a b c
| a . 1 .
| b 1 . 1
| c . 1 .
| >
| > names_(ugM)
| 3 x 3 sparse Matrix of class "dgCMatrix"
|
| [1,] . 1 .
| [2,] 1 . 1
| [3,] . 1 .
|
| Any advice? Thanks in advance
| Søren
| _______________________________________________
| Rcpp-devel mailing list
| Rcpp-devel at lists.r-forge.r-project.org
| https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel
--
Dirk Eddelbuettel | edd at debian.org | http://dirk.eddelbuettel.com
More information about the Rcpp-devel
mailing list