
Using Vectors with the GSL 229
10 const int N = 8;
11 int i;
12 double a[] = {4.45, 1.65, 6.25, 2.48, 1.22, 10.59, 11.22,
8.55};
13 int offset; /
*
starting element of subvector
*
/
14 int nsize; /
*
size of subvector
*
/
15 gsl_vector_view avv = gsl_vector_view_array(a, N); // view
16 printf("Program creates and manipulates a subvector\n\n");
17 printf("Data from vector\n");
18 for (i = 0; i < N; i++)
19 printf ("element %d = %g\n", i, gsl_vector_get
(&avv.vector, i));
20 offset = 2;
21 nsize = 4;
22 gsl_vector_view bvv = gsl_vector_subvector(&avv.vector,
offset, nsize);
23 printf("\nSubvector: \n");
24 for (i = 0; i < nsize; i++)
25 printf ("element %d = %g\n", i,