help-gsl
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Help-gsl] Question about Vector Views


From: Srikant
Subject: [Help-gsl] Question about Vector Views
Date: Tue, 22 Jun 2010 21:59:28 -0700 (PDT)

Hi,

Consider the following code:

size_t rows = 10;
        
gsl_vector *v = gsl_vector_calloc(10);
for(size_t row=0;row<rows;row++)
{
    gsl_vector_set(v,row,row);
}
        
gsl_vector_view v_test = gsl_vector_subvector(v,0,3);

std::cout << gsl_vector_get(&v_test.vector,0) << std::endl;
std::cout << gsl_vector_get(&v_test.vector,1) << std::endl;
std::cout << gsl_vector_get(&v_test.vector,2) << std::endl;
std::cout << gsl_vector_get(&v_test.vector,3) << std::endl;

I was under the impression until now that GSL performs range checking on vector 
views as well and hence the last line should be an error. But, my test program 
suggests otherwise. The output of the above code is:

0
1
2
3 

whereas I was expecting it to be

0
1
2

Is this a bug or is vector view expecting as it is supposed to?

Thanks,
Srikant


      



reply via email to

[Prev in Thread] Current Thread [Next in Thread]