This is really a brilliant observation!
Best,
Yixuan
On Tue, Nov 27, 2018 at 9:34 AM Serguei Sokol wrote:
>
> Le 26/11/2018 à 18:23, Hoffman, Gabriel a écrit :
> > I am developing a statistical model and I have a prototype working in R
> > code. I make extensive use of sparse matrices, so the
I guess BLAS plays a minor role here. The operation here is a sparse
Cholesky decomposition, and if I'm correct the Matrix package in R uses
SuiteSparse/CHOLMOD to do the factorization, while Eigen uses its own
implementation. The underlying numerical algorithms may be different.
Best,
Yixuan
On
e, maybe stick to
>> equivalent functions from a non-R source: C++11, Boost, ...
>>
>> Dirk
>>
>> | (It seems this question must have arisen before, but I haven't found an
>> | answer)
>>
>> --
>> http://dirk.eddelbuettel.com | @eddelbuettel |
pectra
>
> Dirk
>
> --
> http://dirk.eddelbuettel.com | @eddelbuettel | e...@debian.org
> ___
> 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
--
> School of Public Health & Department of Surgery and Cancer
>
> Main address :
>
> Imperial College London | Room 171 | St Mary’s Campus | Norfolk Place | London
> | W2 1PG |
>
>
>
> _______
> Rcpp-devel mailing
m not
sure whether this is a bug of MinGW or of Eigen, but at least we can
start from here to investigate.
(I would suggest future discussions moving to
https://github.com/RcppCore/RcppEigen/issues/34)
Best,
Yixuan
--
Yixuan Qiu
Department of Statistics,
Purdue Unive
| [1] tools_3.2.4 RcppArmadillo_0.6.700.3.0 Rcpp_0.12.4
> | >
> | ___
> | 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
> --
> http://dirk.eddelbuettel.com | @eddelbuettel |
ppreciated.
>
> Regards,
> Charles
>
> ___
> 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
>
--
Yixuan Qiu
Department of Statistics,
project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel
>
--
Yixuan Qiu
Department of Statistics,
Purdue University
___
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
Universidade Federal do Paraná (UFPR)
>
>
> ___
> 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
>
--
Yixuan Qiu
Department of Statistics,
Purdue University
___
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
Oops sorry, should be begin(). data() is the Eigen method.
Best,
Yixuan
On Feb 17, 2015 6:12 PM, "Sparapani, Rodney" wrote:
> On Tue, 2015-02-17 at 15:57 -0500, Yixuan Qiu wrote:
> > Hello Rodney,
> >
> > If all you need is to pass a pointer to the fit() functi
:as(_x); // or a const-correct version
>
>double a=fit(x, n);
>
> ...
>
> Thanks,
>
> Rodney
>
>
> ___
> Rcpp-devel mailing list
> Rcpp-devel@lists.r-forge.r-project.org
> https://li
___
> 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
>
--
Yixuan Qiu
Department of Statistics,
Purdue University
___
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
Hello David,
The general answer to your question is, if the type of your matrix (integer
or numeric) in R is different from the one you declare in Rcpp, Rcpp will
make a copy and cast it to the appropriate type.
For example, 1:12 is of type integer, and if you pass it as a
NumericMatrix, Rcpp will
0.56 NA
> NA
> 2 b2() 10 128.785.762128.38 0.05 NA
> NA
> 3 b3() 10 22.351.000170.68 0.06 NA
> NA
>
>
> On Sun, Nov 23, 2014 at 9:15 PM, Yixuan Qiu wrote:
>
>> Hello Justace,
>> It se
t;
> For information I am running the following:
>
> - Windows 7
> - R: 3.1.2
> - Rtools: 3.1
> - Rcpp: 0.11.3
> - inline: 0.3.13
> - RStudio: 0.99.62
>
>
> ___
> Rcpp-devel mailing list
> Rcpp-devel@lists.r
>>}
>> Con+=-(k/2)*std::log(2*M_PI)
>
>
> it work well.
> I do not know what is wrong with
> *double Con = -(k/2)*std::log(2*M_PI)+rooti.diagonal().log().sum();*
>
>
> ___
> Rcpp-devel mailing list
> Rcpp-de
=CLC_TIME=German_Switzerland.1252
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> other attached packages:
> [1] devtools_1.6
>
> loaded via a namespace (and not attached):
> [1] fortunes_1.5-2 tools_3
nt Solutions
>
> Fiserv
>
> Office: 678-375-5315
>
> www.fiserv.com
>
>
>
>
>
> *From:* rcpp-devel-boun...@r-forge.wu-wien.ac.at [mailto:
> rcpp-devel-boun...@r-forge.wu-wien.ac.at] *On Behalf Of *Yixuan Qiu
> *Sent:* Wednesday, September 03, 2014 11:03 P
giving suggestions.
Happy Rcpping!
Best,
Yixuan
--
Yixuan Qiu
Department of Statistics,
Purdue University
___
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
> | ___
> | 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
>
> --
> http://dirk.eddelbuettel.com | @eddelbuettel | e...@debian.
is too seems daunting (even with
> RcppGSL). I think there may be integration support in the boost headers,
> but my head is too small for this yet.
>
> Any ideas that could help me?
>
> Dave Silkworth
>
>
> ___________
> Rcpp-
t this subsetting procedure was going to be easy...but I
> seem to have got stuck.
>
> thanks
>
> HLM
>
>
> ___
> 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
>
onths back in a threaded mail reader.
>
> --
> Dirk Eddelbuettel | e...@debian.org | http://dirk.eddelbuettel.com
> ___
> Rcpp-devel mailing list
> Rcpp-devel@lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/
rmail/rcpp-devel/2012-May/003792.html
>
> Thanks,
>
> Jianyang
>
> ___
> 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
>
--
Yixu
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 pr
I think that wonderful vignetter by Bates and Eddelbuettel describes
> this :-)
>
> :-)
>
> Dirk
>
> --
> Dirk Eddelbuettel | e...@debian.org | http://dirk.eddelbuettel.com
> _______
> Rcpp-devel mailing list
> Rcpp-devel@lists.r-forge.r-proj
27 matches
Mail list logo
