Re: [Rcpp-devel] Overflow error when creating a large matrix

2015-11-01 Thread Alessandro Mammana
hat I cannot point to the Rcpp code where this happens and propose the fix, I am sitting on a chair with my laptop on my legs, not quite convenient for programming :D. Best, Alessandro On Sun, Nov 1, 2015 at 7:35 PM, Alessandro Mammana wrote: > As Sean O'Riordain correctly pointed out,

Re: [Rcpp-devel] Overflow error when creating a large matrix

2015-11-01 Thread Alessandro Mammana
s utils datasets methods base > > other attached packages: > [1] Rcpp_0.12.1 > > On Sun, Nov 1, 2015 at 3:37 PM, Florian Plaza Oñate > wrote: >> The bug is probably fixed by this commit: >> https://github.com/RcppCore/Rcpp/commit/eb069bf5c20aaa4d38e0ca2897bf04e956

Re: [Rcpp-devel] Overflow error when creating a large matrix

2015-11-01 Thread Alessandro Mammana
ub.com/RcppCore/Rcpp/commit/eb069bf5c20aaa4d38e0ca2897bf04e956cb8c4c > > > Le 31/10/2015 20:08, Alessandro Mammana a écrit : >> >> Dear All, >> When creating a large matrix with the follwing code: >> >> #include >> // [[Rcpp::export]] >> Rcpp::In

[Rcpp-devel] Overflow error when creating a large matrix

2015-10-31 Thread Alessandro Mammana
nt between R and Rcpp in the allocation of a matrix? 2. In a 64 bits machine, what is actually the maximum allowed length of a vector/matrix? does the length of a vector/matrix need to be represented by an int? Thanks a lot and best regards, Alessandro -- Alessandro Mammana, PhD Student Max Plan

[Rcpp-devel] The proper way of wrapping an R object

2015-02-10 Thread Alessandro Mammana
to make sure that nothing is copied (as these vectors can be very long). 2. Even cleaner would be to start a method in this way: // [[Rcpp::export]] Rcpp::List asList(CountSignals x){ do you know if/how I can achieve that? (All of this code is inside an R package). Thanks a lot! --

Re: [Rcpp-devel] Rcpp-devel Digest, Vol 57, Issue 18

2014-07-17 Thread Alessandro Mammana
ng iterators and now the speedup from > parallel is about 5x (more in line with expectations). > -- next part ------ > An HTML attachment was scrubbed... > URL: > <http://lists.r-forge.r-project.org/pipermail/rcpp-devel/attachments/20140714/39077e47/attachment-000

Re: [Rcpp-devel] Sharing Cpp code between packages

2014-05-20 Thread Alessandro Mammana
Makevars :D On Tue, May 20, 2014 at 5:41 PM, Dirk Eddelbuettel wrote: > > On 20 May 2014 at 17:26, Alessandro Mammana wrote: > | Dear All, > | first of all thanks again for Rcpp, I am becoming more and more > | dependent on it and I am filling my directories with .cpp scripts... &

[Rcpp-devel] Sharing Cpp code between packages

2014-05-20 Thread Alessandro Mammana
description file however during compilation I got an error like: undefined symbol: _Z16rcpp_hello_worldv do you know what could be wrong? Thanks a lot and sorry for bothering. Ale -- Alessandro Mammana, PhD Student Max Planck Institute for Molecular Genetics Ihnestraße 63-73 D-14195 Berlin,

[Rcpp-devel] Default empty values for vector and matrices

2014-02-17 Thread Alessandro Mammana
n the future the special value R_NilValue could be converted to an empty vector, it would make things very easy and coherent with R's behaviour. Thx in advance, Ale -- Alessandro Mammana, PhD Student Max Planck Institute for Molecular Genetics Ihnes

Re: [Rcpp-devel] Extending Rcpp with a dummy class

2014-02-14 Thread Alessandro Mammana
when deploying packages... On Fri, Feb 14, 2014 at 1:31 PM, Dirk Eddelbuettel wrote: > > On 14 February 2014 at 11:02, Alessandro Mammana wrote: > | Thanks a lot, I misunderstood the document > | http://cran.r-project.org/web/packages/Rcpp/vignettes/Rcpp-extending.pdf > | . > | But th

Re: [Rcpp-devel] Extending Rcpp with a dummy class

2014-02-14 Thread Alessandro Mammana
/include in the Makevars and a header classdef.h in that directory, but this does not seem to add the desired line. One again I am lost... What should I do? On Thu, Feb 13, 2014 at 4:13 PM, Dirk Eddelbuettel wrote: > > On 13 February 2014 at 15:55, Alessandro Mammana wrote: > | Suppose th

[Rcpp-devel] Extending Rcpp with a dummy class

2014-02-13 Thread Alessandro Mammana
r defining my classes, so how should I do? Thanks in advance! Ale -- Alessandro Mammana, PhD Student Max Planck Institute for Molecular Genetics Ihnestraße 63-73 D-14195 Berlin, Germany ___ 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

Re: [Rcpp-devel] Segfault, is it because of iterators/pointers?

2014-02-12 Thread Alessandro Mammana
I like the "Strict mode" idea, I will use it, thanks! On Wed, Feb 12, 2014 at 2:34 PM, Romain Francois wrote: > > Le 12 févr. 2014 à 13:36, Alessandro Mammana a écrit : > >> Ah wait, my bad (as always T.T), I found a much simpler explanation: >> >> colset &

Re: [Rcpp-devel] Segfault, is it because of iterators/pointers?

2014-02-12 Thread Alessandro Mammana
e you try to cast to is not the one you expect? You might imagine that with large datasets this is important. Sorry for bothering and thanks again, Ale On Wed, Feb 12, 2014 at 1:10 PM, Dirk Eddelbuettel wrote: > > On 12 February 2014 at 11:47, Alessandro Mammana wrote: > | Ok I was able

Re: [Rcpp-devel] Segfault, is it because of iterators/pointers?

2014-02-12 Thread Alessandro Mammana
t? Thanks a lot for your help Ale On Tue, Feb 11, 2014 at 3:44 PM, Dirk Eddelbuettel wrote: > > In essence: "Yes" > > On 11 February 2014 at 15:18, Alessandro Mammana wrote: > | Hi all, > | I got another segfault using Rcpp. It is very difficult to understand > |

[Rcpp-devel] Segfault, is it because of iterators/pointers?

2014-02-11 Thread Alessandro Mammana
ry in the R area of memory with "IntegerVector res(pos.length())" and R decides to move the original values of "vec" to some other place, making the pointer invalid. Is that possible Sorry for being very vague and thx in advance!!! Ale -- Alessandro Mammana, PhD Student Max

[Rcpp-devel] Error when handling default values

2014-02-03 Thread Alessandro Mammana
out of scope it tries to free the memory allocated for the list, and the list tries to free the memory allocated for v1 and v2. But because ucs is not a normal stl container, I guess I need some other rules of thumb to understand how I should manage memory. Thanks a lot in advance, and, don't

[Rcpp-devel] Wierd compilation error

2014-01-29 Thread Alessandro Mammana
hpad.cpp is only 250 lines long... Do you know what could be wrong? Thanks in advance! Ale -- Alessandro Mammana, PhD Student Max Planck Institute for Molecular Genetics Ihnestraße 63-73 D-14195 Berlin, Germany ___ Rcpp-devel mailing list Rcpp-de

Re: [Rcpp-devel] Problems wrapping stl::vector

2014-01-03 Thread Alessandro Mammana
n: > > IntegerVector = as( wrap(a) ) ; > > Romain > > Le 3 janv. 2014 à 17:18, Alessandro Mammana a écrit : > >> Dear all, >> thank you again for this fantastic library. >> Sorry for the very stupid question, I looked before posting, but I >> couldn't

[Rcpp-devel] Problems wrapping stl::vector

2014-01-03 Thread Alessandro Mammana
rn-type] make: *** [scratchpad.o] Error 1 g++ -I/package/mariux64/R/R-3.0.2/lib/R/include -DNDEBUG -I/usr/local/include -I"/package/mariux64/R/R-3.0.2/lib/R/library/Rcpp/include"-fpic -g -O3 -Wall -pedantic -c scratchpad.cpp -o scratchpad.o Error in sourceCpp("scratchpad.cpp") :

[Rcpp-devel] Some beginner questions

2013-11-24 Thread Alessandro Mammana
); so that I am sure that the vector is not copied, the compiler complains saying that as(gr.slot("seqnames")) is an rvalue, and if I want to reference it, the reference should be constant. How do I create a non-constant reference to a slot of a s4 object then? If you made it through th

Re: [Rcpp-devel] Interfacing Rcpp with a C library: memory allocation

2013-11-21 Thread Alessandro Mammana
table for passing objects in and out. Is that right? Thanks again, Ale On Wed, Nov 20, 2013 at 7:25 PM, Dirk Eddelbuettel wrote: > > On 20 November 2013 at 18:54, Alessandro Mammana wrote: > | Thanks a lot for your help, I'll look at the references. I saw Section > | 6.1.2 a

Re: [Rcpp-devel] Interfacing Rcpp with a C library: memory allocation

2013-11-20 Thread Alessandro Mammana
emory errors, but this way they are handled by R? Thanks a lot! Ale On Wed, Nov 20, 2013 at 6:35 PM, Dirk Eddelbuettel wrote: > > On 20 November 2013 at 18:27, Alessandro Mammana wrote: > | Dear all, > | I'm trying to write some efficient code for analyzing sequencing data > |

[Rcpp-devel] Interfacing Rcpp with a C library: memory allocation

2013-11-20 Thread Alessandro Mammana
and much about memory allocation in Rcpp and I couldn't find many resources talking about it. Is there anything R- or Rcpp-specific that I have to keep in mind or should I program as if I were programming in C/C++? Thanks a lot! -- Alessandro Mammana, PhD Stud

[Rcpp-devel] Suggestions for possible optimizations

2013-10-23 Thread Alessandro Mammana
Dear Rcpp developers, I am optimizing some R code by replacing some loops with Rcpp code. After some profiling, it looks like the cpp function below is still the limiting step. Does anybody see some possible source of inefficiency or something easy to improve? Or is it already as fast as it can ge

Re: [Rcpp-devel] Issue with Matrices and Iterators

2013-10-14 Thread Alessandro Mammana
at). > > If instead you select a column mat( _, 0) you have a conversion to a > NumericVector. > I've encountered the same issue a couple of time. > > Matteo > > > On Mon, Oct 14, 2013 at 8:21 PM, Alessandro Mammana > wrote: >> >> Hi all, >> I hav

[Rcpp-devel] Issue with Matrices and Iterators

2013-10-14 Thread Alessandro Mammana
Hi all, I have very little experience with Rcpp and also with c++, so I am sorry if my question might be too stupid. I am trying to use iterators in my code, because I understand them better than vectors, because they should be efficient, and because I should be able to combine different implementa