[Help-gsl] why is gsl_vector_const_view const?

[Brian Austin]

Hi,

I would really like to loop over vector views of a const gsl_matrix. 
Something like this:

void foo( const gsl_matrix *A ){
  gsl_vector_const_view a;
  for( i=0; i<A->size1; i++ ){
    a = gsl_matrix_const_row( A, i );
    //some operations that read elements of a.vector
  }
}//end foo

The compiler seems to dislike this because a is const. I don't want a to 
be const, just the data that a points to. Can anyone suggest a way to do 
this?

Thanks
-Brian