[Rcpp-devel] Vector allocation memory problem

[Cedric Ginestet]

Absolutely right!
How stupid of me.
Cheers,


On 08/04/11 13:18, romain at r-enthusiasts.com wrote:
> Hmmm. I think you just need "i<n" instead of "i<=n" in your for construct.
>
> Romain
>
>
>
> Le 8 avr. 2011 à 10:12 AM, Cedric Ginestet 
> <c.ginestet05 at googlemail.com <mailto:c.ginestet05 at googlemail.com>> a 
> écrit :
>
>> Hi guys,
>>
>> I have encountered a reproducible error in using Rcpp, which I have 
>> reported below. Basically, while allocation of elements in a vector 
>> to a particular value works well for vectors of small sizes, this 
>> returns an error when the size of that vector is increased. Either it 
>> returns a giblc error or it hangs and never resumes. I have 
>> reproduced this error on two different Linux distributions (Ubuntu 
>> 10.10 and 10.04 with g++ 4.4.1 and R 2.12.1).
>>
>> ####################################################
>> library(Rcpp)
>> library(inline)
>>
>> src <- '
>>   Rcpp::IntegerMatrix xA(A);
>>   int n=xA.nrow();
>>   NumericVector dist(n);
>>   for(int i=0;i<=n;i++){dist[i]=1.0;}
>>   return dist;'
>> cxxfun <- 
>> cxxfunction(sig=signature(A="matrix"),body=src,plugin="Rcpp",verbose=TRUE)
>>
>> ## Tests.
>> n <- 20; A <- matrix(0,n,n); dist <- cxxfun(A);
>> n <- 200; A <- matrix(0,n,n); dist <- cxxfun(A);
>> n <- 2000; A <- matrix(0,n,n); dist <- cxxfun(A);
>> *** glibc detected *** /usr/lib/R/bin/exec/R: double free or 
>> corruption (!prev): 0x08cfc278 ***
>>
>> ####################################################
>>
>> Is there something that I doing wrong there?
>> Thank you very much for your help,
>> Ced
>>
>>
>> -- 
>> Cedric Ginestet
>> Centre for Neuroimaging Sciences (L3.04)
>> NIHR Biomedical Research Centre
>> Institute of Psychiatry, Box P089
>> Kings College London
>> De Crespigny Park
>> London
>> SE5 8AF
>> _______________________________________________
>> Rcpp-devel mailing list
>> Rcpp-devel at lists.r-forge.r-project.org 
>> <mailto:Rcpp-devel at lists.r-forge.r-project.org>
>> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/rcpp-devel/attachments/20110408/7d546f0b/attachment.htm>