Re: [Rcpp-devel] Sparse matrix operations

2014-02-04 Thread Dirk Eddelbuettel
On 4 February 2014 at 11:09, Yixuan Qiu wrote: | Sure. What are you planning to do with ARPACK? It seems that the newest version | of Armadillo already supports eigs_sym() and eigs_gen() to retrieve a limited | number of eigenvalues. Recall that we link with (post Rcpp 0.11.0) PKG_LIBS=

Re: [Rcpp-devel] Sparse matrix operations

2014-02-04 Thread Yixuan Qiu
Sure. What are you planning to do with ARPACK? It seems that the newest version of Armadillo already supports eigs_sym() and eigs_gen() to retrieve a limited number of eigenvalues. Best, Yixuan 2014-02-03 Dirk Eddelbuettel : > > On 3 February 2014 at 18:09, Yixuan Qiu wrote: > | A pretty nice

Re: [Rcpp-devel] Sparse matrix operations

2014-02-04 Thread Dirk Eddelbuettel
On 4 February 2014 at 09:30, Søren Højsgaard wrote: | Just to clarify this (for myself): | | In a broader perspective, isn't it so now that there is very limited need for as< >(some_SEXP)? | | When calling from R we can declare input to be NumericMatrices, NumericVectors ets, but we can also

Re: [Rcpp-devel] Sparse matrix operations

2014-02-04 Thread Søren Højsgaard
[mailto:e...@debian.org] Sent: 3. februar 2014 20:20 To: Douglas Bates Cc: Dirk Eddelbuettel; French, Joshua; Søren Højsgaard; rcpp-devel@lists.r-forge.r-project.org Subject: Re: [Rcpp-devel] Sparse matrix operations On 3 February 2014 at 12:59, Douglas Bates wrote: | On Mon, Feb 3, 2014 at 11:4

Re: [Rcpp-devel] Sparse matrix operations

2014-02-03 Thread Dirk Eddelbuettel
On 3 February 2014 at 18:09, Yixuan Qiu wrote: | A pretty nice feature of Eigen is that it allows you to map an existing memory | address to an Eigen object, so the data stored in R can be directly reused by | Eigen. This may be the key point if we want to avoid copying data from R to | Eigen/Arma

Re: [Rcpp-devel] Sparse matrix operations

2014-02-03 Thread Yixuan Qiu
A pretty nice feature of Eigen is that it allows you to map an existing memory address to an Eigen object, so the data stored in R can be directly reused by Eigen. This may be the key point if we want to avoid copying data from R to Eigen/Armadillo. Dirk, does Armadillo has such kind of constructor

Re: [Rcpp-devel] Sparse matrix operations

2014-02-03 Thread Dirk Eddelbuettel
On 3 February 2014 at 12:59, Douglas Bates wrote: | On Mon, Feb 3, 2014 at 11:48 AM, Dirk Eddelbuettel wrote: | | | On 3 February 2014 at 10:23, French, Joshua wrote: | | Soren and Doug, | | | | Thanks for the info about the sparse matrices.  I'll give it a go and | seen wha

Re: [Rcpp-devel] Sparse matrix operations

2014-02-03 Thread Douglas Bates
On Mon, Feb 3, 2014 at 11:48 AM, Dirk Eddelbuettel wrote: > > On 3 February 2014 at 10:23, French, Joshua wrote: > | Soren and Doug, > | > | Thanks for the info about the sparse matrices. I'll give it a go and > seen what > | happens. If it works out, then perhaps I'll be able to provide a nice

Re: [Rcpp-devel] Sparse matrix operations

2014-02-03 Thread Dirk Eddelbuettel
On 3 February 2014 at 10:23, French, Joshua wrote: | Soren and Doug, | | Thanks for the info about the sparse matrices. I'll give it a go and seen what | happens. If it works out, then perhaps I'll be able to provide a nice | RcppGallery example. Yes -- there isn't much yet in RcppArmadillo b

Re: [Rcpp-devel] Sparse matrix operations

2014-02-03 Thread French, Joshua
t;>, Joshua French mailto:joshua.fre...@ucdenver.edu>> Cc: "rcpp-devel@lists.r-forge.r-project.org<mailto:rcpp-devel@lists.r-forge.r-project.org>" mailto:rcpp-de...@r-forge.wu-wien.ac.at>> Subject: RE: [Rcpp-devel] Sparse matrix operations I’ll just add that I have used the sparse

Re: [Rcpp-devel] Sparse matrix operations

2014-01-31 Thread Søren Højsgaard
sts.r-forge.r-project.org] On Behalf Of Douglas Bates Sent: 31. januar 2014 21:45 To: French, Joshua Cc: rcpp-devel@lists.r-forge.r-project.org Subject: Re: [Rcpp-devel] Sparse matrix operations On Fri, Jan 31, 2014 at 1:14 PM, French, Joshua mailto:joshua.fre...@ucdenver.edu>> wrote: Hello

Re: [Rcpp-devel] Sparse matrix operations

2014-01-31 Thread Douglas Bates
On Fri, Jan 31, 2014 at 1:14 PM, French, Joshua wrote: > Hello everyone, > > For those of you who have used the sparse matrix capabilities of > Armadillo/RcppArmadillo, how seamless are the operations in moving between > sparse and dense matrices? I know there are some tricks for getting > sparse

[Rcpp-devel] Sparse matrix operations

2014-01-31 Thread French, Joshua
Hello everyone, For those of you who have used the sparse matrix capabilities of Armadillo/RcppArmadillo, how seamless are the operations in moving between sparse and dense matrices? I know there are some tricks for getting sparseMatrix objects into SpMat objects (Dirk has a nice Rcpp Gallery