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
_______________________________________________
Help-gsl mailing list
Help-gsl@gnu.org
http://lists.gnu.org/mailman/listinfo/help-gsl